InSAR处理doris使用详细步骤.doc

资源描述

InSAR处理doris使用详细步骤.doc

《InSAR处理doris使用详细步骤.doc》由会员分享，可在线阅读，更多相关《InSAR处理doris使用详细步骤.doc（28页珍藏版）》请在冰豆网上搜索。

InSAR处理doris使用详细步骤.doc

doris使用--DORIS干涉处理过程简介。

2009-06-1121:

1，Doris处理可以分步处理，也可以整体处理。

Doris处理要先生成模板文件，doris每步处理都对应一个模板文件，这些模板文件包含了doris每步处理所需要的参数，因此用户可以直接打开模板文件，用以编辑更改然后保存，也可以在linux终端中采用Vi命令，或nedit（相当于window下的写字板）来进行编辑保存。

2，查看版本信息，终端键入：

doris即可查看。

3，生成模板。

命令：

run–g

说明：

会在当前工作目录下生成模板文件，它包含三个文件夹目录：

Inputfiles，Outdata，Outinfo。

Outdata用于存放每步生成的结果数据文件。

Outinfo用于存放每步处理的log文件，记录中间处理过程。

Inputfiles中存放的就是模板文件。

它包含的文件：

对这些模板进行编辑时，文件中每行前面加c和空格，表示程序处理时不读取该行参数，否则都会执行到的，编辑参数选项时直接更改后面的的数值即可。

4，doris处理步骤总体介绍。

Doris的INSAR处理共分八步：

每步又可以分为许多小步，可以通过编辑来选择执行或不执行某些小步骤（小步骤在模板文件中的标识符为process）

第一步run–s1：

读入主影像的头文件及数据体信息以及轨道信息（可以使用粗轨），其对应的模板文件input.m_initial.，这步包含：

读取主影像头文件信息，读取主影像数据体信息，主影像过采样，读取主影像精轨文件。

（如下：

加c表示屏蔽该行，以下都同）

第二步run–s2：

读入辅影像的头文件及数据体信息以及轨道信息（可以使用粗轨），其对应的模板文件：

input.s_initial。

这步包含：

读取辅影像头文件信息，读取辅影像数据体信息，辅影像过采样，读取辅影像精轨文件。

第三步run–s3：

影像配准。

对应的模板文件：

input.coregistration.包含步骤有：

一级轨道级配准，二级像素级粗配准，三级子像素级精配准，主影像方位向滤波，辅影像方位向滤波，以及影像距离向滤波。

第四步run–s4：

影像重采样。

对应的模板文件：

input.resample。

包含步骤有：

估计配准多项式系数，影像重采样，以及距离向滤波。

第五步run–s5：

生成干涉图。

对应的模板文件：

input.products.。

包含步骤有：

干涉图生成，除平地效应，以及生成相干系数图。

其中去平地效应有两种方式：

一种是基于轨道去平地效应（如果在第一二步提供精轨，那么就是基于精轨去平地效应，否则是基于粗轨去平地效应）。

其中comprefpha步骤为用轨道来计算平地相位。

Subtrrefpha步骤为减去平地相位。

第二种是基于粗DEM去平地效应。

其中comprefdem为利用提供的粗DEM来反算求得平地相位，subtrrefdem为减去由粗DEM反算得到的平地相位。

第六步run–s6：

干涉图滤波及解缠。

对应的模板文件：

input.filter_unwrap.包含的步骤有：

干涉图滤波，相位解缠，以及差分干涉处理（针对Dinsar）

注意的是doris没有提供强有力的解缠程序，它只提供一种简单的枝切法解缠，这对于复杂地形下的干涉图解缠会存在很大误差，因此它一般调用比较有名的snaphu解缠算法来完成解缠过程。

Snaphu解缠也可以单独拿出来进行处理，

第七步run–s7：

DEM生成。

对应的模板文件：

input.s2h_geocode.。

它包含的步骤有：

相位高程转换，地理编码

第八步run–s9（整体的步骤包括第三至第七）：

doris分步处理前七步即完成，这里的第八步实际上是整体处理，它直接包含第三至第七步的全过程，对参数设置进行了简化，它对应的模板设置文件：

input.quicklook.。

该模板文件是将前三到七步模板文件中认为对insar处理比较关键的参数进行了一个集中设置。

因此若要进行整体处理时，必须也先要运行第一步和第二步，然后再运行第八步，即可。

doris使用--InSAR处理过程及参数说明

（1）

2009-06-1121:

1，命令介绍.

