遙感數(shù)據(jù)圖像處理實(shí)驗(yàn)三、遙感圖像的幾何校正與裁剪

實(shí)驗(yàn)三、遙感圖像的幾何校正與裁剪實(shí)驗(yàn)內(nèi)容:1.圖像分幅裁剪(SubsetImage2.圖像幾何校正(GeometricCorrection3.圖像拼接處理(MosaicImgaes4.生成三維地形表面(3DSurfacing1.圖像分幅裁剪在實(shí)際工作中,經(jīng)常需要根據(jù)研究工作范圍對(duì)圖像進(jìn)行分幅裁剪,按照ERDASIMAGINE8.4實(shí)現(xiàn)圖像分幅裁剪的過程,可以將圖像分幅裁剪為兩類型:規(guī)則分幅裁剪,不規(guī)則分幅裁剪。1.1規(guī)則分幅裁剪(以c:\ProgramFile\IMAGINE8.4\examples\lanier.img為例規(guī)則分幅裁剪是指裁剪圖像的范圍是一個(gè)矩形,通過左上角和右上角兩點(diǎn)的坐標(biāo)可以確定圖像的裁剪位置,過程如下:方法一:→ERDASIMAGINE8.4圖標(biāo)面板菜單條:Main→DataPreparation(或單擊ERDASIMAGINE8.4圖標(biāo)面板工具條“DataPrep”圖標(biāo)→打開DataPreparation對(duì)話框→單擊SubsetImage按鈕,打開Subset對(duì)話框在Subset對(duì)話框中需要設(shè)置下列參數(shù):→輸入文件名(InputFile:lanier.img→輸出文件名(OutputFile:lanier_sub.img→坐標(biāo)類型(CoordinateType:Map→裁剪范圍(SubsetDefinition:ULX、ULY、LRX、LRY(注:ULX,ULY是指左上角的坐標(biāo),LRX,LRY是指右上角的坐標(biāo),缺省狀態(tài)為整個(gè)圖像范圍→輸出數(shù)據(jù)類型(OutputDataType:Unsigned8Bit→輸出文件類型(OutputLayerType:Continuous→輸出統(tǒng)計(jì)忽略零值:IgnoreZeroInOutputStats→輸出像元波段(SelectLayers:2,3,4→OK(關(guān)閉Subset對(duì)話框,執(zhí)行圖像裁剪方法二:→ERDASIMAGINE8.4圖標(biāo)面板菜單條:Main→StartIMAGINEViewer(或單擊RDASIMAGINE8.4圖標(biāo)面板工具條“Viewer”圖標(biāo)→打開一個(gè)二維視窗→單擊視窗工具條最左端的“打開文件”圖標(biāo)→打開SelectLayerToAdd對(duì)話框在SelectLayerToAdd對(duì)話框完成以下設(shè)置:→LookIn:examples→FileName:lanier.img→Filesoftype:IMAGINEImage→雙擊OK按鈕→在二維視窗中打開lanier.img文件→單擊ERDASIMAGINE8.4圖標(biāo)面板工具條“DataPrep”圖標(biāo)→打開DataPreparation對(duì)話框→單擊SubsetImage按鈕→打開Subset對(duì)話框在Subset對(duì)話框中需要設(shè)置下列參數(shù):→輸入文件名(InputFile:lanier.img→輸出文件名(OutputFile:lanier_sub.img→坐標(biāo)類型(CoordinateType:Map→輸出數(shù)據(jù)類型(OutputDataType:Unsigned8Bit→輸出文件類型(OutputLayerType:Continuous→輸出統(tǒng)計(jì)忽略零值:IgnoreZeroInOutputStats→輸出像元波段(SelectLayers:2,3,4→單擊FromInquireBox按鈕→打開InvalidCoordinateType對(duì)話框→單擊Continue→在顯示圖像文件lanier.img視窗中單擊工具條的“+”按鈕,打開InquireCursor對(duì)話框,在視窗中移動(dòng)十字光標(biāo),確定裁剪范圍左上角和右下角,讀取其坐標(biāo)分別填入SubsetImage對(duì)話框的ULX,ULY中和LRX,LRY中→單擊OK按鈕(關(guān)閉Subset對(duì)話框,執(zhí)行圖像裁剪方法三:首先在視窗中打開lanier.img文件→AOI→Tools打開AOI工具面板→單擊矩形框確定裁剪范圍→File→Save→AOILayerAs→打開SaveAOIAs對(duì)話框,輸入文件名:2→單擊OK(退出SaveAOIAs對(duì)話框→單擊ERDASIMAGINE8.4圖標(biāo)面板工具條“DataPrep”圖標(biāo)→打開DataPreparation對(duì)話框→單擊SubsetImage按鈕→打開Subset對(duì)話框在Subset對(duì)話框中需要設(shè)置下列參數(shù):→輸入文件名(InputFile:lanier.img→輸出文件名(OutputFile:lanier_sub.img→坐標(biāo)類型(CoordinateType:Map→輸出數(shù)據(jù)類型(OutputDataType:Unsigned8Bit→輸出文件類型(OutputLayerType:Continuous→輸出統(tǒng)計(jì)忽略零值:IgnoreZeroInOutputStats→輸出像元波段(SelectLayers:2,3,4→單擊AOI按鈕→打開ChooseAOI對(duì)話框→在ChooseAOI對(duì)話框作如下設(shè)置:→AOISource:File→AOIFile:2→單擊OK(退出ChooseAOI對(duì)話框→單擊OK(退出Subset對(duì)話框,執(zhí)行圖像裁剪→單擊OK(退出Modeler對(duì)話框,完成圖像裁剪1.2不規(guī)則分幅裁剪不規(guī)則分幅裁剪是指裁剪圖像的邊界范圍是個(gè)任意多邊形,無法通過左上角和右下角兩點(diǎn)的坐標(biāo)確定圖像的裁剪位置,而必須事先生成一個(gè)完整的閉合多邊形區(qū)域,可以是一個(gè)AOI多邊形,也可以是ArcInfo的一個(gè)PolygonCoverage,針對(duì)不同的情況采用不同的裁剪過程。(一AOI多邊形裁剪(以c:\ProgramFile\IMAGINE8.4\examples\lanier.img為例→在視窗中打開需要裁剪的圖像→應(yīng)用AOI工具繪制多邊形可以將多邊形AOI文件保存在文件中(*.aoi,也可以暫時(shí)不退出視窗,將圖像與AOI多邊形保留在視窗中,然后進(jìn)行以下操作:→ERDASIMAGINE8.4圖標(biāo)面板菜單條:Main→DataPreparation(或單擊ERDASIMAGINE8.4圖標(biāo)面板工具條“DataPrep”圖標(biāo),打開DataPreparation對(duì)話框→單擊SubsetImage按鈕,打開Subset對(duì)話框在Subset對(duì)話框中需要設(shè)置下列參數(shù):→輸入文件名(InputFile:lanier.img→輸出文件名(OutputFile:lanier_sub.img→CoordinateType:Map→Output:Unsigned8bit→Oupput:Continuous→選定IgnoreZeroinOutputStats→SelectLayers:2,3,4→單擊AOI按鈕,打開ChooseAOI對(duì)話框→選定File→選擇文件:2.aoi→單擊OK(退出ChooseAOI對(duì)話框→單擊OK(退出Subset對(duì)話框,執(zhí)行圖像裁剪→單擊Modeler對(duì)話框中的OK按鈕完成圖像裁剪(二ArcInfo多邊形裁剪(以c:\ProgramFile\IMAGINE8.4\examples\gy_resample.img為例如果是按照行政區(qū)劃邊界或自然區(qū)劃邊界進(jìn)行圖像的分幅裁剪,往往是首先利用ArcInfo或者ERDAS的Vector模塊繪制精確的邊界多邊形(Polygon,然后以ArcInfo的Polygon為邊界進(jìn)行圖像裁剪。對(duì)于這種情況,需要調(diào)用ERDAS其它模塊的功能分兩步完成。第一步:將ArcInfo多邊形轉(zhuǎn)換為柵格圖像文件打開VectorToRaster對(duì)話框有兩種方法:方法一:ERDASIMAGINE8.4圖標(biāo)面板菜單條:Main→ImageInterpreter(或單擊ERDASIMAGINE8.4圖標(biāo)面板工具條“Interpreter”圖標(biāo)→打開ImageInterpreter對(duì)話框→單擊Utilities→打開Utilities對(duì)話框→選擇VectorToRaster→打開VectorToRaster對(duì)話框方法二:ERDASIMAGINE8.4圖標(biāo)面板菜單條:Main→Vector(或單擊ERDASIMAGINE8.4圖標(biāo)面板工具條“Vector”圖標(biāo)→打開VectorUtilities對(duì)話框→選擇VectortoRaster→打開VectortoRaster對(duì)話框本例中選擇前一種方法:在VectortoRaster對(duì)話框中設(shè)置下列參數(shù):→InputVectorFile(輸入矢量文件名稱:boundgy.shp→VectorType(確定矢量文件類型:polygon→UseAttributeasValue(使用矢量屬性值:GYBOUND_ID→OutputFile(輸出柵格文件名稱:raster.img→DataType(柵格數(shù)據(jù)類型:Unsigned8bit→LayerType(柵格文件類型:Thematic→SizeDifinition(轉(zhuǎn)換范圍大小:ULX,ULY,LRX,LRY(缺省條件下為整個(gè)圖像范圍→Units(坐標(biāo)單位:Meters→CellSize(輸出像元大小:X:30/Y:30→選擇SquareCell(正方形像元→單擊OK(關(guān)閉VectorToRaster對(duì)話框,執(zhí)行矢量到柵格的轉(zhuǎn)換第二步:通過掩模運(yùn)算(Mask實(shí)現(xiàn)圖像的不規(guī)則裁剪ERDASIMAGINE8.4圖標(biāo)面板菜單條:Main→ImageInterpreter(或單擊ERDASIMAGINE8.4圖標(biāo)面板工具條“Interpreter”圖標(biāo)→ImageInterpreter→Utilities→Mask→打開Mask對(duì)話框在Mask對(duì)話框中設(shè)置下列參數(shù):→InputFile(輸入需要裁剪的圖像文件名稱:gy_resample.img→InputMaskFile(輸入掩模文件名稱:raster.img→點(diǎn)擊SetupRecode按鈕設(shè)置裁剪區(qū)域內(nèi)新值(NewValue為1,區(qū)域外取0值→選擇Intersection(確定掩模區(qū)域作交集運(yùn)算→OutputFile(輸出圖像文件名稱:mask.img→Output(輸出數(shù)據(jù)類型:Unsigned8bit→單擊OK(關(guān)閉Mask對(duì)話框,執(zhí)行掩模運(yùn)算2.圖像幾何校正(GeometricCorrection幾何校正就是將圖像數(shù)據(jù)投影到平面上,使其符合地圖投影系統(tǒng)的過程;而將地圖坐標(biāo)系統(tǒng)賦予圖像數(shù)據(jù)的過程,稱為地理參考(Georeferencing。由于所有地圖投影系統(tǒng)都遵從于一定的地圖坐標(biāo)系統(tǒng),所以幾何校正過程包含了地理參考過程。在ERDASIMAGINE8.4中進(jìn)行圖像幾何校正,通常有兩種途徑啟動(dòng)幾何校正模塊:(1數(shù)據(jù)預(yù)處理途徑(以c:\ProgramFile\IMAGINE8.4\examples\lanier.img為例方法一:首先在視窗中打開lanier.img文件→ERDASIMAGINE8.4圖標(biāo)面板菜單條:Main→DataPreparation(或單擊ERDASIMAGINE8.4圖標(biāo)面板工具條“DataPrep”圖標(biāo),打開DataPreparation對(duì)話框→選定ImageGeometricCorrection選項(xiàng)→打開SetGeoCorrectionInputFile對(duì)話框場(chǎng)協(xié)調(diào)控制電梯和解救被困人員,同時(shí)將被困電梯的具體方位、電梯編號(hào)、停留樓層通知工程部和電梯公司進(jìn)行解救被困人員。4.2.監(jiān)控中心通過電視監(jiān)控屏觀察電梯內(nèi)人員情況,使用電梯應(yīng)急電話與被困人員聯(lián)系,做好解釋安慰工作。4.3.如遇特殊情況無法消除故障或被困人員中有身體嚴(yán)重不適時(shí),應(yīng)立即報(bào)消防部門或求助有關(guān)部門解決。4.4.當(dāng)值領(lǐng)班管理員主管填寫特發(fā)事件報(bào)告,做好詳細(xì)記錄故障情況及處理經(jīng)過。5.水浸處理程序5.1.當(dāng)值各崗管理員如發(fā)現(xiàn)大廈范圍內(nèi)水浸,應(yīng)立即將出事地點(diǎn)和情況報(bào)告當(dāng)值領(lǐng)班、管理員主管,同時(shí)盡快采用就近的防水設(shè)施保護(hù)好受浸樓層各電梯口,以免電梯受損。5.2.當(dāng)值領(lǐng)班、管理員主管接報(bào)后立即趕往現(xiàn)場(chǎng)查看情況,組織搶險(xiǎn)并及時(shí)向上級(jí)領(lǐng)導(dǎo)(物業(yè)主任、物業(yè)經(jīng)理和公司領(lǐng)導(dǎo)及工程部匯報(bào)情況。5.3.設(shè)法查明浸水原因,并采取有效阻截措施,如水浸是來自大廈外的暴雨洪水,應(yīng)當(dāng)在各低于水位的出入口使用備用攔水閘板和沙包;如水浸是來自市政地下水反溢,應(yīng)當(dāng)暫時(shí)將反溢的地下水道能往大廈的入口封閉,并用排水水泵將大廈的積水抽排出大廈外;如水浸來自大廈機(jī)管設(shè)施的損壞或故障,應(yīng)當(dāng)先關(guān)閉控制有關(guān)故障部位的水閘或供水泵。5.4.組織當(dāng)值人員根據(jù)水浸情況,協(xié)同工程部采取有效措施,如將電梯開高離開受浸范圍、關(guān)閉受浸區(qū)域的電閘、在水漫延的通道上擺放攔水沙包、疏通排水溝、開啟排水泵、用吸水機(jī)等,盡可能減少水浸所致的損失。5.5.水浸排除后,立即通知清潔公司人員清除積水并清理現(xiàn)場(chǎng)環(huán)境,通知工程部查明原因,修復(fù)受損的設(shè)施,盡快恢復(fù)大廈的正常運(yùn)作。5.6.當(dāng)值領(lǐng)班、管理員主管填寫特發(fā)事件報(bào)告,寫出事件詳細(xì)經(jīng)過。6.停電事故的處理程序6.1.接通知大廈將在短時(shí)間內(nèi)停電,應(yīng)在停電前10分鐘,用緊急廣播通知用戶,工程部安排電梯工,提前將電梯停在大廈首層;物業(yè)部派管理員提前10分鐘到達(dá)電梯廳維持秩序,阻止客人使用電梯,以免停電期間有客人困梯。6.2.未預(yù)知的情況下大廈突然發(fā)生停電,立即聯(lián)系工程部,盡快采取措施恢復(fù)供電。6.3.使用緊急照明,保證公共地方及主要通道的照明。6.4.派管理員到各主要電梯廳及大廈出入口維持秩序。6.5.消防監(jiān)控中心和巡樓保安員密切注視大廈各樓層,以防有人趁機(jī)制造治安問題。6.6.當(dāng)值領(lǐng)班、管理員主管填寫特發(fā)事件報(bào)告,寫出事情詳細(xì)經(jīng)過。7.業(yè)主/租戶單元內(nèi)發(fā)生刑事和治安災(zāi)害事故的處理程序。7.1.當(dāng)值管理員應(yīng)迅速向當(dāng)值領(lǐng)班、管理員主管、物業(yè)主任、物業(yè)經(jīng)理匯報(bào)或向公安部門報(bào)案,業(yè)主/租戶受侵害的財(cái)物投購(gòu)保險(xiǎn)的,由被保人通知承保的保險(xiǎn)公司。7.2.當(dāng)值領(lǐng)班、管理員主管應(yīng)根據(jù)具體情況,采取適當(dāng)?shù)姆椒ò颜麄€(gè)現(xiàn)場(chǎng)保護(hù)起來,禁止無關(guān)人員進(jìn)入現(xiàn)場(chǎng),以免破壞現(xiàn)場(chǎng)遺留的痕跡、物證、影響證據(jù)的收集。7.3.抓緊時(shí)機(jī)向發(fā)現(xiàn)人或周圍群眾了解案件、事故發(fā)生發(fā)現(xiàn)經(jīng)過,收集群眾的反映和議論,了解更多的情況并認(rèn)真記錄。7.4.向到達(dá)現(xiàn)場(chǎng)的公安人員認(rèn)真匯報(bào)案件發(fā)生情況,協(xié)助破案。7.5.當(dāng)值領(lǐng)班、管理員主管填寫特發(fā)事件報(bào)告,寫出事情詳細(xì)經(jīng)過。8.搶劫的處理程序在當(dāng)值中如遇有公開使用暴力,威迫或其他手段如打、砸、搶、強(qiáng)行掠取公司或業(yè)主/租戶財(cái)物的犯罪行為時(shí),必須立即按如下程序處理。8.1.迅速制止犯罪,設(shè)法擒獲或召集組織其他人員擒獲罪犯,同時(shí)向當(dāng)值領(lǐng)班、管理員主管、物業(yè)主任、物業(yè)經(jīng)理報(bào)告或公安部門報(bào)案。8.2.如在樓內(nèi)發(fā)生劫案,應(yīng)通知大堂管理員留意,嚴(yán)格檢查,組織人員進(jìn)行清查,如罪犯逃走,追不上時(shí),應(yīng)看清楚犯罪分子的人數(shù)、衣著、面貌明顯特征、身體特點(diǎn)、所用的交通工具及其型號(hào)、品牌、特征等,及時(shí)報(bào)告上級(jí)領(lǐng)導(dǎo)、公安部門。8.3.保護(hù)好現(xiàn)場(chǎng),劫匪遺留兇器、作案工具等不要用手觸摸,不要讓無關(guān)人員進(jìn)入現(xiàn)場(chǎng),如無固定的現(xiàn)場(chǎng)的、無法將劫匪遺留的物品留在原處,要一一收拾起來,交給公安部門處理。8.4.訪問目擊群眾收集發(fā)生劫案情況,提供給公安部門,同時(shí)公安部門未勘查現(xiàn)場(chǎng)或未勘查完畢不能離開。8.5.事主或在場(chǎng)如有受傷的,要立即設(shè)法將傷者送醫(yī)院搶救醫(yī)治,并迅速報(bào)告公安部門。8.6.當(dāng)值領(lǐng)班、管理員主管填寫特發(fā)事件報(bào)告,寫出事件詳細(xì)經(jīng)過。9.對(duì)不執(zhí)行規(guī)定不聽勸阻的處理程序在當(dāng)值過程中如遇到不登記不出示證件強(qiáng)行進(jìn)入,強(qiáng)行在不準(zhǔn)停車的地方停車,在不該堆放垃圾的地方堆放垃圾等,可按如下程序處理:9.1.處理問題時(shí)要有禮貌、態(tài)度和藹、說話和氣、以理服人。9.2.對(duì)不停勸阻者,盡量想辦法查清其姓名或單位,如實(shí)記錄下來以便向領(lǐng)導(dǎo)匯報(bào)。9.3.發(fā)生糾紛,管理人員一定要保持冷靜,避免爭(zhēng)吵。9.4.若遇到個(gè)別人員蠻橫無理,打罵管理人員,可通知上級(jí)領(lǐng)導(dǎo)同此人共同協(xié)商妥善處理,若情節(jié)嚴(yán)重,報(bào)告公安部門依法處理。9.5.當(dāng)值領(lǐng)班、管理員主管填寫特發(fā)事件報(bào)告,寫出事件詳細(xì)經(jīng)過。10.對(duì)醉酒鬧事或精神病人的處理程序醉酒者或精神病人失去了正常的理智,有些處于不能自控的狀態(tài),很容易造成傷害,當(dāng)值人員應(yīng)按如下程序處理:10.1.進(jìn)行勸阻或阻攔,讓其離開大廈管轄范圍。10.2.如果熟悉或認(rèn)識(shí)醉酒者或精神病人,應(yīng)設(shè)法立即通知其家人或工作單位,請(qǐng)他們派人領(lǐng)回,并采取監(jiān)控和監(jiān)護(hù)措施。10.3.若醉酒者或精神病人有危害社會(huì)安全的行為,可將其強(qiáng)制送交公安部門,處理過程中要注意保護(hù)自己。11.在當(dāng)值中拾到遺失物品的處理程序11.1.群眾交來或拾獲群眾遺失或遺留的物品,都要一一登記清楚,并報(bào)告當(dāng)值保安主任或上級(jí)領(lǐng)導(dǎo)。11.2.處理拾遺物品,拾遺物品中有證件提供失主的詳細(xì)地址或電話的,可通知失主前來認(rèn)領(lǐng),但要問明情況,如是貴重物品,必須交由部門領(lǐng)導(dǎo)處理。11.3.如拾遺物品是屬于違禁品的登記后上報(bào)部門、公司領(lǐng)導(dǎo)與公安部門協(xié)調(diào)處理。11.4.對(duì)拾遺物品,必須報(bào)告部門領(lǐng)導(dǎo),不得拾物不報(bào)、少報(bào)或私自挪用,違者重處,情節(jié)嚴(yán)重的以貪污論處。(九、常見的幾種可疑情況及處理措施:加強(qiáng)巡邏檢查是治安防范工作中預(yù)防、發(fā)現(xiàn)和打擊違法犯罪分子的一項(xiàng)行之有效的措施。而對(duì)可疑情況視級(jí)別及采取相應(yīng)的措施是巡邏工作的重點(diǎn),以下幾種常見可疑情況及處理方法:1.可疑情況:在大廈內(nèi)游蕩或借口找人說不出被訪人姓名及所在單元和公司名稱,或接二連三的進(jìn)去單元內(nèi)的。處理方法:密切注意其舉動(dòng),必要時(shí)勸其離開。2.可疑情況:發(fā)現(xiàn)其身上帶有鐵棒、刀具、鉗子、螺絲刀等工具的。處理方法:檢查其扔帶工具的用途,如用途不明的,先約束起來,送公安部門。3.可疑情況:攜帶物品繁多,如電腦、打印機(jī)等貴重物品,又無任何證明的(搬家有管理處放行條的除外。處理方法:暫時(shí)將人、物扣留,待其出具可靠證明后放行,如無任何證明,即通知部門領(lǐng)導(dǎo)或公安部門審查。4.可疑情況:在偏避、隱蔽處清理皮包或錢包的處理方法:立即設(shè)法攔截,詢問證實(shí),如屬盜竊、搶劫財(cái)物的,報(bào)告部門領(lǐng)導(dǎo)或送交公安部門處理。5.可疑情況:單車、摩托車無牌、無證行駛、無鋼印、有耗損痕跡的,或?qū)⑽撮_鎖的單車背走或提走的處理方法:當(dāng)即扣留人、車,待查明后方行。6.可疑情況:遇到保安即轉(zhuǎn)身遠(yuǎn)離或逃跑的人。處理問題:設(shè)法阻截(用對(duì)講機(jī)向各崗位通告及當(dāng)值保安主任報(bào)告擒獲。7.可疑情況:低價(jià)出售物品與實(shí)際價(jià)值相差較大的。視窗菜單條:Raster→GeometricCorrection→打開SetGeometricModel對(duì)話框資源衛(wèi)星圖像校正具體過程:①視窗采點(diǎn)模式:如果已經(jīng)擁有需要校正圖像區(qū)域的數(shù)字圖像、或經(jīng)過校正的圖像、或注記圖層,就可以應(yīng)用視窗采點(diǎn)模式,直接以數(shù)字地圖、或經(jīng)過校正的圖像、或注記圖層作為地理參考,在視窗中打開該圖層,從中采集控制點(diǎn)。(以c:\ProgramFile\IMAGINE8.4\examples下的panAtlanta.img和tmAtlanta.img為例,前者是作為地理參考的圖像,后者為需要校正的圖像第一步:顯示圖像文件→在ERDASIMAGINE8.4圖標(biāo)面板工具條點(diǎn)擊“Viewer”圖標(biāo)兩次,打開兩個(gè)視窗→ERDASIMAGINE8.4圖標(biāo)面板菜單條:Session→TileViewer然后,在Viewer#1中打開需要校正的LandsatTM圖像:tmAtlanta.img在Viewer#2中打開作為地理參考的校正過的SPOT圖像:panAtlanta.img第二步:啟動(dòng)幾何校正模塊Viewer#1菜單條:Raster→GeometricCorrection→打開SetGeometricModel對(duì)話框→選擇多項(xiàng)式幾何校正計(jì)算模型:Polynormial→單擊OK按鈕→同時(shí)打開GeoCorrectionTools對(duì)話框和PolynormialModelProperties對(duì)話框在PolynormialModelProperties對(duì)話框中,定義多項(xiàng)式模型參數(shù)及投影參數(shù)→定義多項(xiàng)式次方(Ploynomial:2,其余取缺省值→單擊Apply按鈕→單擊Close按鈕→打開GCPToolReferenceSetup對(duì)話框說明:該實(shí)例是采用視窗采點(diǎn)模式,作為地理參考的SPOT圖像已經(jīng)含有投影信息,所以這里不需要定義投影參數(shù)。如果不是采用視窗采點(diǎn)模式,或者參考圖像沒有包含投影信息,則必須在這里定義投影信息,包括投影類型及其對(duì)應(yīng)的投影參數(shù)。第三步:啟動(dòng)控制點(diǎn)工具在GCPToolReferenceSetup對(duì)話框中→選擇視窗采點(diǎn)模式:ExistingViewer→單擊OK按鈕(關(guān)閉GCPToolReferenceSetup對(duì)話框→打開ViewerSelectionInstruction指示器→在顯示作為地理參考圖像panAtlanta.img的Viewer#2中點(diǎn)擊左鍵→打開ReferenceMapInformation提示框→單擊OK按鈕(關(guān)閉ReferenceMapInformation提示框此時(shí)整個(gè)屏幕將顯示:兩個(gè)主視窗、兩個(gè)放大窗口、兩個(gè)關(guān)聯(lián)方框(分別位于兩個(gè)視窗中,指示放大窗口與主視窗的關(guān)系、控制點(diǎn)工具窗口、幾何校正工具。表明控制點(diǎn)工具啟動(dòng),進(jìn)入控制點(diǎn)采集狀態(tài)第四步:采集地面控制點(diǎn)控制點(diǎn)工具窗口簡(jiǎn)介:控制點(diǎn)工具窗口(GCPTool有菜單條、工具條、控制點(diǎn)數(shù)據(jù)表和狀態(tài)條四個(gè)部分組成。GCPTool窗口菜單命令及其功能命令功能File:LoadInputSaveInputSaveInputAsLoadReferenceSaveReferenceSaveReferenceASClose文件操作:調(diào)用輸入控制點(diǎn)文件(*.gcc將輸入控制點(diǎn)保存在圖像中保存輸入控制點(diǎn)文件(*.gcc調(diào)用參考控制點(diǎn)文件(*.gcc將參考控制點(diǎn)保存到圖像中保存參考控制點(diǎn)文件(*.gcc關(guān)閉控制點(diǎn)工具View:顯示操作:ViewOnlySelectedGCP’sShowSelectedGCPinTableArrangeFramesonScreenToolsStartChipViewer視窗僅顯示所選擇的控制點(diǎn)在表格顯示所選擇的控制點(diǎn)重新排列屏幕中的組成要素調(diào)出控制點(diǎn)工具圖標(biāo)面板重新打開放大窗口Edit:SetPointType(Control/CheckResetReferenceSourceReferenceProjectionPointPredictionPointMatching編輯操作:設(shè)置采集點(diǎn)的類型(控制/檢查改變參考控制點(diǎn)源文件改變參考文件的投影參數(shù)按照轉(zhuǎn)換方程計(jì)算下一點(diǎn)位置匹配兩幅圖像像元的灰度值Help:HelpforGCPTool聯(lián)機(jī)幫助:關(guān)于GCP工具的聯(lián)機(jī)幫助GCPTool窗口工具圖標(biāo)命令及其功能(按圖標(biāo)從左到右的順序命令功能Toggle自動(dòng)GCP編輯模式開關(guān)鍵Calculate依據(jù)控制點(diǎn)求解幾何校正模型Automatic設(shè)置自動(dòng)轉(zhuǎn)換計(jì)算開關(guān)ComputeError計(jì)算檢查點(diǎn)的誤差,更新RMS誤差SelectGCP激活GCP選擇工具,在視窗中選擇GCPCreateGCP在視窗中選擇定義GCPLock/Unlock鎖住/釋放當(dāng)前命令,以便重復(fù)使用FindinInput選擇尋找輸入圖像中的GCPFindinRefer選擇尋找參考文件中的GCPZValue計(jì)算更新所選GCP的Z值A(chǔ)uto-ZValue自動(dòng)更新所有GCP的Z值GCP窗口數(shù)據(jù)字段及其含義字段含義Point#GCP順序號(hào),系統(tǒng)自動(dòng)產(chǎn)生PointIDGCP標(biāo)識(shí)碼,用戶可以定義>GCP當(dāng)前選擇狀態(tài)提示符號(hào)Color輸入GCP顯示顏色XInput輸入GCP的X坐標(biāo)YInput輸入GCP的Y坐標(biāo)Color參考GCP顯示顏色XReference參考GCP的X坐標(biāo)YReference參考GCP的Y坐標(biāo)TypeGCP的類型(控制點(diǎn)/檢查點(diǎn)XResidual單個(gè)GCP的X殘差YResidual單個(gè)GCP的Y殘差RMSError單個(gè)GCP的RMS誤差Contribution單個(gè)GCP的貢獻(xiàn)率Match兩幅圖像像元灰度值的匹配程度關(guān)于GCPTool窗口的說明:1.輸入控制點(diǎn)(InputGCP是在原始文件視窗中采集的,具有原文件的坐標(biāo)系統(tǒng);而參考控制點(diǎn)(ReferenceGCP是在參考文件視窗中采集的,具有已知的參考坐標(biāo)系統(tǒng),GCP工具將根據(jù)對(duì)應(yīng)點(diǎn)的坐標(biāo)值自動(dòng)生成轉(zhuǎn)換模型。2.在GCP數(shù)據(jù)中,殘差(Residuals、中誤差(RMS、貢獻(xiàn)率(Contribution及匹配程度(Match等參數(shù),是在編輯GCP的過程中自動(dòng)計(jì)算更新的,用戶是不能任意改變的,但可以通過精確GCP位置來調(diào)整3.每個(gè)IMG文件都可以有一個(gè)GCP數(shù)據(jù)集與之相關(guān)聯(lián),GCP數(shù)據(jù)集保存在一個(gè)柵格層數(shù)據(jù)文件中;如果IMG文件有一個(gè)GCP數(shù)據(jù)集存在的話,只要打開GCPTool,GCP點(diǎn)都會(huì)出現(xiàn)在視窗中4.所有的輸入控制點(diǎn)(InputGCP都可以保存在圖像文件中(SaveInput,也可以保存在控制點(diǎn)文件中(SaveInputAs。如果保存在圖像文件中,調(diào)用的方法如3所述,如果是保存在控制點(diǎn)文件中,可以通過加載調(diào)用(LoadInput。5.參考控制點(diǎn)(ReferenceInput也可以類似地保存在參考圖像中(SaveReference或控制點(diǎn)文件中(SaveReferenceAs,可以通過加載調(diào)用(LoadReference采集地面控制點(diǎn)具體過程如下:→在GCPTool窗口中點(diǎn)擊SelectGCP圖標(biāo),進(jìn)入GCP選擇狀態(tài)→在GCP數(shù)據(jù)表中將輸入GCP的顏色(Color設(shè)置為比較明顯的黃色→在Viewer#1中移動(dòng)關(guān)聯(lián)方框位置,尋找明顯的地物特征點(diǎn),作為輸入GCP→在GCPTool窗口中點(diǎn)擊CreateGCP圖標(biāo),并在Viewer#3中點(diǎn)擊左鍵定點(diǎn),GCP數(shù)據(jù)表將記錄一個(gè)輸入GCP,包括其編碼、標(biāo)識(shí)碼。X坐標(biāo)、Y坐標(biāo)→在GCPTool窗口中點(diǎn)擊SelectGCP圖標(biāo),重新進(jìn)入GCP選擇狀態(tài)→在GCP數(shù)據(jù)表中將參考GCP的顏色(Color設(shè)置為比較明顯的紅色→在Viewer#2中移動(dòng)關(guān)聯(lián)方框位置,尋找對(duì)應(yīng)的地物特征點(diǎn),作為參考GCP→在GCPTool窗口中點(diǎn)擊CreateGCP圖標(biāo),并在Viewer#4中點(diǎn)擊左鍵定點(diǎn),系統(tǒng)將自動(dòng)把參考點(diǎn)的坐標(biāo)(XReference,YReference顯示在GCP數(shù)據(jù)表中→在GCPTool窗口中點(diǎn)擊SelectGCP圖標(biāo),重新進(jìn)入GCP選擇狀態(tài),并將光標(biāo)移回到Viewer#1,準(zhǔn)備采集另一個(gè)輸入控制點(diǎn)→不斷重復(fù)以上步驟,采集若干GCP,直到滿足所選定的幾何校正模型為止(多項(xiàng)式變換的次方數(shù)為二時(shí),最少需要六個(gè)控制點(diǎn),而后每采集一個(gè)InputGCP,系統(tǒng)就自動(dòng)產(chǎn)生一個(gè)ReferenceGCP。第五步:采集地面檢查點(diǎn)以上所采集的GCP的類型(Type均為Control(控制點(diǎn),用于控制計(jì)算,建立轉(zhuǎn)換模型及多項(xiàng)式方程.下面所要采集的GCP的類型均是Check(檢查點(diǎn),用于檢驗(yàn)所建立的轉(zhuǎn)換方程的精度和適用性。在GCPTool窗口中→在GCPTool窗口菜單條中確定GCP類型:Edit→SetPointType→Check→在GCPTool窗口菜單條中確定GCP匹配參數(shù)(MatchingParameter:Edit→PointMatching→打開GCPMatching對(duì)話框在GCPMatching對(duì)話框,需要定義下列參數(shù):MatchingParameter(匹配參數(shù):→Max.SearchRadius(最大搜索半徑:3→SearchWindow(搜索窗口大小:X:5;Y:5ThresholdParameters(約束參數(shù):→CorrectionThreshold(相關(guān)閾值:0.8→選定DiscardUnmatchedPoint(刪除不匹配的點(diǎn)→單擊Close(關(guān)閉GCPMatching對(duì)話框→確定地面檢查點(diǎn):在GCPTool工具條中選擇CreateGCP圖標(biāo),并將Lock圖標(biāo)打開鎖住CreateGCP功能,象選擇控制點(diǎn)一樣,分別在Viewer#1和Viewer#2中定義檢查點(diǎn),定義完畢后點(diǎn)擊Unlock圖標(biāo),解除CreateGCP功能→計(jì)算檢查點(diǎn)誤差:在GCPTool工具條中點(diǎn)擊“ComputeError”圖標(biāo),檢查點(diǎn)的誤差值會(huì)顯示在GCPTool的上方,只有所有檢查點(diǎn)的誤差均小于一個(gè)像元,才能繼續(xù)進(jìn)行合理的重采樣第六步:圖像重采樣重采樣過程就是依據(jù)未校正圖像像元值計(jì)算生成一幅校正圖像的過程,原圖像中所有柵格數(shù)據(jù)層都將進(jìn)行重采樣。ERDASIMAGINE8.4提供三種最常用的重采樣方法:NearestNeighbor:鄰近點(diǎn)插值法,將最鄰近像元值直接賦予輸出像元BilinearInterpolation:雙線性插值法,用雙線性方程和2×2窗口計(jì)算輸出像元值CubicConvolution:立方卷積插值法,用立方方程和4×4窗口計(jì)算輸出像元值具體操作如下:在GeoCorrectionTools對(duì)話框中選擇“ImageResample”圖標(biāo)→打開ImageResample對(duì)話框在Resample對(duì)話框中,定義重采樣參數(shù):→OutputFile(輸出圖像文件:rectify.img→ResampleMethod(選擇重采樣方法:NearestNeighbor→OutputCellSize(定義輸出圖像范圍:ULX、ULY;LRX、LRY→選定IgnoreZeroinStats→RecalculateOutputDefault(設(shè)置重新計(jì)算缺省值:SkipFactor:X:1;Y:1→單擊Ok(關(guān)閉Resample對(duì)話框第八步:保存幾何校正模型在GeoCorrectionTools對(duì)話框中點(diǎn)擊Exit按鈕,退出圖像幾何校正過程。按照系統(tǒng)提示將圖像幾何校正模型保存為模式文件resample.gms.第九步:檢驗(yàn)校正結(jié)果檢驗(yàn)校正結(jié)果的基本方法是:同時(shí)在兩個(gè)視窗中打開兩幅圖像,其中一幅是校正后的圖像,一幅是當(dāng)時(shí)的參考圖像,通過視窗地理連接(GeoLink/Unlink功能及查詢光標(biāo)功能進(jìn)行目視定性檢驗(yàn),具體過程如下:→在ERDASIMAGINE8.4圖標(biāo)面板工具條點(diǎn)擊“Viewer”圖標(biāo)兩次,打開兩個(gè)視窗→ERDASIMAGINE8.4圖標(biāo)面板菜單條:Session→TileViewers然后,在Viewer#1中打開校正后的圖像:rectify.img在Viewer#2中打開參考圖像:panAtlanta.img→在Viewer#1中:按住右鍵→快捷菜單→GeoLink/Unlink→在Viewer#2中:點(diǎn)擊左鍵→建立與Viewer#1的連接→在Viewer#1中:按住右鍵→快捷菜單→InquireCursor→打開光標(biāo)查詢對(duì)話框→在Viewer#1中:移動(dòng)查詢光標(biāo),觀測(cè)其在兩屏幕中的位置及匹配程度,并注意光標(biāo)查詢對(duì)話框中數(shù)據(jù)的變化。如果滿意的話,關(guān)閉光標(biāo)查詢窗口。②地圖采點(diǎn)模式:如果只有硬拷貝的地圖或坐標(biāo)紙作為參考,可以采用地圖采點(diǎn)模式。首先在地圖上選點(diǎn)并量算其坐標(biāo),在地圖上尋找六個(gè)以上明顯的地物特征點(diǎn),并量算起坐標(biāo)(多項(xiàng)式變換的次方數(shù)為二時(shí),最少需要六個(gè)控制點(diǎn),然后通過鍵盤輸入坐標(biāo)數(shù)據(jù)。(以c:\ProgramFile\IMAGINE8.4\examples\nj321.img為例首先在一個(gè)二維視窗中打開nj321.img圖像文件視窗菜單條:Raster→GeometricCorrection→打開SetGeometricModel對(duì)話框→選擇多項(xiàng)式幾何校正計(jì)算模型:Polynormial(注:衛(wèi)星圖像校正過程中多采用Polynormial(多項(xiàng)式變換→單擊OK按鈕→同時(shí)打開GeoCorrectionTools對(duì)話框和PolynormialModelProperties對(duì)話框在PolynormialModelProperties對(duì)話框中,定義多項(xiàng)式模型參數(shù)及投影參數(shù):→定義多項(xiàng)式次方(Ploynomial:2,其余取缺省值→單擊Apply按鈕→單擊Close按鈕→打開GCPToolReferenceSetup對(duì)話框→選擇KeyboardOnly(通過鍵盤輸入控制點(diǎn)→單擊Ok→打開ReferenceMapInformation對(duì)話框→在MapUnits(地圖單位中選擇:Other(說明:因?yàn)楸纠胁捎玫牡貓D單位為km,在該下拉菜單中沒有列出→單擊OK→打開GCPTool對(duì)話框在顯示nj321.img圖像的視窗Viewer#1中產(chǎn)生一個(gè)關(guān)聯(lián)方框,同時(shí)產(chǎn)生一個(gè)放大窗口Viewer#2,放大窗口顯示關(guān)聯(lián)方框中圖像的放大效果?！赩iewer#1中用關(guān)聯(lián)方框查找與地圖對(duì)應(yīng)的一個(gè)控制點(diǎn),單擊GCPTool對(duì)話框工具條的“CreateGCP”圖標(biāo),在放大窗口Viewer#2中點(diǎn)擊左鍵定點(diǎn),XInput和YInput值會(huì)自動(dòng)產(chǎn)生,在XRef.和YRef.中輸入控制點(diǎn)在地圖上的坐標(biāo)→按照類似方法輸入其余的五個(gè)控制點(diǎn)的坐標(biāo)→在GCPTool窗口菜單條中確定GCP類型:Edit→SetPointType→Check→當(dāng)在放大窗口Viewer#2中點(diǎn)擊左鍵定點(diǎn)第七個(gè)控制點(diǎn)時(shí),XRef.和YRef.值會(huì)自動(dòng)產(chǎn)生?！^續(xù)在Viewer#1中尋找四個(gè)明顯的點(diǎn),在放大窗口Viewer#2中點(diǎn)擊左鍵定點(diǎn),這五個(gè)點(diǎn)作為檢查點(diǎn)→計(jì)算檢查點(diǎn)誤差:在GCPTool工具條中點(diǎn)擊“ComputeError”圖標(biāo),檢查點(diǎn)的誤差值會(huì)顯示在GCPTool的上方,只有所有檢查點(diǎn)的誤差均小于一個(gè)像元(即<30米,才能繼續(xù)進(jìn)行合理的重采樣,如果誤差大,則必須重復(fù)以上操作直到滿足要求為止?！鷨螕鬐eoCorrectionTools對(duì)話框的“Resample”圖標(biāo)→打開Resample對(duì)話框在Resample對(duì)話框中,定義重采樣參數(shù):→OutputFile(輸出圖像文件名:nj321resample.img→ResampleMethod(重采樣方法:NearestNeighbor(鄰近點(diǎn)插值法→選定IgnoreZeroinStats→其它設(shè)置均取系統(tǒng)默認(rèn)值→單擊OK(執(zhí)行重采樣3.圖像拼接處理(以c:\ProgramFile\IMAGINE8.4\examples下的wasia1_mss