【正文】 k2 g2 %移動曲面擬合法內(nèi)插 3″點坐標 zs=k。 zp=g。 dx=[1 1 1 0 0 1 1 1]39。 %待插值點距離周圍八個數(shù)據(jù)的 X 坐標差 dy=[1 0 1 1 1 1 0 1]39。 %待插值點距離周圍八個數(shù)據(jù)的 Y 坐標差 m1=dx.^2。 m2=dx.*dy。 m3= dy.^2。 m4=dx。 m5=dy。 m6=ones(8,1)。 %生成一個 8*1 元素為 1 的矩陣 M=[m1 m2 m3 m4 m5 m6]。 %系數(shù)矩陣 clear dx dy m1 m2 m3 m4 m5 m6 %清除垃圾變量 P=diag([ 1 1 1 1 ])。 %權陣 (待插點與周圍數(shù)據(jù)點的距離的 平方的倒數(shù) ) for i=4:3:3598 for j=4:3:3598 vk=double([k(i1,j1), k(i1,j), k(i1,j+1), k(i,j1), k(i,j+1), k(i+1,j1), k(i+1,j), k(i+1,j+1)]39。)。 %下面的 *運算要求數(shù)據(jù)為雙精度型 xk=inv(M39。*P*M)*M39。*P*vk。 zs(i,j)=xk(length(xk))。 %xk 中最后一個數(shù)據(jù)便為 3″插值 ， 將其賦值到 zs 中 vg=double([g(i1,j1), g(i1,j), g(i1,j+1), g(i,j1),g(i,j+1), g(i+1,j1), g(i+1,j), g(i+1,j+1)]39。)。 %下面的 *運算要求數(shù)據(jù)為雙精度型 xg= inv(M39。*P*M)*M39。*P*vg。 zp(i,j)=xg(length(xg))。 %xg 中最后一個數(shù)據(jù)便為 3″插值 ， 將其賦值到 zp 中 end end clear i j vk vg xk xg M P %清除垃圾變量 %比較 (aster)擬合前后數(shù)據(jù)之間的差值 shanqucha4=zs(1:3:3601,1:3:3601)k(1:3:3601,1:3:3601)。 %對應 3″數(shù)據(jù)點相減后所得的矩陣 pingyuancha4=zp(1:3:3601,1:3:3601)g(1:3:3601,1:3:3601)。 %對應 3″數(shù)據(jù)點相減后所得的矩陣 sqcmax4= double(max(max(shanqucha4)))。 pycmax4= double(max(max(pingyuancha4)))。 sqcmin4= double(min(min(shanqucha4)))。 pycmin4= double(min(min(pingyuancha4)))。 sqcmean4= double(mean2(shanqucha4))。 pycmean4= double(mean2(pingyuancha4))。 sqcstd4= double(std2(pingyuancha4))。 pycstd4= double(std2(pingyuancha4))。 商丘師范學 院 2021 屆本科畢業(yè)論文（設計） 11 jieguo4=[sqcmax4 sqcmin4 sqcmean4 sqcstd4。pycmax4 pycmin4 pycmean4 pycstd4] jieguo4 = clear shanqucha4 pingyuancha4 sqcmax4 pycmax4 sqcmin4 pycmin4 sqcmean4 pycmean4 sqcstd4 pycstd4 g k %比較擬合后的 aster 與 srtm 數(shù)據(jù) shanqucha5=srtmshanqu(1:1201,1:1201)zs(1:3:3601,1:3:3601)。 %對應 3″數(shù)據(jù)點相減后所得的矩陣 pingyuancha5=srtmpingyuan(1:1201,1:1201)zp(1:3:3601,1:3:3601)。 %對應 3″數(shù)據(jù)點相減后所得 的矩陣 sqcmax5= double(max(max(shanqucha5)))。 pycmax5= double(max(max(pingyuancha5)))。 sqcmin5= double(min(min(shanqucha5)))。 pycmin5= double(min(min(pingyuancha5)))。 sqcmean5= double(mean2(shanqucha5))。 pycmean5= double(mean2(pingyuancha5))。 sqcstd5= double(std2(shanqucha5))。 pycstd5= double(std2(pingyuancha5))。 jieguo5=[sqcmax5 sqcmin5 sqcmean5 sqcstd5。pycmax5 pycmin5 pycmean5 pycstd5] jieguo5 = clear shanqucha5 pingyuancha5 sqcmax5 pycmax5 sqcmin5 pycmin5 sqcmean5 pycmean5 sqcstd5 pycstd5 srtmpingyuan srtmshanqu zp zs