Python投影转换问题

我目前接触到的遥感数据投影类型：

（1）AMSR-E/AMSR2 L3全球亮温数据：EQR(Equidistant Cylindrical) 等经纬度投影（最简单）

（2）哨兵1号雪深数据产品：1 km global cylindrical 1-km Equal-Area Scalable Earth version 2 (EASE-2) 1km的等面积圆柱投影重采样到EASE2网格中

 

 （3）积雪覆盖度产品MOD10A1 006: 正弦投影（带有角点坐标）

我下载了一个小程序计算出了每个像元的中心经纬度(记不太清了)，MRT也可以，不过MRT计算的是顶点坐标

（4）GLDAS LST/AT产品数据：0.25°等经纬度投影

（5）SSM/I、SSMIS亮温数据：EASE2_T

19GHZ用的是6.25km分辨率

 

 37GHZ和85GHZ用的是3.125km分辨率：

 

 （6）其他一些植被覆盖度、DEM、积雪覆盖度等用的都是WGS84经纬度投影

我用到的投影转换方法：

（1）全部变量都转化为0.1°等经纬度投影

方法：循环经纬度，变量在0.1×0.1范围内就进行平均

（2）全部变量都转化为等面积圆柱EASE2网格

方法：用gdal.Warp函数进行转换，转换后的数据由lat/lon变为m,用米数读数据

代码:

option4 = gdal.WarpOptions(format='GTiff', srcSRS='EPSG:4326', dstSRS='EPSG:6933')
FVC = gdal.Warp("E:/mask/Tibetmask_6933.tif", "E:/mask/Tibetmask_4326.tif", options=option4)