無人機(jī)遙感圖像及其三維可視化在汶川地震救災(zāi)中的應(yīng)用

1、物探化探計(jì)算技術(shù)第 32卷 第 2期2010年3月文章編號(hào):1001) 1749( 2010) 02) 0206) 05無人機(jī)遙感圖像及其三維可視化在救災(zāi)中的應(yīng)用1,苗 放1,2,3( 1. 成都理工大學(xué)地球科學(xué)學(xué)院, 四川 成都610059;610041;2. 四川省科技大學(xué)局減災(zāi)救助研究所, 四川 成都3.研究所, 四川 成都610054)地表空間摘 要: 在分析了無人機(jī)低空遙感影像幾何變形原理, 對(duì)無人機(jī)遙感數(shù)據(jù)進(jìn)行圖像拼接與幾何校正, 正校正等處理后與高精度 DEM 融合, 制作出了研究區(qū)景觀的三維遙感虛擬動(dòng)態(tài)模型。通過分析三維可視化圖像, 結(jié)合地形資料, 對(duì)重災(zāi)區(qū)安縣茶坪河流域由山體

2、滑坡和崩塌等地貌變化造成的堰塞湖和相關(guān)災(zāi)害, 進(jìn)行了高分辨率、低空遙感的災(zāi)害獲取和災(zāi)情評(píng)估。實(shí)踐表明, 在災(zāi)害應(yīng)急和復(fù)雜地形條件下, 使用無人機(jī)低空遙感及其三維可視化技術(shù), 進(jìn)行實(shí)地和提取, 能夠?qū)φ鸷蟾鞣N實(shí)時(shí)情景進(jìn)行精確描述, 為決策者提供了準(zhǔn)確、詳細(xì)、及時(shí)的第資料, 從而最大限度地減少損失。: 無人機(jī); 遙感; 三維可視化;號(hào): TP 7文獻(xiàn)標(biāo)識(shí)碼:A經(jīng)過圖像幾何校正、鑲嵌、三維可視化、圖像解譯等快速處理, 實(shí)現(xiàn)了高危堰塞湖, 重災(zāi)區(qū)受災(zāi)狀況等目標(biāo)的實(shí)時(shí)動(dòng)態(tài)監(jiān)測(cè), 為各級(jí)領(lǐng)導(dǎo)和抗震救災(zāi)指揮0前言由于四川盆地多云多雨的特點(diǎn), 造成的決策, 提供了最及時(shí)可靠的數(shù)據(jù)和圖片支大部分被雨云遮蓋, 無

3、法準(zhǔn)確觀測(cè)特定地區(qū), 發(fā)揮持。大范圍以及綜合的優(yōu)勢(shì), 可以在災(zāi)害1數(shù)據(jù)源發(fā)生的第一時(shí)間內(nèi), 獲得即使調(diào)整拍攝角度也無法得到的詳細(xì)災(zāi)區(qū)情形。由于無人機(jī)低空遙感系統(tǒng)具有很高的機(jī)動(dòng)性、靈活性和安全性, 可獲取多角度、高分辨率影像, 不受高度限制和陰云天氣影響, 且系統(tǒng)成本及影像處理費(fèi)用較低。若能發(fā)揮其特有的作用, 就可以為決策者提供準(zhǔn)確、詳細(xì)、及時(shí)的第 資料, 在低空領(lǐng)域、小區(qū)域具有一定的在本文中使用的無人機(jī)低空遙感數(shù)據(jù), 采用了AF- 1000小型無人航攝飛機(jī)獲取。AF- 1000是中短程無人機(jī), 采用常規(guī)空氣動(dòng)力布局, 具有重量輕, 翼載荷小, 起飛降落場(chǎng)地要求外場(chǎng)地即可作業(yè), 特別適合突特點(diǎn)。

