IDL读取MODIS 1B数据中500分辨率下的数据 并显示成伪彩_武汉之光-人机交互小盒...

来源:百度文库 编辑:神马文学网 时间:2024/04/30 08:56:37
查看文章
IDL读取MODIS 1B数据中500分辨率下的数据 并显示成伪彩
2009-04-02 10:55
pro readhdf_cr
; 读取某个HDF文件
FILE_NAME=dialog_pickfile(default_extension= '.hdf',filter='*.hdf')
;读取其中的某个Data Set里面的数据,可以用HDF Explore看一下HDFl里面的文件结构
SDS_NAME="EV_250_Aggr500_RefSB"
AeroWidth=2708 ;影像的宽度
AeroLength=4060 ;影像的高度
;打开HDF文件并初始化HDF接口
sd_id = HDF_SD_START( FILE_NAME, /read )
sds_index = HDF_SD_NAMETOINDEX(sd_id,SDS_NAME)
sds_id = HDF_SD_SELECT( sd_id, sds_index )
HDF_SD_GETDATA, sds_id, data
HDF_SD_ENDACCESS, sds_id
;关闭HDF文件
HDF_SD_END, sd_id
imagedata = bytscl(congrid(data,AeroWidth*0.2,AeroLength*0.2,1));按比例缩小影像
DEVICE,DECOMPOSED=0 ;显示伪彩
LOADCT,36 ;加载第37个调色板
window, 0, Title='Aeromatex',xsize=AeroWidth*0.2, ysize=AeroLength*0.2, retain=2
TVSCL,imagedata
end