Sentinel-2 NDVI影像Resize后全为0值的解决方法求助

你的问题主要出在两个核心点：cv2.resize的尺寸参数顺序错误，以及NDVI数值范围与JPG格式不兼容，下面是具体的修正方案：

问题根源分析 #

1. 尺寸参数顺序颠倒：cv2.resize的目标尺寸参数是(宽度, 高度)，但image.shape返回的是(高度, 宽度)，你直接用new_size[0], new_size[1]会把宽高搞反，可能导致图像变形甚至后续处理异常。
2. NDVI与JPG格式不兼容：NDVI的正常数值范围是[-1, 1]，而JPG仅支持0-255的8位无符号整数，直接保存会把所有负数截断为0，浮点值也会被强制转换，最终出现全0的结果。

修正后的CV2实现代码 #

import cv2
import numpy as np

# 读取原NDVI影像（实际项目建议用Rasterio/GDAL读取浮点型栅格）
image1 = cv2.imread("原NDVI文件路径", cv2.IMREAD_UNCHANGED).astype(np.float32)

# 计算目标尺寸：确保宽高都能被5整除，注意顺序是(宽度, 高度)
height, width = image1.shape
new_width = (width // 5) * 5
new_height = (height // 5) * 5

# 执行Resize，选择适合浮点型的插值方法
resized_ndvi = cv2.resize(image1, (new_width, new_height), interpolation=cv2.INTER_LINEAR)

# 方法1：将NDVI映射到0-255范围后保存为JPG
scaled_ndvi = (resized_ndvi + 1) * 127.5  # 把[-1,1]转成[0,255]
scaled_ndvi = scaled_ndvi.astype(np.uint8)
cv2.imwrite("缩放后的NDVI_8位.jpg", scaled_ndvi)

# 方法2：保存为支持浮点型的TIFF格式（推荐，保留原始精度）
import rasterio
from rasterio.transform import from_origin

# 从原影像读取元数据
src = rasterio.open("原NDVI文件路径.tif")
profile = src.profile
profile.update(width=new_width, height=new_height, dtype='float32')

with rasterio.open("缩放后的NDVI_浮点.tif", 'w', **profile) as dst:
    dst.write(resized_ndvi, 1)

Rasterio专业Resize方案（更适合遥感影像） #

用Rasterio能直接处理栅格的地理信息，同时保留数值精度，是处理遥感影像的首选：

import rasterio
from rasterio.enums import Resampling

with rasterio.open("原NDVI文件路径.tif") as src:
    # 计算缩放比例，确保输出尺寸能被5整除
    scale_width = ((src.width // 5) * 5) / src.width
    scale_height = ((src.height // 5) * 5) / src.height
    
    # 执行Resize
    data = src.read(
        out_shape=(
            src.count,
            int(src.height * scale_height),
            int(src.width * scale_width)
        ),
        resampling=Resampling.bilinear
    )
    
    # 更新元数据（保留地理坐标）
    profile = src.profile
    profile.update(
        height=int(src.height * scale_height),
        width=int(src.width * scale_width),
        transform=src.transform * src.transform.scale(scale_width, scale_height)
    )
    
    # 保存结果
    with rasterio.open("缩放后的NDVI.tif", 'w', **profile) as dst:
        dst.write(data)

关键注意事项 #

• 禁止直接将浮点型NDVI数组保存为JPG，必须先映射到0-255范围，或保存为TIFF/ENVI等支持浮点的栅格格式。
• 使用OpenCV处理栅格时，务必注意(宽度, 高度)的参数顺序，避免尺寸颠倒。
• 处理遥感影像优先选择Rasterio/GDAL这类专业库，能同时保留地理坐标和数值精度，比OpenCV/Pillow更可靠。

内容的提问来源于stack exchange，提问作者gis_helydon