4、在野況下快速部?jī)?yōu)勢(shì) 1 。作者在本文中署, 能夠在五十公里半徑范圍內(nèi)執(zhí)行短時(shí)航拍任務(wù)。無人機(jī)獲取的低空遙感數(shù)據(jù), 具有極高的空間分辨率, 單幅影像輻射范圍小, 影像總量巨大。為了快速實(shí)現(xiàn)影像的/ 災(zāi)情地圖 0作用, 采用了自動(dòng)了無人機(jī)低空遙感數(shù)據(jù)的快速處理, 三維可視化技術(shù)與, 利用無人機(jī)航空平臺(tái)搭載多種遙感傳感器, 在第一時(shí)間內(nèi)獲取大量重點(diǎn)災(zāi)區(qū)的高分辨率航空遙感影像數(shù)據(jù)?；痦?xiàng)目:自然基金項(xiàng)目 ( 60841006) ; 四川省青年基金項(xiàng)目 ( 09ZQ026- 058)收稿日期: 2009- 09- 16改回日期: 2009- 11- 182期等: 無人機(jī)遙感圖像及其三維可視化在207救災(zāi)

5、中的應(yīng)用。圖 1分別是北川縣城和青影像屬于中心投影瞬間一次成像, 即一幅影像上的所有像點(diǎn)共用一個(gè)影像中心和同一個(gè)像片平面, 亦即共用一組外方位元素。因此, 像點(diǎn) a與物點(diǎn) A 之間的空間坐標(biāo) , 可表示為中心投影成像的共線方程式( 1)如下:拼接的圖形處理川東河口堰塞湖, 經(jīng)拼接后的局部影像。2無人機(jī)遙感圖像處理由于發(fā)育在環(huán)境條件惡劣、*xx -= - f a (X - X ) + b ( Y -1As1A地形地貌復(fù)雜等條件下, 所以, 研究任務(wù)要求沿著高山峽谷地帶飛行。因高山峽谷地帶的氣流變化莫測(cè), 而且容易受高空風(fēng)力的影響, 造成飛行航線漂移, 使得拍攝的影像航向重疊度和旁向重疊度都不夠規(guī)

6、則。與傳統(tǒng)的航天影像和航空影像相比, 無Ys ) + c1 (ZA - Zs ) a 3 (XA -X ) + b ( Y - Y ) + c (Z - Z ) s3As3As( 1)*y -y= - f a 2 (XA - X s ) + b2 ( YA -Ys ) + c2 (ZA - Zs ) a 3 (XA -X s ) + b3 ( YA - Ys ) + c3 (ZA - Z s ) 人機(jī)低空遙感影像像幅較小、數(shù)量多、影像的傾角過大且傾斜方向沒有規(guī)律等問題 2 。根據(jù)上述特點(diǎn), 通過分析無人機(jī)低空遙感影像的幾何變形原理, 可以對(duì)無人機(jī)遙感圖像進(jìn)行圖像處理。2. 1 無人機(jī)低空遙感影

7、像幾何變形分析在無人機(jī)遙感中, 大多使用普通相機(jī), 攝影測(cè)量過程從航空攝影開始, 直至獲得像片或模型坐標(biāo)的整個(gè)數(shù)據(jù)獲取過程中, 都會(huì)帶來許多系統(tǒng)誤差。系統(tǒng)誤差主要有:其中a1 = cosU sin k- sinU sin X sink a2 = - cosU sin k- sinUsin X cosk a3 = - sinU cos Xb1 = cos X sin kb2 = cos X cos kb = - sin X3c1 =c2 =c3 =sin U cos k+ cos U sin X sin k- sin U sin k+ cos U sin X cos kcos U cos X(

8、1)攝的畸變差。( 2)攝影材料的系統(tǒng)變形、軟片的壓平誤差。( 3)地球曲率和大氣折光。( 4)量測(cè)儀器的準(zhǔn)系統(tǒng)誤差。( 5)觀測(cè)員的系統(tǒng)誤差。作者在本文采用的無人機(jī)低空遙感影像, 是采x*式中 x、y 是以像主點(diǎn)為原點(diǎn)的像點(diǎn)坐標(biāo);、y*是像點(diǎn)坐標(biāo)偏離真值的誤差; f 為像片主距;XA、YA、ZA 是相應(yīng)地面點(diǎn)的坐標(biāo); a i、bi、ci ( i= 1,2, 3)是三個(gè)外方位角元素 U、X、J 求得的方向余弦用面陣 CCD (耦合器件 )作為感光器件, 所獲圖 1 北川縣城和青川東河口震后影像 (拍攝時(shí)間: 2008- 5- 16)F ig. 1 Be ichuan county and Q

