A Method for Block Adjustment with Airborne InSAR by Aid of Three Dimensional Reconstruction Model
-
摘要: 在机载合成孔径雷达干涉测量（synthetic aperture radar interferometry，InSAR）中，为了获取高精度的数字高程模型（digital elevation model，DEM），研究了机载InSAR视向量正交分解法及三维重建数学模型，分析了中国的机载CASMSAR（synthetic aperture radar system of Chinese Academy of Surveying and Mapping）干涉系统误差来源及观测参数，建立了机载InSAR区域网平差模型。利用国产机载CASMSAR系统获取的X波段干涉数据进行试验，利用高精度的控制点数据进行平差解算，结果表明本文方法能够消除干涉模型间平面和高程差异，DEM成果满足1：1万山地测图精度要求，验证了该模型的正确性和有效性。Abstract: In order to improve the accuracy of DEM (digital elevation model) acquired by airborne InSAR (synthetic aperture radar interferometry), we researched the orthogonal decomposition of visual vector model and three dimensional reconstruction model for airborne InSAR. Error sources and adjustment parameters of domestic airborne InSAR system of CASMSAR (synthetic aperture radar system of Chinese Academy of Surveying and Mapping) were analyzed. A block adjustment method was proposed for airborne InSAR which could be used to derive DEM. An Experiment was carried out with X-band interferometric data acquired by CASMSAR. And we used ground control points to correct this error model. The results show that the proposed algorithm could be used to eliminate the differences in plane and height among multi-strip in the overlapping area. Finally, the results met the requirements of topographic mapping at a scale of 1:10 000, and proved the correctness and effectiveness of the new algorithm with CASMSAR data.
-
机载双天线合成孔径雷达干涉测量(synthetic aperture radar interferometry,InSAR)系统具有机动灵活、按需成像的特点，在地形测图中被广泛应用[1, 2]。在利用InSAR提取数字高程模型(digital elevation model,DEM)时，干涉系统参数误差通过相高转换、地理编码同时影响平面和高程，误差最终传递到DEM上。这些误差主要包括天线位置、速度矢量、基线矢量、斜距(绝对时间延迟)、干涉相位和多普勒频率等。
为了提高DEM精度，需利用地面控制点(ground control point,GCP)对这些参数进行平差改正。国内外学者提出了多种干涉参数改正方法，多采用平面电磁波模型，在研究过程中更偏重于理论研究，与实际应用情况有较大差距[3-5]。对多个单模型影像独立进行平差解算时，由于模型求解精度的不同，使得重叠区域内的点在平面和高程上差异较大；此外，单个模型的平差解算需要一定数量的控制点，在进行大面积测图应用时，不得不面临获取大量控制点的问题[6-9]。本文提出了机载InSAR区域网平差处理方法，该方法能减少控制点需求量、消除模型间平面和高程差异[9-12]，并获取高精度的DEM数据。
中国的机载CASMSAR(synthetic aperture radar system of Chinese Academy of Surveying and Mapping)系统，基线物理长度约2.2 m，在数据获取时间间隔较小的情况下，视为刚性基线，基线物理长度改正量应相同，本文在InSAR区域网平差模型参数设置中，各模型基线长度参量取统一值。由于飞行姿态不一致、雷达产生的随机误差不完全相同，不同的单元模型对应的其它系统参数不同且相互独立。根据敏感度方程[9]，利用连接点和少量控制点进行InSAR区域网平差处理，求得稳定可靠的干涉系统参数和连接点坐标，实现机载SAR数据高精度测图。
1 InSAR三维重建模型与区域网平差解算
1.1 地形三维重建模型
利用InSAR进行地形测绘时，通过距离方程、多普勒方程和干涉相位方程描述InSAR基本测量值与地面目标位置之间的关系，获取地面目标点的三维坐标，称为地形三维重建[5]。表示为：
$$\left\{ \begin{array}{*{35}{l}} {{r}_{A}}_{1}=\left| {{A}_{1}}-P \right| \\ {{f}_{dop}}=\frac{2v\cdot \left( {{A}_{1}}-P \right)}{\lambda \cdot {{r}_{A}}_{1}} \\ \varphi =\phi +2k\pi =\frac{2\pi Q({{r}_{A}}_{1}-{{r}_{A}}_{2})}{\lambda } \\ \end{array} \right.$$ (1) 式中，rA1为主天线到目标点的距离；rA2为副天线到目标点的距离；A1为主天线相位中心；P为目标点坐标；fdop为多普勒频率；v为速度矢量；λ为雷达波长；φ为绝对干涉相位；$\phi $为缠绕干涉相位，k为整周期数，Q为干涉模式。
视向量正交分解的基础是建立在载机移动坐标系上，以O-VNW表示载机移动坐标系，该坐标系以主天线所在相位中心为坐标原点，V轴定义为载机速度矢量方向，W轴定义为载机速度矢量与干涉基线矢量的叉积，N轴由V轴和W轴的叉积共同确定，构成右手坐标系。根据视向量正交分解方法[13]及InSAR三维重建几何关系，在地心空间直角坐标系下描述目标点：
$$\begin{align} & P={{A}_{1}}+{{r}_{1}}={{A}_{1}}+{{r}_{{{A}_{1}}}}\cdot R_{VNW}^{-1}\cdot {{{\hat{r}}}_{VNW}}= \\ & {{A}_{1}}+{{r}_{{{A}_{1}}}}\cdot \text{ }\left[ \begin{matrix} {\hat{v}} & \frac{\left( v\otimes b \right)\otimes v}{\left| \left( v\otimes b \right)\otimes v \right|} & \frac{\left( v\otimes b \right)}{\left| \left( v\otimes b \right) \right|} \\ \end{matrix} \right]\cdot \\ & \left[ \begin{array}{*{35}{l}} \frac{\lambda {{f}_{dop}}}{2\left| v \right|} \\ \frac{{{r}_{A}}_{1}}{2\left| {{b}_{pv}} \right|}\left[ 1+\frac{{{b}^{2}}}{{{r}_{A}}{{_{1}}^{2}}}-{{\left( 1-\frac{\lambda \varphi }{2\pi Q{{r}_{A}}_{1}} \right)}^{2}} \right]-\frac{{{b}_{v}}{{r}_{V}}}{\left| {{b}_{pv}} \right|} \\ \pm 1-{{r}_{V}}^{2}-{{r}_{N}}^{2} \\ \end{array} \right]\text{ } \\ \end{align}$$ (2) 式中，r1为视向量，${\hat{r}}$VNW为单位视向量;rA1=|r1|为视向量的模长，RVNW－1为载机移动坐标系向地心空间直角坐标系的过渡矩阵，b为基线矢量，b为基线矢量的模长，|bpv|为基线矢量垂直于速度方向的分量，bv为基线矢量平行于速度方向分量的模长，${\hat{r}}$VNW=(rV rN rW)T为待求解的单位视向量。
1.2 区域网平差解算
假设区域内有m个干涉单元模型，每个单元模型上分别有n1,n2,…,nm个点(因检查点不参与平差计算，故不统计其数量)，K=n1+n2+…+nm，其中量测点(多个单元模型同名点计为一个点)数量为n，控制点数量为p，连接点数量为q，观测方程个数为3K，未知数个数为(1+4m)+3q，有n＜K，n=p+q。其中，参与平差计算的控制点和连接点均选择相干性较高的点。
解算过程中，控制点有较高的定位精度，不需要进行改正，而连接点的三维重建坐标初值不准确，需要进行平差处理。
在地心空间直角坐标系下，根据式(2)，在目标点X、Y和Z三个方向对干涉参数bl、bα,i、r0,i、φi、fdop,i(i为干涉对编号)及连接点求偏导数，有：
$$\begin{align} & {{F}_{G(i,j)}}({{b}_{l}},{{b}_{\alpha }},{{r}_{0}},\varphi ,{{f}_{\text{dop}}})=0 \\ & {{F}_{Tie(i,k)}}({{b}_{l}},{{b}_{\alpha }},{{r}_{0}},\varphi ,{{f}_{\text{dop}}},{{P}_{X}},{{P}_{Y}},{{P}_{Z}})=0 \\ \end{align}$$ (3) 式中，G(i,j)表示第i个模型中第j个控制点;Tie(i,k)表示第i个模型中第k个连接点。
对全部控制点和连接点按照上式列出3个方程，对未知数求导并线性化得到误差方程：
$$V=A\Delta {{x}_{1}}+B\Delta {{x}_{2}}-L$$ (4) 式中，A为干涉系统参数的系数矩阵；B为控制点和连接点的系数矩阵；均由式(3)计算得到，形式如下：
$$\eqalign{ & V = A\Delta {x_1} + B\Delta {x_2} - L \cr & \mathop V\limits_{3K \times 1} = {\left[ {\matrix{ {{v_{x,G(i,j)}}} & {{v_{y,G(i,j)}}} & {{v_{z,G(i,j)}}} & \cdots & {{v_{x,T(i,k)}}} & {{v_{y,T(i,k)}}} & {{v_{z,T(i,k)}}} & \cdots \cr } } \right]^{\rm{T}}}; \cr & \mathop A\limits_{3K \times \left( {1 + 4m} \right)} = \cr & \left[ {\matrix{ {{{\partial {F_{G(1,j)}}} \over {\partial {b_l}}}} \hfill & {{{\partial {F_{G(1,j)}}} \over {\partial {b_{\alpha l}}}}} \hfill & {{{\partial {F_{G(1,j)}}} \over {\partial {r_{0l}}}}} \hfill & {{{\partial {F_{G(1,j)}}} \over {\partial {\varphi _1}}}} \hfill & {{{\partial {F_{G(1,j)}}} \over {\partial {F_{d1}}}}} \hfill & 0 \hfill & 0 \hfill & 0 \hfill & 0 \hfill & {} \hfill \cr {{{\partial {F_{Tie(1,k)}}} \over {\partial {b_l}}}} \hfill & {{{\partial {F_{Tie(1,k)}}} \over {\partial {b_{\alpha l}}}}} \hfill & {{{\partial {F_{G(1,j)}}} \over {\partial {r_{0l}}}}} \hfill & {{{\partial {F_{Tie(1,k)}}} \over {\partial {\varphi _1}}}} \hfill & {{{\partial {F_{Tie(1,k)}}} \over {\partial {F_{d1}}}}} \hfill & 0 \hfill & 0 \hfill & 0 \hfill & 0 \hfill & {} \hfill \cr {} \hfill & {} \hfill & \vdots \hfill & {} \hfill & {} \hfill & \vdots \hfill & {} \hfill & {} \hfill & {} \hfill & \cdots \hfill \cr {{{\partial {F_{G(2,j)}}} \over {\partial {b_l}}}} \hfill & 0 \hfill & 0 \hfill & 0 \hfill & 0 \hfill & {{{\partial {F_{G(2,j)}}} \over {\partial {b_{\alpha 2}}}}} \hfill & {{{\partial {F_{G(2,j)}}} \over {\partial {r_{02}}}}} \hfill & {{{\partial {F_{G(2,j)}}} \over {\partial {\varphi _2}}}} \hfill & {{{\partial {F_{G(2,j)}}} \over {\partial {F_{d2}}}}} \hfill & {} \hfill \cr \matrix{ {{\partial {F_{Tie(2,k)}}} \over {\partial {b_l}}} \hfill \cr \hfill \cr} \hfill & \matrix{ 0 \hfill \cr \hfill \cr} \hfill & \matrix{ 0 \hfill \cr \vdots \hfill \cr} \hfill & \matrix{ 0 \hfill \cr \hfill \cr} \hfill & \matrix{ 0 \hfill \cr \hfill \cr} \hfill & \matrix{ {{\partial {F_{Tie(2,k)}}} \over {\partial {b_{\alpha 2}}}} \hfill \cr \vdots \hfill \cr} \hfill & \matrix{ {{\partial {F_{Tie(2,k)}}} \over {\partial {r_{02}}}} \hfill \cr \hfill \cr} \hfill & \matrix{ {{\partial {F_{Tie(2,k)}}} \over {\partial {\varphi _2}}} \hfill \cr \hfill \cr} \hfill & \matrix{ {{\partial {F_{Tie(2,k)}}} \over {\partial {F_{d2}}}} \hfill \cr \vdots \hfill \cr} \hfill & \matrix{ \hfill \cr \ddots \hfill \cr} \hfill \cr } } \right] \cr & {{\partial {F_{G(1,j)}}} \over {\partial {b_l}}}{\rm{ = }}\left[ {\matrix{ {{{\partial {F_{xG(1,j)}}} \over {\partial {b_l}}}} \hfill \cr {{{\partial {F_{yG(1,j)}}} \over {\partial {b_l}}}} \hfill \cr {{{\partial {F_{zG(1,j)}}} \over {\partial {b_l}}}} \hfill \cr } } \right];\mathop B\limits_{3K \times 3q} = \cr & \left[ {\matrix{ {{B_{G(1,j)}}} \hfill & 0 \hfill & 0 \hfill & 0 \hfill & {} \hfill \cr 0 \hfill & {{B_{G(1,j + 1)}}} \hfill & 0 \hfill & 0 \hfill & \cdots \hfill \cr {{B_{Tie(1,k)}}} \hfill & 0 \hfill & 0 \hfill & 0 \hfill & {} \hfill \cr 0 \hfill & {{B_{Tie(1,k + 1)}}} \hfill & 0 \hfill & 0 \hfill & {} \hfill \cr {} \hfill & \vdots \hfill & {} \hfill & \vdots \hfill & \cdots \hfill \cr 0 \hfill & 0 \hfill & {{B_{Tie(2,k)}}} \hfill & 0 \hfill & {} \hfill \cr {{B_{G(2,j)}}} \hfill & 0 \hfill & 0 \hfill & 0 \hfill & 0 \hfill \cr 0 \hfill & 0 \hfill & 0 \hfill & {{B_{Tie(2,k + 1)}}} \hfill & {} \hfill \cr 0 \hfill & {{B_{G(1,j + 1)}}} \hfill & 0 \hfill & 0 \hfill & \cdots \hfill \cr {} \hfill & \vdots \hfill & {} \hfill & \vdots \hfill & \ddots \hfill \cr } } \right]; \cr & \left\{ {\matrix{ {{B_{G(i,j)}} = \left[ {\matrix{ 0 \hfill & 0 \hfill & 0 \hfill \cr 0 \hfill & 0 \hfill & 0 \hfill \cr 0 \hfill & 0 \hfill & 0 \hfill \cr } } \right]} \hfill \cr {{B_{Tie(i,k)}} = \left[ {\matrix{ { - 1} \hfill & 0 \hfill & 0 \hfill \cr 0 \hfill & { - 1} \hfill & 0 \hfill \cr 0 \hfill & 0 \hfill & { - 1} \hfill \cr } } \right]} \hfill \cr } } \right. \cr & \mathop L\limits_{3K \times 1} = {\left[ {\matrix{ {{l_{x,G(i,j)}}} \hfill & {{l_{y,G(i,j)}}} \hfill & {{l_{z,G(i,j)}}} \hfill & \cdots \hfill & {{l_{x,Tie(i,k)}}} \hfill & {{l_{y,Tie(i,k)}}} \hfill & {{l_{z,Tie(i,k)}}} \hfill & \cdots \hfill \cr } } \right]^{\rm{T}}}; \cr & \mathop {\Delta {x_1}}\limits_{\left( {1 + 4m} \right) \times 1} = {\left[ {\matrix{ {\Delta {b_l}} \hfill & {\Delta {b_{\alpha ,1}}} \hfill & {\Delta {r_{0,1}}} \hfill & {\Delta {\varphi _1}} \hfill & {\Delta {F_{dop,1}}} \hfill & {\Delta {b_{\alpha ,2}}} \hfill & {\Delta {r_{0,2}}} \hfill & {\Delta {\varphi _2}} \hfill & {\Delta {F_{dop,2}}} \hfill & \cdots \hfill \cr } } \right]^{\rm{T}}}; \cr & \mathop {\Delta {x_2}}\limits_{3q \times 1} = {\left[ {\matrix{ {\Delta {P_{x,k}}} \hfill & {\Delta {P_{y,k}}} \hfill & {\Delta {P_{z,k}}} \hfill & {\Delta {P_{x,k + 1}}} \hfill & {\Delta {P_{y,k + 1}}} \hfill & {\Delta {P_{z,k + 1}}} \hfill & \cdots \hfill \cr } } \right]^{\rm{T}}} \cr} $$ 本文的区域网平差采用迭代计算、逐渐趋近方法[14]，根据最小二乘原理进行区域网平差计算。
2 实验与分析
选取四川省北部若尔盖地区机载InSAR数据，进行区域网平差试验。该区地处青藏高原东部边缘地带，高度2 900～4 200 m，高差1 000 m以上，属于山地地形。飞机飞行方向为由东向西、由西向东两个方向，按照飞行方向、航带号对干涉数据进行编号，试验使用6个航带号，分别为DX_65、DX_66、DX_67、XD_65、XD_66、XD_67，每个航带内选取3个单元模型，编号为50、51和52，总计18对干涉数据，覆盖面积近90 km2，如图 1所示。
试验区范围内控制点数量为7个，检查点数量为19个。本文对比单模型平差方法和InSAR区域网平差方法，从以下三个方面进行评价。
2.1 平差模型精度分析
利用InSAR区域网平差模型改正的干涉系统参数，通过三维重建模型解算控制点的地理坐标，并和控制点已知地理坐标进行比较，统计其中误差，作为评价区域网平差模型是否有效的标准。
试验中，InSAR区域网平差模型使用了7个控制点，两种方法所用控制点残差对比结果(为表述简便，均取绝对值参与统计)见表 1。
表 1 区域网平差模型和单模型平差所用控制点残差统计Table 1. Statistic of Error Using GCPs by This Paper and Single Model控制点残差 统计项 本文区域网平差 单模型平差 最大残差 1.63 4.21 X方向/m 最小残差 0.373 0.65 中误差 0.97 2.04 最大残差 1.33 3.89 Y方向/m 最小残差 0.07 0.54 中误差 0.74 1.78 最大残差 2.91 4.62 Z方向/m 最小残差 0.29 0.68 中误差 1.54 2.93 (1) 采用单模型平差解算，每个单模型至少需要3个控制点，18个模型则至少需要54个控制点，实际使用了67个控制点，而InSAR区域网平差模型使用了7个控制点，有效减少了控制点需求，降低了外业测量成本。
(2) 表 1中，本文InSAR区域网平差模型解算后，控制点X方向的中误差为0.97 m，Y方向的中误差为0.74 m，Z方向的中误差为1.54 m，单模型平差解算后，控制点X、Y和Z三个方向上的残差或中误差较本文区域网平差模型中控制点的残差或中误差明显增大，说明在稀少控制条件下，InSAR区域网平差模型具有明显的优势。区域网平差处理增强了各模型间重叠区域控制点的约束关系，降低了系统噪声和敏感度矩阵病态性的影响，其模型稳健性更强。单模型平差处理无法有效避免这些问题，得到的结果可信度较低。
2.2 干涉系统参数分析
稳定的基线长度、精确的基线倾角、准确的干涉相位、多普勒中心频率和绝对时间延迟等参数是获取高精度DEM的关键因素，分析干涉系统参数的变化规律有利于InSAR区域网平差模型的改进。InSAR区域网平差方法得到的干涉系统参数处理结果见表 2，以其中2个航带(由东向西方向的DX_66和DX_67)的6个单元模型为例，分析两种方法得到的干涉系统参数。
(1) 机载CASMSAR系统，基线物理长度约2.2 m，在时间间隔较小的情况下，视为刚性基线，其基线物理长度改正量应相同。InSAR区域网平差方法得到的基线长度参量为2.192 1 m，而利用单模型平差方法得到的基线长度最大为2.224 0 m，最小为2.155 6 m，互差达到0.068 4 m，与实际测量值比较，本文的区域网平差方法得到的基线长度值接近实测值，更符合物理情况，如表 2所示。
表 2 两种方法得到的干涉系统参数对比Table 2. Interferometric Calibration Parameters Comparison by This Paper and Single Model航带号 统计项 干涉对编号 实际测量值 50 51 52 区域网平差 单模型平差 区域网平差 单模型平差 区域网平差 单模型平差 基线长度/m 2.192 1 2.171 9 2.192 1 2.173 6 2.192 1 2.155 6 2.189 5 基线倾角/(°) 0.342 5 0.748 8 0.481 7 0.811 6 0.540 2 1.251 3 - DX_66 初始斜距改正数/m 1.268 7 0.685 3 -0.172 7 0.470 0 -0.371 9 0.603 6 - 相位改正数/rad -6.465 4 -11.545 7 5.714 9 1.250 5 -0.941 3 -10.033 0 - 多普勒改正数/Hz -0.644 9 -1.174 5 -0.857 6 -0.601 2 -1.278 5 -1.111 5 - 基线长度/m 2.192 1 2.159 9 2.192 1 2.224 0 2.192 1 2.195 5 2.189 5 基线倾角/(°) 0.519 2 1.208 9 0.488 6 -0.103 9 0.5607 0.489 0 - DX_67 初始斜距改正数/m 1.067 7 0.977 7 0.368 4 0.933 6 0.693 8 0.923 2 - 相位改正数/rad -1.633 9 -9.920 2 -1.445 4 6.410 6 4.619 8 5.491 6 - 多普勒改正数/Hz -0.062 8 -0.337 7 -0.623 3 -0.415 7 -0.678 6 -0.420 0 - (2) 单模型平差求解的干涉系统参数，受系统噪声、敏感度矩阵病态性以及控制点精度的影响，其结果只适用于控制点周边区域，无法满足整个测绘带的要求，会导致重叠区域平面和高程信息不一致，引起更大的测量误差。
2.3 DEM精度分析
利用区域范围内19个检查点，分别统计本文方法和单模型方法获取的DEM高程误差，统计结果如表 3所示，本文方法计算得到的DEM成果如图 2所示。
表 3 DEM检查点高程误差统计Table 3. Statistics of Error Using Inspect Points统计项 区域网平差方法/m 单模型平差/m 最大误差 2.31 4.54 最小误差 0.50 0.59 中误差 1.33 3.02 从表 3可以发现：① 本文方法获取的DEM检查点高程中误差为1.33m，而单模型方法得到的中误差为3.02 m，试验结果表明本文方法能够有效减少模型间的高程差异，改善DEM精度；② 在InSAR区域网平差模型参数设置中，对基线长度参量取统一值符合理论和实际情况；③ 由检查点精度可知，满足山地地区1∶1万DEM定位精度要求。
3 结 语
本文提出了一种利用三维重建模型进行视向量正交分解，随后建立机载InSAR区域网平差模型，对干涉系统参数进行模型优化。以四川地区机载X波段干涉数据进行试验，结果表明本文提出的机载InSAR区域网平差方法是有效的，获取的DEM数据满足国家1∶1万地形图测图精度要求。在后续的平差试验中，需要更大范围的数据验证本文提出方法的有效性和实用性。此外，在平差模型的构建中，没有考虑控制点本身的精度，需要剔除误差较大的控制点，进一步提高干涉DEM的精度。
-
表 1 区域网平差模型和单模型平差所用控制点残差统计
Table 1 Statistic of Error Using GCPs by This Paper and Single Model
控制点残差 统计项 本文区域网平差 单模型平差 最大残差 1.63 4.21 X方向/m 最小残差 0.373 0.65 中误差 0.97 2.04 最大残差 1.33 3.89 Y方向/m 最小残差 0.07 0.54 中误差 0.74 1.78 最大残差 2.91 4.62 Z方向/m 最小残差 0.29 0.68 中误差 1.54 2.93 
下载: 导出CSV
表 2 两种方法得到的干涉系统参数对比
Table 2 Interferometric Calibration Parameters Comparison by This Paper and Single Model
航带号 统计项 干涉对编号 实际测量值 50 51 52 区域网平差 单模型平差 区域网平差 单模型平差 区域网平差 单模型平差 基线长度/m 2.192 1 2.171 9 2.192 1 2.173 6 2.192 1 2.155 6 2.189 5 基线倾角/(°) 0.342 5 0.748 8 0.481 7 0.811 6 0.540 2 1.251 3 - DX_66 初始斜距改正数/m 1.268 7 0.685 3 -0.172 7 0.470 0 -0.371 9 0.603 6 - 相位改正数/rad -6.465 4 -11.545 7 5.714 9 1.250 5 -0.941 3 -10.033 0 - 多普勒改正数/Hz -0.644 9 -1.174 5 -0.857 6 -0.601 2 -1.278 5 -1.111 5 - 基线长度/m 2.192 1 2.159 9 2.192 1 2.224 0 2.192 1 2.195 5 2.189 5 基线倾角/(°) 0.519 2 1.208 9 0.488 6 -0.103 9 0.5607 0.489 0 - DX_67 初始斜距改正数/m 1.067 7 0.977 7 0.368 4 0.933 6 0.693 8 0.923 2 - 相位改正数/rad -1.633 9 -9.920 2 -1.445 4 6.410 6 4.619 8 5.491 6 - 多普勒改正数/Hz -0.062 8 -0.337 7 -0.623 3 -0.415 7 -0.678 6 -0.420 0 -
下载: 导出CSV
表 3 DEM检查点高程误差统计
Table 3 Statistics of Error Using Inspect Points
统计项 区域网平差方法/m 单模型平差/m 最大误差 2.31 4.54 最小误差 0.50 0.59 中误差 1.33 3.02
下载: 导出CSV
-
[1] 黄国满,张继贤,赵争. 机载干涉SAR测绘制图应用系统研究[J].测绘学报,2008,37(3):277-279 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB200803003.htm Huang Guoman, Zhang Jixian, Zhao Zheng. Research on Airborne SAR Interferometry Mapping System[J]. Acta Geodaetica et Cartographica Sinica, 2008,37(3):277-279 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB200803003.htm
[2] 刘国祥,丁晓利,李志林,等.使用InSAR建立DEM的试验研究[J].测绘学报,2001,30(4):336-342 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB200104014.htm Liu Guoxiang, Ding Xiaoli, Li Zhilin, et al. Experimental Investigation on DEM Generation Through InSAR[J]. Acta Geodaetica et Cartographica Sinica, 2001,30(4):336-342 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB200104014.htm
[3] Mallorqui J J, Bara M, Broquetas A. Calibration Requirements for Airborne SAR Interferometry[C]. SPIE on SAR Image Analysis, Modeling, and TechniquesⅢ, Barcelona, Spain, 2000
[4] Dall J. Calibration of Airborne Signal-pass Interferometry Data[C]. Proceedings of 2002 EUSAR, Cologne, Germany, 2002
[5] 王彦平.机载干涉SAR定标模型与算法研究[D].北京:中国科学院电子学研究所,2004 Wang Yanping. Studies on Calibration Model and Algorithm for Airborne Interferometric SAR[D]. Beijing:Institute of Electronics, Chinese Academy of Sciences, 2004
[6] 马婧,尤红建,胡东辉. 机载干涉雷达数据联合定位及拼接技术研究[J]. 武汉大学学报·信息科学版,2012,37(4):410-414 http://ch.whu.edu.cn/CN/abstract/abstract182.shtml Ma Jing, You Hongjian, Hu Donghui. A New Method of Block Adjustment of Airborne SAR Image with Few GCPs[J]. Geomatics and Information Science of Wuhan University,2012,37(4):410-414 http://ch.whu.edu.cn/CN/abstract/abstract182.shtml
[7] 云烨,曾琪明,焦健,等.基于参考DEM的机载InSAR定标方法[J]. 测绘学报,2014, 43(1):74-82 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201401013.htm Yun Ye, Zeng Qiming, Jiao Jian,et al. Calibration of Airborne Interferometric SAR Data Based on Reference DEM[J].Acta Geodaetica et Cartographica Sinica, 2014,43(1):74-82 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201401013.htm
[8] 陈刚,汤晓涛,钱方明. 星载InSAR立体基线定标方法[J]. 武汉大学学报·信息科学版,2014,39(1):37-41 http://ch.whu.edu.cn/CN/abstract/abstract2860.shtml Chen Gang, Tang Xiaotao, Qian Fangming. A Three-dimensional Baseline Calibration Method of Spaceborne InSAR[J].Geomatics and Information Science of Wuhan University,2014, 39(1):37-41 http://ch.whu.edu.cn/CN/abstract/abstract2860.shtml
[9] 王萌萌,黄国满,花奋奋,等. 机载双天线InSAR联合定标算法[J]. 测绘学报,2014, 43(12):1259-1265 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201412009.htm Wang Mengmeng, Huang Guoman, Hua Fenfen, et al. Joint Calibration Method of Airborne Dual-antenna Interferometric SAR[J].Acta Geodaetica et Cartographica Sinica, 2014,43(12):1259-1265 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201412009.htm
[10] 黄国满,杨书成,王宁娜,等.稀少控制下机载InSAR区域网联合地理编码方法[J]. 测绘学报,2013,42(3):397-403 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201303014.htm Huang Guoman, Yang Shucheng, Wang Ningna, et al. Block Combined Geocoding of Airborne InSAR with Sparse GCPs[J].Acta Geodaetica et Cartographica Sinica, 2013,42(3):397-403 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201303014.htm
[11] 黄国满,岳昔娟,赵争,等. 基于多项式正射纠正模型的机载SAR影像区域网平差[J]. 武汉大学学报·信息科学版,2008,33(6):569-572 http://ch.whu.edu.cn/CN/abstract/abstract1618.shtml Huang Guoman, Yue Xijuan, Zhao Zheng, et al. Block Adjustment with Airborne SAR Images Based on Polynomial Ortho-retification[J].Geomatics and Information Science of Wuhan University,2008, 33(6):569-572 http://ch.whu.edu.cn/CN/abstract/abstract1618.shtml
[12] 岳昔娟,韩春明,窦长勇,等. 机载InSAR区域网平差数学模型研究[J]. 武汉大学学报·信息科学版,2015,40(1):59-63 Yue Xijuan, Han Chunming, Dou Changyong, et al.Mathematical Model of Airborne InSAR Block Adjustment[J]. Geomatics and Information Science of Wuhan University,2015, 40(1):59-63
[13] Madsen S N, Skou N, Woelders K, et al. EMISAR Single Pass Topographic SAR Interferometer Modes[C]. Proceedings of 1996 IGARSS, Lincoln, USA, 1996
[14] 误差理论与测量平差基础[M]. 武汉:武汉大学出版社,2003 Error Theory and Fundation of Surveying Adjustment[M]. Wuhan:Wuhan University Press,2003
-
期刊类型引用(3)
1. 花奋奋，王萌萌，于洋，卢丽君，黄国满. 融合枝切法和最小费用流的多级相位解缠方法. 测绘科学. 2022(12): 138-146+183 . 
百度学术
2. 汪丙南，向茂生，蒋帅，付希凯. 联合对飞数据的单控制点机载干涉SAR定标算法. 测绘学报. 2018(11): 1495-1505 . 百度学术
3. 王彦军. 关于InSAR沉降监测结果精度比对的方法研究. 北京测绘. 2017(05): 109-111 . 百度学术
其他类型引用(4)
计量
- 文章访问数: 4540
- HTML全文浏览量: 154
- PDF下载量: 460
- 被引次数: 7
