大气遥感实验报告-CUIT - 图文

实验报告

班级： 姓名： 学号：

大气遥感上机实验

实验名称

一、 实验目的

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将结果画图