9、ingchuan estuary ea rthquake im ages ( taken tim e : 2008- 5- 16)20832卷物探化探計(jì)算技術(shù)值; X 、Y、Z 是攝影中心的坐標(biāo)值 3 。為 4個(gè) 6個(gè), 并使其盡可能取在每幅影像中sss式( 1)不僅表示了無人機(jī)遙感影像的構(gòu)像關(guān)系, 同時(shí)也指明了影像幾何變形主要來源于與傳感器本身相 內(nèi)部變形誤差, 以及由傳感器以外因素引入的外部變形誤差, 為影像的幾何糾正提供了理論依據(jù)。2. 2 遙感影像拼接與幾何校正通過上述幾何變形原理的分析, 結(jié)合無人機(jī)影像傾角大而無規(guī)律, 圖像中必然引入了無規(guī)律的幾何變形, 所以對(duì)無人機(jī)遙感圖像進(jìn)行幾何

10、糾正是必須的。通過理論分析, 作者采用多項(xiàng)式法對(duì)無人機(jī)遙感影像中外方位元素引入誤差進(jìn)行幾何糾正。該的基本思想是, 不考慮影像成像過程中的空間幾何, 直接對(duì)影像本身進(jìn)行數(shù)學(xué)模擬。它把地均勻分布, 在地物變化復(fù)雜的地區(qū)適當(dāng)加密。圖 2為安縣茶坪河經(jīng)過幾何糾正后的圖像。圖 2 安縣茶坪河局部幾何糾正圖影像的總體變形看作是、縮放、旋轉(zhuǎn)、仿射、偏F ig. 2 Char iver loca l geom etr ic correction map扭、彎曲, 以及更次基本變形的綜合效果, 因而, 可以用一對(duì)各種傳感糾正前、后影像相應(yīng)點(diǎn)之間的坐標(biāo)個(gè)適當(dāng)?shù)亩囗?xiàng)式進(jìn)行描述。這種2. 3 遙感圖像正校正及三維可視

11、化遙感圖像三維可視化實(shí)現(xiàn)的技術(shù)和思想, 就是依據(jù) DEM 建立表面模型來顯示真實(shí)地形, 然后再將遙感影像或航片等數(shù)字影像進(jìn)行紋理疊加來顯示地表細(xì)節(jié), 充分發(fā)揮計(jì)算機(jī)圖示技術(shù)和虛擬現(xiàn)實(shí)技術(shù)的優(yōu)勢(shì), 利用遙感圖像、地理要素和文字符號(hào)標(biāo)注等多種數(shù)據(jù), 生成三維地形影像 4 。將這些可視化數(shù)據(jù)集成套 三維地形影像, 必須做到不同數(shù)據(jù)間的坐標(biāo)配準(zhǔn), 目的是將不同來源的同一地區(qū)的圖像、地理要素和文字符號(hào)等數(shù)據(jù), 轉(zhuǎn)換到同一坐標(biāo)系中 6 。作者在本文是以地形圖的地理坐標(biāo)作為配準(zhǔn)參考, 其它數(shù)據(jù)均需轉(zhuǎn)換到地形圖所在的坐標(biāo)系中。遙感圖像地質(zhì)解譯三維可視化與影像動(dòng)態(tài)分器形成影像的糾正都是適用的, 但有不同的近似程

12、度。而且, 這種不僅能用于影像對(duì)地面 (地圖)系統(tǒng)的糾正, 還可以用于不同類型影像之間的相互匹配, 能滿足計(jì)算機(jī)理的需要。, 地物變化監(jiān)測(cè)等處假設(shè)影像坐標(biāo)為( u, v), 相應(yīng)像素的地面坐標(biāo)為( x, y), 則二者之間的多項(xiàng)式糾正公式可以表示為公式( 2)。2u= a00 + a 10x + a01 y+ a20 x + a 11xy+a02 y + a 30x + a21 x y+ a 12xy +2322a y3 + ,03( 2)v= b + b x+ b y+ b x2 + b xy+0010012011析系列動(dòng)畫制作工序, 主要以下幾個(gè)內(nèi)容。b y2 + b x3 + b x2

