资源描述:
《大气遥感实验报告-CUIT》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、实验报告班级:姓名:学号:实验名称大气遥感上机实验一、实验目的1、了解常用数据,如MODIS数据等;2、了解MODIS相关产品的特性;3、掌握MATLAB软件对遥感数据的处理过程;4、掌握ENVI软件对遥感数据的处理。二、实验数据FY2D_FDI_ALL_NOM_20080329_0330.HDFMOD03.A2008125.0320.005.2008125132121.hdfMOD05_L2.A2008125.0320.005.2008126075028.hdfout.clrcloud.proMOD35_L2.A2008125.0
2、320.005.2008126050344.hdfFY2D_FDI_ALL_NOM_20080329_0330.HDFFY3A_MERSI_GBAL_L1_20090921_0340_1000M_MS.HDFAVIRIS_1998_scale.txtJasperRidge98av.hdrJasperRidge98av.imgJasperRidge98av_template.txtLandsatTM_JasperRidge_hrf.fstLandsatTM_JasperRidge_b10.fstLandsatTM_JasperRidg
3、e_b20.fstLandsatTM_JasperRidge_b30.fstLandsatTM_JasperRidge_b40.fstLandsatTM_JasperRidge_b50.fstLandsatTM_JasperRidge_b70.fstJasperRidgeTM_template.txt三、实验内容(实验主要内容,软件操作的主要过程截图及实验结果图)1、基于Matlab的hdf格式卫星资料读取/成像1)hdf5的组成信息>>fileinfo=hdf5info('FY2D_FDI_ALL_NOM_20080329_033
4、0.HDF')toplevel=fileinfo.GroupHierarchydataset1=toplevel.Datasets(1)fileinfo=Filename:'FY2D_FDI_ALL_NOM_20080329_0330.HDF'LibVersion:'1.6.2'Offset:0FileSize:220021815GroupHierarchy:[1x1struct]toplevel=Filename:'FY2D_FDI_ALL_NOM_20080329_0330.HDF'Name:'/'Groups:[]Datase
5、ts:[1x19struct]Datatypes:[]Links:[]Attributes:[]dataset1=Filename:'FY2D_FDI_ALL_NOM_20080329_0330.HDF'Name:'/CALChannelIR1'Rank:2Datatype:[1x1struct]Dims:[11024]MaxDims:[11024]Layout:'contiguous'Attributes:[1x6struct]Links:[]Chunksize:[]FillValue:02)成图IR1=hdf5read('FY2
6、D_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');f
7、write(fid,IR3,'uint16');fclose(fid);R=mat2gray(IR1',[01023]);G=mat2gray(IR2',[01023]);B=mat2gray(IR3',[01023]);>>imwrite(R,'IR1.jpg');>>imagesc(IR1',[01023]);colormap(gray);axisequal;>>RGB=cat(3,R,G,B);imwrite(RGB,'rgb.jpg');>>image(RGB);axisequal;灰度图彩色图1、用IDL编写程序读取MOD
8、IS热红外波段反演的大气水汽proextract_data_hdfid=hdf_sd_start('MOD05_L2.A2008125.0320.005.2008126075028.hdf')index=hdf_sd_nametoin