作者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/m.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