MODIS投影转换 Python

mmmmmyyyyyhhhhh

已于 2023-04-25 15:50:56 修改

阅读量787

点赞数

文章标签： python

于 2023-04-19 12:08:00 首次发布

本文链接：https://blog.youkuaiyun.com/qq_42655653/article/details/130240697

版权

本文介绍了如何使用Python处理MODIS遥感数据的投影转换问题。针对无坐标系信息的情况，需要借助MOD03地理定位数据获取1km经纬度，并依据纬度范围定义Albers或Polar Stereographic投影。通过cKDTree构建地理坐标，并利用geotransform进行数据转换。而对于有投影信息的产品，如MOD13A2，可以直接构建geotransform并转化成tif。此外，还提及了gdal.Warp()函数用于转换投影。以MOD10A1的积雪覆盖产品为例，展示了NDSI_Snow_Cover的处理过程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

无坐标系信息

如MOD21_L2

hdf文件有5km的经纬度信息，但是此产品像元大小是1km，需要用到MOD03地理定位数据，读取它的经纬度（1km）。

swath_data = self.read_mxd21(mxd21_file, self.bands,sitelat,sitelon)
lat_mat, lon_mat = self.read_mxd03_latlon(mxd03_file)
swath_data.lats = lat_mat
swath_data.lons = lon_mat

没有投影信息，自己定义一个，纬度<65是Albers，>65是Polar Stereographic。分别为

proj = (f"+proj=aea +datum=WGS84 +ellps=WGS84 "
        f"+lon_0={lon} +lat_0={lat} +lat_1={stdpr1} +lat_2={stdpr2}")

proj = f"+proj=stere +datum=WGS84 +ellps=WGS84 lat_0={lat_0} +lat_ts={lat}"

已知lat，lon可以得到地理坐标࿰