Python读取ETOPO TIFF转NumPy数组获取像素值失败求助

ETOPO2022的15s地表高程影像是32位浮点型GeoTIFF，带有地理空间元数据，普通影像库容易因未处理这些元数据导致读取异常。以下是几种可行的解决方案：

方案1：使用rasterio（推荐） #

rasterio是专门处理地理栅格数据的库，能完美解析GeoTIFF：

1. 安装依赖：

pip install rasterio

2. 读取代码：

import rasterio
import numpy as np

with rasterio.open(filepath) as src:
    # 读取单波段高程数据
    elevation_data = src.read(1)
    # 验证数组形状（应为(3600, 3600)）
    print(elevation_data.shape)
    # 查看具体像素值
    print(elevation_data[0, 0])

方案2：修复PIL读取逻辑 #

如果坚持用PIL，需确保加载完整数据并处理32位浮点模式：

from PIL import Image
import numpy as np

# 打开影像
im = Image.open(filepath)
# 强制加载像素数据（PIL默认懒加载）
im.load()
# 确认模式应为'F'（32位浮点灰度）
print(im.mode)
# 转换为NumPy数组
elevation_data = np.array(im)
print(elevation_data.shape)

若仍异常，升级PIL到最新版本：

pip install --upgrade pillow

方案3：使用GDAL #

GDAL也是地理空间栅格处理的常用工具：

1. 安装依赖：

pip install gdal

2. 读取代码：

from osgeo import gdal
import numpy as np

ds = gdal.Open(filepath)
band = ds.GetRasterBand(1)
elevation_data = band.ReadAsArray()
print(elevation_data.shape)

问题原因说明 #

你遇到的数组大小为1的问题，本质是普通影像库（PIL、imageio、cv2）未正确识别GeoTIFF的结构，误将影像的元数据部分当成了像素数据，或是未加载完整的栅格波段。而rasterio/GDAL专门针对地理空间数据设计，能正确解析这类带元数据的TIFF文件。

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