如何获取MCD12Q1的经纬度信息

寂静山林

本人是用MATLAB的，下载的MCD12Q1的HDF数据不带经纬度信息。根据MODIS说明书的计算方法（如下图），计算出来的经纬度跟MRT里的有差距，差别还是挺大的，有没有人做过？请指教。
我下载了MRT和HCG但也只是看到了四个角的经纬度信息，能导出matlab可用的经纬度信息么？

学习渣渣

MRT 好像是默认计算每个小格子左上角的经纬度，反正计算 就是 读取四个角的经纬度 然后 用这个公式。

寂静山林

学习渣渣 发表于 2020-11-3 22:49
MRT 好像是默认计算每个小格子左上角的经纬度，反正计算 就是 读取四个角的经纬度 然后 用这个公式。

怎么得到每个格点左上角的信息呢？MRT能输出么？

MeteoInfo

寂静山林 发表于 2020-11-4 09:07
怎么得到每个格点左上角的信息呢？MRT能输出么？

可以用MeteoInfoLab脚本将数据Sinusoidal投影的x, y坐标投影转化为经纬度，并输出到文件：

1. fn = 'D:/Temp/hdf/MCD12Q1.A2016001.h26v04.006.2018055072913.hdf'
   
2. f = addfile(fn)
   
3. data = f['LC_Type1'][:]
   
4. x = data.dimvalue(1)
   
5. y = data.dimvalue(0)
   
6. x, y = meshgrid(x, y)
   
7. 
   
8. #Project
   
9. toproj = geolib.projinfo()
   
10. lon, lat = geolib.project(x, y, f.proj, toproj=toproj)
    
11. 
    
12. #Write to file
    
13. binwrite('D:/Temp/lon.bin', lon)
    
14. binwrite('D:/Temp/lat.bin', lat)

fn = 'D:/Temp/hdf/MCD12Q1.A2016001.h26v04.006.2018055072913.hdf'<br /> f = addfile(fn)<br /> data = f['LC_Type1'][:]<br /> x = data.dimvalue(1)<br /> y = data.dimvalue(0)<br /> x, y = meshgrid(x, y)<br /> <br /> #Project<br /> toproj = geolib.projinfo()<br /> lon, lat = geolib.project(x, y, f.proj, toproj=toproj)<br /> <br /> #Write to file<br /> binwrite('D:/Temp/lon.bin', lon)<br /> binwrite('D:/Temp/lat.bin', lat)

寂静山林

MeteoInfo 发表于 2020-11-4 10:39
可以用MeteoInfoLab脚本将数据Sinusoidal投影的x, y坐标投影转化为经纬度，并输出到文件：

我试试，非常感谢

寂静山林

MeteoInfo 发表于 2020-11-4 10:39
可以用MeteoInfoLab脚本将数据Sinusoidal投影的x, y坐标投影转化为经纬度，并输出到文件：

我按照你的写法为什么不对呢？我不会meteoinfo所以还得麻烦你看一下J:\图片1.png

寂静山林

寂静山林 发表于 2020-11-5 09:03
我按照你的写法为什么不对呢？我不会meteoinfo所以还得麻烦你看一下

我改成data = f['LC_Type'][:,:]后错误就只剩module那部分了，这部分我就又不会了，还请王老师指教