基于ERDASIMAGINE的遙感影像裁切方法探討

王海葳　徐少坤

摘 要：我們?cè)贛apGIS生產(chǎn)專(zhuān)題地圖時(shí)，經(jīng)常要用到遙感影像更新地圖要素，由于MapGIS對(duì)大尺寸遙感影像裁切處理功能有限，在此我們借助第三方軟件探討基于ERDAS IMAGINE的遙感影像處理方法。ERDAS IMAGINE是ERDAS公司開(kāi)發(fā)的面向企業(yè)級(jí)的遙感影像處理軟件。以其先進(jìn)的圖像處理技術(shù)、友好的用戶界面、靈活操作方式、以及面向廣闊應(yīng)用領(lǐng)域的產(chǎn)品模塊和3S（遙感圖像處理，地理信息系統(tǒng)和全球定位系統(tǒng)）集成功能，為遙感及相關(guān)應(yīng)用領(lǐng)域的用戶提供了內(nèi)容豐富而功能強(qiáng)大的圖像處理工具。

關(guān)鍵詞：MapGIS ERDAS IMAGINE 影像裁切

中圖分類(lèi)號(hào)：P208 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1672-3791（2015）05（b）-0015-02

在MAPGIS作業(yè)生產(chǎn)中，由于多數(shù)據(jù)源的綜合應(yīng)用，我們需要大量的遙感影像作為基礎(chǔ)數(shù)據(jù)更新，而MAPGIS對(duì)遙感影像的支持有一定的局限性，研究測(cè)試表明MAPGIS對(duì)大尺寸遙感影像TIF文件支持有限，如果文件過(guò)大（約200MB），轉(zhuǎn)換的MSI文件導(dǎo)入MAPGIS中無(wú)法正常顯示。這就需要對(duì)遙感影像數(shù)據(jù)進(jìn)行裁切處理，而MAPGIS對(duì)影像沒(méi)有提供裁切功能，我們只能借助第三方軟件來(lái)進(jìn)行影像裁切，該文應(yīng)用ERDAS IMAGINE軟件對(duì)遙感影像進(jìn)行裁切處理。（該文所用軟件版本MAPGIS6.7、ERDAS IMAGINE 9.2）。

1 遙感影像裁切方法

1.1 AOI感興趣區(qū)域法

AOI可以實(shí)現(xiàn)影像的任意裁切或AOI的批量裁切，其具體方法如下。

第一步，打開(kāi)遙感影像后，選擇AOI菜單下Tools工具欄利用多邊形或矩形選取工具創(chuàng)建要裁切感興趣區(qū)域。

第二步，按Shift鍵選中多個(gè)AOI區(qū)域，單擊Raster菜單下的Mosaic Image項(xiàng)，在Output Image Options，選擇裁切模式為“User-defined AOI”，具體操作如圖所示：

第三步，單擊Output Image Set Output對(duì)話框中Set Output AOI按鈕執(zhí)行Choose AOI對(duì)話框并激活文件Viewer1窗口，在Mosaic Tool對(duì)話框中列表項(xiàng)Output中按Shift鍵選中列表中多個(gè)AOI區(qū)域，選擇文件菜單項(xiàng)Process中Run Mosaic即可批量裁切輸出。

AOI方法可以任意裁切一個(gè)或多個(gè)區(qū)域文件，如果要精確裁切多個(gè)鄰接區(qū)域，我們可借助ERDAS的精確視圖工具，單擊Utility菜單項(xiàng)下的Inquire box出現(xiàn)以下如圖2所示。

在圖2所示的對(duì)話框中，我們可以輸入所要裁切的影像地理坐標(biāo)范圍，其中（ULX，ULY）為圖幅所在的左上角點(diǎn)坐標(biāo)，（LRX，LRY）為圖幅所在的右下角點(diǎn)坐標(biāo)，我們可以根據(jù)這個(gè)顯示區(qū)域創(chuàng)建多個(gè)AOI區(qū)域，來(lái)實(shí)現(xiàn)批量裁剪多個(gè)鄰接區(qū)域。文件輸出的格式選擇tif文件格式，以便我們?cè)贛APGIS中直接導(dǎo)入轉(zhuǎn)換MSI文件。

