[請益] GDAL處理MOD09產品
大家好,最近在用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://www.ptt.cc/bbs/GIS/M.1563873959.A.1EB.html
※ 編輯: fxxkjoe1231 (1.169.142.116 臺灣), 07/24/2019 00:27:45
推
09/01 12:17,
4年前
, 1F
09/01 12:17, 1F
回復一下後來的作法:
用MRT(MODIS REPROJECTION TOOL)重新投影,從MRT所生成的報表內顯示MOD09CMA應為
EPSG3826。用GDAL給定EPSG3826與MRT分別進行重新投影後,兩者pattern是一致的,
但不完全重和,略有偏移。
※ 編輯: fxxkjoe1231 (36.230.38.25 臺灣), 09/13/2019 03:13:13