IDL读取MODIS 1B数据中500分辨率下的数据并显示成伪彩

最新推荐文章于 2024-12-06 20:44:58 发布

原创最新推荐文章于 2024-12-06 20:44:58 发布 · 496 阅读

0 ·

CC 4.0 BY-SA版权

envi/idl 专栏收录该内容

8 篇文章

订阅专栏

本文介绍如何使用IDL语言读取HDF文件中的特定数据集，并进行图像缩放及伪彩色显示。通过具体代码示例，展示了从选择文件到最终图像展示的全过程。

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