數(shù)字微分糾正第1頁(yè)，共37頁(yè)，2023年，2月20日，星期五數(shù)字微分糾正的概念框幅式中心投影影像的數(shù)字微分糾正線性陣列掃描影像的數(shù)字糾正主要內(nèi)容

第2頁(yè)，共37頁(yè)，2023年，2月20日，星期五數(shù)字微分糾正的概念根據(jù)參數(shù)與數(shù)字地面模型，利用相應(yīng)的構(gòu)像方程式，或按一定的數(shù)學(xué)模型用控制點(diǎn)解算，從原始非正射投影的數(shù)字影像獲取正射影像，這種過(guò)程是將影像化為很多微小的區(qū)域逐一進(jìn)行/“測(cè)繪信息網(wǎng)”網(wǎng)友搜集,版權(quán)歸原權(quán)利人所有第3頁(yè)，共37頁(yè)，2023年，2月20日，星期五正射影像（DOM）圖第4頁(yè)，共37頁(yè)，2023年，2月20日，星期五框幅式中心投影影像的數(shù)字微分糾正點(diǎn)元素糾正

線元素糾正

面元素糾正

數(shù)字影像進(jìn)行數(shù)字微分糾正，在原理上最適合點(diǎn)元素微分糾正。但能否真正做到點(diǎn)元素微分糾正，它取決于能否真實(shí)地測(cè)定每個(gè)像元的物方坐標(biāo)X，Y，Z第5頁(yè)，共37頁(yè)，2023年，2月20日，星期五數(shù)字糾正:像素的幾何位置和灰度數(shù)字微分糾正的基本原理

x=fx(X，Y)；

y=fy(X，Y)X=φx(x，y)；Y=φy(x，y)

反解法正解法

數(shù)字微分糾正與光學(xué)微分糾正一樣，是實(shí)現(xiàn)兩個(gè)二維圖像之間的幾何變換

/“測(cè)繪信息網(wǎng)”網(wǎng)友搜集,版權(quán)歸原權(quán)利人所有第6頁(yè)，共37頁(yè)，2023年，2月20日，星期五反解法（間接法）數(shù)字微分糾正

計(jì)算地面點(diǎn)坐標(biāo)X

＝X0

+

M·X’

Y＝Y(jié)0+M·Y’

計(jì)算像點(diǎn)坐標(biāo)DEM(X0,Y0)第7頁(yè)，共37頁(yè)，2023年，2月20日，星期五X，Y，Z直接解求掃描坐標(biāo)行、列號(hào)I，J

第8頁(yè)，共37頁(yè)，2023年，2月20日，星期五數(shù)字?jǐn)z影測(cè)量的內(nèi)定向根據(jù)框標(biāo)解算參數(shù)oxyIJ第9頁(yè)，共37頁(yè)，2023年，2月20日，星期五反解法解算流程XY糾正影像原始影像xy灰度內(nèi)插反算第10頁(yè)，共37頁(yè)，2023年，2月20日，星期五正解法（直接法）數(shù)字微分糾正原始影像正射影像Z?第11頁(yè)，共37頁(yè)，2023年，2月20日，星期五非規(guī)則排列的

二維圖像三維空間（X，Y，Z）正解法缺點(diǎn)第12頁(yè)，共37頁(yè)，2023年，2月20日，星期五正解法解算流程?????????????????????????xy

XY原始影像糾正影像第13頁(yè)，共37頁(yè)，2023年，2月20日，星期五數(shù)字糾正實(shí)際解法及分析

以“面元素”作為“糾正單元”

反算公式計(jì)算該糾正單元4個(gè)“角點(diǎn)”的像點(diǎn)坐標(biāo)，而糾正單元內(nèi)的坐標(biāo)則用雙線性內(nèi)插求得第14頁(yè)，共37頁(yè)，2023年，2月20日，星期五數(shù)字糾正實(shí)際解法及分析

第15頁(yè)，共37頁(yè)，2023年，2月20日，星期五線性陣列掃描影像的數(shù)字糾正由若干條線性陣列掃描影像構(gòu)成像幅第16頁(yè)，共37頁(yè)，2023年，2月20日，星期五SPOT衛(wèi)星（6000條掃描線組成一幅影像）

xy3000條掃描Lp飛行方向/“測(cè)繪信息網(wǎng)”網(wǎng)友搜集,版權(quán)歸原權(quán)利人所有第17頁(yè)，共37頁(yè)，2023年，2月20日，星期五線性陣列掃描影像間接法糾正

Xa2(t)+Yb2(t)+Zc2(t)=A(t)或其中第18頁(yè)，共37頁(yè)，2023年，2月20日，星期五按泰勒級(jí)數(shù)展開為

