使用pyproj和cartopy将SEVIRI卫星图像像素转换为经纬度值

核心逻辑：线性映射 #

SEVIRI静止卫星的该区域图像采用正射圆柱投影（Plate Carrée），在南北纬30°范围内，像素与经纬度呈严格线性对应关系，无需复杂投影公式，直接通过比例换算即可完成转换。

前提确认 #

• 图像尺寸：1024×1024像素
• 地理范围：经度30°E ~ 105°E，纬度30°N ~ -30°S
• 像素原点：默认左上角为(0,0)，x轴向右（经度递增），y轴向下（纬度递减）

转换步骤 #

1. 计算经度

• 经度总跨度：105° - 30° = 75°
• 每像素经度增量：75° / 1024 ≈ 0.0732421875°/px
• 转换公式：

  经度 = 30° + (x像素值 × 0.0732421875)
  

• 示例（x=800px）：
  30 + 800×0.0732421875 = 88.59375°E

2. 计算纬度

• 纬度总跨度：30° - (-30°) = 60°
• 每像素纬度增量：60° / 1024 ≈ 0.05859375°/px
• 转换公式（因y轴向下对应纬度降低，用初始纬度减去像素增量）：

  纬度 = 30° - (y像素值 × 0.05859375)
  

• 示例（y=305px）：
  30 - 305×0.05859375 = 12.12890625°N

特殊情况说明 #

如果图像像素原点为左下角（y轴向上为纬度递增），纬度公式需调整为：

纬度 = -30° + (y像素值 × 0.05859375)

实际应用中，SEVIRI发布的常规图像基本以左上角为原点，优先使用前文公式。

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