一種面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法及裝置制造方法
【專利摘要】本發(fā)明提供一種面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法及裝置，該方法包括：將一區(qū)域雷達(dá)影像的所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)及其頂點進行編號；對每個子網(wǎng)的四個邊界點，查找在預(yù)置初始步長L范圍內(nèi)的所有點目標(biāo)及其所屬的子網(wǎng)；計算所有點目標(biāo)與當(dāng)前子網(wǎng)邊界點的距離，并將距離由小到大排序，依次計算得到的每條邊的最大模型相干系數(shù)，大于給定閾值時即認(rèn)為該邊為有效邊，否則繼續(xù)計算下一條邊；當(dāng)遍歷完所有的子網(wǎng)后，重新統(tǒng)計子網(wǎng)數(shù)目，若該子網(wǎng)數(shù)目只有1個，則無需擴展連接計算；若該子網(wǎng)數(shù)目超過1，再依次以2L，3L，4L，5L，……nL為步長，重復(fù)上述步驟進行擴展連接，直至循環(huán)停止。其能逐步減少子網(wǎng)總數(shù)，實現(xiàn)不連續(xù)子網(wǎng)的快速連接，有利于大區(qū)域的地表形變監(jiān)測。
【專利說明】-種面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法及裝置

【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及InSAR(SyntheticApertureRadarInterferometry，合成孔徑雷達(dá) 干涉測量）【技術(shù)領(lǐng)域】，尤其涉及一種面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法及裝置。

【背景技術(shù)】
[0002] 合成孔徑雷達(dá)（SAR，SyntheticApertureRadar)是20世紀(jì)50年代發(fā)展起來的最 重要的對地觀測技術(shù)。但是SAR技術(shù)只能獲取地表目標(biāo)物的二維信息，缺乏獲取目標(biāo)點高 程信息和監(jiān)測目標(biāo)微小形變的能力。將干涉測量技術(shù)與SAR技術(shù)結(jié)合而形成的合成孔徑雷 達(dá)干涉測量技術(shù)（InSAR，SyntheticApertureRadarInterferometry)提供了獲取地面三 維信息的全新方法，它通過兩副天線同時觀測或通過一副天線兩次平行觀測，獲取地面同 一景觀的復(fù)圖像對，根據(jù)地面各點在兩幅復(fù)圖像中的相位差，得出各點在兩次成像中微波 的路程差，從而獲得地面目標(biāo)的高程信息。基于InSAR技術(shù)的發(fā)展，合成孔徑雷達(dá)差分干涉 測量（DifferentialSyntheticApertureRadarInterferometry,DInSAR)技術(shù)是對兩幅 以上的干涉圖或?qū)σ环缮鎴D加一幅地面數(shù)字高程模型（DEM,DigitalElevationModel) 進行再處理的一種技術(shù)，它通過去除地形引起的干涉相位，獲得關(guān)于地表形變的信息，在火 山監(jiān)測、地震位移測量、地面沉降等領(lǐng)域具有極大的應(yīng)用潛力。
[0003] 基于重復(fù)軌道的DInSAR技術(shù)容易受空間失相干、時間失相干和大氣干擾等因素 的影響，難以進行常態(tài)化的實際應(yīng)用。引起空間失相干的原因包括大的垂直基線和大的形 變梯度。時間失相干則是由于兩次成像時刻環(huán)境的不一致導(dǎo)致同一地面像元中散射體的散 射特性發(fā)生改變而引起的。易變的大氣條件可能會導(dǎo)致在兩幅干涉影像上不一致的相位延 遲，從而引起形變測量誤差。
[0004] 為了克服傳統(tǒng)DInSAR技術(shù)的這些限制，自上世紀(jì)90年代末，一些新的InSAR處理 技術(shù)被提出。這些技術(shù)的共同特點是：基于時間序列SAR影像進行處理，處理的對象不是整 幅影像的全部像元，而是其中具有穩(wěn)定散射特性從而能在較長時間間隔內(nèi)保持高相干的像 元子集，也就是高相干點。這些技術(shù)總體上可以概括為兩類：以永久散射體干涉（Permanent Scatterer或者PersistentScattererInterferometry，或PS-InSAR)為代表的單一主影 像時間序列InSAR技術(shù)和以小基線集技術(shù)（Smallbaselinesubsetinterferometry，或 SBASInSAR)為代表的多主影像時間序列InSAR技術(shù)。為敘述方便，將這兩種技術(shù)統(tǒng)稱為時 間序列InSAR技術(shù)。時間序列InSAR技術(shù)對上述三種制約因素均有良好的免疫力，目前已 經(jīng)取代傳統(tǒng)的DInSAR技術(shù)在火山、地震、滑坡、地面沉降等領(lǐng)域得到大量應(yīng)用。由于城市地 區(qū)擁有密集的天然點目標(biāo)，如路燈、屋頂?shù)龋虼藭r間序列InSAR技術(shù)在城市地面沉降測量 方面應(yīng)用最為廣泛。
[0005] 在時間序列InSAR技術(shù)中，當(dāng)提取完所有點目標(biāo)后，通常利用Delaunay三角網(wǎng)連 接所有的點目標(biāo)，通過相鄰點相位差分并建立線性相位模型求解形變參數(shù)。但受大氣相關(guān) 距離限制以及模型相干系數(shù)閾值影響，Delaunay三角網(wǎng)常常不能完整地連接所有點目標(biāo)。 當(dāng)點目標(biāo)質(zhì)量好且密度較大時，如城區(qū)，通過Delaunay建網(wǎng)能夠較好的連接各點目標(biāo)，實 現(xiàn)形變參數(shù)的解算；但是當(dāng)點目標(biāo)質(zhì)量較差或相隔較遠(yuǎn)時，如非城區(qū)，容易出現(xiàn)大量孤立的 子網(wǎng)，這時就無法解算形變信息。
[0006] 對于上述問題，目前研究人員提出了兩種解決方案，一種是點目標(biāo)自由連接方法， 將一定距離內(nèi)的所有點兩兩連接，從而極大增加了三角網(wǎng)的邊數(shù)；一種是復(fù)雜網(wǎng)絡(luò)方法，它 首先利用Delaunay三角網(wǎng)對所有點目標(biāo)快速建網(wǎng)，提取城區(qū)密集點目標(biāo)的形變參數(shù)，然后 在大氣相關(guān)距離范圍內(nèi)連接所有沒有被估計的點目標(biāo)，建立復(fù)雜網(wǎng)絡(luò)，重新計算形變參數(shù)。 上述兩種方法均能提高點目標(biāo)之間的連接關(guān)系，但顯而易見的是，它們所帶來的計算量非 常巨大，不適合大區(qū)域海量數(shù)據(jù)處理。如何快速且穩(wěn)健地連接點目標(biāo)子網(wǎng)，以實現(xiàn)大區(qū)域地 表形變監(jiān)測，這是本領(lǐng)域技術(shù)人員亟待解決的關(guān)鍵問題。


