作者fxxkjoe1231 (鲁到入味)
看板GIS
标题[请益] GDAL处理MOD09产品
时间Tue Jul 23 17:25:52 2019
大家好，最近在用GDAL处理MOD09CMA的产品。
用gdal_translate将HDF格式中要的资料转为GeoTiff後，用gdalinfo去查看他的座标资讯如下：
Coordinate System is:
GEOGCS["Unknown datum based upon the Clarke 1866 ellipsoid",
DATUM["Not_specified_based_on_Clarke_1866_spheroid",
SPHEROID["Clarke 1866",6378206.4,294.9786982139006,
AUTHORITY["EPSG","7008"]]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]]
Origin = (-180.000000000000000,90.000000000000000)
Pixel Size = (0.050000000000000,-0.050000000000000)
找不到对应的EPSG码，去看产品的手册，上面只写到：
its projection is in linear latitude and longitude (Plate Carre), and its resolution is 0.05°
http://modis-sr.ltdri.org/guide/MOD09_UserGuide_v1_3.pdf
用关键字去找，找到了EPSG:32662，却发现已经被弃用，要用32663，但我的GDAL不支援，所以索性用32662跑，但进行重新投影时都会产生错误如下
gdalwarp.exe -overwrite -s_srs EPSG:32662 -t_srs EPSG:4326 -r bilinear -tr 0.05 0.05 AOD550_CROP.tif AOD550_CROP_RP.tif
ERROR 1: Too many points (441 out of 441) failed to transform,
unable to compute output bounds.
上网查了一下，很大的机率是来源的投影给错了。
GDAL版本是GDAL 1.10.1, released 2013/08/26
想请问是否有人知道MOD09产品对应的EPSG码或是我的作法哪部分出了问题，
谢谢大家
-----
Sent from JPTT on my Sony G8441.
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 101.15.194.98 (台湾)
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/GIS/M.1563873959.A.1EB.html
※ 编辑: fxxkjoe1231 (1.169.142.116 台湾), 07/24/2019 00:27:45
1F：推 psilotum: 参数可以用 prj2epsg.org 查看看，看起来像EPSG:4008 09/01 12:17
回复一下後来的作法：
用MRT(MODIS REPROJECTION TOOL)重新投影，从MRT所生成的报表内显示MOD09CMA应为
EPSG3826。用GDAL给定EPSG3826与MRT分别进行重新投影後，两者pattern是一致的，
但不完全重和，略有偏移。
※ 编辑: fxxkjoe1231 (36.230.38.25 台湾), 09/13/2019 03:13:13