| 查看: 839 | 回复: 1 | ||
ZZM_ZZM铁虫 (初入文坛)
|
[求助]
trmm数据经纬度问题
|
|
下载的是TRMM的3b43数据，用idl怎么校正经纬度，下面的程序怎么用，哪里有问题 SATE_TRMM_3B43_CASE.pro程序文件 pro SATE_TRMM_3B43_CASE ; 1*1440*400 FILENAME='SATE_L3_TRM_MUTDS_MWB_3B43_GLB_V6-199807.HDF' print,'*',FILENAME SD_ID=HDF_SD_START(FILENAME,/READ) P='precipitation' ; range:0-10, unit:mm hr-1 SDS_INDEX_P=HDF_SD_NAMETOINDEX(SD_ID,P) SDS_ID_P=HDF_SD_SELECT(SD_ID,SDS_INDEX_P) HDF_SD_GETDATA,SDS_ID_P,DATA_P HDF_SD_ENDACCESS,SDS_ID_P HDF_SD_END,SD_ID PRINT,DATA_P(0,0) print,size(DATA_P) OPENW,LUN,'0_50N_0_180E_199807_output.txt',/GET_LUN ; 0-50N,0-180E east asia FOR J=0,399 DO BEGIN FOR I=0,1439 DO BEGIN IF J gt 200 THEN BEGIN ;0-50N (J>200) IF I gt 720 THEN BEGIN ;0-180E (1440<I<720), asia, Australia PRINTF,LUN,(J+1)/4.0-50.125,(I+1)/4.0-180.125,DATA_P(J,I), $ format='(2(f7.3,2x),f7.3,2x)' ; lat, lon, precipitation ENDIF ENDIF ;ELSE BEGIN ; 180-360E --> -180~0,, America RINTF,LUN,(J+1)/4.0-50.125,(I+1)/4.0-180.125,DATA_P(J,I,0), $ ;format='(2(f6.2,2x),f6.4,2x)' ;ENDELSE ;1-1440, I, --> 0.125-359.875, Longitude ;1-720 -> (Europea, America), 0.125-179.875W, Lon=I/4.0-180.125<0; 721-1440 -> (Asia, Australia)0.125-179.875E, Lon=I/4.0-180.125<180 ; ;1-400, J, -->49.875N-49.875S, Latitude ;1-200 ->Lat=J/4.0-50.125<0, 49.875S-0.125S; 201-400 -> Lat=J/4.0-50.125>0, 0.125N-49.875N ENDFOR ENDFOR CLOSE,LUN print,'end' END |
回复此楼» 猜你喜欢
SCIENCE ADVANCE投稿状态咨询
已经有5人回复
-
地质学复试调剂
已经有1人回复
-
地质学论文润色/翻译怎么收费?
已经有103人回复
-
明年散金币
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助各省边界的经纬度数据？？
已经有5人回复
- 知道遥感数据的四个顶点的经纬度，如何不做几何校正直接提取想要的经纬度位置的值
已经有5人回复
- 如何将经纬度数据生成Vector文件啊？
已经有12人回复
- 经纬度数据如何转换为80坐标？
已经有4人回复
- 经纬度数据如何转换为北京六分带投影坐标系
已经有3人回复
- 能否基于DEM数据提取经纬度？
已经有9人回复
2楼2017-12-20 22:19:14
投票:
30