【發(fā)明內(nèi)容】

[0007] 本發(fā)明實施例提供一種面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法及裝置，以解 決不連續(xù)的Delaunay三角網(wǎng)穩(wěn)健連接問題，從而更好的實現(xiàn)大區(qū)域地表形變監(jiān)測。
[0008] -方面，本發(fā)明實施例提供了一種面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法， 所述面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法包括：
[0009] 將一區(qū)域雷達(dá)影像的所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)及所述所有不連續(xù)的 Delaunay三角網(wǎng)子網(wǎng)頂點進行編號；
[0010] 對所述所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)中的每個子網(wǎng)的四個邊界點，查找在 預(yù)置初始步長L范圍內(nèi)的所有點目標(biāo)及所述所有點目標(biāo)所屬的子網(wǎng)；
[0011] 計算所述所有點目標(biāo)與當(dāng)前子網(wǎng)邊界點的距離，并將所述距離由小到大排序，依 次計算每條邊的最大模型相干系數(shù)，當(dāng)該最大模型相干系數(shù)大于給定閾值時即認(rèn)為該邊為 將該點目標(biāo)所屬的子網(wǎng)與當(dāng)前子網(wǎng)連接的有效邊，否則繼續(xù)計算下一條邊；
[0012] 當(dāng)遍歷完所有的子網(wǎng)后，重新統(tǒng)計子網(wǎng)數(shù)目，若該子網(wǎng)數(shù)目只有1個，說明全部子 網(wǎng)均已連接成功，則無需擴展連接計算；若該子網(wǎng)數(shù)目超過1，說明還有未被連接的孤立子 網(wǎng)，再依次以2L，3L，4L，5L，……nL為步長，重復(fù)上述步驟進行擴展連接，直至循環(huán)停止，此 時所有子網(wǎng)均已完成擴展連接；若還存在若干與主要網(wǎng)絡(luò)相隔nL以上的子網(wǎng)，將所述與主 要網(wǎng)絡(luò)相隔nL以上的子網(wǎng)拋棄，其中η為大于1的自然數(shù)。
[0013] 本發(fā)明實施例通過這種多層級擴展方式，能夠逐步減少子網(wǎng)總數(shù)，從而實現(xiàn)不連 續(xù)子網(wǎng)的快速連接，特別有利于大區(qū)域的地表形變監(jiān)測。
[0014] 優(yōu)選的，在本發(fā)明一實施例中，所述將一區(qū)域雷達(dá)影像的所有不連續(xù)的Delaunay 三角網(wǎng)子網(wǎng)及所述所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)頂點進行編號，包括：將所述所有不 連續(xù)的Delaunay三角網(wǎng)子網(wǎng)編號按照點目標(biāo)的順序從小到大，依次進行編號，分別以子網(wǎng) 1、子網(wǎng)2、子網(wǎng)3.......標(biāo)識，直至完成所有點目標(biāo)的遍歷；同時記錄每個子網(wǎng)中含有的點 目標(biāo)及該點目標(biāo)編號。本發(fā)明實施例通過上述編號方式，進一步實現(xiàn)了不連續(xù)子網(wǎng)的快速 連接，特別有利于大區(qū)域的地表形變監(jiān)測。
[0015] 優(yōu)選的，在本發(fā)明一實施例中，所述預(yù)置初始步長L根據(jù)SAR影像的分辨率和所 有點目標(biāo)的密度進行調(diào)整，但最大步長nL小于或等于大氣相關(guān)距離，即3000m;所述預(yù)置初 始步長L可設(shè)為分辨率的50倍：如當(dāng)SAR像兀分辨率為IOm時，該預(yù)置初始步長L可設(shè)為 500m;當(dāng)SAR像元分辨率為20m時，該預(yù)置初始步長L可設(shè)為IOOOm;當(dāng)分辨率為3m時，該 預(yù)置初始步長L可設(shè)為150m。本發(fā)明實施例通過上述預(yù)置初始步長L的設(shè)置方式，通過這 種多層級擴展方式，進一步能夠逐步減少子網(wǎng)總數(shù)，從而實現(xiàn)不連續(xù)子網(wǎng)的快速連接，特別 有利于大區(qū)域的地表形變監(jiān)測。
[0016] 優(yōu)選的，在本發(fā)明一實施例中，所述最大模型相干系數(shù)，反映了相對形變速率和相 對高程誤差對相鄰點目標(biāo)的擬合程度，其最大值為1，最小值為0 ;將所述最大模型相干系 數(shù)的給定閾值設(shè)為〇. 7,當(dāng)某條邊的最大模型相干系數(shù)大于0. 7時，認(rèn)為該條邊是可靠的連 接關(guān)系，將兩個子網(wǎng)相連；否則該條邊不能作為兩個子網(wǎng)的連接邊，需要繼續(xù)計算下一條邊 的最大模型相干系數(shù)，直至完成所有邊循環(huán)。本發(fā)明實施例通過上述最大模型相干系數(shù)的 設(shè)置方式，進一步實現(xiàn)不連續(xù)子網(wǎng)的快速連接，特別有利于大區(qū)域的地表形變監(jiān)測。
[0017] 優(yōu)選的，在本發(fā)明一實施例中，當(dāng)預(yù)置初始步長設(shè)為L時，所述當(dāng)遍歷完所有的 子網(wǎng)后，重新統(tǒng)計子網(wǎng)數(shù)目，若該子網(wǎng)數(shù)目只有1個，說明全部子網(wǎng)均已連接成功，則無需 擴展連接計算；若該子網(wǎng)數(shù)目超過1，說明還有未被連接的孤立子網(wǎng)，再依次以2L，3L，4L， 5L，……nL為步長，重復(fù)上述步驟進行擴展連接，直至循環(huán)停止，此時所有子網(wǎng)均已連接 成功；若還存在若干與主要網(wǎng)絡(luò)相隔nL以上的子網(wǎng)，將所述與主要網(wǎng)絡(luò)相隔nL以上的子 網(wǎng)拋棄，包括：在完成L內(nèi)的所有子網(wǎng)搜索后，重新統(tǒng)計不連續(xù)子網(wǎng)數(shù)目，并對所述不連續(xù) 子網(wǎng)頂點進行編號；若子網(wǎng)數(shù)目只有1個，說明全部子網(wǎng)均已連接成功，則無需擴展連接計 算；若子網(wǎng)數(shù)目超過1，再以2L為步長，重復(fù)上述步驟，對每個子網(wǎng)進行相鄰子網(wǎng)搜索；以此 類推，當(dāng)新的子網(wǎng)數(shù)目仍然超過1時，再進一步擴大搜索步長距離，逐漸以3L，4L，5L，…… 3000m為步長循環(huán)計算，進行子網(wǎng)鏈接，直至循環(huán)停止，此時所有子網(wǎng)均已完成擴展連接； 若還存在若干與主要網(wǎng)絡(luò)相隔3000m以上的子網(wǎng)，將所述與主要網(wǎng)絡(luò)相隔3000m以上的子 網(wǎng)拋棄。本發(fā)明實施例通過上述擴展連接計算的方式，通過這種多層級擴展方式，能夠逐步 減少子網(wǎng)總數(shù)，從而實現(xiàn)不連續(xù)子網(wǎng)的快速連接，特別有利于大區(qū)域的地表形變監(jiān)測。
[0018] 另一方面，本發(fā)明實施例提供了一種面向時間序列InSAR的不連續(xù)子網(wǎng)連接裝 置，所述面向時間序列InSAR的不連續(xù)子網(wǎng)連接裝置包括：
[0019] 編號單元，用于將一區(qū)域雷達(dá)影像的所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)及所述 所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)頂點進行編號；
[0020] 查找單元，用于對所述所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)中的每個子網(wǎng)的四個 邊界點，查找在預(yù)置初始步長L范圍內(nèi)的所有點目標(biāo)及所述所有點目標(biāo)所屬的子網(wǎng)；
[0021] 計算單元，用于計算所述所有點目標(biāo)與當(dāng)前子網(wǎng)邊界點的距離，并將所述距離由 小到大排序，依次計算每條邊的最大模型相干系數(shù)，當(dāng)該最大模型相干系數(shù)大于給定閾值 時即認(rèn)為該邊為將該點目標(biāo)所屬的子網(wǎng)與當(dāng)前子網(wǎng)連接的有效邊，否則繼續(xù)計算下一條 邊；
[0022] 遍歷單元，用于當(dāng)遍歷完所有的子網(wǎng)后，重新統(tǒng)計子網(wǎng)數(shù)目，若該子網(wǎng)數(shù)目只有1 個，說明全部子網(wǎng)均已連接成功，則無需擴展連接計算；若該子網(wǎng)數(shù)目超過1，說明還有未 被連接的孤立子網(wǎng)，再依次以2L，3L，4L，5L，……nL為步長，重復(fù)上述步驟進行擴展連接， 直至循環(huán)停止，此時所有子網(wǎng)均已完成擴展連接；若還存在若干與主要網(wǎng)絡(luò)相隔nL以上的 子網(wǎng)，將所述與主要網(wǎng)絡(luò)相隔nL以上的子網(wǎng)拋棄，其中η為大于1的自然數(shù)。
[0023] 本發(fā)明實施例通過這種多層級擴展方式，能夠逐步減少子網(wǎng)總數(shù)，從而實現(xiàn)不連 續(xù)子網(wǎng)的快速連接，特別有利于大區(qū)域的地表形變監(jiān)測。
[0024] 優(yōu)選的，在本發(fā)明一實施例中，所述編號單元，進一步用于將所述所有不連續(xù)的 Delaunay三角網(wǎng)子網(wǎng)編號按照點目標(biāo)的順序從小到大，依次進行編號，分別以子網(wǎng)1、子網(wǎng) 2、子網(wǎng)3.......標(biāo)識，直至完成所有點目標(biāo)的遍歷；同時記錄每個子網(wǎng)中含有的點目標(biāo)及 該點目標(biāo)編號。本發(fā)明實施例通過上述編號方式，進一步實現(xiàn)了不連續(xù)子網(wǎng)的快速連接，特 別有利于大區(qū)域的地表形變監(jiān)測。
[0025] 優(yōu)選的，在本發(fā)明一實施例中，所述預(yù)置初始步長L根據(jù)SAR影像的分辨率和所 有點目標(biāo)的密度進行調(diào)整，但最大步長nL小于或等于大氣相關(guān)距離，即3000m;所述預(yù)置初 始步長L可設(shè)為分辨率的50倍：如當(dāng)SAR像兀分辨率為IOm時，該預(yù)置初始步長L可設(shè)為 500m;當(dāng)SAR像元分辨率為20m時，該預(yù)置初始步長L可設(shè)為IOOOm;當(dāng)分辨率為3m時，該 預(yù)置初始步長L可設(shè)為150m。本發(fā)明實施例通過上述預(yù)置初始步長L的設(shè)置方式，通過這 種多層級擴展方式，能夠逐步減少子網(wǎng)總數(shù)，從而實現(xiàn)不連續(xù)子網(wǎng)的快速連接，特別有利于 大區(qū)域的地表形變監(jiān)測。
[0026] 優(yōu)選的，在本發(fā)明一實施例中，所述最大模型相干系數(shù)，反映了相對形變速率和相 對高程誤差對相鄰點目標(biāo)的擬合程度，其最大值為1，最小值為〇 ;將所述最大模型相干系 數(shù)的給定閾值設(shè)為〇. 7,當(dāng)某條邊的最大模型相干系數(shù)大于0. 7時，認(rèn)為該條邊是可靠的連 接關(guān)系，將兩個子網(wǎng)相連；否則該條邊不能作為兩個子網(wǎng)的連接邊，需要繼續(xù)計算下一條邊 的最大模型相干系數(shù)，直至完成所有邊循環(huán)。本發(fā)明實施例通過上述最大模型相干系數(shù)的 設(shè)置方式，進一步實現(xiàn)不連續(xù)子網(wǎng)的快速連接，特別有利于大區(qū)域的地表形變監(jiān)測。
[0027] 優(yōu)選的，在本發(fā)明一實施例中，當(dāng)預(yù)置初始步長設(shè)為L時，所述遍歷單元，進一步 用于在完成L內(nèi)的所有子網(wǎng)搜索后，重新統(tǒng)計不連續(xù)子網(wǎng)數(shù)目，并對所述不連續(xù)子網(wǎng)頂點 進行編號；若子網(wǎng)數(shù)目只有1個，說明全部子網(wǎng)均已連接成功，則無需擴展連接計算；若子 網(wǎng)數(shù)目超過1，再以2L為步長，重復(fù)上述步驟，對每個子網(wǎng)進行相鄰子網(wǎng)搜索；以此類推，當(dāng) 新的子網(wǎng)數(shù)目仍然超過1時，再進一步擴大搜索步長距離，逐漸以3L，4L，5L，……3000m為 步長循環(huán)計算，進行子網(wǎng)鏈接，直至循環(huán)停止，此時所有子網(wǎng)均已連接成功；若還存在若干 與主要網(wǎng)絡(luò)相隔3000m以上的子網(wǎng)，將所述與主要網(wǎng)絡(luò)相隔3000m以上的子網(wǎng)拋棄。本發(fā)明 實施例通過上述擴展連接計算的方式，通過這種多層級擴展方式，能夠逐步減少子網(wǎng)總數(shù)， 從而實現(xiàn)不連續(xù)子網(wǎng)的快速連接，特別有利于大區(qū)域的地表形變監(jiān)測。
[0028] 上述技術(shù)方案具有如下有益效果：本發(fā)明實施例所述面向時間序列InSAR的不連 續(xù)子網(wǎng)連接方法包括：將一區(qū)域雷達(dá)影像的所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)及所述所 有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)頂點進行編號；對所述所有不連續(xù)的Delaunay三角網(wǎng)子 網(wǎng)中的每個子網(wǎng)的四個邊界點，查找在預(yù)置初始步長L范圍內(nèi)的所有點目標(biāo)及所述所有點 目標(biāo)所屬的子網(wǎng)；計算所述所有點目標(biāo)與當(dāng)前子網(wǎng)邊界點的距離，并將所述距離由小到大 排序，依次計算每條邊的最大模型相干系數(shù)，當(dāng)該最大模型相干系數(shù)大于給定閾值時即認(rèn) 為該邊為將該點目標(biāo)所屬的子網(wǎng)與當(dāng)前子網(wǎng)連接的有效邊，否則繼續(xù)計算下一條邊；當(dāng)遍 歷完所有的子網(wǎng)后，重新統(tǒng)計子網(wǎng)數(shù)目，若該子網(wǎng)數(shù)目只有1個，說明全部子網(wǎng)均已連接成 功，則無需擴展連接計算；若該子網(wǎng)數(shù)目超過1，說明還有未被連接的孤立子網(wǎng)，再依次以 2L，3L，4L，5L，……nL為步長，重復(fù)上述步驟進行擴展連接，直至循環(huán)停止，此時所有子網(wǎng) 均已連接成功；若還存在若干與主要網(wǎng)絡(luò)相隔nL以上的子網(wǎng)，將所述與主要網(wǎng)絡(luò)相隔nL以 上的子網(wǎng)拋棄，其中η為大于1的自然數(shù)。本發(fā)明實施例達(dá)到了如下技術(shù)效果：通過這種多 層級擴展方式，能夠逐步減少子網(wǎng)總數(shù)，從而實現(xiàn)不連續(xù)子網(wǎng)的快速連接，特別有利于大區(qū) 域的地表形變監(jiān)測。

