如何用Python或Java将3波段Float32类型GeoTIFF转为JPG？

我来帮你搞定这个GeoTIFF格式转换的问题～你遇到的情况确实大概率是Float32类型的波段数据范围超出了普通图像格式的有效值区间（比如8位图像仅支持0-255，16位是0-65535），直接转图像自然会失败。下面给你一套可行的处理流程：

步骤1：完整读取所有波段数据 #

先把三个波段的Float32数据完整读出来，注意保留原始类型：

from osgeo import gdal
import numpy as np

# 读取影像
ds = gdal.Open('你的影像路径.tif', gdal.GA_ReadOnly)
if not ds:
    raise Exception("无法打开影像文件！")

# 读取三个波段的数组
band1 = ds.GetRasterBand(1).ReadAsArray().astype(np.float32)
band2 = ds.GetRasterBand(2).ReadAsArray().astype(np.float32)
band3 = ds.GetRasterBand(3).ReadAsArray().astype(np.float32)

# 关闭数据集释放资源
ds = None

步骤2：对影像波段（1、2）做归一化处理 #

因为Float32的数值范围可能跨度极大（比如包含负数、远大于255的值），需要把它们缩放到目标图像格式的有效区间。这里以转8位PNG/JPG为例，用线性归一化把数据映射到0-255：

def normalize_band(band_array):
    # 计算波段的最小/最大值（如果有异常值，可先过滤再计算）
    min_val = np.nanmin(band_array)
    max_val = np.nanmax(band_array)
    
    # 避免除以0的极端情况
    if max_val == min_val:
        return np.zeros_like(band_array, dtype=np.uint8)
    
    # 线性缩放到0-255区间
    normalized = ((band_array - min_val) / (max_val - min_val)) * 255
    return normalized.astype(np.uint8)

# 对波段1、2执行归一化
norm_band1 = normalize_band(band1)
norm_band2 = normalize_band(band2)

步骤3：处理入射角波段（波段3） #

入射角通常也是Float32类型，处理方式分两种情况：

• 如果入射角范围是0-90度，可直接转成8位整数（因为90<255）
• 如果入射角范围更大，同样用上面的归一化函数处理

# 情况1：入射角范围0-90，直接转uint8
angle_band = band3.astype(np.uint8)

# 情况2：入射角范围更大，用归一化函数
# angle_band = normalize_band(band3)

步骤4：合并波段并保存为目标格式 #

用PIL合并波段后保存成普通图像格式，示例如下：

from PIL import Image

# 合并三个波段（注意PIL默认通道顺序是RGB，按需调整顺序）
merged_array = np.dstack((norm_band1, norm_band2, angle_band))

# 转成PIL图像并保存
img = Image.fromarray(merged_array)
img.save('转换后的图像.png')

额外提示 #

• 如果影像存在NoData值，记得在归一化前先处理（比如用np.nan替换，或者过滤后再计算min/max）
• 如果要保存回GeoTIFF格式而非普通图像，不需要归一化，直接用GDAL创建新数据集写入即可：

# 示例：保存回Float32格式的GeoTIFF
driver = gdal.GetDriverByName('GTiff')
out_ds = driver.Create('输出路径.tif', band1.shape[1], band1.shape[0], 3, gdal.GDT_Float32)
# 复制原影像的投影和地理变换信息
out_ds.SetProjection(ds.GetProjection())
out_ds.SetGeoTransform(ds.GetGeoTransform())
# 写入三个波段的数据
out_ds.GetRasterBand(1).WriteArray(band1)
out_ds.GetRasterBand(2).WriteArray(band2)
out_ds.GetRasterBand(3).WriteArray(band3)
# 关闭数据集
out_ds = None

内容的提问来源于stack exchange，提问作者송준석