《envi影像镶嵌,拼接线色彩平衡.doc》由会员分享，可在线阅读，更多相关《envi影像镶嵌,拼接线色彩平衡.doc（8页珍藏版）》请在金锄头文库上搜索。
1、envi影像镶嵌,拼接线/色彩平衡2010-02-04 14:34（包括：将以像元为基础的图像或应用地理坐标系的图像镶嵌起来两种方式）1）基于pixel：将以像元为基础的图像镶嵌起来，MapMosaickingPixel based。envidataavmosaicdv06_2.imgdv06_3.img操作步骤：2 选择MapMosaickingPixel based，打开Pixel based Mosaic镶嵌界面。2 ImportImport Files选择要镶嵌的影像，定义镶嵌范围（行列号、影像拉框选择、地理坐标，引入其他文件范围等），键入镶嵌输出的维数（或尺寸），选择镶嵌所用的波段O
2、K。2 选择或键入镶嵌尺寸的大小OK。2 调整、确定影像的相对位置。2 分别选择影像右键点击Edit Entry：对背景值、羽化值、匀光等进行设置。2 所有参数设置完成以后FileSave Template：生成临时效果影像，方便检查、调整。2 生成最终结果：FileApply确定输出象素大小、重采样方式、输出路径及文件名、背景值等OK。2）基于地理坐标：MapMosaickingGeoreferenced。envidataavmosaiclch_olw.img；lch_02w.img操作步骤：2 选择MapMosaickingGeoreferenced，打开Map based Mosaic镶
3、嵌界面。2 ImportImport Files选择要镶嵌的影像，定义镶嵌范围（行列号、影像拉框选择、地理坐标，引入其他文件范围等），键入镶嵌输出的维数（或尺寸），选择镶嵌所用的波段。2 分别选择影像右键点击Edit Entry：对背景值、羽化值、匀光等进行设置。2 所有参数设置完成以后FileSave Template：生成临时效果影像，方便检查、调整。2 生成最终结果：FileApply确定输出象素大小、重采样方式、输出路径及文件名、背景值等OK。3）折线、曲线等其他方式确定镶嵌影像范围：2 主影像窗口OverlyAnnotation打开添加注记对话框选择折线、曲线等线形；2 画线：选择在
4、Images、Scoll或Zoom窗口进行画线，定义范围，注重选择线的粗细和颜色（画线注重：一定要将线的头尾画到影像的边缘，左键开始，右键结束，右键确定）；2 添加标记：Annotation：Rectangle 窗口objectSymbol选择符号更改颜色和大小将标记加在不要的影像一边；2 画好线、添加完标记以后，进行存储：Annotation：Rectangle 窗口FileSave Annotation确定路径和文件名，保存为：*.ann文件；2 ImportImport Files引入要镶嵌的影像；2 分别选择影像右键点击Edit Entry，Cutline FeatheringSele
5、ct Cutline Annotation File引入*.ann文件，给定Cutline羽化值，OK。4）匀光：envidataavmosicmosaic1_equal.dat；mosaic_2.datEntryColor Balancing分别确定基准影像和调整影像Ok。Envi调用MODIS Reprojection Tool（MRT）对MODIS产品进行批处理拼接、重投影、裁切Science & TechnologyGeoscience,remote sensing,climate,Crop【转】Envi调用MODIS Reprojection Tool（MRT）对MODIS产品进行批
6、处理拼接、重投影、裁切 1 熟悉MRTMODIS产品的类型不同，一景HDF格式的影像包含的波段也各不相同。MRT处理时需要选择处理波段，0表示不作处理，1表示处理，首先要确定影像的波段数。1、拼接Mrtmosaic.exe程序用来拼接影像。调用方式为：mrtmosaic -i g:n%1.txt -s 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 -o g:temp.hdf命令行帮助如下Usage: mrtmosaic -i input_filenames_file -t -h -o output_fil
7、ename -s spectral_subset b1 b2 . bN -g filename for the log file需要指定三个参数：-I指定输入文件，可以采取两种方式1、1.hdf 2.hdf 3.hdf2、把影像的完整路径保存到txt文件中，作为input参数。如-I MOD092008001.txt。-s 指定需要处理的波段，同样可以采用两种方式1、直接给出，如 s “0 0 0 0 0 0 0 0 0 0 0 0”，注意英文引号；2、指定一个txt路径，让程序读取；-o 指定输出路径，一般直接给出1、直接给出，如 o g:tmp.hdf。注意直接存为HDF格式，便于后续处理
8、。2、重投影、裁切Resample.exe用来重投影、裁切是MRT程序的核心。调用方式为：命令行帮助如下：Usage: resample -p parameter_file optionsOptions that override parameter file specifications: -i input_file_name -o output_file_name -r resampling_type NN BI CC NONE -t projection_type AEA ER GEO HAM IGH ISIN LA LCC MERCAT MOL PS SIN TM UTM -j pro
9、jection_parameter_list p1 p2 . p15 -s spectral_subset b1 b2 . bN If using the -s switch, the SDSs should be represented as an array of 0s and 1s. A 1 specifies to process that SDS; 0 specifies to skip that SDS. Unspecified SDSs will not be processed. If the -s switch is not specified, then all SDSs
10、will be processed. -a spatial_subset_type INPUT_LAT_LONG INPUT_LINE_SAMPLE OUTPUT_PROJ_COORDS -l spatial_subset ULlat ULlong LRlat LRlong -or-ULline ULsample LRline LRsample (0-based) -or-ULprojx ULprojy LRprojx LRprojy NOTE: line/sample must be specified for the highest resolutionof all SDSs specif
11、ied to be processed in the product. -u UTM_zone -x pixel_size -g filename for the log file可以只指定1个参数：用 p 读入prm参数文件，进行处理。例如resample -p g:prrmMOD092008001.prm，prm文件如下：INPUT_FILENAME = g:tmp_%1.hdf输入文件SPECTRAL_SUBSET = ( 1 1 1 1) 前面mrtmosaic拼接结果有4个波段SPATIAL_SUBSET_TYPE = INPUT_LAT_LONG 经纬度裁切方式SPATIAL_SU
12、BSET_UL_CORNER = ( 33.0 108.0 ) 左上纬经度SPATIAL_SUBSET_LR_CORNER = ( 28.0 117.0 ) 右下纬经度OUTPUT_FILENAME = F:MRT_outmyd%1.tif 输出路径，不同波段自动区分RESAMPLING_TYPE = NEAREST_NEIGHBOR 最近邻采样方法OUTPUT_PROJECTION_TYPE = UTM 输出文件投影方式utmOUTPUT_PROJECTION_PARAMETERS = (投影参数0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
13、0.0 0.0 0.0 ) DATUM = WGS84 大地水准面UTM_ZONE = 50投影分带带号如果不会设置，可以先在图形界面里设置一次，把设置的结果保存下来。打开刚刚保存的参数文件去掉以#号打头的注释，文件显示如下：INPUT_FILENAME = G:TmpMosaic.hdfSPECTRAL_SUBSET = ( 1 1 1 1 1 1 1 1 1 1 1 1 )SPATIAL_SUBSET_TYPE = INPUT_LAT_LONGSPATIAL_SUBSET_UL_CORNER = ( 32.0 107.0 )SPATIAL_SUBSET_LR_CORNER = ( 29.0
14、 116.0 )OUTPUT_FILENAME = G:mod092001001.tifRESAMPLING_TYPE = NEAREST_NEIGHBOROUTPUT_PROJECTION_TYPE = UTMOUTPUT_PROJECTION_PARAMETERS = ( 0.0 0.0 0.00.0 0.0 0.00.0 0.0 0.00.0 0.0 0.00.0 0.0 0.0 )DATUM = WGS84UTM_ZONE = 50如果研究区的经纬度范围不同，可以修改SPATIAL_SUBSET_UL_CORNER = ( 32.0 107.0 )SPATIAL_SUBSET_LR_CORNER = ( 29.0 116.0 )两行如果需要不同的投影，如经纬度投影，可以修改OUTPUT_PROJECTION_TYPE = UTMOUTPUT_PROJECTION_PARAMETERS = ( 0.0 0.0 0.00.0 0.0 0.00.0 0.0 0.00.0 0.0 0.00.0 0.0 0.0 )DATUM = WGS84UTM_ZONE = 50如果拼接时只选择了1个波段，可以修改SPECTRAL_SUBSET = ( 1 )当然这句INPUT_FILENAME =
