实验报告
班级： 姓名： 学号：
大气遥感上机实验
实验名称
一、 实验目的
1、 2、 3、 4、
了解常用数据，如MODIS数据等； 了解MODIS相关产品的特性；
掌握MATLAB软件对遥感数据的处理过程； 掌握ENVI软件对遥感数据的处理。
二、 实验数据
FY2D_FDI_ALL_NOM_20080329_0330.HDF
MOD03.A2008125.0320.005.2008125132121.hdf MOD05_L2.A2008125.0320.005.2008126075028.hdf out.clr cloud.pro
MOD35_L2.A2008125.0320.005.2008126050344.hdf FY2D_FDI_ALL_NOM_20080329_0330.HDF
FY3A_MERSI_GBAL_L1_20090921_0340_1000M_MS.HDF AVIRIS_1998_scale.txt JasperRidge98av.hdr JasperRidge98av.img
JasperRidge98av_template.txt LandsatTM_JasperRidge_hrf.fst LandsatTM_JasperRidge_b10.fst LandsatTM_JasperRidge_b20.fst LandsatTM_JasperRidge_b30.fst LandsatTM_JasperRidge_b40.fst LandsatTM_JasperRidge_b50.fst LandsatTM_JasperRidge_b70.fst JasperRidgeTM_template.txt
三、 实验内容（实验主要内容，软件操作的主要过程截图及实验结果图）
1、 基于Matlab的hdf格式卫星资料读取/成像
1）hdf5的组成信息
>> fileinfo = hdf5info('FY2D_FDI_ALL_NOM_20080329_0330.HDF') toplevel = fileinfo.GroupHierarchy dataset1 = toplevel.Datasets(1) fileinfo =
Filename: 'FY2D_FDI_ALL_NOM_20080329_0330.HDF' LibVersion: '1.6.2'
Offset: 0
FileSize: 220021815 GroupHierarchy: [1x1 struct] toplevel =
Filename: 'FY2D_FDI_ALL_NOM_20080329_0330.HDF' Name: '/' Groups: []
Datasets: [1x19 struct] Datatypes: [] Links: [] Attributes: [] dataset1 =
Filename: 'FY2D_FDI_ALL_NOM_20080329_0330.HDF' Name: '/CALChannelIR1' Rank: 2
Datatype: [1x1 struct] Dims: [1 1024] MaxDims: [1 1024] Layout: 'contiguous' Attributes: [1x6 struct] Links: [] Chunksize: [] FillValue: 0 2）成图
IR1=hdf5read('FY2D_FDI_ALL_NOM_20080329_0330.HDF','/NOMChannelIR1'); IR2=hdf5read('FY2D_FDI_ALL_NOM_20080329_0330.HDF','/NOMChannelIR2'); IR3=hdf5read('FY2D_FDI_ALL_NOM_20080329_0330.HDF','/NOMChannelIR3'); >> fid = fopen('fy2d_IR.raw','wb') fwrite(fid, IR1,'uint16'); fwrite(fid, IR2,'uint16'); fwrite(fid, IR3,'uint16'); fclose(fid);
R=mat2gray(IR1',[0 1023]); G=mat2gray(IR2',[0 1023]); B=mat2gray(IR3',[0 1023]); >> imwrite(R,'IR1.jpg'); >> imagesc(IR1',[0 1023]); colormap(gray); axis equal;
>> RGB=cat(3,R,G,B); imwrite(RGB,'rgb.jpg'); >> image(RGB); axis equal;
灰度图
彩色图
2、 用IDL编写程序读取MODIS热红外波段反演的大气水汽
pro extract_data_
hdfid=hdf_sd_start('MOD05_L2.A2008125.0320.005.2008126075028.hdf') index=hdf_sd_nametoindex(hdfid,'Water_Vapor_Near_Infrared') varid=hdf_sd_select(hdfid,index)
hdf_sd_getdata,varid,jin hdf_sd_endaccess,varid hdf_sd_end,hdfid
hdfid=hdf_sd_start('MOD03.A2008125.0320.005.2008125132121.hdf') index=hdf_sd_nametoindex(hdfid,'Longitude') varid=hdf_sd_select(hdfid,index) hdf_sd_getdata,varid,lon hdf_sd_endaccess,varid hdf_sd_end,hdfid
hdfid=hdf_sd_start('MOD03.A2008125.0320.005.2008125132121.hdf') index=hdf_sd_nametoindex(hdfid,'Latitude') varid=hdf_sd_select(hdfid,index) hdf_sd_getdata,varid,lat hdf_sd_endaccess,varid hdf_sd_end,hdfid a=size(jin)
openw,1,'F:\\Temp\\data.txt' for ii=0,a[1]-1 do begin for jj=0,a[2]-1 do begin
printf,1,lon[ii,jj],lat[ii,jj],jin[ii,jj]*0.001 endfor endfor close,1 end
完成后用Surfer将结果画图
大气遥感实验报告-CUIT - 图文