13、y+ b xy2 +02302112( 1)遙感圖像數(shù)字處理。遙感圖像數(shù)字處理b y3 + ,03是數(shù)字區(qū)調(diào)工作中提取的重要環(huán)節(jié)。圖像處多項(xiàng)式的項(xiàng)數(shù)(即系數(shù)的個(gè)數(shù) )N, 與其階數(shù) n有固定的:理的目的, 是對(duì)原始遙感圖像進(jìn)行輻射校正、幾何校正和投影差改正等, 最終制作出統(tǒng)一規(guī)格標(biāo)準(zhǔn)的高質(zhì)量遙感圖像, 以提高數(shù)字區(qū)調(diào)地質(zhì)解譯應(yīng)用效果 5 。N = ( n + 1) ( n + 2) /2( 3)因此反求 n 階糾正多項(xiàng)式的系數(shù), 至少需要N個(gè)點(diǎn),其在影像和地面的坐標(biāo)值。若已知( 2)高精度 DEM 生成。高精度三維圖像的點(diǎn)數(shù)大于 N, 則可以利用最小二乘法原理, 求解多項(xiàng)式的系數(shù)。作者在本文中

14、, 以數(shù)字化地形圖的地理坐標(biāo)作為配準(zhǔn)參考, 圖像數(shù)據(jù)轉(zhuǎn)換到地形圖在公里網(wǎng)坐標(biāo)需要高精度的 DEM 來支撐, 在遙感圖像的正射處理中, 各像點(diǎn)的投影差改正也需要對(duì)應(yīng)點(diǎn)的高精度DEM 6 。( 3)三維飛行路線選取。根據(jù)工作區(qū)的地質(zhì)構(gòu)造復(fù)雜程度和工作需要, 按照一定的規(guī)則進(jìn)行飛系中。利用一階多項(xiàng)式進(jìn)行擬合糾正,點(diǎn)的選2期等: 無人機(jī)遙感圖像及其三維可視化在209救災(zāi)中的應(yīng)用行勘察路線部署。( 4)三維可視化系列動(dòng)畫制作。首先進(jìn)行三維飛行的參數(shù)設(shè)置, 如航高、時(shí)速、夸大系數(shù)、屏幕大小、視角設(shè)置及背景效果等; 然后再根據(jù)布置的飛行路線完成三維動(dòng)畫制作。( 5)三維可視化輸出。根據(jù)選擇的飛行路線, 逐條

15、生成遙感圖像地質(zhì)解譯三維可視化及影像動(dòng)態(tài)分析系列動(dòng)畫, 并把這些轉(zhuǎn)化成了通用動(dòng)畫所支持的格式打包, 如 Qu ickT ime M ovie格式等, 刻盤并提交給有使用。(滑坡、堰塞湖( 6)提取。圖 3 安縣茶坪河三維影像 (局部區(qū)域 )等)、地質(zhì)構(gòu)造、巖溶地貌解譯和影像判讀等。遙感圖像三維可視化, 是在高程表面模型 ( DEM)上覆蓋遙感圖像、地理要素和文字符號(hào)標(biāo)注等多種數(shù)據(jù), 生成的三維地形影像。不同類型數(shù)據(jù)的集成套合, 是以地理坐標(biāo)為組織的, 因此在套三維影像時(shí), 必須做到不同數(shù)據(jù)間的坐標(biāo)配準(zhǔn), 將同一地區(qū)不同來源的影像、地理要素和文字符號(hào)轉(zhuǎn)換到同一坐標(biāo)系中 7 。作者在本文中以地形圖

16、的地理坐標(biāo)作為配準(zhǔn)參考, 其它數(shù)據(jù)均需轉(zhuǎn)換到地形圖所在的坐標(biāo)系中。其中 DEM由地形圖上數(shù)字化得來的等高線或高程點(diǎn)生成, 因此已實(shí)現(xiàn)了與地形圖的坐標(biāo)配準(zhǔn)。無人機(jī)影像在進(jìn)行幾何糾正時(shí), 已實(shí)現(xiàn)了與地形圖間的配準(zhǔn)。這樣 DEM 和影像都是依據(jù)地形圖內(nèi)容而進(jìn)行的特征數(shù)字化或F ig. 3 Char ive r three- dmi ensiona l mi age ( loca l area)無人機(jī)影像具有數(shù)量多、相幅小的特點(diǎn), 但是由于高度和分辨率的特殊性, 要適應(yīng)抗震救災(zāi)的實(shí)時(shí)性要求, 必須在較短時(shí)間內(nèi)完成影像的快速拼接、糾正和解譯等相關(guān)工作。低空遙感能提供完整影像, 且分辨率較高, 經(jīng)過校正和