生成模板：

run–g

编辑模板：

编辑第一步模板：

run–e1

编辑第二步模板：

run–e2

编辑第三步模板：

run–e3

编辑第四步模板：

run–e4

编辑第五步模板：

run–e5

编辑第六步模板：

run–e6

编辑第七步模板：

run–e7

编辑第八步模板：

run–e9

（每行前加C表示屏蔽该行。

打开模板文件后，按“i“表示可以对进行插入操作，如果删去时，直接用delete或后退键。

按ESC键时则退出当前的编辑状态。

若要退出文件时，则先按ESC，然后：

q！

，或：

wq，其中：

q！

表示退出不保存修改，：

wq为保存修改。

）

运行处理：

运行第一步：

run–s1

运行第二步：

run–s2

运行第三步：

run–s3

运行第四步：

run–s4

运行第五步：

run–s5

运行第六步：

run–s6

运行第七步：

run–s7

运行第八步：

run–s9

2，参数解说及设置。

2.1，第一步：

读取主影像信息。

1）在文件input.m_initial中处理

（1）需要处理步骤。

必需：

m_readfiles（读取主影像头文件信息）.

必需：

m_crop（读取主影像数据体信息）.

不需：

m_ovs（主影像过采样步骤）：

作干涉处理时，一般不运行。

可选：

m_prorbits（读取主影像精轨文件）:

若要使用提供的精密轨道数据时，就执行该行，否则直接C掉。

（2）小步骤参数说明及设置。

############################################################

INSAR通用参数设置，每步都有的。

###########################################################

名称：

LOGFILE

用处：

用以记载DORIS处理过程中所有信息。

文件名：

doris.log.

名称：

M_RESFILE

用处：

用以记载主影像文件所有信息。

文件名：

master.res

名称：

S_RESFILE

用处：

用以记载辅影像文件所有信息。

文件名：

slave.res

名称：

I_RESFILE

用处：

用以记载主辅影像配准及以后步骤处理的所有信息。

文件名：

master_slave.res

############################################################

分步骤：

M_readfiles（读取主影像头文件信息）

###########################################################

名称：

M_IN_METHOD

用处：

使用的原始数据的传感器类型。

有（ERS（CEOS格式），ASAR（ENVISAT格式），JERS（CEOS格式），RSAT（RSAT格式），ATLANTIS（RSAT格式））。

说明：

必需

针对CEOS和RSAT格式的参数文件｛

名称：

M_IN_VOL SLC的VOLUME文件

名称：

M_IN_LEA SLC的LEADER文件

名称：

M_IN_DAT SLC的DAT文件

名称：

M_IN_NULL SLC的null文件，默认为dummy

｝

针对ENVISAT格式的参数文件｛

名称：

M_IN_DAT SLC的DAT文件（ASAR只有一个文件）

｝

############################################################

分步骤：

M_crops（读取主影像数据体信息）

###########################################################

名称：

M_CROP_IN

用处：

输入数据体的路径。

说明：

必需.dat文件。

名称：

M_CROP_OUT

用处：

输出的raw数据体路径。

说明：

必需

名称：

M_DBOW

用处：

要处理的数据区的大小，直接按影像大小来计数。

说明：

可选，C掉时表示默认为整块影像大小。

名称：

M_DBOW_GEO

用处：

要处理的数据区的大小，按按经纬度来计数。

说明：

一般不用这个，它与M_DBOW不能同选。

############################################################

分步骤：

M_ovr（主影像过采样，针对PSINSAR处理的）

###########################################################

############################################################

分步骤：

M_porbits（读取精密轨道信息）

###########################################################

名称：

M_ORBDIR

用处：

精密轨道文件路径。

说明：

针对精轨处理才给定路径。

2.2，第二步：

读辅影像信息。

2）在文件input.s_initial中处理

设置基本与第一步类同。

这里不作详细介绍。

2.3，第三步：

影像配准。

3）在文件input.coregistration中处理

（1）需要处理步骤

必需：

coarseorb一级轨道级配准。

必需：

coarsecorr二级像素级粗配准

必需：

m_filtazi主影像方位向滤波

必需：

s_filtazi辅影像方位向滤波

不需：

filtrange距离向滤波

必需：

fine 三级子像素精配准

（2）小步骤参数设置说明

############################################################

分步骤：

coarseorb（一级轨道级配准）

###########################################################

############################################################

分步骤：

coarsecorr（二级像素级粗配准）

###########################################################

CC_METHOD magfft|magspace