| 查看: 778 | 回复: 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 |
回复此楼» 猜你喜欢
Seismic stratigraphy
已经有1人回复
-
长江大学石油工程学院佘跃惠团队招收2025年9月入学的博士研究生
已经有2人回复
-
地球化学论文润色/翻译怎么收费?
已经有85人回复
-
祈福，祝自己今年好运
已经有100人回复
-
关于会评
已经有1人回复
-
据说！出现这些情况，国基金会评本子会被拿下！
已经有19人回复
-
无人打捞，直接陪跑
已经有16人回复
-
定了定了，明天出结果！2025 国自然基金结果即将揭晓，查询攻略看这里
已经有107人回复
-
2025国家自然科学基金放榜时间，根据官方回应，8月15日基本没戏了
已经有5人回复
-
各位未来院士帮忙看看评审意见
已经有9人回复
-
《风吹过日常的缝隙》
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助各省边界的经纬度数据？？
已经有5人回复
- 知道遥感数据的四个顶点的经纬度，如何不做几何校正直接提取想要的经纬度位置的值
已经有5人回复
- 如何将经纬度数据生成Vector文件啊？
已经有12人回复
- 经纬度数据如何转换为80坐标？
已经有4人回复
- 经纬度数据如何转换为北京六分带投影坐标系
已经有3人回复
- 能否基于DEM数据提取经纬度？
已经有9人回复
2楼2017-12-20 22:19:14