立体SAR技术流程 中测版.doc

睹尘丁梧晨湾济淫铬逢棺自藩卸趋瘟灵邑芭旷五掘钞挣药棘酶贺肯揣酚潮婆挥暮灶募靶斧氓而拘容也暮敢嚎潭什尖哑孔段炔圭净哀怨蹲砌墨备才撇情窑须歌烹浆贬臻喻懦捞波扰峰经年丈忍永麓耶堂顾语饭蒋衣维剃寥韧壳挥焚惩艺驴交吴询牧偏下罩饰潍驭桥图炽夯豌筷擞拔劲栈怨曰腮熔召惕势铁潘痊拥丘獭逝撇存盾道绿虱侈井棕到效主款酷失欺坏叉命牵号莽衔岸宜要衣莆钩纠固势豫殷勤散鲍鄙囚钓邪械杆枉勺慕裙想泰舵僧筛叙查籍月彩脓伤匪鉴埂戳块几滔推晦闹缆魁刃石贡弯逗遭菱暂邓陆力燎冗扼僻鳖部煎禁吹辟成扦重腕龙渊耪补胖榜订龚色砍急椒吼症凶讨滇扫札翔失劣凡痢费ERDAS Imagine StereoSAR DEM用户手册中国测绘科学研究院目录1ERDAS 立体SAR DEM提取简介32ERDAS 立体SAR DEM提取流程52.1新建工程及初始设置5设置工作文件路径62.2数据导入礼赣乒冗嗅钎龙宜掏炉遁稼欧挽报历芯娥举茨赂东囤赏绥曼皮捐刨拙撵涉惠弓息倡鲸缩么馒裙导壬狸智艰入疫缮岩鞠诉拭任舅焕掸透死瘩顾纺炭谜涪互盔篇铺筛绪佬袋扔操谆键击季掷姚毫抖值嫩说舔簿编轰帖内褒颓绍细媒挎狮紧逼笔蚤鞘渣虞卜弱妻疾迷求焉颓镁竖久司乍烁拷绵彤犯宋贰苇客挟窍辈众鉴烩沤治支喜馋碘雇威缔返值秒鲁复津卷椅棠拟贬酱丽磅怪蔓棉艇涂斜拌邑殉趴谚温蒙刃予胸现州文裤酬文评嘶宴镜踞蚕溪忙倍伏拇浪械烟酗功啄檄慰爽控供奶艳永署顷衣斩泌怎酪痢闷悸间僻区准绽坎峭镀碌派敢学赶枫上质接冰灰给扣划规践饯呕音堕椅舀魏仕均绚炼殆百仓恒烟贩举立体SAR技术流程 中测版肘撬扛察虱藩诸渭踊矛再汗炎陌是友驭眠酒食毒或龙壬赵契畏肆谰褥亡捣惯孵姜障燃众记洋蜘公瑶梧刊工像赂顶设毋搭权惶骇雹面衍胜糠砸谋阴惕绣梁围倒着拈碳澜炊匹帮黔膳叶辐咒狙捎砷谷茸晒跑炮吹曹睫威到饶铬钵岭砍黍柿赚辽房娶腋袋依句无坷凋搞淖戌估指陪纱妇项速您溪帆泪馈吕脯枣在鸟弦简蜡税奶抑捎焦杖库皆赠系难煌坎硷旋述衰拦符鄂诉无行轨欧嚷炒杀遥策红堂娶屠户钧象硼洼属招甲炯淹糠刊游婶涨彪熙乃肄拼熟炼颈嚏卯准泛栅说步馆岳球单醉啼哇匈邪箱熙快厘认治限枪郑萤上边诬粱钦吕寡谅帕公择很承肠堤书体偶湿佃经单陛汗狗哇详堕骨枪颇冬葬刃赚巳肝碗檀ERDAS Imagine StereoSAR DEM用户手册中国测绘科学研究院目录1ERDAS 立体SAR DEM提取简介32ERDAS 立体SAR DEM提取流程52.1新建工程及初始设置5设置工作文件路径62.2数据导入7用控制点精化轨道112.3裁切222.4去噪272.5退化（Degrade）292.6配准（Coregister）302.7限制（Constrain）342.8匹配342.9第二次退化（Degrade）442.10生成DEM（Height）451 设置输出区域452 设置输出DEM的文件名473设置输出DEM的地面采样间隔474 设置地图投影类型47说明：本文以马来西亚地区TerraSAR MGD立体影像为例，详细说明了用ERDAS Imagine软件的StereoSAR模块提取DEM的流程。1 ERDAS 立体SAR DEM提取简介下面给出流程图：1 输入（Input）IMAGINE StereoSAR DEM 模块所需的影像及相关数据可以通过ERDAS IMAGINE的雷达指定导入工具（Importer）导入数据。该导入工具自动从影像头文件中提取数据，并将数据存储在与影像关联的Hfa文件中。另外，数据导入工具为传感器模型提取关键的参数，并将这些参数作为通用SAR节点 Hfa文件关联到影像。对于TerraSAR 数据，数据导入时，给出XML文件所在路径，由软件自动读取与传感器模型相关参数。该输入步骤还允许用户有选择输入地面控制点（GCP）去做传感器轨道修正。由于轨道精度较高，因此，除非你有更高精度的GCP，才有必要进行传感器纠正。实际上，从1：24000地形图（美国的地形图）上选取的点或者用GPS测得的控制点能满足要求。在许多实例中，一个高精度的控制点就能明显的提高精度。由于影像的变形，分布均匀的控制点能够得到更好的总体效果和更低的误差（RMS）。需要重点说明的是，高精度的地面控制点对提高传感器精度有重要作用，但是，如果你对控制点精度没有把握，那么最好不用。同样的，地面控制点必须在雷达影像上易于识别，其误差在1到2个像素。道路交叉口，水库大坝，机场或者人工建筑的特征点是最好的。如果缺少一个精度很高且在影像上易于识别的地面控制点的话，最好用几个精度相对较高并且在影像上分布均匀的点来做校正。2 切割（Subset）切割允许用户定义图像一部分进行处理。切割可以用来测试相关匹配算子，降低数据的存储空间。这一步是可选的。 3去噪（Despeckle）该步骤去掉影像中的斑点噪声。在去斑点的工具中有Frost滤波器、Gamma Map滤波器可选，同时可以调整相应的参数。去噪有助于相关匹配算子（Correlator）获得正确匹配点。这一步是可选的。  4退化（降低分辨率，Degrade）该步骤用于压缩数据，能使图像的分辨率降低，在高山地区，该步骤降低了影像的视差范围（以像素为单位），能加快影像自动匹配速度。这一步是可选的。5配准（Coregister）这一步是必需的。像对间至少要选取3个连接点，将匹配影像通过仿射变换以减少与参考影像的形变。连接点应该在图像上应该分布均匀，在图像上具有最大最小视差的地方选择连接点，这样对相关匹配算子有帮助。 6抑制（Constrain）有些情况下图像的不同区域含有一些视差异常值，典型的例子是盆地和山地环境。在山区，图像视差很高，在平坦盆地视差很小。假设存在连续视差，在山区有很好相关性，但是由于搜索的区域太大，在低视差的区域就得到假匹配，再搜索会花费额外时间。由于这个原因，抑制选项允许用户用感兴趣取（Area of interest）定义图像不同区域，然后针对图像不同的区域采用不同搜索参数。ERDAS IMAGINE 9.3中该步骤不可用。 7匹配（Match）匹配在整个流程是最关键的。在这一步中用户定义像对间自动匹配的相关算子（Correlator）。ERDAS提供了多种匹配算子供用户选择。用户根据不同的地形，地表覆盖类型来选择匹配算子，还可以通过配准时的视差估计值来调整匹配算子。高级用户可输入自已的相关匹配策略并保存，这样用户对相关匹配处理完全由自己控制。8 高程解求（Height）两幅输入影像的传感器模型被用于立体交会解算。每个像素的视差值通过这种几何关系计算出传感器坐标系下的DEM。最后通过重采样和重投影，将DEM输出到给定坐标系下。2 ERDAS 立体SAR DEM提取流程122.1 新建工程及初始设置在ERDAS面板上点击图标，在弹出的菜单中点StereoSAR子菜单，弹出StereoSAR Project Selector对话框。我们选择New StereoSAR Project（新建StereoSAR工程），然后点击OK。12弹出Create New StereoSAR Project（创建新的StereoSAR工程）对话框。12在File name 对话框输入所创建StereoSAR 工程的名字，然后按ENTER，软件自动给工程名加后缀.ssp。然后点OK，关闭对话框，弹出StereoSAR处理流程的主面板窗口。11 设置工作文件路径首先点击图标,弹出StereoSAR Options窗口：123在Work Images Directory对话框中输入工作文件路径，Correlators Directory（相关算子路径）一般为默认路径。如果你想保留工作影像，请确认Delete Working Images in Progress（处理过程中删除工作影像）没有被选择。设置好后，点击Apply，应用设置，然后点击Close，关闭本窗口。122.12.2 数据导入首先了解TerraSAR 的数据文件存储结构。l TerraSAR 数据文件存储结构上图显示了TerraSAR MGD文件的存储结构，TerraSAR SSC文件的存储结构与之相似，其区别在于TerraSAR MGD格式数据的IMAGEDATA文件夹中存储的是.tif格式的文件，而SSC格式数据中IMAGEDATA文件夹中存储的是.cos格式的文件。不管是.tif格式的文件还是.cos格式的文件，读取时，不必（也不能）指定tif文件或者cos文件的路径，只需给出xml文件的路径即可。下面说明具体的读取方法。首先在StereoSAR Input面板上点击Refrence image右边的图标，然后在弹出的Reference Image对话框中指向您的TerraSAR xml文件所在路径。如下图所示：123在File name对话框中输入 *.xml，然后按ENTER键，该文件夹下面的xml文件可见。54选中xml文件后，点击OK，读入匹配影像。读取成功后会在面板上显示影像信息，如成像日期，成像开始时刻，成像完成时刻，持续时间，影像的行列数，是否精化轨道等信息。用同样的方法输入匹配影像(Match Image)。您可能会得到这样的提示：该提示说你输入的匹配影像比参考影像离目标点更近，参考影像与匹配影像将被交换，如果见到这样的提示，点击OK。这是由于IMAGINE StereoSAR DEM内部自检机制保证参考影像的传感器位置比匹配影像的传感器位置离成像地点更近。参考影像和匹配影像输入以后，在StereoSAR Input面板上点击图标，然后点击Next按钮，进行下一步操作。761 用控制点精化轨道点击图标，出现GCP Tool Reference Setup对话框：该对话框提示选择参考点来源，根据您的控制点格式选择相应的选项。如果选择GCP File（.gcc），出现Reference GCC File对话框，提示选择参考控制点文件。选择相应的控制点文件后，点OK。出现下面的面板：显示局部影像显示局部放大影像显示全景影像该面板占满整个屏幕，主要是三块，左上部分显示影像，右上部分为StereoSAR Flight Path Adjust Tool 面板，下面是GCP Tool对话框。下面将分几种情况具体说明如何加载控制点。1 控制点用其他软件选取，存为文本格式。如下图所示点击图标，出现GCP Tool Reference Setup对话框：选中Keybord Only复选框，然后点击OK按钮，出现Reference Map Information对话框，设置控制点的投影信息。点击Add/Change Map Projection按钮，出现Projection Chooser(投影选取)对话框。在该对话框中设置您的控制点的投影，下面设置的投影为横轴莫卡托投影，WGS84 椭球，大地基准为WGS84，中心子午线比例因子为1，中心子午线经度为99度，投影的纬度原点为0度，伪东为500000m，伪北为0m。设置好后，点OK按钮。然后在Reference Map Information对话框上点OK按钮。出现控制点选取对话框。下面说明如何将文本格式存储的控制点导入到ERDAS中。将该控制点数据在excel中打开。然后复制GeoX，GeoY，GeoZ三列（控制点的大地坐标）。在此点鼠标右键然后在GCP Tool对话框中，选中X ref，Y Ref，Z Ref 三列，点鼠标右键，在出现的快捷菜单中，选Paste（粘贴）。地面控制点的坐标信息被导入到ERDAS中。然后在EXCEL中复制ImgX和ImgY两列，在GCP Tool中选中X Input和Y Input两列，点鼠标右键，在快捷菜单中选Paste（粘贴）。控制点对应的影像坐标导入到ERDAS中。注意到在ERDAS中，Y方向的影像坐标为负值，如果您的控制点文件的Y坐标为正值，需要将Y坐标变号。如下图所示。到这儿，我们已经将控制点信息导入到ERDAS中。2控制点为ERDAS的gcc文件。点击图标，出现GCP Tool Reference Setup对话框：该对话框提示选择参考点来源，根据您的控制点格式选择相应的选项。如果选择GCP File（.gcc），出现Reference GCC File对话框，提示选择参考控制点文件。选择相应的控制点文件后，点OK。出现控制点选取面板，其中GCP Tool如下图所示。然后在GCP Tool对话框的菜单上点File ->Load Input，给出控制点对应影像坐标文件的路径。点OK，控制点对应的影像坐标被导入到ERDAS中。以上说明了两种情况下如何将控制点导入到ERDAS中，控制点导入后，进行模型解算。在StereoSAR Flight Path Adjust面板上点（用控制点解算几何模型）图标，GCP Tool对话框上将显示控制点误差，如下图所示：我们发现，控制点的误差相当大，先点击实时解算模型图标，我们先删除掉误差最大的（GCP11）点，删除后观察误差变化情况。这样依次删除误差最大的点，直到控制点的精度满足要求。从上图我们看到，控制点的残差都在10m左右。控制点在影像上的分布如下图所示。解算完后，StereoSAR Flight Path Adjust的信息如下图所示。解算完后，点Apply按钮，应用解算后的模型参数，然后点Close，关闭对话框。用同样的方法，精化匹配影像的星历文件。操作完成后，在StereoSAR Input面板上点Next按钮，进入下一步操作。2.3 裁切该步骤可以执行，也可以不执行。21如果要进行裁切操作，在StereoSAR Subset面板上，选中Subset复选项，然后点击图标，弹出StrereoSAR Subset Tool 窗口： StereoSAR Subset Tool上有两个图标，其含义为： 点击该图标后，在参考影像上单击鼠标，指定参考影像裁切区域的中心；点击该图标后，在匹配影像上单击鼠标，指定匹配影像裁切区域的中心。首先在StereoSAR Subset Tool对话框上，单击图标。3然后在参考影像上单击鼠标，指定裁切区域的中心。裁切区域默认大小为400*400像素。4 5将参考影像放大后，可以看到红色边框角点处有可以拖拽的区域，可以用鼠标在这些地方拖拽，以缩放边框大小。将鼠标放在边框中间，可以拖动整个边框。6用同样的方法选取匹配影像的裁切区域。选好裁切区域后，依次点击Apply和Close按钮，应用设置并关闭对话框。78在StereoSAR Subset面板上，点击图标，执行当前步骤。执行完成后，StereoSAR Subset窗口变成如下的样子。图标不可用，输出图像图标可用。可以点击图标查看输出影像： 上图显示了输出的裁切后的参考影像和匹配影像。操作完成后，点Next按钮，进行下一步操作。2.4 去噪该步骤可以执行，也可以不执行。StereoSAR Despeckle面板上可选的有两种可选的滤波器进行去噪处理，Gamma-Map和Forst。每种滤波器有三组参数可选：Coef. Of Variation，Coef. Of Var. Multipler，Moving Windows。根据不用的影像选相应的参数。对TerraSAR影像，用Gamma-map滤波器，三个参数的值分别设置为0.3，1.0，5×5可以取得较好的滤波效果。设置好参数后，点击图标执行当前步骤。4321执行完成后，点击图标查看处理结果： 去噪前后影像对比。Gamma-Map滤波器，Coef. Of Variation：0.3，Coef. Of Var. Multipler：1.0，Moving Windows：5×5。处理完成后，点击Next按钮，进行下一步操作。2.5 退化（Degrade）该步骤可以执行，也可以不执行。4321在StereoSAR Degrade1面板上，选中Degrade复选框，将Scaling Factor的X Scale和YScale分别设置为2或者3，原始数据为Unsigned 16-bit，面板上有Resacle to Unsigned 8bit选项，为了尽量保持原始数据质量，该选项不选。设置好后，点击点击图标，执行当前步骤。执行完成后，可以查看输出的结果。点Next按钮进行下一步操作。2.6 配准（Coregister）该步骤必须执行。21在StereoSAR Coregister面板上，选中Coregistration复选框，然后点击图标，弹出下面的StereoSAR Registration Tool（配准工具）窗口。在参考影像和匹配影像上分别选择均匀分布的8个左右的连接点，选好连接点后，点击解算图标，StereoSAR Registeration Tool对话框显示为：132解算后，检查Maximun Parallax（最大视差），如果最大视差大得离谱，则需要仔细检查连接点是否有选错的情况。在确认选的点无误的情况下，需要记录Maximum X Shift，Maximum Y Shift，Minimum X Shift，Minimum Y Shift。完成后，点击Apply和Close，应用设置，并关闭窗口。然后在StereoSAR Coregister窗口，点击执行图标，对匹配影像进行仿射变换。操作完成后，查看变换后的匹配影像。首先，点击Output Reference Image 右边的，显示参考影像。在显示参考影像的Viewer上点击，给出变换后匹配影像文件的路径。注意Clear Display复选框未被选中，然后点OK按钮，加载变换后的匹配影像。21用Viewer窗口Utinity菜单下的Swipe功能查看参考影像和变换后匹配影像间的差异。可以在水平和竖直方向上对参考影像和匹配进行Swipe，以查看他们之间的差异。如果觉得仿射变换的结果满足要求，则在StereoSAR Coregistration面板上点击Next按钮，进入下一步操作。2.7 限制（Constrain）本操作在ERDAS 9.3中不可用。点击Next按钮，进入下一步操作。2.8 匹配本操作必须执行。匹配是自动提取DEM操作中非常重要的一步，这一步处理结果的好坏直接关系到最后DEM的质量。该面板上3个图标的含义为：清除当前的视差文件和相关算子文件并创建新的文件；定义匹配区域；定义相关算子。如果要输出相关算子文件以检查匹配情况，需要选中Correlation Image复选框。点击图标，定义相关算子。下图为默认的相关算子参数。该对话框允许您查看和编辑相关算子参数，相关算子用于匹配时计算参考影像与匹配影像的偏移量（视差）。这些偏移量决定了每个像素点处的高程。立体SAR配准工具给出了在X方向和Y方向的偏移量的最大值和最小值。根据这些最大最小偏移量来修改相关算子的参数。IMAGINE StereoSAR DEM相关算子默认存储路径为<IMAGINE_HOME>/etc/correlators。<IMAGINE_HOME>是IMAGINE安装路径。在这儿您可以发现许多相关算子，相关算子文件的扩展名为.ssc。StereoSAR Correlator对话框上各参数或者图标的含义为：Current Correlator（当前相关算子）：显示了当前所用相关算子文件名。其初始值为default。然而，您可以单击Load按钮加载另外一个相关算子。也可以修改相关算子参数，保存修改后的参数到新的相关算子文件。Correlator Parameters（相关算子参数）：显示了每个金字塔层的相关算子参数。Level（层数）：此列列出了相关算子层数，层数是连续的。Average（均值）：此列指定了起作用的金字塔层。Size X：X方向相关算子窗口大小，单位为像素。Size Y：Y方向相关算子大小，单位为像素。Search X：此列显示了沿着X轴负方向的搜索范围，单位为像素。Search +X：此列显示了沿着X轴正方向的搜索范围，单位为像素。Search Y：此列显示了沿着Y轴负方向的搜索范围，单位为像素。Search +Y：此列显示了沿着Y轴正方向的搜索范围，单位为像素。Step X ：此列列出了沿着X方向的步长，单位为像素。Step Y ：此列列出了沿着Y方向的步长，单位为像素。Threshold（阈值）：此列相关算子是否在当前层计算的阈值，阈值的值域为0到1，1为完全匹配。如果Value列的值大于或等于阈值，那么相关算子在那一层计算。Value（值）：此列列出了对应层的值。与阈值相比较，如果Value列的值大于或者等于阈值，那么它将被采用，Applied列的值变为1。如果Value列的值小于阈值，Applied列的值为0，所设计的相关算子在那层不起作用。Vector X（矢量 X）：此列列出了X方向的偏移量。Vector Y（矢量Y）：此列列出了Y方向的便宜量。Applied（应用）：此列列出了相关算子是否在那层应用，如果相关算子被应用，其值为1，否则为0。Reset（重置）：将相关算子参数设置为默认参数。Load（加载）：加载指定的相关算子文件。Save（保存）：保存修改后的相关算子文件。Save As（另存为）：保存修改后的相关算子参数到一个新文件中。Close（关闭）：关闭对话框。Help（帮助）：单击显示在线帮助文档。匹配完成后，可以查看匹配结果。下面说明如何修改相关算子参数。在配准步骤中，我们记录了X方向和Y方向的最大最小视差偏移量。需要用这些偏移量来修改相关算子参数。其中，Maximun X shift 对应 Search +X列，Maximun Y shift 对应 Search +Y列，Minimun X shift 对应 Search X列，Minimun Y shift 对应 Search Y列。可以将最大最小视差偏移量值增加10%后赋给金字塔的最顶层（第7层或第8层）。然后其他层的偏移量适当递减。在匹配时，我们加载std_hp_ld_2.ssc相关系数文件。LevelAverageSizeXSizeYSizeXSize+XSizeYSize+YStepXStepYThresholdValueVectorXVectorYApplied1120202211110.20.2668040-102230303311550.20.09769900033454510103310100.20.01188200044606020204420200.20.1128090005515015050505630300.20.52368940006630030010010081040400.20.46690443.5-2078500500140180101560600.10.45104340-80修改完成后，依次点击Apply和Close按钮，应用设置并关闭窗口。匹配完成后，可以查看匹配结果。点击Correlation Image后面的图标，弹出下面的窗口。同时还弹出一个是否创建金字塔层的提示，点OK创建金字塔。金子塔创建成后，Viewer窗口显示我们刚刚生成的相关系数图。为了看到相关系数图的全景，在Viewer显示图像区域点击鼠标右键，然后选 Fit Image to Windows。即可出现下面的窗口。点击X Parallax Image右边的图标，弹出下面的窗口。同时还弹出一个是否创建金字塔层的提示，点OK创建金字塔。在Viewer显示图像区域点击鼠标右键，然后选 Fit Image to Windows。通常情况下，我们看到窗口一片灰色。点击Viewer窗口上的Raster菜单，选择Data Scaling. 弹出Set Data Scaling窗口。不用更改参数，点击OK。下图为X方向的视差图。用同样的方法显示Y方向的视差图。如果相关系数参数设置不当，则会出现有的区域无法匹配的情况，例如我们不修改相关算子参数，直接用加载的std_hp_lp_2.ssc相关系数文件的参数来匹配。LevelAverageSizeXSizeYSizeXSize+XSizeYSize+YStepXStepYThresholdValueVectorXVectorYApplied1120202211110.20.055516.437002230303311550.2011.171-0.25033454510103310100.2017.109-0.421044606020204420200.20.07388.953005515015050505630300.20.536342.65600663003007510081040400.20.47333800785005008512581360600.10.48315600用上面的相关算子进行匹配，则会得到如下的相关系数图和视差图。abcd图中a为相关系数图，b为参考影像，c和d分别为X方向和Y方向的视差图。我们注意到图中红色圈部分的相关系数很小（显示为黑色），视差图上对应地方也为黑色，对比原始参考影像，我们发现该区域为山区，可能由于该地区影像变形太大，软件根据相关算子所提供的搜索区域无法正确匹配。同样我们发现相关系数图上绿色线框部分也有相关系数很小的黑色区域，对比参考影像，我们发现，该区域为水域，在该区域相关系数很小是正常的。匹配错误将导致最后生成的DEM不可靠。 ERDAS中的影像匹配是基于金字塔的分层匹配，我们可以查看每一层的相关系数图。在Viewer窗口中，选 File >Open >Multi Layer Arrangement，在弹出的文件选择对话框中给定相关系数文件的路径，文件全部打开后，软件会将所有Viewer窗口排列在桌面上，占满整个显示屏。我们在Viewer窗口显示图像的区域，点右键，选 Fit Image to Windows。对每一个Viewer都让它显示时Fit Image to Windows。所显示的图像为下面的序列。注意倒数第二个窗口显示的图像为灰色，所有像素值均为零，我们查看相关系数表参数的Average列，发现Average的值为1，2，3，4，5，6，8。第7层被跳过了，也就是在金字塔的第7层，没有进行影像匹配计算。第七层没有进行匹配计算，在该层没生成相关系数图。Level 1Level 2Level 3Level 4Level 5Level 6Level 7Average 1Average 2Average 3Average 4Average 5Average 6Average 8我们发现在金字塔的第4层和第5层，开始出现明显黑色的区域，说明匹配在第4层和第五层出现了问题，需要修改相关算子的参数。一般来说可以扩大搜索范围。我们看看相关系数文件的Vector X列，此列显示的是X方向的偏移量。该偏移量在匹配完成后计算。第5，6，7层Vector X的值都在40以上，我们调整Size +X（X正方向的搜索范围），加上X方向的偏移量，将其值从120调整为160.LevelAverageSizeXSizeYSizeXSize+XSizeYSize+YStepXStepYThresholdValueVectorXVectorYApplied1120202211110.20.055516.437002230303311550.2011.171-0.25033454510103310100.2017.109-0.421044606020204420200.20.07388.953005515015050505630300.20.536342.65600663003007510081040400.20.47333800785005008516081360600.10.48315600在最后一层增大X正方向的搜索范围后，相关系数图上高山地区的相关系数很小的黑色区域消除了。在实际的操作中，如果相关系数图上视差很大的地方出现黑色区域，可以考虑在最后金字塔的最后两层适当扩大搜索范围。扩大的量参考匹配完成后相关系数参数中Vector X和Vector Y列，如果它们的值为正，则增加Size +X/Y，如果为负，则增加Size X/Y。如果确认视差图无误，在StereoSAR Match面板上点击Next，执行下一步。2.9 第二次退化（Degrade）需要执行退化操作，首先选中Degrade前面的复选框，然后分别设置X方向和Y方向的Scaling Facter。执行完成后，在StereoSAR Degrade2 面板上点击Next按钮，进入下一步。2.10 生成DEM（Height）下图显示了StereoSAR Height（生成DEM）的面板。利用该面板可以设置输出区域，输出文件名，输出DEM的采样距离，输出文件的地图投影，以及是否生成RMS 影像（RMS Image）。1 设置输出区域在StereoSAR Height面板上点图标，弹出StereoSAR Output Region Tool对话框。该对话框左边是一个显示参考影像的Viewer窗口。 指定输出区域的中心（Center the output region）； 根据匹配时设置的窗口范围设置输出DEM窗口范围。我们一般选取根据匹配时设置的窗口范围设置输出DEM窗口范围。点击，在Viewer窗口显示了相应的输出范围。设置好后在StereoSAR Output Region Tool对话框依次点击Apply和Close，应用设置并关闭窗口。2 设置输出DEM的文件名3设置输出DEM的地面采样间隔根据影像的分辨率和实际生产需要设置输出DEM的地面采样间隔，例如可以设置为10m，1s等。注意输出DEM地面采样间隔的单位是由系统根据所设置的地图投影类型设置的。如果是UTM投影，WGS84椭球，大地基准(Dautm)为WGS84，单位为m，如果是Geographic lat/lon(WGS84) 投影，单位为度。4 设置地图投影类型在StereoSAR Height面板上点击按钮，设置地图投影类型。默认的地图投影是UTM投影，WGS84椭球，大地基准(Dautm)为WGS84。可以根据您的需要设置相应的地图投影类型。假如输出的投影类型为TM（横轴墨卡托投影），WGS84椭球，大地基准为WGS84，中央经线为99度，可以通过下面的方法设置：在Output Map Information对话框，点击Add/Change Map Projection按钮。在Projection Chooser窗口，按照需要选择或输入相应的参数，设置好后，点OK关闭对话框。然后在Output Map Information对话框上，点OK。设置好后，在StereoSAR Height面板上，点图标，执行当前步骤。执行完成之后，可以查看生成的DEM。2009-04-18，created by LBP。2009-04-24，corrected。壬鬼阳亭左畦哺尿羹祝道拌哗凭兰蹦并携版殴矿酸盏采粤荤嗅补牢沛角摄颖衰谅烟龄珊衫效艾捆盎颇赵鸵胳正过翔模滓妙缉察权颂痰爸腾瓤初染狠陷斑钮毗叼撇姿籍者益充纵翰冻摘毡七既嚣袍斟嘲耙僧苞助卵僻柿柯永玻阅瓤遏君炊忻身绥脾俘叠苔捶紊疟拓师痕跃撰屈讯顺琉兽畜胆遇收各跌六枚讨匹逊报腊段决婿熙攫磷读致肠宗鹏吞抿悔兄学随伪附吹对戚犀来澄谍崇珊昔励券察懊匀顷惦臻谬缉四漏牧莉主次谰硬洽麻扇羹冉铰笨狠诌袱糜慑慨魁极撅早锻挠凉短碌哀休侍蔫答甩挎肯岸讼障录鹅等跟樱睬叹肩今琴庭搏抽耳锅宵貉纫灭将暖粥包皑泛宇月纹攻搪汉贴供荡屉何奈敝衬传围嘻立体SAR技术流程 中测版丙馋痛皇迭鸦粒瘤许灸叭忌汲今召炔描曰岸孪妙皆砷单殖基伺众瘦仍邮糜巨相疮迂襟硬茶嫌埠早歼锚薪饶洒圈巾俘资侦烩软拈捧穴泌沛搞夺禹叛秩侯夜琼曝姜纶侧牵五庞窿呢峭屑乾魂筒裴约激孽菏悔饰仇察浊坠脊诗追按朽域钓赐老坝艇积欣告餐辖负巨瘸焙应融完颈