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选择要镶嵌的影像，定义镶嵌范围（行列号、影像拉框选择、地理坐标，引入其他文件范围等），键入镶嵌输出的维数（或尺寸），选择镶嵌所
2、用的波段OK。2 选择或键入镶嵌尺寸的大小OK。2 调整、确定影像的相对位置。2 分别选择影像右键点击Edit Entry：对背景值、羽化值、匀光等进行设置。2 所有参数设置完成以后FileSave Template：生成临时效果影像，方便检查、调整。2 生成最终结果：FileApply确定输出象素大小、重采样方式、输出路径及文件名、背景值等OK。2）基于地理坐标：MapMosaickingGeoreferenced。envidataavmosaiclch_olw.img；lch_02w.img操作步骤：2 选择 MapMosaickingGeoreferenced，打开 Map based
3、Mosaic 镶嵌界面。2 ImportImport Files选择要镶嵌的影像，定义镶嵌范围（行列号、影像拉框选择、地理坐标，引入其他文件范围等），键入镶嵌输出的维数（或尺寸），选择镶嵌所用的波段。2 分别选择影像右键点击Edit Entry：对背景值、羽化值、匀光等进行设置。2 所有参数设置完成以后FileSave Template：生成临时效果影像，方便检查、调整。2 生成最终结果：FileApply确定输出象素大小、重采样方式、输出路径及文件名、背景值等OK。3）折线、曲线等其他方式确定镶嵌影像范围：2 主影像窗口OverlyAnnotation 打开添加注记对话框选择折线、曲线等线形
4、；2 画线：选择在 Images、Scoll 或 Zoom 窗口进行画线，定义范围，注重选择线的粗细和颜色（画线注重：一定要将线的头尾画到影像的边缘，左键开始，右键结束，右键确定）；2 添加标记：Annotation：Rectangle 窗口objectSymbol选择符号更改颜色和大小将标记加在不要的影像一边；2 画好线、添加完标记以后，进行存储：Annotation：Rectangle 窗口FileSave Annotation确定路径和文件名，保存为：*.ann 文件；2 ImportImport Files引入要镶嵌的影像；2 分别选择影像右键点击Edit Entry，Cutline
5、FeatheringSelect Cutline Annotation File 引入*.ann 文件，给定 Cutline 羽化值，OK。4）匀光：envidataavmosicmosaic1_equal.dat；mosaic_2.datEntryColor Balancing分别确定基准影像和调整影像Ok。Envi 调用 MODIS Reprojection Tool（MRT）对 MODIS 产品进行批处理拼接、重投影、裁切Science 0 specifies to skip that SDS. Unspecified SDSs will not be processed.If the
6、-s switch is not specified, then all SDSs 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
7、 for the highest resolutionof all SDSs specified 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 个波
8、段SPATIAL_SUBSET_TYPE = INPUT_LAT_LONG 经纬度裁切方式SPATIAL_SUBSET_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
9、 = (投影参数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.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
10、_CORNER = ( 32.0 107.0 )SPATIAL_SUBSET_LR_CORNER = ( 29.0 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如果研究区的经纬度范围不同，可以修改SPA
11、TIAL_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_F
12、ILENAME = G:TmpMosaic.hdf 批处理是必须改掉的了关于像素分辨率，一般留空，也就是默认不改变影像的像素大小。2 关于 DOS 批处理裁切 MOD 产品单次处理Dir g:mod11*A%1*.hdf /s /b g:n%1.txtf:mrtmodisbinmrtmosaic -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_%1.hdf“f:mrtmodisbinresample -p “g:prm%1.prm“del g
13、:prm%1.prm /q /fdel g:n%1.txt /q /fdel g:temp_%1.hdf /q /f批量处理for /f %i in (e:lis.txt) do 单次处理文件名 %i如果双击运行，用%i 号表示变量。在命令行中%i 表示变量。其中 e:lis.txt 是年份及天数2008001200800920080172008025200803320080412008049G:n%1.txt 在运行时为 g:n2008001.txt，其中包括 2008 年第 1 天的多景 MODIS 产品影像名。Dir g:mod11*A%1*.hdf /s /b g:n%1.txt用来列
14、出 g 盘（数据存储盘）以 mod11 打头的指定日期的的 hdf 格式文件，即 MODIS 影像。这里要说明一下，MODIS 产品命名的规则MYD11A2.A2002185.h27v05.005.2007221162948.hdf从左至右，前 7 位表示产品类型 MYD11 指下午星(MYD)第 11 号产品（ land surface temperature）， A2 表示处理级别。A2002185 表示过境时间是 2002 年第 185 天，h27v05 是 MODIS 产品在全球的轨道行列号，中国一般是 h27v05、h27v06 、h28v05、h28v06，005 表示 HDF5
15、文件格式，早期 MODIS 采用 HDF4 格式存储，所以有时也可见到 004 的产品。2007 等一串数字表示影像处理的年月日时分秒，hdf 表示存储格式为 EOS-HDF。了解了 MODIS 产品命名规则之后，可以用 DOS dir 命令来挑选符合条件的 MODIS 产品路径了。 g:n%1.txt 表示把结果转存为 g:n%1.txt%1 表示接收的第一个参数，即 20080013Envi 生成时间列表2008001200800920080172008025这样的列表文件可以用在 Envi 中定义一个 pro 来生成：pro daylist,startdate,enddate,resul
16、t_txt_path;example daylist,2008001,2008013,example.txt;a=(startdate mod 1000) /8 *8+(startdate /1000)*1000+1b=(enddate mod 1000) /8 +1)*8+(enddate /1000)*1000+1help,a,b step=8openw,lun,result_txt_path,/get_lunres=a; for i=a,b,step do beginprintf,lun,string(res,format=(%“%7d“)res=res+stependforfree_lun,lunprint,done;close,/allend在命令行窗口敲入daylist,2009001,2009033,f:g4.txt打开 f:g4.txt 如下：2009001200900920090172009025200903320090414Envi 平台下整合前面工作留给读者思考提示：spawn 函数
