2010-02-04 14:34
（包括：将以像元为基础的图像或应用地理坐标系的图像镶嵌起来两种方式） 1）基于pixel：将以像元为基础的图像镶嵌起来，Map—Mosaicking—Pixel based。
envidata\avmosaic\dv06_2.imgdv06_3.img
操作步骤：
一种关注
^2 选择Map—Mosaicking—Pixel based，打开Pixel based Mosaic镶嵌界面。 ^2 Import—Import Files—选择要镶嵌的影像，定义镶嵌范围（行列号、影像拉框选择、地理坐标，引入其他文件范围等），键入镶嵌输出的维数（或尺寸），选择镶嵌所用的波段—OK。 ^2 选择或键入镶嵌尺寸的大小—OK。
^2 调整、确定影像的相对位置。
^2 分别选择影像—右键点击—Edit Entry：对背景值、羽化值、匀光等进行设置。
^2 所有参数设置完成以后—File—Save Template：生成临时效果影像，方便检查、调整。
^2 生成最终结果：File—Apply—确定输出象素大小、重采样方式、输出路径及文件名、背景值等—OK。 2）基于地理坐标：Map—Mosaicking—Georeferenced。
envidata\avmosaic\lch_olw.img；lch_02w.img
操作步骤：
^2 选择Map—Mosaicking—Georeferenced，打开Map based Mosaic镶嵌界面。
^2 Import—Import Files—选择要镶嵌的影像，定义镶嵌范围（行列号、影像拉框选择、地理坐标，引入其他文件范围等），键入镶嵌输出的维数（或尺寸），选择镶嵌所用的波段。
^2 分别选择影像—右键点击—Edit Entry：对背景值、羽化值、匀光等进行设置。
万次火柴^2 所有参数设置完成以后—File—Save Template：生成临时效果影像，方便检查、调整。
^2 生成最终结果：File—Apply—确定输出象素大小、重采样方式、输出路径及文件名、背景值等—OK。
3）折线、曲线等其他方式确定镶嵌影像范围：
^2 主影像窗口—Overly—Annotation打开添加注记对话框—选择折线、曲线等线形；
^2 画线：选择在Images、Scoll或Zoom窗口进行画线，定义范围，注重选择线的粗细和颜（画线注重：一定要将线的头尾画到影像的边缘，左键开始，右键结束，右键确定）；
^2 添加标记：Annotation：Rectangle 窗口—object—Symbol—选择符号—更改颜和大小—将标记加在不要的影像一边；
^2 画好线、添加完标记以后，进行存储：Annotation：Rectangle 窗口—File—Save Annotation—确定路径和文件名，保存为：*.ann文件；
^2 Import—Import Files—引入要镶嵌的影像；
^2 分别选择影像—右键点击—Edit Entry，Cutline Feathering—Select Cutline Annotation File引入*.ann文件，给定Cutline羽化值，OK。
4）匀光：
envidata\avmosic\mosaic1_equal.dat；mosaic_2.dat
Entry—Color Balancing—分别确定基准影像和调整影像—Ok。
Envi调用MODIS Reprojection Tool（MRT）对MODIS产品进行批处理拼接、重投影、裁切
Science & Technology
Geoscience,remote sensing,climate,Crop
【转】Envi调用MODIS Reprojection Tool（MRT）对MODIS产品进行批处理拼接、重投影、裁切
1 熟悉MRT
MODIS产品的类型不同，一景HDF格式的影像包含的波段也各不相同。MRT处理时需要选择处理波段，0表示不作处理，1表示处理，首先要确定影像的波段数。
1、拼接
程序用来拼接影像。调用方式为：
mrtmosaic -i "g:\" -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_filename
-s spectral_subset "b1 b2 ... bN"
-g filename for the log file
需要指定三个参数：
-I 指定输入文件，可以采取两种方式
1、1.hdf 2.hdf 3.hdf
2、把影像的完整路径保存到txt文件中，作为input参数。如-。
-s 指定需要处理的波段，同样可以采用两种方式
1、直接给出，如 –s “0 0 0 0 0 0 0 0 0 0 0 0”，注意英文引号；
2、指定一个txt路径，让程序读取；
-o 指定输出路径，一般直接给出
1、直接给出，如 –o g:\tmp.hdf。注意直接存为HDF格式，便于后续处理。
2、重投影、裁切
用来重投影、裁切是MRT程序的核心。调用方式为：
命令行帮助如下：
Usage: resample -p parameter_file [options]
Options 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 projection_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 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 resolution
of 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个波段
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_out\myd%1.tif 输出路径，不同波段自动区分
支脚
RESAMPLING_TYPE = NEAREST_NEIGHBOR 最近邻采样方法
OUTPUT_PROJECTION_TYPE = UTM 输出文件投影方式utm
OUTPUT_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
0.0 0.0 0.0 )
DATUM = WGS84 自动化机械手臂大地水准面
UTM_ZONE = 50投影分带带号
如果不会设置，可以先在图形界面里设置一次，把设置的结果保存下来。
打开刚刚保存的参数文件
去掉以#号打头的注释，文件显示如下：
INPUT_FILENAME = G:\TmpMosaic.hdf
SPECTRAL_SUBSET = ( 1 1 1 1 1 1 1 1 1 1 1 1 )
SPATIAL_SUBSET_TYPE = INPUT_LAT_LONG
SPATIAL_SUBSET_UL_CORNER = ( 32.0 107.0 )
SPATIAL_SUBSET_LR_CORNER = ( 29.0 116.0 )
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议