摄影测量学数字微分纠正

数字微分纠正《摄影测量学》（下）第五章武汉大学遥感信息工程学院摄影测量教研室数字微分纠正的概念框幅式中心投影影像的数字微分纠正线性阵列扫描影像的数字纠正主要内容数字微分纠正的概念根据参数与数字地面模型，利用相应的构像方程式，或按一定的数学模型用控制点解算，从原始非正射投影的数字影像获取正射影像，这种过程是将影像化为很多微小的区域逐一进行正射影像（DOM）图框幅式中心投影影像的数字微分纠正点元素纠正线元素纠正面元素纠正数字影像进行数字微分纠正，在原理上最适合点元素微分纠正。但能否真正做到点元素微分纠正，它取决于能否真实地测定每个像元的物方坐标X，Y，Z数字纠正:像素的几何位置和灰度数字微分纠正的基本原理x=fx(X，Y)；y=fy(X，Y)X=φx(x，y)；Y=φy(x，y)反解法正解法数字微分纠正与光学微分纠正一样，是实现两个二维图像之间的几何变换反解法（间接法）数字微分纠正计算地面点坐标X＝X0+M·X’Y＝Y0+M·Y’计算像点坐标)()()()()()()()()()()()(33322203331110ssssssssssssZZcYYbXXaZZcYYbXXafyyZZcYYbXXaZZcYYbXXafxxDEM(X0,Y0)fJJIInnmmZZYYXXcbacbacbafyyxxSSS0021213332221110001000011111098765111094321LYLXLLZLYLXLJLYLXLLZLYLXLIX，Y，Z直接解求扫描坐标行、列号I，JSSSZZYYXXcbacbacbannmmfJJII3332221112121001000''0''数字摄影测量的内定向fJJIInnmmfyyxx0021210010000根据框标解算参数oxyIJ反解法解算流程XY纠正影像原始影像xy灰度内插反算正解法（直接法）数字微分纠正fcycxcfbybxbZYfcycxcfayaxaZX321321321321原始影像正射影像Z?非规则排列的二维图像三维空间（X，Y，Z）正解法缺点正解法解算流程•••••••••••••••••••••••••xyXY原始影像纠正影像数字纠正实际解法及分析以“面元素”作为“纠正单元”])()())([(1),(])()())([(1),(3421234212ijyjyinyjniyjninnjiyijxjxinxjnixjninnjix反算公式计算该纠正单元4个“角点”的像点坐标，而纠正单元内的坐标则用双线性内插求得数字纠正实际解法及分析线性阵列扫描影像的数字纠正由若干条线性阵列扫描影像构成像幅)()()()()()()()()()()()(10333222111tZZtYYtXXtctbtatctbtatctbtafxSSSSPOT卫星（6000条扫描线组成一幅影像）xy3000条扫描Lp飞行方向线性阵列扫描影像间接法纠正))]()()()()()(()()()([10222222tctZtbtYtatXtZctYbtXaSSSXa2(t)+Yb2(t)+Zc2(t)=A(t))()()()()()()(222tctZtbtYtatXtASSS或其中按泰勒级数展开为2222222222222222)2()1()0()2()1()0()2()1()0()2()1()0()()()()(tAtAAtAtctcctctbtbbtbtataata)(])()[()1(2222)2(222)0(222)1()1()1()2()2()2()0()0()0(AZcYbXatAZcYbXaAZcYbXat取二次项0][][][2)2(222)1(222)0(222)2()2()2()1()1()1()0()0()0(tAZcYbXatAZcYbXaAZcYbXa线性阵列扫描影像间接法纠正求出各元素对应的a1(t)，a2(t)，…c3(t)，XS(t)，YS(t)，ZS(t)，然后求出该相应像点的y及t，或y及x。211112111121111)2()1()0()2()1()0()2()1()0()()()(tctcctctbtbbtbtataata对a1（t），b1（t），…用多项式表达tllyP)(0t值实际上表达了上图坐标系中像点p在时刻t的y坐标lp，l0扫描线行数；δ为CCD一个探测像元的宽度μ为扫描线的时间间隔)]())(()())(()())([(1111tctZZtbtYYtatXXxSSS)]())(()())(()())([(1333tctZZtbtYYtatXXfSSS再求像点p的x坐标)())(()())(()())(()())(()())(()())((333111tctZZtbtYYtatXXtctZZtbtYYtatXXfxSSSSSS或tZZtZtYYtYtXXtXttttttSSSSSSSSS)0()()0()()0()()0()()0()()0()(为影像中心行外方位元素通常可认为各参数是t的线性函数为变化率参数直接法纠正))(()()()()()())(()()()()()(31313131tZZftcxtcftbxtbtYYtZZftcxtcftaxtatXXSSSSa1(t)，a2(t)，…c3(t)，XS(t)，YS(t)，ZS(t)为像点（x，y）对应的外方位元素，可由其行号li计算：SiSSiSiSSiSiSSiiiiiiiZllZZYllYYXllXXllllll)()()()()()(000000000000直接法纠正给定高程初始值Z0计算出地面平面坐标近似值（X1，Y1）用DEM与（X1，Y1）内插出高程Z1,对应的地面点（X，Y，Z，）直接法与间接法相结合的纠正方案影像规则格网点对应的地面坐标的解算得到地面一非规则格网直接法))(()()()()()())(()()()()()(31313131tZZftcxtcftbxtbtYYtZZftcxtcftaxtatXXSSSS内插出地面规则格网点该四点的地面坐标（X’11，Y’11），（X’12，Y’12），（X’21，Y’21），（X’22，Y’22）及其影像坐标（x’11，y’11），（x’12，y’12），（x’21，y’21），（x’22，y’22）拟合两平面：p11p12p21p22P’11P’12P’13P’14间接法进行纠正''''''210210YbXbbyYaXaax由Pll点的地面坐标（Xll，Y11）,计算其相应的影像坐标（x11，y11）。多项式纠正影像变形规律可近似地看作为平移、缩放、旋转、仿射、偏扭、弯曲等基本形变的合成)()(25423210iiiiiiiYcYXcXcYcXccx)(29282736iiiiiiYcYXcYXcXc)()(25423210iiiiiiiYdYXdXdYdXddy)(29282736iiiiiiYdYXdYXdXd反解法的多项式对每个控制点，已知其地面坐标Xi，Yi，可以列出上面两个方程Δxi=xi（近似计算值）-xi’（量测值）Δyi=yi（近似计算值）-yi’（量测值）iSiiiiSiZZYfyXZZfxarctg正解法的多项式)()(25423210iiiiiiiyayxaxayaxaaX)(39282736iiiiiiyayxayxaxa)()(25423210iiiiiiiybyxbxbybxbbY)(39282736iiiiiiybyxbyxbxbΔXi＝Xi（已知值）-Xi（近似计算值）ΔYi＝Yi（已知值）-Yi（近似计算值正射影像精度的检查与质量控制野外检测：检查正射影像的绝对精度，与等高线图或线划地图套合后进行目视检查左影像和右影像制作同一地区的两幅正射影像，量测两幅正射影像上同名点的视差。正射影像的影像质量一般采用目视检查，有合适的反差，均匀的色调接边不仅涉及几何方面的精度问题，还涉及不同影像之间色调的不一致正射影像精度的检查与质量控制无缝镶嵌Seamline