17、災(zāi)害區(qū)域、對(duì)象增強(qiáng)處理, 結(jié)合部份地理, 能滿足對(duì)堰塞湖災(zāi)害的快速定量評(píng)估與解譯要求。:堰塞湖災(zāi)害定量評(píng)估與解譯( 1)滑坡。根據(jù)遙感影像的比例尺和方向, 判別滑坡的長(zhǎng)度和寬度, 以及滑動(dòng)方向。由于影像有30% 以上的航向和旁向重疊, 通過像對(duì), 可快速提取滑坡體高度和厚度, 再通過地質(zhì)圖的資料對(duì)比, 可準(zhǔn)確得出的主要巖性。文字符號(hào)注記, 與地形圖于同一坐標(biāo)系中, 無( 2)崩塌。崩塌體長(zhǎng)度、寬度的提取, 與需進(jìn)行再次配準(zhǔn) 8 ?？紤]到模擬飛行觀察效果, 計(jì)算機(jī)處理能力, 以及編輯操作簡(jiǎn)便易行等因素, 疊合后生成的三維影像, 實(shí)現(xiàn)了三維地形影像模擬飛行的動(dòng)態(tài)觀測(cè)。對(duì)安縣茶坪河圖像進(jìn)行三維可視化

18、及實(shí)現(xiàn)飛行觀察, 示例見圖 3。在圖 3中, 河道位置清晰, 堰塞湖、滑坡、泥石流等地質(zhì)災(zāi)害可清楚解譯?；禄疽粯?。災(zāi)害區(qū)崩塌體的平面規(guī)模與其厚度有相關(guān)性, 發(fā)生崩塌所形成的倒石堆坡度37b左右, 可通過崩塌體的長(zhǎng)度推算其高度。( 3)泥石流。泥石流溝的判讀主要是通過對(duì)在溝道內(nèi)松散固體物質(zhì)的辨識(shí)獲得。通過知識(shí)庫及相關(guān)經(jīng)驗(yàn),具備爆發(fā)泥石流所需要的地形條件(溝道比降)。震后所的大量滑坡和3崩塌固體物, 提供了形成泥石流的必備條件。( 4)堰塞湖。堰塞湖是由于河道二岸滑坡 (崩塌)阻塞河道所致, 堰塞壩為阻塞河道的滑坡體。確定堰塞壩體積, 需判讀出堰塞壩的平面規(guī)模, 結(jié)合其高度進(jìn)行。建立震后災(zāi)區(qū)的

19、三維動(dòng)態(tài)影像, 可以根據(jù)堰塞體的回水位置, 確定其高程, 結(jié)合地形圖數(shù)據(jù), 計(jì)算災(zāi)害遙感取與分析的快速提在發(fā)生后, 災(zāi)區(qū)河谷地帶災(zāi)情嚴(yán)重, 形成多處堰塞湖, 后期隱患。由于災(zāi)區(qū)道路損毀情況嚴(yán)重, 交通線中斷, 無人機(jī)作為能深入了解災(zāi)情的有效, 在抗震救災(zāi)要求快速、及時(shí)、準(zhǔn)確反映受災(zāi)地區(qū)的真實(shí)情況時(shí), 可以在較短的時(shí)間內(nèi)壩置水深。通過壩體前、后的有水和無水區(qū)位提供給決策層詳細(xì)的數(shù)據(jù), 方便定 。決策的制置, 可以確定壩體的高度, 從而計(jì)算出堰塞壩的體積。堰塞湖的流域面積可通過地形圖量算, 再配合 921032卷物探化探計(jì)算技術(shù)水文以及氣象資料, 計(jì)算出有關(guān)匯流以及水位上漲。壩體的性評(píng)估, 除了對(duì)

20、堰塞體的組成物質(zhì)進(jìn)行外, 還需要深入到現(xiàn)場(chǎng)考察。另外, 由參考文獻(xiàn): 1, 等. 微型無人機(jī)遙感影像的糾偏與定位 J. 南京氣象學(xué)院學(xué)報(bào), 2005, 28( 5): 632., 在圖像與震前地形不符合的地引起的地形變化, 例如滑坡、堰于包含了高程方, 可以顯示出 2馬,. 微型無人機(jī)遙感應(yīng)用初探 J. 廣西氣象, 2005, 26: 180.塞湖, 以及由此引起的河道改變等?？梢詾楹Ａ窟b感數(shù)據(jù)的快速解譯與定量分析提供依據(jù), 可以實(shí)現(xiàn)災(zāi)害特征及感目標(biāo)點(diǎn)的解譯與空間統(tǒng)計(jì)分析, 3,與指揮林劍,. 微小型無人機(jī)遙感圖像應(yīng)用 J. 火力, 2009, 34( 7): 158., 陳紹求, 等. GI