【專利附圖】

【附圖說明】
[0029] 為了更清楚地說明本發(fā)明實施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對實施例或現(xiàn) 有技術(shù)描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本 發(fā)明的一些實施例，對于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動的前提下，還可以 根據(jù)這些附圖獲得其他的附圖。
[0030] 圖1為本發(fā)明實施例一種面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法流程圖；
[0031] 圖2為本發(fā)明實施例一種面向時間序列InSAR的不連續(xù)子網(wǎng)連接裝置結(jié)構(gòu)示意 圖；
[0032] 圖3是本發(fā)明應(yīng)用實例使用的河北省某試驗區(qū)24景歐洲空間局ENVISATASAR影 像的平均幅度圖：
[0033] 圖4是本發(fā)明應(yīng)用實例一種面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法流程圖；
[0034] 圖5為本發(fā)明應(yīng)用實例使用的試驗區(qū)不連續(xù)Delaunay三角網(wǎng)子網(wǎng)示意圖；
[0035] 圖6是本發(fā)明應(yīng)用實例以500米為步長，最大搜索距離為3000米，利用多層級快 速連接方法得到的網(wǎng)絡(luò)示意圖；
[0036] 圖7為本發(fā)明應(yīng)用實例基于多層級快速連接方法得到的5179個點目標(biāo)的線性地 表形變速率示意圖（單位：毫米/年）；
[0037] 圖8是利用現(xiàn)有技術(shù)復(fù)雜網(wǎng)絡(luò)連接方法得到的網(wǎng)絡(luò)示意圖；
[0038] 圖9是現(xiàn)有技術(shù)基于復(fù)雜網(wǎng)絡(luò)連接方法得到的5141個點目標(biāo)的線性地表形變速 率示意圖（單位：毫米/年）。