第19頁(yè)，共37頁(yè)，2023年，2月20日，星期五取二次項(xiàng)線性陣列掃描影像間接法糾正

第20頁(yè)，共37頁(yè)，2023年，2月20日，星期五求出各元素對(duì)應(yīng)的a1(t)，a2(t)，…c3(t)，XS(t)，YS(t)，ZS(t)，然后求出該相應(yīng)像點(diǎn)的y及t，或y及x。對(duì)a1（t），b1（t），…用多項(xiàng)式表達(dá)第21頁(yè)，共37頁(yè)，2023年，2月20日，星期五t值實(shí)際上表達(dá)了上圖坐標(biāo)系中像點(diǎn)p在時(shí)刻t的y坐標(biāo)lp，l0掃描線行數(shù)；δ為CCD一個(gè)探測(cè)像元的寬度μ為掃描線的時(shí)間間隔第22頁(yè)，共37頁(yè)，2023年，2月20日，星期五再求像點(diǎn)p的

x坐標(biāo)

或/“測(cè)繪信息網(wǎng)”網(wǎng)友搜集,版權(quán)歸原權(quán)利人所有第23頁(yè)，共37頁(yè)，2023年，2月20日，星期五為影像中心行外方位元素通?？烧J(rèn)為各參數(shù)是t的線性函數(shù)

為變化率參數(shù)第24頁(yè)，共37頁(yè)，2023年，2月20日，星期五直接法糾正a1(t)，a2(t)，…c3(t)，XS(t)，YS(t)，ZS(t)為像點(diǎn)（x，y）對(duì)應(yīng)的外方位元素，可由其行號(hào)li計(jì)算：第25頁(yè)，共37頁(yè)，2023年，2月20日，星期五直接法糾正第26頁(yè)，共37頁(yè)，2023年，2月20日，星期五給定高程初始值Z0計(jì)算出地面平面坐標(biāo)近似值（X1，Y1）

用DEM與（X1，Y1）內(nèi)插出高程Z1,對(duì)應(yīng)的地面點(diǎn)（X，Y，Z，）第27頁(yè)，共37頁(yè)，2023年，2月20日，星期五直接法與間接法相結(jié)合的糾正方案影像規(guī)則格網(wǎng)點(diǎn)對(duì)應(yīng)的地面坐標(biāo)的解算得到地面一非規(guī)則格網(wǎng)

直接法第28頁(yè)，共37頁(yè)，2023年，2月20日，星期五內(nèi)插出地面規(guī)則格網(wǎng)點(diǎn)

該四點(diǎn)的地面坐標(biāo)（X’11，Y’11），（X’12，Y’12），（X’21，Y’21），（X’22，Y’22）及其影像坐標(biāo)（x’11，y’11），（x’12，y’12），（x’21，y’21），（x’22，y’22）擬合兩平面：p11p12p21p22P’11P’12P’13P’14第29頁(yè)，共37頁(yè)，2023年，2月20日，星期五間接法進(jìn)行糾正

由Pll點(diǎn)的地面坐標(biāo)（Xll，Y11）,計(jì)算其相應(yīng)的影像坐標(biāo)（x11，y11）。第30頁(yè)，共37頁(yè)，2023年，2月20日，星期五多項(xiàng)式糾正影像變形規(guī)律可近似地看作為平移、縮放、旋轉(zhuǎn)、仿射、偏扭、彎曲等基本形變的合成

反解法的多項(xiàng)式

第31頁(yè)，共37頁(yè)，2023年，2月20日，星期五對(duì)每個(gè)控制點(diǎn)，已知其地面坐標(biāo)Xi，Yi，可以列出上面兩個(gè)方程

Δxi=xi（近似計(jì)算值）-xi’（量測(cè)值）

Δyi=yi（近似計(jì)算值）

-yi’（量測(cè)值）

第32頁(yè)，共37頁(yè)，2023年，2月20日，星期五正解法的多項(xiàng)式

ΔXi

＝Xi（已知值）-Xi（近似計(jì)算值）

ΔYi

＝Y(jié)i（已知值）-Yi（近似計(jì)算值第33頁(yè)，共37頁(yè)，2023年，2月20日，星期五正射影像精度的檢查與質(zhì)量控制野外檢測(cè)：檢查正射影像的絕對(duì)精度，與等高線圖或線劃地圖套合后進(jìn)行目視檢查左影像和右影像制作同一地區(qū)的兩幅正射影像，量測(cè)兩幅正射影像上同名點(diǎn)的視差。第34頁(yè)，共37頁(yè)，2023年，2月20日，星期五