《武汉大学数字摄影测量学ppt课件.ppt》由会员分享,可在线阅读,更多相关《武汉大学数字摄影测量学ppt课件.ppt(37页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、数字微分纠正摄影测量学摄影测量学(下)(下)第五第五章章遥感信息工程学院遥感信息工程学院 摄影测量教研室摄影测量教研室n数字微分纠正的概念数字微分纠正的概念n框幅式中心投影影像的数字微分纠框幅式中心投影影像的数字微分纠正正n线性阵列扫描影像的数字纠正线性阵列扫描影像的数字纠正主要内容主要内容 数字微分纠正的概念根据参数与数字地面模型,利用相应的构像方程式,或按一定的数学模型用控制点解算,从原始非正射投影的数字影像获取正射影像,这种过程是将影像化为很多微小的区域逐一进行正射影像(正射影像(DOM)图)图框幅式中心投影影像的数字微分纠正框幅式中心投影影像的数字微分纠正点元素纠正 线元素纠正 面元素
2、纠正 数字影像进行数字微分纠正,在原理上最适合点元素微分纠正。但能否真正做到点元素微分纠正,它取决于能否真实地测定每个像元的物方坐标X,Y,Z数字纠正: 像素的几何位置和灰度数字微分纠正的基本原理数字微分纠正的基本原理 x = fx (X,Y) ; y = fy (X,Y) X = x (x,y); Y = y (x,y) 反解法正解法 数字微分纠正与光学微分纠正一样,是实现两个二维图像之间的几何变换 反解法(间接法)数字微分纠正反解法(间接法)数字微分纠正l 计算地面点坐标X X0 + M X Y Y0 + M Y 计算像点坐标)()()()()()()()()()()()(33322203
3、331110ssssssssssssZZcYYbXXaZZcYYbXXafyyZZcYYbXXaZZcYYbXXafxxDEM(X0,Y0)fJJIInnmmZZYYXXcbacbacbafyyxxSSS0021213332221110001000011111098765111094321LYLXLLZLYLXLJLYLXLLZLYLXLIX,Y,Z直接解求扫描坐标行、列号I,J SSSZZYYXXcbacbacbannmmfJJII33322211121210010000数字摄影测量的内定向fJJIInnmmfyyxx0021210010000l根据框标解算参数oxyIJ反解法解算流程反解法
4、解算流程XY纠正影像原始影像xy灰度内插反算正解法(直接法)数字微分纠正正解法(直接法)数字微分纠正fcycxcfbybxbZYfcycxcfayaxaZX321321321321原始影像正射影像Z? 非规则排列的 二维图像 三维空间(X,Y,Z)正解法正解法缺点正解法解算流程正解法解算流程xy XY原始影像纠正影像数字纠正实际解法及分析数字纠正实际解法及分析 以“面元素”作为“纠正单元” )()()(1),()()()(1),(3421234212ijyjyinyjniyjninnjiyijxjxinxjnixjninnjix反算公式计算该纠正单元4个“角点”的像点坐标,而纠正单元内的坐标则
5、用双线性内插求得 数字纠正实际解法及分析数字纠正实际解法及分析 线性阵列扫描影像的数字纠正线性阵列扫描影像的数字纠正由若干条线性阵列扫描影像构成像幅 )()()()()()()()()()()()(10333222111tZZtYYtXXtctbtatctbtatctbtafxSSSSPOT卫星(6000条扫描线组成一幅影像) xy3000条扫描Lp飞行方向线性阵列扫描影像线性阵列扫描影像间接法纠正间接法纠正 )()()()()()()()()(10222222tctZtbtYtatXtZctYbtXaSSS Xa2(t) + Yb2(t) + Zc2(t) =A(t) )()()()()()
6、()(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取二次项 02) 2(222) 1 (222) 0(222)2()2()2() 1 () 1 () 1 ()0()0()0(tAZcYbXatAZcYbXaAZcYbX
7、a线性阵列扫描影像线性阵列扫描影像间接法纠正间接法纠正 求出各元素对应的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)()()()()()(1333tct
8、ZZtbtYYtatXXfSSS再求像点p的 x坐标 )()()()()()()()()()()()(333111tctZZtbtYYtatXXtctZZtbtYYtatXXfxSSSSSS或tZZtZtYYtYtXXtXttttttSSSSSSSSS)0()()0()()0()()0()()0()()0()(为影像中心行外方位元素通常可认为各参数是t的线性函数 为变化率参数直接法纠正直接法纠正)()()()()()()()()()()()(31313131tZZftcxtcftbxtbtYYtZZftcxtcftaxtatXXSSSSa1(t),a2(t),c3(t),XS(t),YS(t)
9、,ZS(t)为像点(x,y)对应的外方位元素,可由其行号li计算:SiSSiSiSSiSiSSiiiiiiiZllZZYllYYXllXXllllll)()()()()()(000000000000直接法纠正直接法纠正给定高程初始值Z0计算出地面平面坐标近似值(X1,Y1) 用DEM与(X1,Y1)内插出高程Z1,对应的地面点(X,Y,Z,)直接法与间接法相结合的纠正方案直接法与间接法相结合的纠正方案影像规则格网点对应的地面坐标的解算得到地面一非规则格网 直接法直接法)()()()()()()()()()()()(31313131tZZftcxtcftbxtbtYYtZZftcxtcftaxt
10、atXXSSSS内插出地面规则格网点 该四点的地面坐标(X11,Y11),(X12,Y12),(X21,Y21),(X22,Y 2 2) 及 其 影 像 坐 标(x11,y11),(x12,y12),(x21,y21),(x22,y22)拟合两平面:p11p12p21p22P11P12P13P14间接法进行纠正 210210YbXbbyYaXaax由Pll点的地面坐标(Xll,Y11),计算其相应的影像坐标(x11,y11)。多项式纠正多项式纠正影像变形规律可近似地看作为平移、缩放、旋转、仿射、偏扭、弯曲等基本形变的合成 )()(25423210iiiiiiiYcYXcXcYcXccx)(29
11、282736iiiiiiYcYXcYXcXc)()(25423210iiiiiiiYdYXdXdYdXddy)(29282736iiiiiiYdYXdYXdXd反解法的多项式 对每个控制点,已知其地面坐标Xi,Yi,可以列出上面两个方程 xi = xi(近似计算值) - xi(量测值) yi = yi(近似计算值) - yi(量测值) iSiiiiSiZZYfyXZZfxarctg正解法的多项式 )()(25423210iiiiiiiyayxaxayaxaaX)(39282736iiiiiiyayxayxaxa)()(25423210iiiiiiiybyxbxbybxbbY)(39282736iiiiiiybyxbyxbxb Xi Xi(已知值)- Xi(近似计算值) Yi Yi(已知值)- Yi(近似计算值 正射影像精度的检查与质量控制正射影像精度的检查与质量控制l野外检测野外检测:检查正射影像的绝对精度,l与等高线图或线划地图套合后进行目视检查 l左影像和右影像制作同一地区的两幅正射影像,量测两幅正射影像上同名点的视差。正射影像的影像质量一般采用目视检查,有合适的反差,均匀的色调 接边不仅涉及几何方面的精度问题,还涉及不同影像之间色调的不一致正射影像精度的检查与质量控制正射影像精度的检查与质量控制无缝镶嵌无缝镶嵌Seam line
限制150内