【具體實施方式】
[0039] 下面將結(jié)合本發(fā)明實施例中的附圖，對本發(fā)明實施例中的技術(shù)方案進行清楚、完 整地描述，顯然，所描述的實施例僅僅是本發(fā)明一部分實施例，而不是全部的實施例?；?本發(fā)明中的實施例，本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動前提下所獲得的所有其他 實施例，都屬于本發(fā)明保護的范圍。
[0040] 如圖1所示，為本發(fā)明實施例一種面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法流 程圖，所述面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法包括：
[0041] 101、將一區(qū)域雷達(dá)影像的所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)及所述所有不連續(xù) 的Delaunay三角網(wǎng)子網(wǎng)頂點進行編號；
[0042] 102、對所述所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)中的每個子網(wǎng)的四個邊界點，查 找在預(yù)置初始步長L范圍內(nèi)的所有點目標(biāo)及所述所有點目標(biāo)所屬的子網(wǎng)；
[0043] 103、計算所述所有點目標(biāo)與當(dāng)前子網(wǎng)邊界點的距離，并將所述距離由小到大排 序，依次計算每條邊的最大模型相干系數(shù)，當(dāng)該最大模型相干系數(shù)大于給定閾值時即認(rèn)為 該邊為將該點目標(biāo)所屬的子網(wǎng)與當(dāng)前子網(wǎng)連接的有效邊，否則繼續(xù)計算下一條邊；
[0044] 104、當(dāng)遍歷完所有的子網(wǎng)后，重新統(tǒng)計子網(wǎng)數(shù)目，若該子網(wǎng)數(shù)目只有1個，說明全 部子網(wǎng)均已連接成功，則無需擴展連接計算；若該子網(wǎng)數(shù)目超過1，說明還有未被連接的孤 立子網(wǎng)，再依次以2L，3L，4L，5L，……nL為步長，重復(fù)上述步驟進行擴展連接，直至循環(huán)停 止，此時所有子網(wǎng)均已連接成功；若還存在若干與主要網(wǎng)絡(luò)相隔nL以上的子網(wǎng)，將所述與 主要網(wǎng)絡(luò)相隔nL以上的子網(wǎng)拋棄，其中η為大于1的自然數(shù)。
[0045] 本發(fā)明實施例通過這種多層級擴展方式，能夠逐步減少子網(wǎng)總數(shù)，從而實現(xiàn)不連 續(xù)子網(wǎng)的快速連接，特別有利于大區(qū)域的地表形變監(jiān)測。
[0046] 優(yōu)選的，所述將一區(qū)域雷達(dá)影像的所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)及所述 所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)頂點進行編號，包括：將所述所有不連續(xù)的Delaunay 三角網(wǎng)子網(wǎng)編號按照點目標(biāo)的順序從小到大，依次進行編號，分別以子網(wǎng)1、子網(wǎng)2、子網(wǎng) 3.......標(biāo)識，直至完成所有點目標(biāo)的遍歷；同時記錄每個子網(wǎng)中含有的點目標(biāo)及該點目 標(biāo)編號。本發(fā)明實施例通過上述編號方式，進一步實現(xiàn)了不連續(xù)子網(wǎng)的快速連接，特別有利 于大區(qū)域的地表形變監(jiān)測。
[0047] 優(yōu)選的，所述預(yù)置初始步長L根據(jù)SAR影像的分辨率和所有點目標(biāo)的密度進行調(diào) 整，但最大步長nL小于或等于大氣相關(guān)距離，即3000m;所述預(yù)置初始步長L可設(shè)為分辨率 的50倍：如當(dāng)SAR像元分辨率為IOm時，該預(yù)置初始步長L可設(shè)為500m;當(dāng)SAR像元分辨 率為20m時,該預(yù)置初始步長L可設(shè)為IOOOm;當(dāng)分辨率為3m時,該預(yù)置初始步長L可設(shè)為 150m。本發(fā)明實施例通過上述預(yù)置初始步長L的設(shè)置方式，通過這種多層級擴展方式，能夠 逐步減少子網(wǎng)總數(shù)，從而實現(xiàn)不連續(xù)子網(wǎng)的快速連接，特別有利于大區(qū)域的地表形變監(jiān)測。 [0048] 優(yōu)選的，所述最大模型相干系數(shù)，反映了相對形變速率和相對高程誤差對相鄰點 目標(biāo)的擬合程度，其最大值為1，最小值為〇;將所述最大模型相干系數(shù)的給定閾值設(shè)為 0. 7,當(dāng)某條邊的最大模型相干系數(shù)大于0. 7時，認(rèn)為該條邊是可靠的連接關(guān)系，將兩個子 網(wǎng)相連；否則該條邊不能作為兩個子網(wǎng)的連接邊，需要繼續(xù)計算下一條邊的最大模型相干 系數(shù)，直至完成所有邊循環(huán)。本發(fā)明實施例通過上述最大模型相干系數(shù)的設(shè)置方式，進一步 實現(xiàn)不連續(xù)子網(wǎng)的快速連接，特別有利于大區(qū)域的地表形變監(jiān)測。
[0049] 優(yōu)選的，當(dāng)預(yù)置初始步長設(shè)為L時，所述當(dāng)遍歷完所有的子網(wǎng)后，重新統(tǒng)計子網(wǎng)數(shù) 目，若該子網(wǎng)數(shù)目只有1個，說明全部子網(wǎng)均已連接成功，則無需擴展連接計算；若該子網(wǎng) 數(shù)目超過1，說明還有未被連接的孤立子網(wǎng)，再依次以2L，3L，4L，5L，……nL為步長，重復(fù) 上述步驟進行擴展連接，直至循環(huán)停止，此時所有子網(wǎng)均已連接成功；若還存在若干與主要 網(wǎng)絡(luò)相隔nL以上的子網(wǎng)，將所述與主要網(wǎng)絡(luò)相隔nL以上的子網(wǎng)拋棄，包括：在完成L內(nèi)的 所有子網(wǎng)搜索后，重新統(tǒng)計不連續(xù)子網(wǎng)數(shù)目，并對所述不連續(xù)子網(wǎng)頂點進行編號；若子網(wǎng)數(shù) 目只有1個，說明全部子網(wǎng)均已連接成功，則無需擴展連接計算；若子網(wǎng)數(shù)目超過1，再以2L 為步長，重復(fù)上述步驟，對每個子網(wǎng)進行相鄰子網(wǎng)搜索；以此類推，當(dāng)新的子網(wǎng)數(shù)目仍然超 過1時，再進一步擴大搜索步長距離，逐漸以3L，4L，5L，……3000m為步長循環(huán)計算，進行子 網(wǎng)鏈接，直至循環(huán)停止，此時所有子網(wǎng)均已連接成功；若還存在若干與主要網(wǎng)絡(luò)相隔3000m 以上的子網(wǎng)，將所述與主要網(wǎng)絡(luò)相隔3000m以上的子網(wǎng)拋棄。本發(fā)明實施例通過上述擴展 連接計算的方式，通過這種多層級擴展方式，能夠逐步減少子網(wǎng)總數(shù)，從而實現(xiàn)不連續(xù)子網(wǎng) 的快速連接，特別有利于大區(qū)域的地表形變監(jiān)測。
[0050] 如圖2所示，為本發(fā)明實施例一種面向時間序列InSAR的不連續(xù)子網(wǎng)連接裝置結(jié) 構(gòu)示意圖，所述面向時間序列InSAR的不連續(xù)子網(wǎng)連接裝置包括：
[0051] 編號單元21，用于將一區(qū)域雷達(dá)影像的所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)及所 述所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)頂點進行編號；
[0052] 查找單元22,用于對所述所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)中的每個子網(wǎng)的四 個邊界點，查找在預(yù)置初始步長L范圍內(nèi)的所有點目標(biāo)及所述所有點目標(biāo)所屬的子網(wǎng)；
[0053] 計算單元23,用于計算所述所有點目標(biāo)與當(dāng)前子網(wǎng)邊界點的距離，并將所述距離 由小到大排序，依次計算每條邊的最大模型相干系數(shù)，當(dāng)該最大模型相干系數(shù)大于給定閾 值時即認(rèn)為該邊為將該點目標(biāo)所屬的子網(wǎng)與當(dāng)前子網(wǎng)連接的有效邊，否則繼續(xù)計算下一條 邊；
[0054] 遍歷單元24,用于當(dāng)遍歷完所有的子網(wǎng)后，重新統(tǒng)計子網(wǎng)數(shù)目，若該子網(wǎng)數(shù)目只有 1個，說明全部子網(wǎng)均已連接成功，則無需擴展連接計算；若該子網(wǎng)數(shù)目超過1，說明還有未 被連接的孤立子網(wǎng)，再依次以2L，3L，4L，5L，……nL為步長，重復(fù)上述步驟進行擴展連接，直 至循環(huán)停止，此時所有子網(wǎng)均已連接成功；若還存在若干與主要網(wǎng)絡(luò)相隔nL以上的子網(wǎng)， 將所述與主要網(wǎng)絡(luò)相隔nL以上的子網(wǎng)拋棄，其中η為大于1的自然數(shù)。
[0055] 本發(fā)明實施例通過這種多層級擴展方式，能夠逐步減少子網(wǎng)總數(shù)，從而實現(xiàn)不連 續(xù)子網(wǎng)的快速連接，特別有利于大區(qū)域的地表形變監(jiān)測。
[0056] 優(yōu)選的，所述編號單元21，進一步用于將所述所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng) 編號按照點目標(biāo)的順序從小到大，依次進行編號，分別以子網(wǎng)1、子網(wǎng)2、子網(wǎng)3.......標(biāo) 識，直至完成所有點目標(biāo)的遍歷；同時記錄每個子網(wǎng)中含有的點目標(biāo)及該點目標(biāo)編號。本發(fā) 明實施例通過上述編號方式，進一步實現(xiàn)了不連續(xù)子網(wǎng)的快速連接，特別有利于大區(qū)域的 地表形變監(jiān)測。
[0057] 優(yōu)選的，所述預(yù)置初始步長L根據(jù)SAR影像的分辨率和所有點目標(biāo)的密度進行調(diào) 整，但最大步長nL小于或等于大氣相關(guān)距離，即3000m;所述預(yù)置初始步長L可設(shè)為分辨率 的50倍：如當(dāng)SAR像元分辨率為IOm時，該預(yù)置初始步長L可設(shè)為500m;當(dāng)SAR像元分辨 率為20m時,該預(yù)置初始步長L可設(shè)為IOOOm;當(dāng)分辨率為3m時,該預(yù)置初始步長L可設(shè)為 150m。發(fā)明實施例通過上述預(yù)置初始步長L的設(shè)置方式，通過這種多層級擴展方式，能夠逐 步減少子網(wǎng)總數(shù)，從而實現(xiàn)不連續(xù)子網(wǎng)的快速連接，特別有利于大區(qū)域的地表形變監(jiān)測。
[0058] 優(yōu)選的，所述最大模型相干系數(shù)，反映了相對形變速率和相對高程誤差對相鄰點 目標(biāo)的擬合程度，其最大值為1，最小值為〇;將所述最大模型相干系數(shù)的給定閾值設(shè)為 0. 7,當(dāng)某條邊的最大模型相干系數(shù)大于0. 7時，認(rèn)為該條邊是可靠的連接關(guān)系，將兩個子 網(wǎng)相連；否則該條邊不能作為兩個子網(wǎng)的連接邊，需要繼續(xù)計算下一條邊的最大模型相干 系數(shù)，直至完成所有邊循環(huán)。本發(fā)明實施例通過上述最大模型相干系數(shù)的設(shè)置方式，進一步 實現(xiàn)不連續(xù)子網(wǎng)的快速連接，特別有利于大區(qū)域的地表形變監(jiān)測。
[0059] 優(yōu)選的，當(dāng)預(yù)置初始步長設(shè)為L時，所述遍歷單元24,進一步用于在完成L內(nèi)的所 有子網(wǎng)搜索后，重新統(tǒng)計不連續(xù)子網(wǎng)數(shù)目，并對所述不連續(xù)子網(wǎng)頂點進行編號；若子網(wǎng)數(shù)目 只有1個，說明全部子網(wǎng)均已連接成功，則無需擴展連接計算；若子網(wǎng)數(shù)目超過1，再以2L 為步長，重復(fù)上述步驟，對每個子網(wǎng)進行相鄰子網(wǎng)搜索；以此類推，當(dāng)新的子網(wǎng)數(shù)目仍然超 過1時，再進一步擴大搜索步長距離，逐漸以3L，4L，5L，……3000m為步長循環(huán)計算，進行子 網(wǎng)鏈接，直至循環(huán)停止，此時所有子網(wǎng)均已連接成功；若還存在若干與主要網(wǎng)絡(luò)相隔3000m 以上的子網(wǎng)，將所述與主要網(wǎng)絡(luò)相隔3000m以上的子網(wǎng)拋棄。本發(fā)明實施例通過上述擴展 連接計算的方式，通過這種多層級擴展方式，能夠逐步減少子網(wǎng)總數(shù)，從而實現(xiàn)不連續(xù)子網(wǎng) 的快速連接，特別有利于大區(qū)域的地表形變監(jiān)測。
[0060] 上述技術(shù)方案具有如下有益效果：本發(fā)明實施例所述面向時間序列InSAR的不連 續(xù)子網(wǎng)連接方法包括：將一區(qū)域雷達(dá)影像的所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)及所述所 有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)頂點進行編號；對所述所有不連續(xù)的Delaunay三角網(wǎng)子 網(wǎng)中的每個子網(wǎng)的四個邊界點，查找在預(yù)置初始步長L范圍內(nèi)的所有點目標(biāo)及所述所有點 目標(biāo)所屬的子網(wǎng)；計算所述所有點目標(biāo)與當(dāng)前子網(wǎng)邊界點的距離，并將所述距離由小到大 排序，依次計算每條邊的最大模型相干系數(shù)，當(dāng)該最大模型相干系數(shù)大于給定閾值時即認(rèn) 為該邊為將該點目標(biāo)所屬的子網(wǎng)與當(dāng)前子網(wǎng)連接的有效邊，否則繼續(xù)計算下一條邊；當(dāng)遍 歷完所有的子網(wǎng)后，重新統(tǒng)計子網(wǎng)數(shù)目，若該子網(wǎng)數(shù)目只有1個，說明全部子網(wǎng)均已連接成 功，則無需擴展連接計算；若該子網(wǎng)數(shù)目超過1，說明還有未被連接的孤立子網(wǎng)，再依次以 2L，3L，4L，5L，……nL為步長，重復(fù)上述步驟進行擴展連接，直至循環(huán)停止，此時所有子網(wǎng) 均已連接成功；若還存在若干與主要網(wǎng)絡(luò)相隔nL以上的子網(wǎng)，將所述與主要網(wǎng)絡(luò)相隔nL以 上的子網(wǎng)拋棄，其中η為大于1的自然數(shù)。本發(fā)明實施例達(dá)到了如下技術(shù)效果：通過這種多 層級擴展方式，能夠逐步減少子網(wǎng)總數(shù)，從而實現(xiàn)不連續(xù)子網(wǎng)的快速連接，特別有利于大區(qū) 域的地表形變監(jiān)測。
[0061] 以下舉應(yīng)用實例進行說明：
[0062] 如圖3所示，是本發(fā)明應(yīng)用實例使用的河北省某試驗區(qū)24景歐洲空間局ENVISAT ASAR影像的平均幅度圖，分辨率約為20米，行數(shù)為740個像元，列數(shù)為600個像元，以下進 行詳細(xì)說明：
[0063] 如圖4所示，是本發(fā)明應(yīng)用實例一種面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法 流程圖，所述方法包括：
[0064] 401、將一區(qū)域雷達(dá)影像的所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)及所述所有不連續(xù) 的Delaunay三角網(wǎng)子網(wǎng)頂點進行編號。
[0065] 首先對所有的Delaunay三角網(wǎng)子網(wǎng)連接的點目標(biāo)進行編號，按照先列后行的順 序排列，可得到點1、點2、點3.......、點N個點目標(biāo)；
[0066] 其次按照點目標(biāo)編號順序從小到大，依次搜索與該點目標(biāo)相連的子網(wǎng)，對其編號， 分別以子網(wǎng)1、子網(wǎng)2、子網(wǎng)3........子網(wǎng)M標(biāo)識，直至完成所有點目標(biāo)的遍歷；
[0067] 在對子網(wǎng)搜索編號的過程中，同時記錄該子網(wǎng)含有的點目標(biāo)及其編號。
[0068] 402、對所述所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)中的每個子網(wǎng)的四個邊界點，查 找在以500m為半徑范圍內(nèi)的所有點目標(biāo)及所述所有點目標(biāo)所屬的子網(wǎng)。
[0069]從子網(wǎng)1開始循環(huán)計算，對于當(dāng)前子網(wǎng)，以最大外接矩形找到該子網(wǎng)的邊界點，即 最大外接矩形與子網(wǎng)的交點，通常情況下，會得到4個邊界點；但是當(dāng)該子網(wǎng)只包含兩個或 三個頂點時，則該子網(wǎng)的邊界點數(shù)目將小于4。
[0070] 然后以各邊界點為中心，以500m為半徑，在1000m*1000m的空間搜索窗口內(nèi)，查找 非當(dāng)前子網(wǎng)的點目標(biāo)，并記錄各點目標(biāo)所屬的子網(wǎng)編號。
[0071] 可選的，這里搜索半徑可根據(jù)不同分辨率SAR影像適當(dāng)調(diào)整，增大或減小均可。
[0072] 403、計算所述所有點目標(biāo)與當(dāng)前子網(wǎng)邊界點的距離，并將所述距離由小到大排 序，依次計算每條邊的最大模型相干系數(shù)，當(dāng)該最大模型相干系數(shù)大于給定閾值時即認(rèn)為 該邊為將該點目標(biāo)所屬的子網(wǎng)與當(dāng)前子網(wǎng)連接的有效邊，否則繼續(xù)計算下一條邊。
[0073] 通過各點目標(biāo)的坐標(biāo)，計算搜索到的點目標(biāo)與當(dāng)前邊界點距離，并按由小到大排 序，依次計算每條邊的最大模型相干系數(shù)。模型相干系數(shù)具體計算方法如下： IM
[0074] Υη^Αχ,η^,"^,〇ν") =-· /-〇 (I)
[0075] 式中，Y為模型相干系數(shù)；M為干涉像對個數(shù)；δ 為相鄰點目標(biāo)的二次差分相 位；SφΠΜ?ε1為相鄰高相干點的線性模型相位；（xm，ym)、（xn，yn)為兩相鄰高相干點的位置坐 標(biāo)Ji為第i幅差分相位圖的時間基線，Ti =tm-ts。
[0076] 式⑴中，相鄰點目標(biāo)（xm，ym)，（xn，yn)間的二次差分相位可表示為：