21、S 數(shù)據(jù)融入遙感圖像理 4提供災(zāi)后恢復(fù)與重建詳細(xì)統(tǒng)計(jì)數(shù)據(jù)的撐 10 。支解的模型初探 J. 物探化探計(jì)算技術(shù), 2002, 24( 1):62., FRANO IS CAUNEAU, TH RRY RANCHN, 54結(jié)論與探討等. 應(yīng)用多種數(shù)據(jù)融合對(duì)多時(shí)相 SPOT圖像處理提取三峽工程大壩區(qū)地面動(dòng)態(tài)變化探計(jì)算技術(shù), 2005, 27( 3): 267. J . 物探化作者在本文無人機(jī)低空遙感數(shù)據(jù)在災(zāi)害監(jiān)測(cè)中的應(yīng)急應(yīng)用, 能較完整地處理震后無人 6,三維建模,. 大重疊度無人機(jī)遙感影像的機(jī)遙感影像, 處理結(jié)果清晰、生動(dòng), 并在研究 J. 測(cè)繪科學(xué), 2005, 30( 2): 36.救 7,

22、等. 地質(zhì)構(gòu)造三維可視化模災(zāi)工作中得以實(shí)現(xiàn)。數(shù)據(jù)成果為勘察工作者和決策者提供了三維的背景資料, 為其它數(shù)字化成果提型探討 J. 地質(zhì)與勘探, 2001, 37( 4): 60. 8,湖, 等. 遙感圖像三維可視化技供了可視化載體, 使遙感影像呈化、動(dòng)態(tài)化模術(shù)在西部高原區(qū)機(jī)場(chǎng)建設(shè)工程中的應(yīng)用) 以小擬顯示,可身臨其境, 高空俯視, 景觀漫哨機(jī)場(chǎng)為例 J. 物探化探計(jì)算技術(shù), 2007, 29 ( 3):260.游, 對(duì)決策者具有十分重要的意義。無人機(jī)低空遙感數(shù)據(jù)在/ 5# 120抗震救災(zāi)中, 發(fā)揮了不可替代的科學(xué)作用。同時(shí), 在特定災(zāi)害環(huán)境和復(fù)雜地理?xiàng)l件 9李冰海,. 遙感圖像數(shù)字處理新) 模擬

23、反射率圖像 J. 物探化探計(jì)算技術(shù), 2006, 29( 1): 58.下, 無人機(jī)低空遙感數(shù)據(jù)在應(yīng)用中也術(shù)問題, 例如, 數(shù)字影像內(nèi)部誤差糾正, 無出許多技點(diǎn)糾 10易美華,縮中的應(yīng)用及, 黃國祥, 等. 小波變換在遙感圖像壓實(shí)現(xiàn) J . 物探化探計(jì)算技術(shù),正, 快速自動(dòng)配準(zhǔn)與拼接,災(zāi)害的自動(dòng)提2003, 25( 3): 270.取等。在后續(xù)研究中, 建立無人機(jī)高分辨率低空遙感數(shù)據(jù)的快速處理系統(tǒng), 及災(zāi)情快速解譯和定量分析技術(shù), 可以提高無人機(jī)獲取海量遙感數(shù)據(jù)的處理速度, 更進(jìn)一步的提高數(shù)據(jù)成果的實(shí)時(shí)性和應(yīng)用性。作者簡(jiǎn)介:圖學(xué)與地理( 1978- ), 男,系統(tǒng)。, 主要專業(yè)為地# 5#Ma

24、r2010COMPUT ING TECHN IQUES FOR GEOPHYS ICAL AND GEOCHEM ICAL EXPLORAT IONFU We i2b ing, YANG Bi n, L I Q iang, et al ( State K ey Laboratory ofO il and Gas Reservoir Geology and Exp loration " Chengdu Un iversity of Technology ",opt im ization filter mand the calcu lation formu laare give

25、n. F inally, the processing of actua lmagnetic m easurement data using the above methods shows the good resu lts.K ey wor ds: magne tic anomalies; three d irections ex2 pansion; interfering factors; filter optmi izat ionChengdu 610059,). COMP UTI NG TECHNIQUES FOR GEOPH YSICAL AND GEOCHEM ICAL EXP L

