envi影像镶嵌,拼接线色彩平衡

envi影像镶嵌,拼接线/色彩平衡2010-02-0414:34（包括：将以像元为基础的图像或应用地理坐标系的图像镶嵌起来两种方式）1）基于pixel：将以像元为基础的图像镶嵌起来，Map—Mosaicking—Pixelbased。envidata\avmosaic\dv06_2.imgdv06_3.img操作步骤：^2选择Map—Mosaicking—Pixelbased，打开PixelbasedMosaic镶嵌界面。^2Import—ImportFiles—选择要镶嵌的影像，定义镶嵌范围（行列号、影像拉框选择、地理坐标，引入其他文件范围等），键入镶嵌输出的维数（或尺寸），选择镶嵌所用的波段—OK。^2选择或键入镶嵌尺寸的大小—OK。^2调整、确定影像的相对位置。^2分别选择影像—右键点击—EditEntry：对背景值、羽化值、匀光等进行设置。^2所有参数设置完成以后—File—SaveTemplate：生成临时效果影像，方便检查、调整。^2生成最终结果：File—Apply—确定输出象素大小、重采样方式、输出路径及文件名、背景值等—OK。2）基于地理坐标：Map—Mosaicking—Georeferenced。envidata\avmosaic\lch_olw.img；lch_02w.img操作步骤：^2选择Map—Mosaicking—Georeferenced，打开MapbasedMosaic镶嵌界面。^2Import—ImportFiles—选择要镶嵌的影像，定义镶嵌范围（行列号、影像拉框选择、地理坐标，引入其他文件范围等），键入镶嵌输出的维数（或尺寸），选择镶嵌所用的波段。^2分别选择影像—右键点击—EditEntry：对背景值、羽化值、匀光等进行设置。^2所有参数设置完成以后—File—SaveTemplate：生成临时效果影像，方便检查、调整。^2生成最终结果：File—Apply—确定输出象素大小、重采样方式、输出路径及文件名、背景值等—OK。3）折线、曲线等其他方式确定镶嵌影像范围：^2主影像窗口—Overly—Annotation打开添加注记对话框—选择折线、曲线等线形；^2画线：选择在Images、Scoll或Zoom窗口进行画线，定义范围，注重选择线的粗细和颜色（画线注重：一定要将线的头尾画到影像的边缘，左键开始，右键结束，右键确定）；^2添加标记：Annotation：Rectangle窗口—object—Symbol—选择符号—更改颜色和大小—将标记加在不要的影像一边；^2画好线、添加完标记以后，进行存储：Annotation：Rectangle窗口—File—SaveAnnotation—确定路径和文件名，保存为：*.ann文件；^2Import—ImportFiles—引入要镶嵌的影像；^2分别选择影像—右键点击—EditEntry，CutlineFeathering—SelectCutlineAnnotationFile引入*.ann文件，给定Cutline羽化值，OK。4）匀光：envidata\avmosic\mosaic1_equal.dat；mosaic_2.datEntry—ColorBalancing—分别确定基准影像和调整影像—Ok。Envi调用MODISReprojectionTool（MRT）对MODIS产品进行批处理拼接、重投影、裁切Science&TechnologyGeoscience,remotesensing,climate,Crop【转】Envi调用MODISReprojectionTool（MRT）对MODIS产品进行批处理拼接、重投影、裁切1熟悉MRT MODIS产品的类型不同，一景HDF格式的影像包含的波段也各不相同。MRT处理时需要选择处理波段，0示不作处理，1表示处理，首先要确定影像的波段数。1、拼接Mrtmosaic.exe程序用来拼接影像。调用方式为：mrtmosaic-i"g:\n%1.txt"-s"0000000000000000000000000001110000"-o"g:\temp.hdf"命令行帮助如下Usage:mrtmosaic-iinput_filenames_file-t-h-ooutput_filename                -sspectral_subset"b1b2...bN"                -gfilenameforthelogfile需要指定三个参数：-I 指定输入文件，可以采取两种方式1、1.hdf2.hdf3.hdf2、把影像的完整路径保存到txt文件中，作为input参数。如-IMOD092008001.txt。-s指定需要处理的波段，同样可以采用两种方式1、直接给出，如–s“000000000000”，注意英文引号；2、指定一个txt路径，让程序读取；-o指定输出路径，一般直接给出1、直接给出，如–og:\tmp.hdf。注意直接存为HDF格式，便于后续处理。2、重投影、裁切Resample.exe用来重投影、裁切是MRT程序的核心。调用方式为：命令行帮助如下：Usage:resample-pparameter_file[options]Optionsthatoverrideparameterfilespecifications:  -iinput_file_name  -ooutput_file_name  -rresampling_type[NNBICCNONE]  -tprojection_type[AEAERGEOHAMIGHISINLALCCMERCATMOLPSSINTMUTM]  -jprojection_parameter_list"p1p2...p15"  -sspectral_subset"b1b2...bN"  Ifusingthe-sswitch,theSDSsshouldberepresentedasan  arrayof0sand1s.A'1'specifiestoprocessthatSDS;  '0'specifiestoskipthatSDS.UnspecifiedSDSswillnotbeprocessed.  Ifthe-sswitchisnotspecified,thenallSDSswillbeprocessed.  -aspatial_subset_type[INPUT_LAT_LONGINPUT_LINE_SAMPLEOUTPUT_PROJ_COORDS]  -lspatial_subset"ULlatULlongLRlatLRlong"              -or- "ULlineULsampleLRlineLRsample(0-based)"              -or- "ULprojxULprojyLRprojxLRprojy"     NOTE:line/samplemustbespecifiedforthe            highestresolution ofallSDSsspecified            tobeprocessedintheproduct.  -uUTM_zone  -xpixel_size  -gfilenameforthelogfile可以只指定1个参数：用–p读入prm参数文件，进行处理。例如resample-p"g:\prrmMOD092008001.prm"，prm文件如下：INPUT_FILENAME=g:\tmp_%1.hdf输入文件SPECTRAL_SUBSET=(1111)前面mrtmosaic拼接结果有4个波段SPATIAL_SUBSET_TYPE=INPUT_LAT_LONG经纬度裁切方式SPATIAL_SUBSET_UL_CORNER=(33.0108.0)左上纬经度SPATIAL_SUBSET_LR_CORNER=(28.0117.0)右下纬经度OUTPUT_FILENAME=F:\MRT_out\myd%1.tif输出路径，不同波段自动区分RESAMPLING_TYPE=NEAREST_NEIGHBOR最近邻采样方法OUTPUT_PROJECTION_TYPE=UTM输出文件投影方式utmOUTPUT_PROJECTION_PARAMETERS=(投影参数 0.00.00.0 0.00.00.0 0.00.00.0 0.00.00.0 0.00.00.0)DATUM=WGS84大地水准面UTM_ZONE=50投影分带带号如果不会设置，可以先在图形界面里设置一次，把设置的结果保存下来。打开刚刚保存的参数文件去掉以#号打头的注释，文件显示如下：INPUT_FILENAME=G:\TmpMosaic.hdfSPECTRAL_SUBSET=(111111111111)SPATIAL_SUBSET_TYPE=INPUT_LAT_LONGSPATIAL_SUBSET_UL_CORNER=(32.0107.0)SPATIAL_SUBSET_LR_CORNER=(29.0116.0)OUTPUT_FILENAME=G:\mod092001001.tifRESAMPLING_TYPE=NEAREST_NEIGHBOROUTPUT_PROJECTION_TYPE=UTMOUTPUT_PROJECTION_PARAMETERS=( 0.00.00.0 0.00.00.0 0.00.00.0 0.00.00.0 0.00.00.0)DATUM=WGS84UTM_ZONE=50如果研究区的经纬度范围不同，可以修改SPATIAL_SUBSET_UL_CORNER=(32.0107.0)SPATIAL_SUBSET_LR_CORNER=(29.0116.0)两行如果需要不同的投影，如经纬度投影，可以修改OUTPUT_PROJECTION_TYPE=UTMOUTPUT_PROJECTION_PARAMETERS=( 0.00.00.0 0.00.00.0 0.00.00.0 0.00.00.0 0.00.00.0)DATUM=WGS84UTM_ZONE=50如果拼接时只选择了1个波段，可以修改SPECTRAL_SUBSET=(1)当然这句INPUT_FILENAME=G:\TmpMosaic.hdf批处理是必须改掉的了关于像素分辨率，一般留空，也就是默认不改变影像的像素大小。2关于DOS批处理裁切MOD产品单次处理Dirg:\mod11*A%1*.hdf/s/b>>g:\n%1.txtf:\mrt\modis\bin\mrtmosaic-i"g:\n%1.txt"-s"0000000000000000000000000001110000"-o"g:\temp_%1.hdf"f:\mrt\modis\bin\resample-p"g:\prm%1.prm"delg:\prm%1.prm/q/fdelg:\n%1.txt/q/fdelg:\temp_%1.hdf/q/f批量处理for/f%%iin(e:\lis.txt)do@单次处理文件名%%i如果双击运行，用%%i号表示变量。在命令行中%i表示变量。其中e:\lis.txt是年份及天数2008001200800920080172008025200803320080412008049G:\n%1.txt在运行时为g:\n2008001.txt，其中包括2008年第1天的多景MODIS产品影像名。Dirg:\mod11*A%1*.hdf/s/b>>g:\n%1.txt用来列出g盘（数据存储盘）以mod11打头的指定日期的的hdf格式文件，即MODIS影像。这里要说明一下，MODIS产品命名的规则MYD11A2.A2002185.h27v05.005.2007221162948.hdf从左至右，前7位表示产品类型MYD11指下午星(MYD)第11号产品（landsurfacetemperature），A2表示处理级别。A2002185表示过境时间是2002年第185天，h27v05是MODIS产品在全球的轨道行列号，中国一般是h27v05、h27v06、h28v05、h28v06，005表示HDF5文件格式，早期MODIS采用HDF4格式存储，所以有时也可见到004的产品。2007等一串数字表示影像处理的年月日时分秒，hdf表示存储格式为EOS-HDF。了解了MODIS产品命名规则之后，可以用DOSdir命令来挑选符合条件的MODIS产品路径了。>>g:\n%1.txt表示把结果转存为g:\n%1.txt%1表示接收的第一个参数，即20080013Envi生成时间列表2008001200800920080172008025这样的列表文件可以用在Envi中定义一个pro来生成：prodaylist,startdate,enddate,result_txt_path       ;exampledaylist,2008001,2008013,'example.txt'       ;       ;       a=(startdatemod1000)/8*8+(startdate/1000)*1000+1              b=((enddatemod1000)/8+1)*8+(enddate/1000)*1000+1       help,a,b              step=8       openw,lun,result_txt_path,/get_lun       res=a;            fori=a,b,stepdobegin                      printf,lun,string(res,format='(%"%7d")')                      res=res+step       endfor       free_lun,lun       print,'done';       close,/all       end在命令行窗口敲入daylist,2009001,2009033,'f:\g4.txt'打开f:\g4.txt如下：2009001200900920090172009025200903320090414Envi平台下整合前面工作留给读者思考 提示：spawn函数