【權(quán)利要求】
1. 一種面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法，其特征在于，所述面向時間序列 InSAR的不連續(xù)子網(wǎng)連接方法包括： 將一區(qū)域雷達(dá)影像的所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)及所述所有不連續(xù)的 Delaunay三角網(wǎng)子網(wǎng)頂點進行編號； 對所述所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)中的每個子網(wǎng)的四個邊界點，查找在預(yù)置 初始步長L范圍內(nèi)的所有點目標(biāo)及所述所有點目標(biāo)所屬的子網(wǎng)； 計算所述所有點目標(biāo)與當(dāng)前子網(wǎng)邊界點的距離，并將所述距離由小到大排序，依次計 算每條邊的最大模型相干系數(shù)，當(dāng)該最大模型相干系數(shù)大于給定閾值時即認(rèn)為該邊為將該 點目標(biāo)所屬的子網(wǎng)與當(dāng)前子網(wǎng)連接的有效邊，否則繼續(xù)計算下一條邊； 當(dāng)遍歷完所有的子網(wǎng)后，重新統(tǒng)計子網(wǎng)數(shù)目，若該子網(wǎng)數(shù)目只有1個，說明全部子網(wǎng)均 已連接成功，則無需擴展連接計算；若該子網(wǎng)數(shù)目超過1，說明還有未被連接的孤立子網(wǎng)， 再依次以2L，3L，4L，5L，……nL為步長，重復(fù)上述步驟進行擴展連接，直至循環(huán)停止，此時 所有子網(wǎng)均已連接成功；若還存在若干與主要網(wǎng)絡(luò)相隔nL以上的子網(wǎng)，將所述與主要網(wǎng)絡(luò) 相隔nL以上的子網(wǎng)拋棄，其中η為大于1的自然數(shù)。
2. 如權(quán)利要求1所述面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法，其特征在于，所述將 一區(qū)域雷達(dá)影像的所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)及所述所有不連續(xù)的Delaunay三角 網(wǎng)子網(wǎng)頂點進行編號，包括： 將所述所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)編號按照點目標(biāo)的順序從小到大，依次進 行編號，分別以子網(wǎng)1、子網(wǎng)2、子網(wǎng)3.......子網(wǎng)Μ標(biāo)識，直至完成所有點目標(biāo)的遍歷；同 時記錄每個子網(wǎng)中含有的點目標(biāo)及該點目標(biāo)編號，其中Μ為大于1的自然數(shù)。
3. 如權(quán)利要求1所述面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法，其特征在于，所述預(yù) 置初始步長L根據(jù)SAR影像的分辨率和所有點目標(biāo)的密度進行調(diào)整，但最大步長nL小于或 等于大氣相關(guān)距離，即3000m ;所述預(yù)置初始步長L設(shè)為影像分辨率的50倍：當(dāng)SAR像元分 辨率為l〇m時，該預(yù)置初始步長L設(shè)為500m ;當(dāng)SAR像元分辨率為20m時，該預(yù)置初始步長 L設(shè)為1000m ;當(dāng)分辨率為3m時，該預(yù)置初始步長L設(shè)為150m。
4. 如權(quán)利要求1所述面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法，其特征在于，所述最 大模型相干系數(shù)，反映了相對形變速率和相對高程誤差對相鄰點目標(biāo)的擬合程度，其最大 值為1，最小值為〇 ;將所述最大模型相干系數(shù)的給定閾值設(shè)為〇. 7,當(dāng)某條邊的最大模型相 干系數(shù)大于0. 7時，認(rèn)為該條邊是可靠的連接關(guān)系，將兩個子網(wǎng)相連；否則該條邊不能作為 兩個子網(wǎng)的連接邊，需要繼續(xù)計算下一條邊的最大模型相干系數(shù)，直至完成所有邊循環(huán)。
5. 如權(quán)利要求1所述面向時間序列InSAR的不連續(xù)子網(wǎng)連接方法，其特征在于，當(dāng)預(yù) 置初始步長設(shè)為L時，所述當(dāng)遍歷完所有的子網(wǎng)后，重新統(tǒng)計子網(wǎng)數(shù)目，若該子網(wǎng)數(shù)目只有 1個，說明全部子網(wǎng)均已連接成功，則無需擴展連接計算；若該子網(wǎng)數(shù)目超過1，說明還有未 被連接的孤立子網(wǎng)，再依次以2L，3L，4L，5L，……nL為步長，重復(fù)上述步驟進行擴展連接， 直至循環(huán)停止，此時所有子網(wǎng)均已完成擴展連接；若還存在若干與主要網(wǎng)絡(luò)相隔nL以上的 子網(wǎng)，將所述與主要網(wǎng)絡(luò)相隔nL以上的子網(wǎng)拋棄，包括： 在完成L內(nèi)的所有子網(wǎng)搜索后，重新統(tǒng)計不連續(xù)子網(wǎng)數(shù)目，并對所述不連續(xù)子網(wǎng)頂點 進行編號；若子網(wǎng)數(shù)目只有1個，說明全部子網(wǎng)均已連接成功，則無需擴展連接計算；若子 網(wǎng)數(shù)目超過1，再以2L為步長，重復(fù)上述步驟，對每個子網(wǎng)進行相鄰子網(wǎng)搜索；以此類推，當(dāng) 新的子網(wǎng)數(shù)目仍然超過1時，再進一步擴大搜索步長距離，逐漸以3L，4L，5L，……3000m為 步長循環(huán)計算，進行子網(wǎng)鏈接，直至循環(huán)停止，此時所有子網(wǎng)均已完成擴展連接；若還存在 若干與主要網(wǎng)絡(luò)相隔3000m以上的子網(wǎng)，將所述與主要網(wǎng)絡(luò)相隔3000m以上的子網(wǎng)拋棄。
6. -種面向時間序列InSAR的不連續(xù)子網(wǎng)連接裝置，其特征在于，所述面向時間序列 InSAR的不連續(xù)子網(wǎng)連接裝置包括： 編號單元，用于將一區(qū)域雷達(dá)影像的所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)及所述所有 不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)頂點進行編號； 查找單元，用于對所述所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)中的每個子網(wǎng)的四個邊界 點，查找在預(yù)置初始步長L范圍內(nèi)的所有點目標(biāo)及所述所有點目標(biāo)所屬的子網(wǎng)； 計算單元，用于計算所述所有點目標(biāo)與當(dāng)前子網(wǎng)邊界點的距離，并將所述距離由小到 大排序，依次計算每條邊的最大模型相干系數(shù)，當(dāng)該最大模型相干系數(shù)大于給定閾值時即 認(rèn)為該邊為將該點目標(biāo)所屬的子網(wǎng)與當(dāng)前子網(wǎng)連接的有效邊，否則繼續(xù)計算下一條邊； 遍歷單元，用于當(dāng)遍歷完所有的子網(wǎng)后，重新統(tǒng)計子網(wǎng)數(shù)目，若該子網(wǎng)數(shù)目只有1個， 說明全部子網(wǎng)均已連接成功，則無需擴展連接計算；若該子網(wǎng)數(shù)目超過1，說明還有未被連 接的孤立子網(wǎng)，再依次以2L，3L，4L，5L，……nL為步長，重復(fù)上述步驟進行擴展連接，直至循 環(huán)停止，此時所有子網(wǎng)均已完成擴展連接；若還存在若干與主要網(wǎng)絡(luò)相隔nL以上的子網(wǎng)， 將所述與主要網(wǎng)絡(luò)相隔nL以上的子網(wǎng)拋棄，其中η為大于1的自然數(shù)。
7. 如權(quán)利要求6所述面向時間序列InSAR的不連續(xù)子網(wǎng)連接裝置，其特征在于，所述 編號單元，進一步用于將所述所有不連續(xù)的Delaunay三角網(wǎng)子網(wǎng)編號按照點目標(biāo)的順序 從小到大，依次進行編號，分別以子網(wǎng)1、子網(wǎng)2、子網(wǎng)3.......子網(wǎng)Μ標(biāo)識，直至完成所有 點目標(biāo)的遍歷；同時記錄每個子網(wǎng)中含有的點目標(biāo)及該點目標(biāo)編號，其中Μ為大于1的自然 數(shù)。
8. 如權(quán)利要求6所述面向時間序列InSAR的不連續(xù)子網(wǎng)連接裝置，其特征在于，所述預(yù) 置初始步長L根據(jù)SAR影像的分辨率和所有點目標(biāo)的密度進行調(diào)整，但最大步長nL小于或 等于大氣相關(guān)距離，即3000m ;所述預(yù)置初始步長L設(shè)為影像分辨率的50倍：當(dāng)SAR像元分 辨率為l〇m時，該預(yù)置初始步長L設(shè)為500m ;當(dāng)SAR像元分辨率為20m時，該預(yù)置初始步長 L設(shè)為1000m ;當(dāng)分辨率為3m時，該預(yù)置初始步長L設(shè)為150m。
9. 如權(quán)利要求6所述面向時間序列InSAR的不連續(xù)子網(wǎng)連接裝置，其特征在于，所述最 大模型相干系數(shù)，反映了相對形變速率和相對高程誤差對相鄰點目標(biāo)的擬合程度，其最大 值為1，最小值為〇 ;將所述最大模型相干系數(shù)的給定閾值設(shè)為〇. 7,當(dāng)某條邊的最大模型相 干系數(shù)大于0. 7時，認(rèn)為該條邊是可靠的連接關(guān)系，將兩個子網(wǎng)相連；否則該條邊不能作為 兩個子網(wǎng)的連接邊，需要繼續(xù)計算下一條邊的最大模型相干系數(shù)，直至完成所有邊循環(huán)。
10. 如權(quán)利要求6所述面向時間序列InSAR的不連續(xù)子網(wǎng)連接裝置，其特征在于，當(dāng)預(yù) 置初始步長設(shè)為L時，所述遍歷單元，進一步用于在完成L內(nèi)的所有子網(wǎng)搜索后，重新統(tǒng)計 不連續(xù)子網(wǎng)數(shù)目，并對所述不連續(xù)子網(wǎng)頂點進行編號；若子網(wǎng)數(shù)目只有1個，說明全部子網(wǎng) 均已連接成功，則無需擴展連接計算；若子網(wǎng)數(shù)目超過1，再以2L為步長，重復(fù)上述步驟，對 每個子網(wǎng)進行相鄰子網(wǎng)搜索；以此類推，當(dāng)新的子網(wǎng)數(shù)目仍然超過1時，再進一步擴大搜索 步長距離，逐漸以3L，4L，5L，……3000m為步長循環(huán)計算，進行子網(wǎng)鏈接，直至循環(huán)停止，此 時所有子網(wǎng)均已完成擴展連接；若還存在若干與主要網(wǎng)絡(luò)相隔3000m以上的子網(wǎng)，將所述 與主要網(wǎng)絡(luò)相隔3000m以上的子網(wǎng)拋棄。
【文檔編號】G06F17/30GK104239419SQ201410410540
【公開日】2014年12月24日 申請日期:2014年8月20日 優(yōu)先權(quán)日:2014年8月20日
【發(fā)明者】吳宏安, 張永紅 申請人:中國測繪科學(xué)研究院