【正文】 ?？傊诋厴I(yè)設(shè)計(jì)這個(gè)環(huán)節(jié)中，X老師嚴(yán)謹(jǐn)治學(xué)的態(tài)度對(duì)我有很深的影響，希望今后對(duì)待學(xué)業(yè)事業(yè)也能夠如X老師般嚴(yán)謹(jǐn)認(rèn)真。[x,y]=size(A)。 P=eye(a)。 %中心化后的x坐標(biāo)矩陣n=B(:,2)a2。*P*V/(a1))C=load(39。a5=ones(s,1)。T=eye(s)。)ylabel(39。) end endcount。)title(39。)。r1=ones(a,1)。 P=eye(a)。D=inv(S39。E=sum(abs(l1))。./(E.*H)。D1=inv(T39。 r4(:,1)=a1。 k2=[m2 n2 l2]。BB=AA(:,g)。)xlabel(39。 text(i,T(i,1),39。)ylabel(39。) endendcount1。 subplot(1,1,1)%figure(2) [c,h] =contour(x,y,z,15)。)共 38 頁(yè) 第 34 頁(yè)。v439。 text(i,V(i,1),39。)xlabel(39。)count=0。subplot(2,1,2)plot(t0,T,39。Sdj]。 %中心化后的x坐標(biāo)矩陣n2=C(:,2)r5。[s,t]=size(C)。T=S(:,g)。endH=reshape(G,a,1)。K=l139。Sij]。 %中心化后的y坐標(biāo)矩陣l1=10.*B(:,3)。)。) endendcount1附錄 2 多面函數(shù)內(nèi)插MATLAB代碼基于正交最小二乘多面函數(shù)擬合代碼A=load(39。)ylabel(39。 text(i,Q(i,1),39。)xlabel(39。M=[ones(s,1) z u z.^2 z.*u u.^2 z.^3 (z.^2).*u z.*(u.^2) u.^3 z.^4 (z.^3).*u (z.^2).*(u.^2) z.*(u.^3) u.^4 z.^5 (z.^4).*u (z.^3).*(u.^2) (z.^2).*(u.^3) z.*(u.^4) u.^5]。a4=ones(s,1)。V=N*X+a310.*B(:,3)。 a3(:,1)=k3。)。F:\matlab\toolbox\nmm\data\dat\39。首先我要感謝XXX老師，感謝張老師在畢業(yè)設(shè)計(jì)上對(duì)我進(jìn)行的悉心的指導(dǎo)。另外，在擬合前未對(duì)趨勢(shì)面進(jìn)行分析。本文嘗試采用多項(xiàng)式擬合法與多面函數(shù)擬合法進(jìn)行空間數(shù)據(jù)插值，以此填補(bǔ)缺失數(shù)據(jù)。另外，四種核函數(shù)相比較而言，正雙曲核函數(shù)擬合、檢核效果較好，這與擬合點(diǎn)以及待插點(diǎn)屬性密切相關(guān)。通過(guò)本實(shí)例截取的西安市地面沉降監(jiān)測(cè)數(shù)據(jù)分析可以發(fā)現(xiàn)，實(shí)例分析所用的地形情況相對(duì)而言數(shù)據(jù)變化不是很劇烈，分布也還適中，另外由于本實(shí)例選用的數(shù)據(jù)比較集中，平緩所以擬合、檢核的精度也相對(duì)比較好。 本實(shí)例分析基于MATLAB數(shù)據(jù)實(shí)驗(yàn)軟件，采用了基于DInSAR的西安市地面沉降監(jiān)測(cè)數(shù)據(jù)，應(yīng)用多項(xiàng)式曲面擬合方法對(duì)缺失點(diǎn)的地面沉降量進(jìn)行了擬合內(nèi)插。各沉降槽一般是北深南淺．地裂縫南側(cè)沉降量大，形成地形變陡變帶，地形上多呈陡坎或陡坡。內(nèi)符合精度數(shù)是隨著次數(shù)的增加會(huì)提高，但是外符合精度在一定次數(shù)范圍內(nèi)先增加后降低。由于平滑因子的選取與和函數(shù)的選擇關(guān)系密切，而且選取的隨機(jī)性較強(qiáng)，一定情況下對(duì)內(nèi)插結(jié)果影響并不太大，所以在本文中暫不予以討論。本文在對(duì)核函數(shù)的選擇上比較謹(jǐn)慎，不同核函數(shù)的多面函數(shù)內(nèi)插結(jié)果的精度會(huì)有較大的差別。在上節(jié)中解多面函數(shù)系數(shù)時(shí)，由于A中每一列向量的貢獻(xiàn)是不同的，這意味著每一節(jié)點(diǎn)對(duì)內(nèi)插值的影響不同，因此可以根據(jù)各節(jié)點(diǎn)貢獻(xiàn)的大小來(lái)選擇核函數(shù)的節(jié)點(diǎn)。()式～()式列出了空間數(shù)據(jù)擬合中常用三種核函數(shù)具體表達(dá)形式。它的基本思想是任何一個(gè)規(guī)則或不規(guī)則的連續(xù)曲面均可以由若干簡(jiǎn)單面(或稱(chēng)單值數(shù)學(xué)面)來(lái)疊加逼近。誤差方程可表示為： ()根據(jù)最小二乘原理，解得： ()將上式解得的代入()式，便可獲得確定的曲面函數(shù)模型，從而插求出其它未知點(diǎn)的地面沉降量。由于雷達(dá)衛(wèi)星本身的姿態(tài)偏轉(zhuǎn)以及衛(wèi)星軌道的誤差，ERS系列雷達(dá)圖像如果不借助外部控制點(diǎn)，其圖像的定位精度大約在一個(gè)像元以?xún)?nèi)(20米左右)。}為數(shù)學(xué)期望，上標(biāo)*為復(fù)共軛,是復(fù)相關(guān)系數(shù)，M為主圖像的復(fù)值，S為從圖像的復(fù)值。b、重采樣獲取了從圖像相對(duì)于主圖像的配準(zhǔn)參數(shù)之后，需要對(duì)從圖像進(jìn)行重采樣，使從圖像上的每一個(gè)像元值精確的對(duì)應(yīng)于主圖像上的每一個(gè)像元。 以下簡(jiǎn)要說(shuō)明這五個(gè)主要步驟所采用的處理方法： SAR數(shù)據(jù)的讀取 從單視復(fù)數(shù)圖SLC數(shù)據(jù)中讀取干涉測(cè)量所必需的參數(shù)以及SAR圖像的基本參數(shù)是差分干涉測(cè)量的第一步。　　以上給出的解釋模型是針對(duì)兩部天線分別獨(dú)立的發(fā)射與接收回波的系統(tǒng)推導(dǎo)而來(lái)的，對(duì)于一部天線發(fā)射、兩部天線同時(shí)接收回波的系統(tǒng)，只要把公式中的代之以就可以了。如果對(duì)某一地物兩次成像時(shí)天線都處于相同的空間位置，雷達(dá)信號(hào)往返的路徑之差應(yīng)該是天線視線方向的變化，具體可以表現(xiàn)為相位差。論文結(jié)合西安市地面沉降InSAR監(jiān)測(cè)數(shù)據(jù)，探討多項(xiàng)式擬合和多面函數(shù)擬合方法。而時(shí)間失相干則是衛(wèi)星兩次飛行期間，由于地表擾動(dòng)(如農(nóng)業(yè)活動(dòng))、植被生長(zhǎng)、風(fēng)的影響或地物濕度的變化引起地面物體散射特性的變化，導(dǎo)致失相干?？傊珼InSAR技術(shù)作為新興的用于地面沉降監(jiān)測(cè)的手段已突顯其優(yōu)勢(shì)。DInSAR, Differential InSAR,即重復(fù)軌道差分干涉測(cè)量技術(shù)，是InSAR應(yīng)用的一個(gè)拓展，雷達(dá)干涉圖的差分可用于監(jiān)測(cè)雷達(dá)視線方向厘米級(jí)或更微小的地球表面形變，以揭示許多地球物理現(xiàn)象，如地震形變、火山運(yùn)動(dòng)、冰川漂移、地面沉降以及山體滑坡等。水準(zhǔn)測(cè)量是通過(guò)布設(shè)水準(zhǔn)路線，定期高精度地測(cè)定若干水準(zhǔn)點(diǎn)，以此獲取該區(qū)域沉降變化規(guī)律的一種常用作業(yè)手段分層標(biāo)是通過(guò)對(duì)各含水層地下水位變化、土層變形動(dòng)態(tài)規(guī)律的分析，確定影響地面沉降的主要因素及其權(quán)重，提出具體的地面沉降防治方案。s production and life seriously affected. Based on land subsidence damage to the environment and the gravity of economic development, the current number of regions and cities have land subsidence monitoring harm reduction as an important prevention work. In recent years, extensive research and practice shows that, InSAR (Interferometric Synthetic Aperture Radar) technology to highprecision monitoring of the surface of the small surface deformation, allweather, alltime, wide coverage, high degree of automation and high precision ground deformation monitoring of the advantages of effective to make up for leveling and GPS measurements the problem of inadequate spatial resolution, is the subsidence is an important side lesson plan techniques. However, InSAR monitoring tools also have their own deficiencies, InSAR technology, while high spatial resolution, but the operation of the radar satellite because of its inherent cycle can not meet monitoring requirements. In addition, the time lost coherence or mountain top and bottom inverted SAR image of the interferogram so easily lead to data loss occurs, and thus difficult to achieve regional monitoring purposes. In this paper, the phenomenon of missing data, research is based on polynomial fitting and multifaceted function fitting of the missing data fitting method, focusing on the multifaceted function fitting the selection of nodes, and finally to ground subsidence monitoring data of Xi39。┊┊┊┊┊┊┊┊┊┊┊┊┊裝┊┊┊┊┊訂┊┊┊┊┊線┊┊┊┊┊┊┊┊┊┊┊┊┊二 ○ 一 ○屆 畢 業(yè) 論 文InSAR數(shù)據(jù)后處理研究— 缺失數(shù)據(jù)擬合學(xué) 院: 地質(zhì)工程與測(cè)繪專(zhuān) 業(yè)：地理信息系統(tǒng)姓 名：XXX學(xué) 號(hào)：XXX指導(dǎo)教師：XXX完成時(shí)間：二〇一〇年六月共 32 頁(yè) 第 2 頁(yè)摘 要我國(guó)的地面沉降現(xiàn)象一直比較嚴(yán)重。關(guān)鍵詞：地面沉降，InSAR，DInSAR，多項(xiàng)式擬合，多面函數(shù)擬合ABSTRACTSerious land subsidence phenomenon happen in china for a long time. Subsidence damage to buildings and production facilities, damage to underground line facilities to the peo