26、ORATION, 2010, 32( 2): 0188The further h igh e fic iency exploitat ion of gasreservoir is very d ifficu lty because of the complexity of the spat ia l d istribut ion of single sand body in S area, Sulige G as F ield. Based on the regiona l geology, sed2 imentary facies, well logging, well drilling a

27、nd core data in S area,the study of high2resolut ion sequence stratigraphy and electrofacies in the low Sh ihezi For2 mat ion ofHE28m ember in the area was carried out by adopting the scheme of multistage sequence cyc le d i2 vision and assoc iates analysis2technique of electrofa2 cies and gu ided b

28、y the theory, technique and method of Cro sps h igh2resolution sequence stratigraphy. This paper not on ly recogn izes 2 long2term, 4 m id2term and 11 short2term base leve l cycle sequences after fine c lassification and corre lation but also estab lishes stratigraph ic corre lation framework w ith

29、H igh2resolu2 t ion in th is area. The division of sed imentary fac ies and establishm ent of subsequence correlat ion frame2 works above m entioned enhance the understanding for the disc ipline of the vertical2horizonta l d istribut ion of sand body in S area, wh ich make geologica l founda2 t ions

30、 for the p lanning of reasonab le cycling exp loita2 t ion2 development scheme in next step.K ey w ord s: we ll logging;h igh2resolution sequence stratigraphy; Su lige gas fie ld;electrofac ies ana lysis; base level cyc le; strat igraphy corre lationGEOCHEM ICALCHARACTER ISTICSAND M ETALLOGEN ICPROG

31、NO SIS O F M ENG2 J IAN LEAD2Z INC M INE AREA IN GENGMA, YUNNAN PROV INCELI W en2ku i1, YANG Gu i2bing2, TANG W en2chun2 ( 1. H unan Provinc ial Bureau of Geological Explora2tion Team 417, H engyang 421001,2. G eo2chem ica l Explorat ion Team , Sichuan Provincia l Bu2 reau of Geology and M ineral re

32、sources ,Deyang618000,). COMP UTING TECHN IQUES FORGEOPH YSICAL AND GEOCHEMICAL EXP LORA2 TION, 2010, 32( 2): 0200The research on 1: 50000 sedmi ent ofM engjian Lead2Z inc m ine area in G engma shows that the con2 tent of Au, Ag, Pb, Sb ou tclass the average content of country and Three R iver Regio

33、n. Study of ore2foenergy m inera lizat ion show that e lem ent a2noma lies of Pb, Zn, Cu, Ag, Au, Sb distribute as a Northeast2Southwest line and the h igher e lement a2 noma lies is in the centre and Northeast of the area. Itcan be divided into wt o parts of ore2foprospectarea, one is in the N orth

34、east for Pb2Cu exp lorer area, the other is in the southwest for Cu, Pb, Zn etc. ,and better area of the ore2foit cou ld be a industrial depositis in southwes,t andK ey w ord s: r iver sed iments; method of ore2for energy m inera lizat ion; m etallogen ic prognosisTH E STUDY O F EDGE EXPANS ION AND

35、OP 2 T IM IZAT ION M ED IAN FI LTER I NG M ETH 2 OD I N M AGNET IC DATA PRO CESSINGMA Guo2qing, MENG L ing2shun, DU X iao2juan, et al ( College of Geo2Exploration Science and T echnol2UNMANNED AER IAL VEH IC LE REMO TE SEN SI NG IM AGE AND TH REE2D IM EN SION2 AL V ISUAL IZAT ION APPL IED I N W EN2

36、CHUAN EARTH QUAK Eogy, Jilin Un iversity, Changchun 130026,).112COMP UTI NG TECHN IQUES FOR GEOPH YSICAL AND GEOCHEM ICAL EXP LORATION, 2010, 32( 2): 0194Based on the characteristics of magnetic anoma2 lies in the ground caused by subsurface geological bod ies, A n edge expansion m ethod of three di