創(chuàng)建AOI除運(yùn)用Moscaic工具進(jìn)行裁切外，我們還可以直接編寫(xiě)AOI裁切影像的批處理文件*.bcf，文件編寫(xiě)格式如下：

modeler -nq subset.pmdl -meter -state"$（Input）"Float 1"$（Output）" Float_Single Float '$（Input.ULX）''$（Input.ULY）''$（Input.LRX）''$（Input.LRY）' Map ignore Continuous"Continuous" "Continuous" "$（aoi）" ''

編寫(xiě)好的*.bcf我們可以在菜單項(xiàng)Session中打開(kāi)Open Batch command file文件調(diào)用所編寫(xiě)的bcf文件即可實(shí)現(xiàn)批量裁切。

1.2 建立ASCII文件實(shí)現(xiàn)影像的批量裁切

建立一個(gè)txt分幅文本，可以自己手動(dòng)輸入，也可以用編程實(shí)現(xiàn)。第一行是分幅的幅號(hào)，任意設(shè)定，后面是每幅的左上角和右下角坐標(biāo)。其格式為（“圖幅號(hào)”“ULX，ULY”“LRX，LRY”）中間用空格隔開(kāi)。如圖3所示。

分幅文件編寫(xiě)好后，在圖1中的Mosaic界面上，選中Output Image Options，選擇裁切模式為“ASII Sheet File”，文件的路徑指向所建立的txt文件進(jìn)行輸出即可。

1.3 矢量數(shù)據(jù)分幅裁切

矢量數(shù)據(jù)裁切主要是用*.shp文件的Polygon文件類(lèi)型裁切，這里如果是經(jīng)緯度裁切（方里網(wǎng)裁切或其他面域裁切與此相同），我們可以在MAPGIS中將生成的經(jīng)緯度文件拓?fù)錇槊嫖募⑥D(zhuǎn)換生成shp格式，在Output Image Options中，選擇裁切模式為Polygon Vector File，文件路徑指向所要裁切的shp文件，輸出選項(xiàng)選中“Muliple Files”文件即可實(shí)現(xiàn)批量裁切多個(gè)文件。具體操作如圖4所示。

2 幾點(diǎn)注意事項(xiàng)

（1）由于Erdas Imagine軟件對(duì)中文字符支持不是很好（最新的Edas Imagine2010也是如此），所以在文件夾及文件命名時(shí)盡量避免使用中文字符。

（2）當(dāng)做大量影像裁切時(shí)，有時(shí)會(huì)出現(xiàn)“efio_Ftruncate failed”錯(cuò)誤提示，這是由于缺少臨時(shí)空間來(lái)處理臨時(shí)文件，裁切時(shí)需至少所有裁切影像的2.5倍的臨時(shí)處理空間來(lái)處理，臨時(shí)空間位置的確定可以在Session-Preferences-User-Interface Session-Default Output Directory進(jìn)行設(shè)置，比如設(shè)置為D：＼temp。這樣才會(huì)保證文件正常輸出。

（3）在影像做分幅裁切前，我們可以先將多幅影像拼接為一幅大的影像（Eerdas支持遙感影像TIF格式大小4G，金字塔重建Image格式48TB以下），利用圖1的Masaic Tool面板一次添加多幅影像，拼接輸出Image文件格式，然后在運(yùn)用以上三種方法批量分幅裁切。這樣可以大大提高分幅輸出的效率。

參考文獻(xiàn)

[1] 孫家柄.遙感原理與應(yīng)用[M].武漢：武漢大學(xué)出版社，2009.

[2] 賈永紅.數(shù)字圖像處理[M].武漢：武漢大學(xué)出版社，2010

[3] 黨安榮.ERDAS IMAGINE遙感圖像處理方法[M].北京：清華大學(xué)出版社，2003.