沒有原始坐標(biāo)的圖片的幾何校正方法

如果你有一副照片想給其賦予坐標(biāo)并想求其中的面積、周長等等屬性，那么下面的步驟值得你細(xì)看，希望對你們有所幫助。一、準(zhǔn)備步驟1.準(zhǔn)備一張方格紙，最好面積是照片面積的2倍以上2.把照片放在方格紙的正中央，然后利用數(shù)碼相機拍下照片，用于幾何校正注意事項：如果有很多照片需要比對的，最好焦距放到一定的大小，也就是能放大到同樣的倍數(shù)。二、選擇相應(yīng)的軟件現(xiàn)在市場上流行的有三大主流軟件：PCI、Erdas和Envi。這三款軟件都能滿足此項校正的需要，但是由于現(xiàn)在Envi是免費軟件，在其官網(wǎng)上能夠下載最新版本的軟件壓縮包，所以現(xiàn)在還是以Envi作為例子來說明。1.由于獲得的原始照片沒有坐標(biāo)系，現(xiàn)在我們先給其強行配備一個坐標(biāo)系統(tǒng)。當(dāng)然這個坐標(biāo)系統(tǒng)的參數(shù)只要設(shè)置的統(tǒng)一，怎樣設(shè)置可以根據(jù)自己的興趣。利用ARCGIS9.3打開arcmap下的工具箱ArcToolbox，找到里面的DataManagementTools項，然后子目錄下的ProjectionsandTransformations，接著選擇Raster下的DefineProjection,打開對話框選擇自己比較熟悉的投影坐標(biāo)系統(tǒng)（因為地理坐標(biāo)系統(tǒng)一般為經(jīng)緯度，而投影坐標(biāo)系統(tǒng)為米），對要進(jìn)行幾何校正的照片賦予其同樣的坐標(biāo)系統(tǒng)。2.打開Envi軟件，并打開要校正的照片（要在窗口中顯示出來，僅在Availablebandlists中看到不行）3.選擇Envi菜單欄中的Map下的Registration：——選擇selectGCPs：imagetomap，打開對話框?qū)υ捒蛑幸x擇自己最為熟悉的投影系統(tǒng)，進(jìn)行統(tǒng)一的參數(shù)選擇和設(shè)置，單位設(shè)置成米，而且分辨率最好設(shè)置的小一點（例如我現(xiàn)在設(shè)置其為1米）4.點擊ok就會出現(xiàn)GroundControlPointsSelection對話框來選擇自己的“地面控制點”。在圖像窗口點擊相應(yīng)的位置，然后輸入自己想要的坐標(biāo)（注意這時候我們最好根據(jù)分辨率把實際單位長度在輸入坐標(biāo)時放大成整數(shù)倍，例如實際照片上的1cm我可以輸入為200m，這時候?qū)?yīng)于1cm里就有200個小的像元，因為分辨率為1m，使得我們的照片在校正之后看上去大小適度，且細(xì)節(jié)比較清楚。還有輸入的坐標(biāo)最好不要跳躍的輸入否者照片變形較大，我建議最好每隔1cm輸入一個坐標(biāo)這樣幾何校正的精度會在一定程度上滿足你的需要）5.坐標(biāo)輸入完成后，一定要保存其為.pts格式的文件，以便連續(xù)的調(diào)用。步驟如下：在GroundControlPointsSelection對話框的菜單欄選擇File的下拉菜單的第一項“saveGCPSw/mapcoodinations”打開OutputRegistrationPoints對話框，選擇保存路徑，點擊OK即可。6.如果沒有關(guān)閉GroundControlPointsSelection對話框那么選擇菜單欄選擇File的下拉菜單的restoreGCPsfromASCII選項，打開EnterGroundConntrolPointsFilename對話框，選擇剛保存的.pts文件。然后選擇菜單欄的Options下的Warpdisplayedband選項就會打開RegistrationParameters對話框，現(xiàn)在只需要設(shè)置右邊的參數(shù)metod下有多項式，三角測量方法等，一般選擇Polynomial即多項式方法，后面的degree是多項式的最高此項，Resampling：最近鄰法、雙線性內(nèi)插和三次卷積（這個方法的優(yōu)劣這里不做解釋）；background一搬給出0值，即是背景色；然后選擇是保存文件還是在內(nèi)存中。最后點擊OK即可7.如果6中所說的GroundControlPointsSelection對話框已經(jīng)關(guān)閉，那么可以選擇Envi菜單欄中的Map下的Registration：——選擇warpfromGCPs：imagetomap，打開對話框，打開EnterGroundConntrolPointsFilename對話框，選擇剛保存的.pts文件，點擊Ok出現(xiàn)imagetomapRegistration對話框，設(shè)置的參數(shù)和第3步設(shè)置的參數(shù)要一致，點擊OK后出現(xiàn)Inputwarpimage對話框，在此窗口中選擇要進(jìn)行幾何校正的照片后點擊OK,就會打開RegistrationParameters對話框，現(xiàn)在只需要設(shè)置右邊的參數(shù)metod下有多項式，三角測量方法等，一般選擇Polynomial即多項式方法，后面的degree是多項式的最高此項，Resampling：最近鄰法、雙線性內(nèi)插和三次卷積（這個方法的優(yōu)劣這里不做解釋）；background一搬給出0值，即是背景色；然后選擇是保存文件還是在內(nèi)存中。最后點擊OK即可。8.通過上面的七步已經(jīng)完成了照片的幾何校正，那么如何求取照片中所關(guān)注部分的面積、周長、半徑或者其他屬性呢？利用上面的三個軟件的先關(guān)功能也可以求部分屬性，但是精度不如利用ARCGIS中的相關(guān)功能所求取的。9.打開ArcCatalog，在自己感興趣的文件中創(chuàng)建一個人數(shù)據(jù)庫“personGeodatabase.mdb”接著新建一個Featuredatabase，然后在其下新建一個FeatureClass，然后在特征類基礎(chǔ)上新建一topology。（說明如果要求面積周長等屬性，特征類應(yīng)是多邊形ploygon）10.在ArcGis中打開剛新建的特征類和已經(jīng)經(jīng)過幾何校正的照片。11.在編輯Editor工具條下選擇Startediting然后選擇后面的編輯工具，對照片感興趣的區(qū)域進(jìn)行編輯（注意這個編輯Target應(yīng)是新建的polygon上進(jìn)行）12.編輯完成后，選擇stopediting13.打開ploygon，右擊選擇openattributetable，打開屬性表，自動生成了編輯區(qū)域的面積和周長，如果要計算其他的屬性，可以點擊屬性表右下方的options的向下的三角號選擇AddField增加字段（設(shè)置相應(yīng)的字段名，數(shù)據(jù)類型，及其精度），然后選擇要進(jìn)行計算的字段列右擊選擇FieldCalculator，打開FieldCalcu