37、rect ions is proposed. This method can reduce the boundary effects and make the results more reliable. In the ac2 tualmeasuremen,t magnetic data is often interfered by var ious w ith factors and the factors makes the data processing and data interpretat ion more d ifficu lt It is the first t ime tha

38、t we introduce the med ian filtering technology to magnetic data processing, and we use the new filtering method to solve th is prob lem. TheH E Le i, M IAO, TANG Shu2ya ,et a l ( 1.College of Earth Science ,Chengdu Un iversity ofTechnology, Chengdu 610059,Re lief D isaster R eduction, Sichuan m ini

39、strat ion, Chengdu 610041,2. Institute of Earthquake Ad2). COMP UT2ING TECHN IQUES FOR GEOPH YSICAL AND GEO2 CHEMICAL EXP LORATION, 2010, 32( 2): 0206Th is paper analyzes the geometric distortion prin2 c ip les of unmanned aerial veh icle ( UAV ) remote sensing image. A fter mak ing the treatment of

40、 mi age stitching, geometric correct ion and ejaculat ion correc2 tion, the remote sensing data of UAV is integrated w ith high2precision DEM, then the three2d mi ensionalvirtual dynam ic mof research landscape is bu ilt# 6#ABSTRACTSVol 32 No 2Through ana lyzing three2d mi ensiona l visualization im

41、2 age, and combining w ith terrain data, it atta ins the high2precision, low a lt itude remote sensing disaster in formation and eva luates the disaster situa tion aboutthe dammed lakes and assoc iated hazards which caused by changes in topography such as landslidesproposedfeasib le and it can impro

42、ve the testeffic iency as well as the leve l of automation.K ey words: au tomatic test software; m aster stat ion software; slave station software; IP instrument; ether2 netand ava lanches in Cha2river of An county inTHE M ETHOD OF GRD GR ID I N PACK I NG PROF ILE BY GRADED COLOR BASED ON SECONDARY

43、DEVELO PM ENT O F M APG ISLIAO Gu i2xiang1, LI Zhen2hu 2i ( 1. Institute of G eo2physica l and G eochem ical Exploration, CAGS, Lang2W en2chuan earthquake2h it area. The fact shows that in the disaster emergency and the comp licated terrain situation, data collect ion and inform ation extraction can

44、 dep ict the true environmen t exactly by using the technology of UAV low2a lt itude remote sensing and three2dmi ensional visua lizat ion. Then it can provideH ebe i 065000, tute of CAPF, Lang2. Gold G eological Inst i2H ebei 065000,). COM 2decision2m akers w ith accurate,ed,tmi e ly andP UTING TEC

45、HN IQUES FOR GEOPH YSICAL AND GEOCHEMICAL EXP LORATION,2010,32 ( 2 ):0217first2hand informat ion, and m inim ize the losses.K ey words: unmanned aeria l vehicle;remote sens2 W enchuaning;three2d mi ensionalvisualization;The paper introducesera l pr inciple ofearthquakeGRD gridd ing method in draw in

46、g profiles map by gra2ded color filling. Itspprogram SaveDemToSurfG rd( ) functionprocess is utilizing d byMAPGISDEVELOPM ENTOFAUTOM AT ICTESTSO FTW ARE FOR M ULT I2FREQUENC Y IP RE2 CE IVER BASED ON TH E ETHERNETYAO H ong2chun1, CHEN Ru2jun1, 2, 3, LIU Sheng2li1( 1. School of Info2physics and Geom

47、atic Engineering, Cen tra l South Un iversity, Changsha 410083,2. Post2doctor Station of Computer Sc ience and Tech2 nology, Centra l South U niversity, Changsha 410083,3. Post2doctor Research Station of BGP,on secondary development p latform based on V isual Basic to rea lize the converting survey

48、data to normal GRD grid da ta, and converting GRD grid data toMSI image file in image processing module of MAPG IS. The advantage and lim itat ion of the method is also discussed in th is paper.K ey wor ds: MAPGIS; secondary developmen;t V is2ua l Basic; GRD grid data; profiles mapCPNC,Zhuozhou 072751,). COMP UTINGTECHN IQUES FOR GEOPH YSICAL AND GEO 2 CHEMICAL EXPLORATION, 2010, 32( 2): 0211In order to improve test e ficiency of pseudo2ran2 dom multi2frequency IP receiver, we designed testing software based on E thernet and the characteristics of the receiver.