主权项 |
一种自适应几何关系强弱的有理函数模型区域网平差方法,它包括以下步骤:1)读取卫星遥感影像的有理函数模型参数、归一化因子,以及连接点、控制点和检查点的影像坐标;2)根据读取的有理函数模型参数、归一化因子和连接点坐标,计算每一个连接点的所有交会角,并取每一个连接点的最大的交会角作为立体几何关系强弱的量化指标;3)根据每一个连接点最大交会角大小,将连接点分为第一类连接点和第二类连接点;4)根据连接点的最大交会角大小,算法自适应确定连接点的物方大地坐标初始值,具体过程为:针对第一类连接点:首先获取影像范围内的高差dh=Height<sub>max</sub>‑Height<sub>min</sub>,式中的Height<sub>max</sub>表示最大高程;Height<sub>min</sub>表示最小高程;然后计算因高差引起的最大投影差ds=dh×tanθ<sub>max</sub>,式中的θ<sub>max</sub>表示连接点的最大交会角,当ds≤k×r时,高程初值取值为参数HEIGHT_OFF,经纬度初值用有理函数模型和高程初值直接求解,式中的k表示精度要求;r表示影像分辨率;当ds>k×r时,利用初始高程HEIGHT_OFF和有理函数模型参数,直接求解出初始的大地坐标,然后利用该大地坐标经纬度在辅助数字高程模型上内插出新的高程值,比较新的高程值与原始高程的差值,如果差值较小,则该经纬度与高程作为解算初值,如果差值较大,则利用新的高程替代HEIGHT_OFF为初值迭代求解;针对第二类连接点:利用有理函数模型空间前方交会的方式直接得到初始值;5)建立误差方程,进行区域网平差,其中:对于控制点,其线性化误差方程为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>F</mi><mi>x</mi></msub><mo>=</mo><msub><mi>F</mi><mrow><mi>x</mi><mn>0</mn></mrow></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>x</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>e</mi><mn>0</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δe</mi><mn>0</mn></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>x</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>e</mi><mn>1</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δe</mi><mn>1</mn></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>x</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>e</mi><mn>2</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δe</mi><mn>2</mn></msub><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA0001081283810000011.GIF" wi="798" he="119" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>F</mi><mi>y</mi></msub><mo>=</mo><msub><mi>F</mi><mrow><mi>y</mi><mn>0</mn></mrow></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>y</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>f</mi><mn>0</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δf</mi><mn>0</mn></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>y</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>f</mi><mn>1</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δf</mi><mn>1</mn></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>y</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>f</mi><mn>2</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δf</mi><mn>2</mn></msub><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA0001081283810000012.GIF" wi="790" he="122" /></maths>对于第一类连接点,其线性化误差方程为<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>F</mi><mi>x</mi></msub><mo>=</mo><msub><mi>F</mi><mrow><mi>x</mi><mn>0</mn></mrow></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>x</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>e</mi><mn>0</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δe</mi><mn>0</mn></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>x</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>e</mi><mn>1</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δe</mi><mn>1</mn></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>x</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>e</mi><mn>2</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δe</mi><mn>2</mn></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>x</mi></msub></mrow><mrow><mo>∂</mo><mi>l</mi><mi>a</mi><mi>t</mi></mrow></mfrac><mo>·</mo><mi>Δ</mi><mi>l</mi><mi>a</mi><mi>t</mi><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>x</mi></msub></mrow><mrow><mo>∂</mo><mi>l</mi><mi>o</mi><mi>n</mi></mrow></mfrac><mo>·</mo><mi>Δ</mi><mi>l</mi><mi>o</mi><mi>n</mi><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA0001081283810000013.GIF" wi="1293" he="124" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>F</mi><mi>y</mi></msub><mo>=</mo><msub><mi>F</mi><mrow><mi>y</mi><mn>0</mn></mrow></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>y</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>f</mi><mn>0</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δf</mi><mn>0</mn></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>y</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>f</mi><mn>1</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δf</mi><mn>1</mn></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>y</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>f</mi><mn>2</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δf</mi><mn>2</mn></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>y</mi></msub></mrow><mrow><mo>∂</mo><mi>l</mi><mi>a</mi><mi>t</mi></mrow></mfrac><mo>·</mo><mi>Δ</mi><mi>l</mi><mi>a</mi><mi>t</mi><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>y</mi></msub></mrow><mrow><mo>∂</mo><mi>l</mi><mi>o</mi><mi>n</mi></mrow></mfrac><mo>·</mo><mi>Δ</mi><mi>l</mi><mi>o</mi><mi>n</mi><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA0001081283810000014.GIF" wi="1285" he="135" /></maths>对于第二类连接点,其线性化误差方程为<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>F</mi><mi>x</mi></msub><mo>=</mo><msub><mi>F</mi><mrow><mi>x</mi><mn>0</mn></mrow></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>x</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>e</mi><mn>0</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δe</mi><mn>0</mn></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>x</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>e</mi><mn>1</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δe</mi><mn>1</mn></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>x</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>e</mi><mn>2</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δe</mi><mn>2</mn></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>x</mi></msub></mrow><mrow><mo>∂</mo><mi>l</mi><mi>a</mi><mi>t</mi></mrow></mfrac><mo>·</mo><mi>Δ</mi><mi>l</mi><mi>a</mi><mi>t</mi><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>x</mi></msub></mrow><mrow><mo>∂</mo><mi>l</mi><mi>o</mi><mi>n</mi></mrow></mfrac><mo>·</mo><mi>Δ</mi><mi>l</mi><mi>o</mi><mi>n</mi><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>x</mi></msub></mrow><mrow><mo>∂</mo><mi>h</mi></mrow></mfrac><mo>·</mo><mi>Δ</mi><mi>h</mi><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA0001081283810000021.GIF" wi="1478" he="127" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>F</mi><mi>y</mi></msub><mo>=</mo><msub><mi>F</mi><mrow><mi>y</mi><mn>0</mn></mrow></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>y</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>f</mi><mn>0</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δf</mi><mn>0</mn></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>y</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>f</mi><mn>1</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δf</mi><mn>1</mn></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>y</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>f</mi><mn>2</mn></msub></mrow></mfrac><mo>·</mo><msub><mi>Δf</mi><mn>2</mn></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>y</mi></msub></mrow><mrow><mo>∂</mo><mi>l</mi><mi>a</mi><mi>t</mi></mrow></mfrac><mo>·</mo><mi>Δ</mi><mi>l</mi><mi>a</mi><mi>t</mi><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>y</mi></msub></mrow><mrow><mo>∂</mo><mi>l</mi><mi>o</mi><mi>n</mi></mrow></mfrac><mo>·</mo><mi>Δ</mi><mi>l</mi><mi>o</mi><mi>n</mi><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>F</mi><mi>y</mi></msub></mrow><mrow><mo>∂</mo><mi>h</mi></mrow></mfrac><mo>·</mo><mi>Δ</mi><mi>h</mi><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA0001081283810000022.GIF" wi="1478" he="127" /></maths>上式中,F<sub>x</sub>表示方程在列方向的误差,F<sub>y</sub>表示方程在行方向的误差;F<sub>x0</sub>、F<sub>y0</sub>表示泰勒级数0次项,e<sub>0</sub>,e<sub>1</sub>,e<sub>2</sub>和f<sub>0</sub>,f<sub>1</sub>,f<sub>2</sub>为仿射变换参数;x、y为精化后的影像行列坐标;<img file="FDA0001081283810000023.GIF" wi="652" he="63" />表示经度、纬度和高程的一阶偏导数;lat表示经度;lon表示纬度;h表示高程;6)输出有理函数模型精化参数,完成基于有理函数模型的区域网平差。 |