ICS35.240.70基于InSAR技術的地殼形變監測規范國家市場監督管理總局國家標準化管理委員會IGB/T44146—2024 Ⅲ 1 1 1 13.2縮略語 3 4 4 5 5 5 6 6 6 7 9質量保障 9.1質量檢查 9.2真實性檢驗 10.1數據成果 10.2圖件成果 10.3文檔成果 附錄A(資料性)常用星載SAR數據格式及參數列表 附錄B(資料性)典型形變場的特征及處理方法 附錄C(資料性)不同波段SAR數據的地殼形變InSAR監測精度 附錄D(資料性)誤差類型與誤差校正 附錄E(資料性)InSAR地殼形變監測結果的對比驗證方法 附錄F(資料性)成果圖示例 附錄G(資料性)地殼形變模型反演方法 Ⅲ本文件按照GB/T1.1—2020《標準化工作導則第1部分：標準化文件的結構和起草規則》的規定起草。請注意本文件的某些內容可能涉及專利。本文件的發布機構不承擔識別專利的責任。本文件由中國科學院提出。本文件由全國遙感技術標準化技術委員會(SAC/TC327)歸口。本文件主要起草人：張景發、曾琪明、李振洪、單新建、王超、劉照言、廖明生、李志偉、葛大慶、GB/T44146—2024發展迅速，應用領域越來越廣，特別是在地殼形變測量工作中發揮著越來越重要的作用。目前，SAR數據來源多樣，InSAR技術策略各有側重、處理流程復雜，不同InSAR處理軟件要求也存在一定差異，且InSAR技術應用人員專業背景和知識水平不盡相同。因此，把現階段業已成熟的經驗總結出來，從InSAR處理各環節上進行規范，以保證InSAR處理成果質量并提高規范性與工作效率，既必要1GB/T44146—2024基于InSAR技術的地殼形變監測規范本文件適用于基于星載SAR數據開展的地殼形變監測。利用星載SAR數據開展的其他地表固體形變監測參照使用。2規范性引用文件本文件沒有規范性引用文件。3.1.13.1.2發震斷層在地表造成的永久靜態位移。3.1.3同一地域內兩次大震之間產生的地殼形變。3.1.4利用雷達與目標的相對運動，把尺寸較小真實天線孔徑用信號處理方法合成為一個較大等效天線3.1.5單視復數圖像single-lookcomplexim未作多視處理的復數形式的SAR圖像。3.1.6視線向lineofsight;LOS雷達天線相位中心指向地面目標的方向。23.1.7主圖像primaryimage在兩幅單視復數圖像進行干涉處理時，作為配準參照的圖像。3.1.8在兩幅單視復數圖像進行干涉處理時，參照主圖像進行配準的圖像。3.1.9空間基線spatialbaseline對同一目標進行兩次SAR成像時，雷達天線相位中心間的空間連線向量。3.1.10垂直基線perpendicularbaseline空間基線在垂直于主圖像視線向的投影分量。3.1.11兩幅SAR圖像獲取的時間間隔。3.1.12在無其他失相干因素影響下，使地面散射體的兩次觀測SAR回波信號之間相干性降為零的空間垂直基線長度。3.1.13在無其他失相干因素影響下，使地面散射體的兩次觀測SAR回波信號之間保持相干性的最長時間間隔。3.1.14對同一地區不同時相所獲得的雷達復圖像，經幾何變換使其同名點在空間位置對應同一地面目標的處理方法。3.1.15合成孔徑雷達干涉測量interferometricsyntheticapertureradar;InSAR利用兩幅經配準的SAR復圖像共軛相乘得到的干涉相位以獲取圖像覆蓋區域地面高程或地表形變信息的一種空間對地觀測技術。3.1.163.1.17取值為兩幅SAR復圖像配準區域中對應像素的相位相關系數，用于衡量兩幅SAR復圖像對應像素的相位相關程度和干涉圖的質量。[來源：GB/T32874—2016,3.1.10,有修改]3GB/T44146—20243.1.18失相干decorrelation由于空間基線、時間基線、熱噪聲、多普勒中心頻率漂移以及地面散射特性的改變等導致的相干性降低的現象。3.1.19因SAR視線向距離的變化而在地球參考橢球面上產生的干涉相位。[來源：GB/T32874—2016,3.1.13,有修改]3.1.20相位解纏phaseunwrapping將干涉相位由主值(模為2π)恢復為全值的過程。[來源：GB/T32874—2016,3.1.14,有修改]3.1.21合成孔徑雷達差分干涉測量differentialinterferometricsyntheticapertureradar;D-InSAR對覆蓋同一地區的兩期及以上SAR圖像進行干涉處理，得到圖像間同名像元點的的相位差，并聯合外部DEM等參考數據去除地形、平地相位貢獻，以得到衛星雷達視線方向上距離變化信息的一種InSAR處理方法。3.1.22同一顆或編隊SAR衛星對同一地區以相同的成像模式、入射角和極化方式獲取的多次重復觀測SAR圖像序列，進行干涉處理，通過時序分析，估計并削弱軌道誤差、地形數據誤差、大氣延遲等影響，以得到長時序地表形變信息的一種InSAR處理方法。3.1.23通過對在較長時間基線和較大空間基線，甚至超過臨界空間基線情況下，仍保持較高相干性的點目地表形變信息的一種InSAR處理方法。3.1.24對同一空間范圍的長時序SAR圖像集，設定較短時空基線閾值以確保干涉圖的相干性，形成一組/多組干涉對集，利用最小二乘或奇異值分解開展時序分析，估計并削弱軌道誤差、地形數據誤差、大氣延遲等影響，以得到長時序地表形變信息的一種InSAR處理方法。3.1.25地殼時序形變圖timeseriescrustaldeformationmap利用TS-InSAR方法獲取的地殼形變信息制作的體現地殼形變與時間關聯性的序列圖組。下列縮略語適用于本文件。DEM:數字高程模型(digitalelevationmodel)D-InSAR:合成孔徑雷達差分干涉測量(differentialinterferometricsyntheticapertureradar)GNSS:全球導航衛星系統(globalnavigationsatellitesystem)InSAR:合成孔徑雷達干涉測量(interferometricsyntheticapertureradar)PS-InSAR:永久散射體合成孔徑雷達干涉測量(persistentscattererInSAR)4GB/T44146—2024SAR:合成孔徑雷達(syntheticapertureradar)SBAS-InSAR:小基線集合成孔徑雷達干涉測量(smallbaselinesubsetInSAR)SLC:單視復數圖像(single-lookcompleximage)SRTM:航天飛機雷達地形測繪任務(shuttleradartopographymission)TOPS:步進掃描地形觀測(terrainobservationbyprogressivescans)TS-InSAR:合成孔徑雷達時序干涉測量(timeseriesInSAR)4基本要求采用InSAR技術開展地殼形變監測，應滿足以下基本要求。a)星載SAR數據成像后統一采用SLC數據格式。b)測量單次地殼形變，至少獲取≥50mm地殼形變發生前后各1期SAR數據；測量時序地殼形變，宜積累不少于15期SAR數據。d)不宜在如下地區開展基于InSAR技術的地殼形變監測：2)完全被濃密植被覆蓋地區。基于InSAR技術的地殼形變監測總體流程見圖1,主要內容如下。a)資料與數據準備。根據工作目標及地殼形變特性確定工作區范圍，收集工作區的基礎資料和SAR數據。常用星載SAR數據格式及參數列表見附錄A。b)形變起始參考點選擇。宜選擇無形變區高相干點作為形變計算起始參考點。c)InSAR數據處理。根據SAR數據累積量等因素選擇D-InSAR或TS-InSAR(包括PS-InSAR和SBAS-InSAR)方法進行干涉處理。典型形變場的特征及處理方法見附錄B。d)質量保障。檢查處理后獲得的地殼形變圖質量，并與其他獨立的地殼形變測量成果進行不低于觀測成果5%樣本量的對比驗證。優先采用直接對比驗證，適宜時可采用交叉對比驗證和間接對比驗證。5GB/T44146—2024資料與j數據準備PS-1nSARSBAS-InSAR地殼形變圖質量檢查真實性檢驗數據成果圖件成果文檔成果成果編制工作報告質量保障圖1基于InSAR技術的地殼形變監測總體流程圖6資料與數據準備6.1基礎資料收集應收集以下資料。a)與地殼應力和形變測量有關的各類資料。例如：工作區內水準網、GNSS網的布設情況及測量況資料。b)與干涉失相干有關的資料。例如厚植被、冰雪覆蓋、強降雨導致洪澇等。c)地面覆蓋資料。包括：工作區主要地面覆蓋類型；對SAR觀測有影響的各種地物分布情況，包e)DEM數據。其覆蓋范圍應略大于相應SAR數據的覆蓋范圍；格網間隔宜接近SAR數據分辨率。6.2SAR數據獲取應收集工作區范圍內SAR數據以及衛星所提供的精密軌道數據。組成數據像對或數據集的SAR數據獲取要求如下。6GB/T44146—2024干波源為SAR數據源。b)軌道。應獲取同一軌道且覆蓋相同工作區的SAR數據；衛星飛行姿態穩定，L、S、C波段的衛星軌道定位精度宜優于0.05m,X波段的衛星軌道定位精度宜優于0.03mm。數據。不同波段SAR數據的地殼形變InSAR監測精度見附錄C。d)極化方式。優先選用同一極化方式的SAR數據，如垂直極化或者水平極化的SAR數據。e)圖像重疊率。當大尺度工作區需要同軌或跨軌的多景SAR圖像進行拼接時，同軌相鄰兩景SAR圖像的重疊率應超過沿軌圖像長度的15%;跨軌數據相鄰兩景圖像的重疊率應超過圖像幅寬的15%。f)時間基線。應選擇時間間隔較短的SAR數據組成像對或數據集。臨界時間基線與SAR波長、信號帶寬、氣候條件、植被覆蓋程度等有關。干旱、少植被地區，即使時間基線較長(如10年甚至更長),干涉圖依舊可保持較好相干性；植被較密的地區，即使時間基線較短(如幾天或更短),干涉失相干現象依舊可能嚴重。g)空間基線。應選擇小于臨界空間基線的SAR數據組成像對或數據集。臨界空間基線與SAR保證較高的相干性，較理想的空間基線在X波段小于或等于100m,C波段小于或等于300m,L波段小于或等于1000m。7形變起始參考點選擇宜在無形變區內選擇易辨別的高相干地物作為形變計算起始參考點。不宜在礦區塌陷區、油氣開域選擇形變起始參考點。典型形變的形變起始參考點選取要求如下：a)同震形變，宜選擇距震中較遠的無形變區內的易辨別、高相干地物；b)震間形變，宜選擇垂直于斷層方向，距活動斷層較遠的無形變區內的易辨別、高相干地物；c)震后形變，宜選擇距發震斷層較遠的無形變區內的易辨別、高相干地物；d)火山形變，宜選擇距火山口較遠的無形變區內的易辨別、高相干地物；f)地質災害形變，宜選擇遠離地質災害區的無形變區內的易辨別、高相干地物。形變起始參考點亦可選取其他高精度技術手段持續監測的易辨別的高相干地物。8InSAR數據處理8.1處理方法選擇根據地殼形變量(典型形變場的形變量見附錄B)、持續時間、SAR數據累積量、數據時間間隔分布、工作目標等選擇InSAR處理方法。a)當形變量大于或等于0.02m時，宜選擇D-InSAR方法；如有多期SAR數據，亦可選擇TS-In-SAR方法。b)當地表形變量小于0.02m,宜選擇TS-InSAR方法。c)當選擇TS-InSAR方法時，可根據累積SAR數據集內各數據像對的垂直基線情況，選擇PS-InSAR或SBAS-InSAR方法：根據垂直基線分成幾個集合組，每個集合組內垂直基線均滿足6.2g)中有關空間基線的要求，且單個集合組內SAR數據累積量至少15期時，應選擇SBAS-7GB/T44146—2024InSAR方法；當集合組只有1個時，應選擇PS-InSAR方法。8.2D-InSAR數據處理D-InSAR數據處理流程見圖2。SAR數據成像處理精密軌道數據精密軌道數據原始干涉相位圖生成去除平地相位模擬地形相位去除地形相位模擬地形相位差分T涉相位圖相位噪聲濾波相干系數計算和位解纏省達視線向形變計算地理編碼地殼形變圖圖2D-InSAR處理流程圖原始數據可選擇下列一種成像算法生成SLC數據： 距離-多普勒算法(Range-DopplerAlgorithm,RDA);——頻調變標算法(ChirpScalingAlgorithm,CSA);-—波數域算法(w-kAlgorithm,w-k)。同一組SLC數據集生成應采用單一成像算法。8GB/T44146—2024按以下步驟對主圖像、輔圖像和DEM進行圖像配準和重采樣，可輸入精密軌道數據參與圖像配準a)宜使用與SAR數據分辨率接近、源于SRTM的DEM數據，應依據主圖像的衛星軌道參數將DEM轉換成主圖像格網的SAR模擬圖像。b)宜使用互相關和幾何配準方法計算主圖像和輔圖像之間特征像元點的偏移量，實現粗配準。c)按以下方式實現精配準：——非TOPS模式下宜使用最大相關系數法、最大頻譜法等精配準方法實現SAR圖像精配準，配準精度應達到1/10個SAR像元； TOPS模式下宜使用增強譜分集法(enhancedspectraldiversity,ESD)實現SAR圖像精配準，配準精度應達到1/1000個SAR像元。按以下步驟生成原始干涉相位圖。直基線和斜視角不同導致的不重疊頻譜。φin=angle(S?·S?)=angle(|S?|·|S?|e(?-k?))…S?輔圖像S?經過復共軛運算后的復數值； 注2:基于InSAR技術的地殼形變監測是去除原始干涉相位圖中，除地殼形變相位之外的其他相位成分的過程。原始干涉相位圖的干涉相位構成及誤差類型與誤差校正見附錄D。宜使用以下方法估計干涉空間基線：按以下步驟去除原始干涉相位圖中的平地相位。a)獲取SAR數據中軌道參數，計算主圖像和輔圖像沿方位向每一行的平行基線。9GB/T44146—2024式中：…………λ———SAR波長，單位為米(m);B/——空間平行基線，單位為米(m)。c)將步驟b)的計算結果進行二維多項式擬合，得到整幅圖像各點的平地相位，也可以直接按公式(3)計算整幅圖像的二維平地相位，再將其轉化為平地參考相位復數值。f(x,y)——二維平地相位；x,y——圖像點位置坐標；ai-;——多項式系數；d——多項式系數的階，d=1,2,3。…………注：平地參考相位計算原理：選擇參考橢球基準面，在同一參照坐標系下，確定基準面中任一點的坐標和該點所對應的兩次成像衛星的位置，計算場景內基準面上的點到衛星的距離差。對多個點的計算結果進行最小二乘擬合，就得到基準面的參考相位。從干涉相位中減去參考相位就得到由地形高差以及形變引起的相位圖。d)按公式(4),將主圖像、輔圖像及平地參考相位復數值進行復共軛相乘即得到去除平地相位后的干涉相位圖。φflatening=angle(S?·S?·R*)…………(4)——去除平地相位后的干涉相位圖，單位為弧度(rad);S?--—主圖像對應像元的復數值；S?———輔圖像S?經過復共軛運算后的復數值；R*———平地參考相位圖像R對應像元經過復共軛運算后的復數值。8.2.2.6模擬地形相位宜使用與SAR圖像分辨率接近的源于SRTM的高精度DEM數據，按以下步驟模擬地形相位：a)建立SAR成像參數的圖像坐標與DEM大地坐標之間映射關系，形成DEM數據的模擬地形相位圖像；b)將模擬地形圖像與真實SAR圖像進行精確配準，獲得真實SAR圖像坐標系的高程數據；c)利用主圖像和輔圖像的精密軌道數據計算真實SAR圖像中每個像元的模擬地形相位。8.2.2.7去除地形相位從8.2.2.5中獲得的去除平地相位后的干涉相位中減去8.2.2.6獲取的模擬地形相位，得到去除地形相位后的干涉相位圖(差分干涉相位圖)。8.2.2.8相干系數計算按公式(5)計算去除地形相位后的干涉相位圖中各個像素點的相干系數，得到相干性圖。…………GB/T44146—2024式中：S,、S,分別為滑動窗口中主圖像和輔圖像對應像元的復數值；S?輔圖像S?經過復共軛運算后的復數值。注：計算相關系數的滑動窗口一般設為5×5。8.2.2.9相位噪聲濾波可采用下列濾波方法去除8.2.2.7得到的干涉相位圖中的噪聲影響：——空間域濾波方法：等值線濾波、中值濾波、均值濾波、圓周期濾波等；對于濾波后的干涉相位圖，應檢查工作區內相干性圖的相干系數。如果區域的大部分相干系數小于0.3,則需要變換濾波方法進行相位噪聲濾波。8.2.2.10相位解纏查看8.2.2.9得到的差分干涉相位圖中的相干性分布，應掩膜低相干區域，使用周邊可靠數值填充或插值；應采用一定的解纏方法處理較高相干性的區域；對同震形變梯度較大區域，可采用彈性半空間位錯模型輔助解纏，從干涉圖中扣除源于模型模擬的一階地殼形變，再解纏殘差相位。按公式(6)計算解纏后的真實相位。Φi,;=φi,;+2kπ………式中：Φ差分干涉相位圖解纏后的真實相位；k——正整數，為相位纏繞的整周期數，可通過判別相鄰像元的相位差值跳躍變化而確定。8.2.2.11雷達視線向形變計算按公式(7)計算解纏后的真實相位對應的雷達視線向形變。式中：…Φ,;——解纏后的真實相位。8.2.2.12地理編碼利用8.2.2.6a)建立的SAR圖像坐標與大地坐標之間映射關系，基于反向查找表，對8.2.2.11得到的雷達視線向形變從SAR圖像坐標轉換到大地坐標，得到地殼形變圖。8.2.2.13地殼形變圖地理編碼后得到的地殼形變圖，在確定InSAR-LOS向形變量主要由垂直向或水平向形變引起的GB/T44146—2024Du=Dlos/cosθine DH=Dlos/sinθine 8.3PS-InSAR數據處理PS-InSAR數據處理流程見圖3。去誤差干涉圖時間及空間域形變和位估算雷達視線向形變計算形變計算地殼時序形變圖離散點相位解纏線性形變相位和DEM誤差相位估算人氣相位估算形變時序相位估算地理編碼GB/T44146—2024a)主圖像選擇：從時序SAR數據中選擇1幅主圖像，與時序SAR數據其他所有圖像組成多組像對，主圖像宜處于時間中點，與其他圖像組成像對的空間基線不大于臨界空間基線；b)復圖像配準、重采樣和裁切：按8.2.2.2對主圖像與各輔圖像組成的像對進行復圖像配準、重采樣和裁切，形成多組干涉像對；c)DEM與主圖像配準和裁切：按8.2.2.2對DEM數據與主圖像進行配準和裁切；d)空間基線估計：按8.2.2.4執行；e)模擬地形相位：按8.2.2.6計算多組干涉像對對應的模擬地形相位。按8.2.2.3～8.2.2.9的步驟對多組干涉像對分別執行差分干涉數據處理，獲得去除平地相位、去除地形相位后的多幅干涉相位圖和相干性圖。時間及空間域形變相位估算步驟如下。a)PS點選擇：宜選擇穩定的強散射性點目標，其對應圖像像元的相干系數、幅度等數值較大，可通過計算多幅圖像同名點的相干系數、幅度的離差等，設置閾值以提取PS點；應提取一定密度的PS點(如≥100個/km2),并去除在雪面、水體等區域出現的偽PS點。b)離散點相位解纏：可按8.2.2.10對PS點進行相位解纏，也可使用周邊可靠數值插值或模擬數據進行相位解纏。c)線性形變相位和DEM誤差相位估算：宜采用最小二乘法估計PS點的線性形變相位和DEM誤差相位。d)殘余相位低通濾波：可采用最小二乘法求解每個PS點的殘余相位(包含非線性形變相位、大氣相位和其他噪聲相位);對此殘余相位進行空間低通濾波處理，以去除空間高頻噪聲相位。e)大氣相位估算：應基于SAR數據自身或外部數據計算大氣相位。可用時序方法分析SAR數據自身相位時空特性，估計并分離出大氣相位；也可使用外部數據計算大氣相位，計算方法見附錄D中D.3。f)非線性形變相位估算：從低通濾波后的殘余相位中去除高頻大氣相位，得到初步的非線性形變g)時間維度低通濾波：對各成像時刻的殘余相位進行時間維度低通濾波，將時間高頻的大氣相位、其他噪聲相位與時間低頻的非線性形變相位進一步分離，得到非線性時序形變相位。h)形變時序相位估算：疊加前面步驟得到的線性形變相位和非線性形變相位。形變計算步驟如下：a)雷達視線向形變計算：按8.2.2.11執行；c)地殼時序形變圖：按8.2.2.13執行，多個時相組合成地殼時序形變圖。SBAS-InSAR數據處理流程見圖4。GB/T44146—2024時序時序SAR數據DEM數據小基線干涉像對組合DEM與上圖像配準和裁切組內主圖像選擇數據預處理復圖像配準、重采樣和裁切十涉數據集生成線性形變相位和DJM誤差相位估算殘余相位低通濾波奇異值分解殘余相位和線性形變相位大氣相位估算線性形變相位估算非線性形變相位估算形變時序相位估算占達視線向形變計算地理編碼地殼時序形變圖空問基線估計差分T涉數據處理模擬地形相位去誤差干涉圖去誤差干涉圖1時間及空間域形變和位估算形變計算數據預處理應按以下步驟執行。a)小基線干涉像對組合：綜合考慮空間垂直基線與時間基線以及各圖像對之間多普勒中心頻率差異，對時序SAR數據進行分組，形成若干個小基線干涉像對組合，使得組內的SAR圖像干b)組內主圖像選擇：每個組合中選擇1幅圖像作為主圖像，與組內其他圖像組成像對，主圖像宜處于時間中點，與其他圖像組成像對的空間基線遠小于臨界空間基線、時間基線不高于臨界時間基線。c)復圖像配準、重采樣和裁切：按8.2.2.2對主圖像與各輔圖像組成的像對進行配準、重采樣和裁d)DEM與主圖像配準和裁切：按8.2.2.2對DEM數據與主圖像進行配準和裁切。GB/T44146—2024地形相位后的多幅干涉相位圖和相干性圖。時間及空間域形變相位估算步驟如下：a)線性形變相位和DEM誤差相位估算：宜采用最小二乘方法，估算小基線組內線性形變相位和DEM誤差相位，并從干涉相位中扣除由于DEM配準不精確造成的DEM誤差相位；e)奇異值分解(SingularValueDecoposition,SVD):將由于長基線限制而分開的獨立原始數據變時序相位。形變計算步驟如下：9質量保障InSAR地殼形變監測結果的質量檢查應至少包括表1所列內容。表1InSAR地殼形變監測結果質量檢查內容檢查內容檢查方法主、輔復圖像配準精度是否符合8.2.2.2要求干涉圖相位與地形的相關性檢查干涉圖相位變化是否與地形呈高相關性。相關性弱，則質量較好；相關性強，則質量較差差分干涉相位統計值統計值直方圖是否大致服從正態分布差分干涉圖與相干性圖對比應符合8.2.2.9有關相干性閾值的要求。小于0.3,質量差；0.3～0.5,質量良好；大于0.5,質量很好差分干涉圖紋路檢查檢查紋路是否清晰可辨。如清晰可辨，對應區域干涉圖質量較好；反之，則干涉圖質量較差GB/T44146—2024表1InSAR地殼形變監測結果質量檢查內容(續)檢查內容檢查方法解纏誤差殘差點分布圖采用殘差點空間分布的統計方法核驗分布圖。零星分布，質量較好；不均勻分布或密集分布，質量較差。構建時空閉合環，檢查解纏相位的閉合差，探測解纏誤差PS點分布是否符合8.3.2.3a)的要求。PS點目標分布密度較均勻，多于100個/km2。采用檢查PS點地物屬性的方法。多與反射能力較強的實際地物對應，水體中極少出現，則質量較好，否則質量較差小基線干涉像對組合是否符合8.4.2.1的要求。干涉像對組內空間基線如果遠小于臨界空間基線，時間基線也不大于臨界時間基線，主圖像處于時間中點位置，則質量較好；反之，質量較差時間序列圖的離散度和曲線采用檢查離散度收斂情況和跨斷層時序曲線方法。離散度低，質量好；跨斷層時序曲線，與活動斷層的升降對應較好，如相對上升盤對應于正斷層的下盤或逆斷層的上盤，則質量較好，否則質量較差9.2真實性檢驗9.2.1直接對比驗證應利用同時期地殼形變的地面監測結果，對InSAR地殼形變監測結果進行直接對比驗證。InSAR地殼形變監測結果宜與水準測量、GNSS測量等高精度監測結果進行對比。a)與水準測量結果對比。宜采用時序InSAR監測結果的多點均值與對應時段的水準測量結果對比，將所查找到的一定空間范圍內高相干目標InSAR監測結果的均值與水準測量點結果對比。b)與GNSS測量結果對比。使用研究區內相同時段高等級GNSS測站點的監測結果進行檢驗。選取距離高等級GNSS站點一定范圍內的InSAR高相干點目標參與驗證。獲取擬合數據參與驗證。直接對比驗證時應考慮InSAR、GNSS和水準測量在測量原理、測量時間、空間尺度及基準方面的差異。InSAR地殼形變監測結果的直接對比驗證方法見附錄E中E.1。使用具有重疊時空的SAR數據得到的InSAR監測結果進行交叉驗證。這些SAR數據可來源于相同衛星不同軌道、不同成像模式，或來源于不同衛星、將不同來源的InSAR地殼形變監測結果投影到同一方向，選出同名點統計差值及均方根誤差。InSAR地殼形變監測結果的交叉對比驗證方法見E.2。使用InSAR地殼形變監測結果統計數據，根據地震學原理或災害情況進行對比驗證。a)同震形變，可通過兩種間接對比方法開展驗證。1)利用地震波和/或其他大地觀測數據反演得到地震斷層參數，采用彈性半空間的位錯模型模擬同震形變，進而與InSAR監測結果進行相關性分析，求取兩者差值中誤差。2)基于InSAR獲得的形變結果統計數據，根據地震學原理或震害情況定性檢驗InSAR測GB/T44146—2024b)地質災害形變，可對照地質災害形變區是否與采礦、地下水開采等活動有關，再加以評估和驗證。a)原始數據與資料，包括參與InSAR處理的原始SAR觀測數據，處理所需的各類輔助資料；a)制圖范圍：圖幅內圖廓線范圍涵蓋實際b)坐標系：制圖坐標系應采用2000國家大地坐標系(CGCS2000,ChinaGeodeticCoordinateSys-tem2000);c)比例尺：根據不同應用目的選擇1:250000、1:100000、1:50000、1:25000或1:10000InSAR處理生成地殼形變分布成果圖的成果格式、示例圖及制圖說明表見附錄F。文檔成果宜包括但不限于如下內容：d)數據處理方法及關鍵參數選取；e)地殼形變信息提取結果及其時空分布，也可給出地殼形變模型反演方法獲得的地殼運動模型(見附錄G);SAR處理成果滿足要求。GB/T44146—2024(資料性)常用星載SAR數據格式及參數列表表A.1給出了常用星載SAR數據的格式。SAR衛星名稱數據格式圖像文件頭文件ERS-1/2CEOS(*.IMG,*.001,*.dat)CEOS(*.001)CEOS(*.IMG,*.001,*.dat)CEOS(*.001)ENVISATASAREnvisat(*.N1)Envisat(*.N1)ALOS-1ALOSL1.1(*H1.1_A)Radarsat-1CEOS(×.IMG,*.001,*.dat)CEOS(*.001,*.TRL,*.TRA)Radarsat-2RADARSAT2(*.tif)Terrasar-XTERRASAR-X(*.cos)Cosmo-SkyMedStripmapHimageCosmo-skymed(*.h5)ALOS-2ALOSL1.1Sentinel-1IWSentinel-1file(*.tif)SAOCOMSAOCOMfile(*.tif)CapellaSARCapellafile(*.tif)海絲一號TY-1file(*.tif)Lutan-1LT-1file(*.tif)A.2星載SAR數據參數列表星載SAR數據一般包括以下幾種獲取模式：a)聚束模式(spotlight),分辨率一般為幾分米至幾米；b)(超級)條帶模式(stripmap),分辨率一般為幾米至幾十米；c)寬幅模式(burst、ScanSAR、TOPS等),分辨率一般為十幾米至百米。表A.2給出了可用于InSAR處理的常用星載SAR數據的主要參數。名稱所屬國家或機構運行時間軌道高度/km波長/m極化方式入射角/最短間隔/d地面分辨率(m)及模式RAW數據存檔情況幅寬/可否編程ERS-1歐洲空間局1991年—2000年C(0.056)VV是全球覆蓋20次以上否ERS-2歐洲空間局1995年—2012年C(0.056)VV是全球覆蓋20次以上否JERS-1日本1992年—1998年L(0.235)HH否全球覆蓋5次以上否Radarsat-1加拿大1995年—2013年C(0.056)HH8～30是部分地區覆蓋50～500否ENVISAT-ASAR歐洲空間局2002年—2012年C(0.056)雙極化25～100是全球覆蓋25次以上否ALOS-PAISAR日本2006年—2011年L(0.236)全極化是全球覆蓋15次以上20～350否Radarsat-2加拿大2007年月至今C(0.056)單極化、雙極化、全極化3(超級條帶)5(條帶)16～100(其他)否2013年間多期存檔20(超級條帶)50(條帶)50～500(其他)可TerraSAR-X德國2007年月至今X(0.031)單極化、雙極化、全極化20～5511(單)1(聚束)3(條帶)18.5(掃描)40(寬掃描)否大部分地區需要編程觀測20(聚束)30(條帶)150(掃描)270(寬掃描)可TanDEM-X星座(2)德國2010年6月至今X單極化、雙極化、全極化11(單星)1(聚束)3(條帶)18.5(掃描)40(寬掃描)否大部分地區需要編程觀測20(聚束)30(條帶)150(掃描)270(寬掃描)可COSMO-SkyMed星座(6)意大利2007年6月、2007年12月、2008年10月、2010年11月、2019年12月、2022年2月X單極化、雙極化16.36~52.0616(單星)1(雙星/31(聚束)1、3、5(條帶)30(掃描)100(寬掃描)否大部分地區需要編程觀測7～10(聚束)30～40(條帶)100～00(掃描)可表A.2常用星載SAR數據主要參數(續)名稱所屬國家或機構運行時間軌道高度/km波長/m極化方式入射角/最短間隔/d地面分辨率(m)及模式RAW數據存檔情況幅寬/可否編程Sentinel-1星座歐洲空間局2014年4月至今CHH十VV、VV+VH12(單星)6(雙星)5(聚束)5×20(條帶)20(掃描)20(寬掃描)是全球觀測計劃，每6天/12天1次同軌觀測20(聚束)80(條帶)250(掃描)400(加寬掃描)否ALOS-2(PAL-SAR-2)日本2014年5月至今L全極化8~1×3(聚束)3、6、10(條帶)100(掃描)否全球觀測計劃25(聚束)30～50(條帶)350～490(掃描)可Lutan-1A中國2022年1月至今L單極化、雙極化、全極化8(單)4(雙)3(條帶1)12(條帶2)20～30(寬幅)30(掃描)是全球觀測計劃50(條帶1)100(條帶2)150～500(寬幅)400(掃描)可Lutan-1B中國2022年2月至今L單極化、雙極化、全極化8(單)4(雙)3(條帶1)12(條帶2)20～30(寬幅)30(掃描)是全球觀測計劃50(條帶1)100(條帶2)150～500(寬幅)400(掃描)可阿根廷2018年10月、2020年8月至今L單極化、雙極化、全極化16(單星)8(雙星)10(條帶1)10(條帶2)30(掃描窄幅1)50(掃描窄幅2)50(掃描寬幅1)100(掃描寬幅2)是全球觀測計劃20(條帶1)40(條帶2)100(掃描窄幅1)150(掃描窄幅2)220(掃描寬幅1)350(掃描寬幅2)可品SAR衛星名稱所屬國家或機構運行時間軌道高度/km波長/m極化方式入射角/最短間隔d地面分辨率(m)及模式數據存檔情況幅寬/可否編程星座(27)芬蘭2018年1月至今~580X單極化22(單星)1(27星)0.5(SPOTFINE)3(STRIP)是全球觀測計劃5(聚束)30(條帶)120(掃描)可Capella星座(10)美國2018年12月至今X單極化15(單星)0.35(聚束)0.6(滑塊聚束)0.8(條帶)是全球觀測計劃5(聚束)200(滑塊聚束)10～50(條帶)否海絲一號中國至今C單極化1(聚束)3(條帶)10(掃描NS)20(掃描ES)是全球觀測計劃5(聚束)20(條帶)50(掃描)100(加寬掃描)可巢湖一號中國2022年2月至今C單極化1(聚束)2(滑聚)3(條帶)12(掃描NS)20(掃描ES)10(TOPS模式)是全球觀測計劃7(聚束)25(條帶)100(掃描NS)170(掃描ES)可高景二號星座中國2022年7月至今X單極化0.5(滑動聚束)1(滑動聚束)2(條帶)是全球觀測計劃10(滑動聚束)20(條帶)可宏圖一號星座中國2023年3月至今X單極化0.5(聚束1)1(聚束2)3(條帶)5(TOPSAR)3(多星條帶干涉)是全球觀測計劃5(聚束1)10(聚束2)20(條帶)50(TOPSAR)20(多星條帶干可GB/T44146—2024(資料性)典型形變場的特征及處理方法表B.1給出了典型形變場的特征及處理方法。表B.1典型形變場的特征及處理方法形變類型形變場特征時間過程分類處理方法空間分布范圍形變量或形變速率瞬時緩慢D-InSARTS-InSAR地震√√√震間形變<20mm/年√√震后形變√√火山火山隆起/下沉<20mm/年√√火山爆發√√滑坡滑坡前兆0.5～5√√地面沉降地下水開采0.5～100(20～100)mm/年√√√礦產開采0.1～10(20～100)mm/年√√√√的地殼形變時加以分辨并剔除。“事件過程持續數秒至數分鐘。事件過程持續數天至數月，甚至數年。“包括兩次特征地震期間與長期應變積累(震間)、地震發生(同震)等過程相關的地殼表面變形。d根據統計分析，利用D-InSAR方法能探測到的可靠形變量50mm(對應最小震級地震Ms6),利用TS-InSAR方法能探測到的可靠形變量20mm(對應最小震級地震Ms4.5)。GB/T44146—2024(資料性)不同波段SAR數據的地殼形變InSAR監測精度在植被覆蓋密集區域表現優于C波段和X波段；在形變量大于或等于1m的地區，如同震形變、滑坡時序變形和礦區開采形變等地區，有較好的性能；在有利條件下，D-InSAR地殼形變測量精度可達10mm,多時相數據TS-InSAR測量精度可達5mm/年～8mm/年。C.2C波段數據在植被相對稀少區域表現優于X波段數據；在有利條件下，D-InSAR地殼形變測量精度可達8mm,多時相數據TS-InSAR測量精度可達2mm/年～5mm/年。C.3X波段數據在裸露區及重大人工建筑目標地區具有較好的探測效果，如城區精細形變測量等；在有利條件下D-InSAR地殼形變測量精度可達5mm,多時相數據TS-InSAR測量精度可達2mm/年～3mm/年。GB/T44146—2024(資料性)誤差類型與誤差校正D.1原始干涉相位圖構成原始干涉相位圖的干涉相位主要包含6種相位成分，按公式(D.1)計算。…(D.1)式中：pfat——平地相位；φde——SAR兩次觀測期間目標沿衛星視線方向移動引起的形變相位；pam——兩次觀測間大氣環境差異引起的延遲相位；在差分InSAR數據處理中，要獲取精確的地殼形變相位φde,就要去除原始干涉相位圖的干涉相D.2誤差類型成因分析InSAR形變處理結果主要包含以下誤差。a)失相干誤差。主要包括時間失相干誤差、空間失相干誤差以及熱噪聲。時間失相干誤差與地面覆蓋物散射特性變化有關；空間失相干誤差由SAR天線視角差導致的微波波數差異造成。b)大氣延遲誤差。指由于雷達信號在傳輸過程中，大氣折射的非均質性引起的雷達信號傳播路徑彎曲和雷達信號傳播延遲；大氣對雷達信號傳播的影響包含受對流層和電離層的影響。c)衛星軌道誤差。屬于系統誤差，通常以殘差條紋形式存在于干涉圖中。d)DEM誤差。在利用DEM對地形相位進行估計時，所用DEM誤差會導致干涉相位中殘留地形相位誤差項，其大小與垂直基線長度和雷達入射角有關。D.3誤差校正常見誤差校正方法如下。a)失相干誤差。濾波方法和時序分析方法相結合方法去除失相干地物，保留高相干地物。b)大氣延遲誤差，校正方法如下。1)單景干涉圖像。利用外部數據校正大氣誤差；外部數據包括GNSS大氣延遲產品、MER-IS和MODIS等的水汽產品，以及大氣數值模型大氣產品(ECMWF、MERRA2和GA-COS);利用線性或者冪律模型削弱與地形相關的對流層延遲影響。2)多時相干涉圖像。利用時序方法分離出大氣誤差相位分量，用最小二乘方法或者奇異值分解(SVD分解)估計每景SAR圖像對應的大氣相位；可用外部數據，計算出一階大氣誤差相位，并從殘余相位中減去，然后用主成分分析法估計殘余大氣誤差相位。3)電離層誤差。使用圖像本身的電離層統計分量、電離層模型，或者split-spectrum算法(方位向)進行校正。c)衛星軌道誤差。基于地面控制點估算一個相位趨勢面，從原始相位中減去，實現軌道誤差校正。GB/T44146—2024(資料性)InSAR地殼形變監測結果的對比驗證方法E