遙感地學(xué)分析-圖像解譯（PCI幾何糾正 ArcInfoWorkStation目視解譯 PCI監(jiān)督分類）【Beta10】.doc

幾何校正圖像解譯圖像監(jiān)督分類附錄By fhys（不少錯(cuò)漏，多指正）大致流程：幾何校正幾何粗校正重投影裁剪部分區(qū)域幾何精校正圖像解譯統(tǒng)一投影信息創(chuàng)建主控層創(chuàng)建解譯圖層添加解譯屬性字段解譯各土地利用類型建立拓?fù)湫薷膽覓旎『蚻ab點(diǎn)AML批處理實(shí)用Arcmap中成圖ArcMap中coverage文件與shp文件的轉(zhuǎn)換圖像監(jiān)督分類附錄相關(guān)概念coverage格式文件概念如何在win7打開舊版help文件Arc命令&watchArcedit命令de【drawenvironment】Arc命令createArcedit命令sel【select】 Arcedit命令putArc命令buildArc命令ab【label】Arc命令additem屬性表的相關(guān)操作命令A(yù)rc命令CLEANCLEAN與build區(qū)別AML批處理實(shí)用nice的命令一、 幾何校正【Geometric correction】 本次所用軟件：PCI Geomatica V9.0.0。 操作過程中有不懂的，可以查看每個(gè)操作框的“help”獲取本機(jī)幫助，或者在PCI的安裝目錄的“manuals”文件夾找到對應(yīng)的PDF幫助文檔，如果我的PCI是裝在 “C:Geomatica_V90”下，那么對應(yīng)的PDF幫助文檔在“C:Geomatica_V90manuals”中，看幫助文檔很直接，比在網(wǎng)上搜索教程詳細(xì)很多。 這里把“配準(zhǔn)”和“校正”認(rèn)為是同一個(gè)概念，不做區(qū)分。 這里的“分辨率”指的是地物分辨率（如一個(gè)像元代表地物30*30m的分辨率就為30m），而不是指圖像的行列數(shù)的相乘 先列出原圖與粗校正后的圖像比較（左原，右后），便于后面清楚的步驟： #幾何粗校正利用PCI Geomatica中的GCPWorks模塊，進(jìn)行幾何粗校正。進(jìn)入GCPWorks模塊，在該菜單下，選擇適當(dāng)?shù)倪x項(xiàng)，然后Accept繼續(xù):幾何粗校正菜單應(yīng)選備注Processing RequirementsFull Processing可選的校正過程（分別是）：1、 包括校正所有步驟2、 只是選擇或查看GCP（地面控制點(diǎn)）3、 只是拼接圖像Mathematical ModelPolynomial配準(zhǔn)方式：1、 多項(xiàng)式的方式配準(zhǔn)2、 單純的圖像拼接？3、 薄板樣條的方式配準(zhǔn)Source of GCPsUser Entered Coordinates地面控制點(diǎn)（GCP）的來源：1、 用一幅已配準(zhǔn)的圖像來對比選點(diǎn)2、 使用數(shù)字化儀（有點(diǎn)古老）3、 用矢量圖中的點(diǎn)4、 用戶自定義輸入控制點(diǎn)（輸入點(diǎn)的行列位置及其對應(yīng)的投影坐標(biāo)位置）5、 ？在下一個(gè)對話框中，點(diǎn)擊Select Uncorrected Image，選擇待校正的圖像，即給定的gz.pix，選定需要顯示的三個(gè)波段（R5,G4,B3），合成后即為標(biāo)準(zhǔn)假彩色圖像。 接著，點(diǎn)擊Define Georeferencing Units，進(jìn)入對話框，在原本顯示“Meter”的下拉列表框中選擇“Other”，進(jìn)入對話框，選擇TM-Transverse Hercator橫軸墨卡托投影。在Earth Model中，選擇TM-Transverse Hercator橫軸墨卡托投影方式。 接著，設(shè)定投影參數(shù)，在Longitude中設(shè)定True origin的值為111d0000.0000“E（TM投影中，中央經(jīng)線的位置，因?yàn)閺V州所跨經(jīng)度約在112-113E之間，其最接近的中央經(jīng)線是111E）,在False中設(shè)定Easting的值為500000.000（在TM投影中，坐標(biāo)原點(diǎn)西移500km，為了使其x、y坐標(biāo)均為正），如圖：點(diǎn)擊Accept進(jìn)入下一個(gè)對話框，定義橢球類型，選擇E015-Krassovsky1940（E015的E代表Ellipsoids，橢球體的意思），如圖： 在投影信息顯示框點(diǎn)擊Accept，確定設(shè)置，回到GCPWorks模塊的對話框。接著，點(diǎn)擊Collect GCPs，進(jìn)行控制點(diǎn)選擇。由于是粗校正，所以只確定四個(gè)角點(diǎn)的大地坐標(biāo)以及其相對應(yīng)的行列號。方法可以是:l 直接在主窗口中輸入E、N和對應(yīng)的P、L然后點(diǎn)Accept as GCPl (用這個(gè)比較方便)也可以先寫一個(gè)txt文件，然后一次性導(dǎo)入（Fileload GCP file） 具體是先把TM圖像對應(yīng)的頭文件的四個(gè)邊角點(diǎn)的投影的X、Y值按照容易理解的方式放在Excel表格中然后復(fù)制到txt文本中，txt中的格式可以參見PCI圖框右下角的Help獲取最詳細(xì)的幫助文檔，這里用的是IPLXY格式（即序號 列數(shù) 行數(shù) 投影的X值 投影的Y值，中間用空格或tab符隔開）接著導(dǎo)入到PCI中選擇IPLXY格式，點(diǎn)Apply Format點(diǎn)Accept，返回到GCP主窗口再點(diǎn)Close，返回到GCPWorks主窗口最后是輸出保存粗校正的圖像Perform Registration to Disk，但要先進(jìn)行下面步驟：由于在PCI Geomatica V9.0.0中，需要新建一個(gè)含相同投影的（TM、E015、111E中央經(jīng)線、500km西移）、有足夠圖幅大小的圖像來保存這個(gè)粗校正后的圖像，才能保存粗校正后的圖像（BigBigBug）。所以打開ImageWorks，隨便選一個(gè)圖像或者User Defined，波段任意（主要是為了新建圖像，不是使用那個(gè)打開的圖），然后accept & load（如果是User Defined就點(diǎn)Accept）接著，打開菜單欄FileNew，進(jìn)入新建文件框，選擇pix格式，OK，進(jìn)入新建文件的詳細(xì)設(shè)置（這里有些需要特別注意，最好按后面的順序，便于按照這里進(jìn)行說明）：7個(gè)波段（因?yàn)楸恍Ｕ膱D有7個(gè)波段）選擇“Use bounds and resolution”（用邊界【左上角和右上角】和分辨率（30*30m代表一個(gè)像元）來決定行列數(shù)，所以不用設(shè)置行列數(shù)，它會隨著下面的設(shè)置而自定改變）選擇投影，跟前面一樣（TM、E015、111E中央經(jīng)線、500km西移）設(shè)置圖像的左上角和右下角（因?yàn)橹恍枰@兩個(gè)點(diǎn)就可以確定圖幅的大?。?，左上角對應(yīng)TM圖像對應(yīng)頭文件中Xmin、Ymax，右下角對應(yīng)頭文件的Xmax、Ymin到這步，Pixel Size會變成其他數(shù)字，把它設(shè)為30*30（這時(shí)會發(fā)現(xiàn)右下角的X、Y值會有輕微改變，這是因?yàn)橄裨钦麄€(gè)的，不能只有半個(gè)像元，為了剛好能把圖像“切成”整數(shù)個(gè)像元，軟件會自動設(shè)置UL和LR的與對應(yīng)E和N的差為分辨率的倍數(shù)(這里為30)，原則是至少保留邊角點(diǎn)在生成的圖像中）。到這里，可以看到圖像的行列數(shù)已經(jīng)變成比未校正的圖像的行列數(shù)大一些（這是因?yàn)榇中Ｕ?，圖像會按照對應(yīng)投影的X、Y值變得傾斜，導(dǎo)致圖幅要更大才【行列數(shù)更多】能容得下校正后的圖像，所以不用改成與原圖一樣的行列數(shù)）最后點(diǎn)Create生成新建文件新建文件完畢后，回到粗校正的輸入文件的窗口中，選“Default”（自動按順序?qū)?yīng)）重采樣（Resampling Mode）默認(rèn)最近鄰域插值（Nearest）Model默認(rèn)background默認(rèn)Memory（緩存），可以設(shè)成最大，可以少加快生成速度最后Perform Registration輸出最終結(jié)果到文件中打開看看是否正確、#重投影 之所以要重投影是因?yàn)?，精校正中的參考圖（俗稱“準(zhǔn)確的圖”）的投影不是TM E015 （橫軸墨卡托投影【也即高斯克里格投影】，克拉索夫斯基1940橢球體），而是ACEA E015（Albers投影【正軸等積割圓錐投影】），要統(tǒng)一地理基準(zhǔn)面和投影系統(tǒng)才能準(zhǔn)確進(jìn)行精校正。其實(shí)也可以把“準(zhǔn)確的圖”轉(zhuǎn)為TM E015，只要統(tǒng)一就行。 在我國，一般Albers投影的中心經(jīng)度是105E（我國區(qū)域的中心大致在115E），標(biāo)準(zhǔn)緯線是25N和47N，原點(diǎn)經(jīng)緯度（105E，0），偏移量為0【？】。中國地圖一般采用這個(gè)方式的投影，但在各省里，可以根據(jù)各自地理位置設(shè)置自定的參數(shù)（自定義中央經(jīng)線、標(biāo)準(zhǔn)緯線、坐標(biāo)原點(diǎn)），這里采用對廣州地區(qū)范圍的參數(shù)設(shè)置會更準(zhǔn)確，但統(tǒng)一起見，直接采用全國地圖的設(shè)置方式，也就是下面截圖的設(shè)置。 這時(shí)的Destination File保存的文件可以直接指定，而不像之前那樣要新建了才能用來保存。 最后在Reprojection窗口中設(shè)置投影和分辨率時(shí)，步驟跟前面基本一樣：選投影選擇Use bounds and resolution方式設(shè)分辨率為30*30，提示更改時(shí)與前面一樣，也是更改change lower right 不用修改左上角點(diǎn)和右下角點(diǎn)，會自動隨分辨率而變。全選到目標(biāo)圖層destination layersselectallreproject輸出圖像。 得到結(jié)果會發(fā)現(xiàn)丟失了部分圖像數(shù)據(jù)，如下面的圖，這大概是因?yàn)镻CI只把原圖像四個(gè)邊角點(diǎn)（不是衛(wèi)星圖的輪廓邊角點(diǎn)，而是原圖像顯示范圍的四個(gè)邊角點(diǎn)）轉(zhuǎn)換投影后的坐標(biāo)作為轉(zhuǎn)投影后的圖像的界限，而當(dāng)轉(zhuǎn)投影后圖像上下邊界的一些地區(qū)超出了該范圍后就變成丟失了。 但這里需要精校正的只是其中部分區(qū)域，所以丟失一點(diǎn)邊界的信息也不影響效果。若一定要保留所有圖像信息，可以把對應(yīng)的UL和LR設(shè)置到合適的值。 另外，要查看一幅圖像的坐標(biāo)到底有沒有像TM（橫軸墨卡托）投影那樣西移了500km的話，要在上面的Reprojection中選other，接著選一個(gè)投影，接著accept才會看到這幅圖默認(rèn)是否偏移了。 #裁剪部分區(qū)域 之所以要裁剪只是因?yàn)椤皽?zhǔn)確的圖”只有一小部分，所以為了保證圖像的精校正精度，最好是：“準(zhǔn)確的圖”的范圍=“不準(zhǔn)確的圖”的范圍。先看看效果（左為“準(zhǔn)確的圖”，右為對重投影后的大圖“裁剪后的圖”）： 步驟：依然在Utility的tools中選擇Subset，進(jìn)入Subset主窗口輸入輸出文件的位置和名字（可以直接指定文件名，而不用新建文件才能保存）點(diǎn)preview預(yù)覽，拉伸框框到合適區(qū)域（這個(gè)預(yù)覽圖看起來太小了）然后select all，并add到右邊去最后輸出圖像#幾何精校正 類似于粗校正，也是在GCPWorks中，進(jìn)入GCPWorks 主窗口，在Progressing Requirements中選擇Full Processing，在Mathematical Model選擇Thin Plate Spline（薄板樣條函數(shù)）【也可以選擇Polynomial多項(xiàng)式方式，各有優(yōu)缺點(diǎn)】，accept進(jìn)入下一步。選擇未精校正圖像，這里是gz_new1reprojSubset.pix，選5 4 3作為顯示波段，load and close回到校正流程的主窗口，接著選擇“準(zhǔn)確的圖”m440101.pix【Select Georeferenced Database】然后開始選GCP點(diǎn)（地面控制點(diǎn)），步驟與粗校正一樣【粗與精校正主要不同就是控制點(diǎn)數(shù)量不同】，保存結(jié)果也是要用新建的圖像來保存，也可以保存CGP點(diǎn)位TXT便于以后修改。具體精校正結(jié)果這里不再列出來。說明：在校正流程的主窗口中還可以發(fā)現(xiàn)Colour Matching等的選項(xiàng)，意思見下表：幾何精校正菜單備注Colour Matching自動或手動設(shè)置兩個(gè)圖像的顏色大致匹配Pre-registration Check預(yù)覽當(dāng)前所選GCP點(diǎn)校正的結(jié)果 最終得到文件-精校正后的圖像【這才是準(zhǔn)確的圖像】-二、目視解譯【visual interpretation】 本次所用軟件：ArcInfo的Workstation【若要安裝，安裝的路徑和目錄最好是沒有空格和中文，因?yàn)樗cCMD類似，會對空格識別為下一個(gè)命令或參數(shù)】【ArcMap這種可視化操作出來之前ArcGIS的主要軟件就是ArcInfo Workstation】，詳細(xì)概念可見后面附錄。 Arc命令窗口出了點(diǎn)擊【Arc】那個(gè)圖標(biāo)外，還可以在【運(yùn)行】【運(yùn)行的快捷鍵是Super鍵+R鍵】中輸入【Arc】就可快速打開Arc命令窗口。 workstation自帶的幫助文檔可以很好幫助學(xué)習(xí)，甚至是至關(guān)重要，因?yàn)樽钤敿?xì)的每個(gè)命令的使用說明和意義都在里面。,若是win7系統(tǒng)下,在arc下輸入help打不開hlp舊版幫助文檔,可以下載一個(gè)更新來解決問題,詳細(xì)見后面附錄。 這里說Arc命令是指：【Arc：】開頭的命令。同樣ArcEdit命令也即【ArcEdit：】開頭的命令，其他類似。 Arc命令【Arc：describe A】 可查看文件A的所有信息（包括投影信息） Arcedit命令【Arcedit：list 選中的要素】 可查看該對象的type類型等信息 【&Watch】這是一個(gè)比較重要的命令最好在操作之前知道，用與記錄命令窗口輸入和顯示的所有命令到一個(gè)文本文件中。詳細(xì)用法見后面附錄。一、 統(tǒng)一投影信息 無論做任何GIS分析或操作，涉及到兩個(gè)及以上圖層時(shí)都要統(tǒng)一這些圖層的投影、橢球體、坐標(biāo)偏移等信息，否則會出現(xiàn)操作的結(jié)果會是偏移和變形。這里coverage矢量文件與柵格圖的統(tǒng)一投影信息為:克拉索夫斯基橢球體（Krasovsky）Albers投影（等積正軸割圓錐投影）中央經(jīng)線：105E標(biāo)準(zhǔn)緯線：25N和47N（全國地圖常用）坐標(biāo)原點(diǎn)為：（105E，0N）？二、 創(chuàng)建主控層 ArcInfo Workstation中圖層的概念見后面附錄 命令行下面一行的紅色文字是對該行命令的解釋 命令后面的括號內(nèi)容如ae【Arcedit】，表示ae命令可用Arcedit代替，ae是縮寫了的命令 workstation中的命令不區(qū)分大小寫，就如CMD命令行一樣 可將在同一個(gè)coverage圖像中的arc、node、polygon等要素理解為這個(gè)coverage的圖層，可相互獨(dú)立編輯Arc: w D:InterpretationW【W(wǎng)orkspace】，設(shè)置工作目錄到“D:Interpretation”Arc: create gz_master gz_area 新建主控圖層【暫為空圖層】，所謂主圖層在這里只是一個(gè)大輪廓，采用gz_area的投影和圖像大小【boundary】，可用Arc命令describe gz_master來查看新建的主控圖層的信息Arc:ae進(jìn)入ae【Arcedit】工作模塊，命令行前面變成【Arcedit：】Arcedit:disp 9999disp【display】，設(shè)置顯示終端為電腦屏幕【9999代表電腦屏幕】，也可以設(shè)為0000【關(guān)閉顯示，圖像顯示框消失】Arcedit:ec gz_area ec【edit coverage】，設(shè)置要編輯的coverage文件名是gz_areaArcedit:de arcde【drawenvironment】，arc表示要畫的圖層是arc ，【de】的詳細(xì)用法見后面附錄。Arcedit:draw把當(dāng)前文件【gz_area】指定的要素【arc】顯示在屏幕上【9999】,draw之前的Arcedit命令都是為了設(shè)置顯示環(huán)境Arcedit：ef arcef【editfeature】，arc表示要編輯的圖層是arcArcedit：sel outlinesel【select】，outline表示要要選的是該段的最小輪廓，選中后的要素默認(rèn)會變成黃綠色見圖1【如果變黃色的部分不是完整的輪廓，可以輸入draw刷新一下】，有select當(dāng)然也有unselect，詳細(xì)見后面附錄Arcedit：put gz_masterput，存儲所選的要素及其attribute到gz_master中Arcedit：save保存上面put所做的修改圖1三、 創(chuàng)建解譯圖層Arc: copy gz_master gz_a從剛才剛得到arc類型輪廓的gz_master文件中復(fù)制為gz_a作為面圖層。但注意的是，這里發(fā)生的事情只是復(fù)制，也就是說gz_a中的內(nèi)容與gz_master一樣，只是后面建屬性表后再增加的字段是面類型【pat：polygon attribute】而已Arc: copy gz_master gz_l類似上面，從剛才剛得到arc類型輪廓的gz_master文件中復(fù)制為gz_l作為線圖層。四、 添加解譯屬性字段 創(chuàng)建線、面圖層Arc: ae進(jìn)入Arcedit模塊Arcedit: disp 9999Arcedit: ec gz_a編輯gz_a文件，面圖層，coverage格式Arcedit: de arc lab設(shè)置要顯示的要素是arc和lab【label：標(biāo)志點(diǎn)】 添加0 lab點(diǎn)【先要了解label點(diǎn)作用，見后面附錄】【主要是為了設(shè)置輪廓外的區(qū)域的值，這里的0 lab點(diǎn)的意思是該label點(diǎn)的值為0，且應(yīng)該點(diǎn)在大輪廓的外面，但由于還未建立build對應(yīng)的屬性表來存這個(gè)0值，所以這一步可以放在建立了屬性表后再點(diǎn)0 lab點(diǎn)】Arcedit: ef labEditfeature，要編輯的要素是labelArcedit: add添加label點(diǎn)，輸入了這個(gè)add命令后，把鼠標(biāo)移到顯示圖像的窗口，會變成十字絲，表示可以點(diǎn)label點(diǎn)，這里點(diǎn)在輪廓邊界的外面，只點(diǎn)一個(gè)點(diǎn)Arcedit: save保存這個(gè)label點(diǎn)Arcedit: qq【quit】，退出Arcedit模塊，回到Arc模塊 創(chuàng)建屬性表Arc: build gz_abuild，是給指定的coverage圖像【這里是gz_a】創(chuàng)建對應(yīng)要素類型的屬性表,其實(shí)可以用【build gz_a poly】來代替這個(gè)命令，因?yàn)閎uild默認(rèn)的最后一個(gè)參數(shù)是poly【polygon】，也就是建立一個(gè)對應(yīng)于多邊形【polygon】要素的屬性表Arc: build gz_l line類似上面，這里后面的參數(shù)是【line】，表示建立一個(gè)對應(yīng)于線狀要素的屬性表?！静煌豧eature會對應(yīng)不同類型的屬性表】 增加類型屬性字段Arc: additem gz_a.pat gz_a.pat LUtype 4 4 I為面圖層的屬性表添加屬性【列】，各參數(shù)含義見下面附錄【additem 輸入的作為基礎(chǔ)屬性表 輸出的屬性表可以不同名字，存為另外一個(gè)表 新建列的名字 新建列的顯示寬度 輸出新建列列的寬度 新建列的類型 新建列的其實(shí)位置】pat代表polygon attribute更詳細(xì)具體參數(shù)含義見后面附錄Arc: additem gz_l.aat gz_l.aat LUtype 4 4 I為線圖層屬性表添加屬性【列】，aat代表arc attribute 查看剛才的字段Arc: items gz_a.patitems 列出該屬性表的所有記錄lists and describes the items for all records of the specified INFO data file.Arc: items gz_l.aat同上，【如果不小心新建錯(cuò)了名字，要刪除可以用Arc：killinfo 屬性表命令來刪除】五、 解譯各土地利用類型 解譯到面圖層Arc: ae進(jìn)入Arcedit模塊Arcedit: disp 9999設(shè)置顯示Arcedit: ec gz_a設(shè)置要編輯的coverage文件是gz_aArcedit: de arc設(shè)置要顯示的是arc要素Arcedit: intersectarc all設(shè)置自動剪斷【相交則剪斷為分段】，這個(gè)很重要，如果忘記了的話，即使圈了再多圈圈也要重來，因?yàn)椴蛔詣哟驍鄷购竺鏌o法刪掉懸掛弧，如果不想打斷可以用【Arcedit：intersectarc off】Arcedit: image gz_Ok.tif composite 1 2 3Image，設(shè)置顯示gz_OK.tif圖像為背景圖層，其RGB取值分別與tif文件的1 2 3 波段對應(yīng)?！綼dds an image or a group of images from an image catalog to the background draw environment.】【這里要注意的是，選作背景的圖的地理投影和橢球體一定要與gz_a這個(gè)coverage的投影和橢球體相同，否則解譯出來的圖像不符合實(shí)際，這里統(tǒng)一的投影見上面的第一小節(jié)統(tǒng)一投影信息】【由于gz_OK.tif這個(gè)文件的投影信息放在了與其同目錄下的gz_OK.aux文件中，這是PCI軟件對tif格式文件產(chǎn)生的補(bǔ)充文件，而workstation無法識別這個(gè)文件so does Arcmap，只能識別tif文件中所帶的坐標(biāo)信息和邊界信息，詳細(xì)可用Arc：describe gz_OK.tif查看，見下面附圖。所以即使workstation無法識別tif的投影信息，但只要先在Arcedit中打開了gz_a這個(gè)coverage文件ec gz_a，有了gz_a的投影坐標(biāo)作為基礎(chǔ)，也可以準(zhǔn)確加載gz_OK.tif到顯示的背景中】Arcedit: draw顯示coverage文件和gz_OK.tif背景Arcedit: ef arc設(shè)置要編輯的是coverage文件中的arc要素圖層Arcedit: add開始圈圈【圈出各種土地利用類型】，這時(shí)會有提示，見下面的附圖，其中按數(shù)字鍵2【node】是代表開始和結(jié)束，1【vertex】是添加轉(zhuǎn)折點(diǎn)，3【curve】使直線變成曲線圓滑,4【delete】刪除一個(gè)直線的節(jié)點(diǎn)，5【delete arc】刪除連續(xù)的整條曲線不包括直線部分，6【spline on/off】？點(diǎn)了奇數(shù)個(gè)點(diǎn)時(shí)會自動變成曲線，7【square on/off】使后面畫的直線都與前面一條直線的角度成90不能與spline 一起用，8是其他一些參數(shù)的設(shè)置，可暫不用細(xì)究，9【quit】退出。每次圈到的多邊形都應(yīng)該首尾有交叉【這樣才能閉合多邊形】，但要注意，不能交叉的出來的線【稱懸掛弧】不能太長，會影響后面對這些“尾巴”的剪掉【因?yàn)閏lean命令可自動按照指定最長伸出范圍來裁剪，如果尾巴太長，要剪掉就需要設(shè)比較大的閾值，可能導(dǎo)致有些不是“尾巴”的線都被自動裁減掉】這里圈圈時(shí)候要注意，工作量會很大，所以很有必要知道上面這些參數(shù)的意義【起碼知道圈錯(cuò)了怎樣才能后退一步】。還有，圈了不少地物后最好用【save】命令保存一下，否則萬一死機(jī)或操作出錯(cuò)也不至于要重頭再來。這樣圈出來的都是arc類型要素【因?yàn)锳rcedit: ef arc】，所以還要進(jìn)行clean【建立拓?fù)涞拿?，后面用到】?？梢韵劝阉钥梢苑蛛x出來的地物都圈出來，圈完后再統(tǒng)一按類別點(diǎn)上label點(diǎn)【用于標(biāo)識多邊形屬性的點(diǎn)】，再統(tǒng)一賦類型值【calculate命令】Arcedit: ef lab編輯label點(diǎn)，即給每個(gè)圈出來的圈圈點(diǎn)上唯一的一個(gè)label點(diǎn)，最好是統(tǒng)一一次性把所有的相同地物類型都點(diǎn)上，然后一次性用【calculate】命令賦上屬性值A(chǔ)rcedit: add開始添加label點(diǎn)Arcedit: calculate LUtype = 2給label點(diǎn)賦值，【注意】【type = 2】 之間要有空格隔開，否則報(bào)錯(cuò)Arcedit:list選中對應(yīng)要素后查看其屬性相關(guān)信息Arcedit: save保存所有的修改：圈圈和label點(diǎn)，其實(shí)可以在上一步圈圈時(shí)就save一次 解譯到線圖層Arc: ae進(jìn)入Arcedit模塊，后面基本都與【解譯到面圖層】類似，只是畫完線圖層后，不用clean來建拓?fù)洹疽驗(yàn)閏lean會把交叉弧段打斷，而實(shí)際的線圖層不需打斷】Arcedit: disp 9999Arcedit: ec gz_l設(shè)置要編輯的coverage文件是gz_lArcedit: de arcArcedit: intersectarc all設(shè)置自動剪斷【相交則剪斷為分段】，這個(gè)很重要，如果忘記了的話，即使圈了再多圈圈也要重來，因?yàn)椴蛔詣哟驍鄷购竺鏌o法刪掉懸掛弧，如果不想打斷可以用【Arcedit：intersectarc off】Arcedit: image gz_Ok.tif composite 1 2 3- 至于線圖層的解譯，類似。一般小道路和小河流等線狀物應(yīng)該解譯到在線圖層中，而像珠江這些大河最好用解譯到面圖層。六、 建立拓?fù)?建拓?fù)浯致哉f就是把點(diǎn)線面這三大類要素分清楚，并且符合實(shí)際【舉例來說，不能有重疊的兩個(gè)多邊形，易導(dǎo)致統(tǒng)計(jì)分析時(shí)候出錯(cuò)，但如果只是一般成圖，可適當(dāng)省去】。在Arcmap和Arccatalog中操作時(shí)，其實(shí)剛畫出來的多邊形等東西都也是未建立拓?fù)潢P(guān)系的，在ArcMap和ArcCatalog中具體建拓?fù)湟奱rccatalog拓?fù)浜痛司W(wǎng)址。 而workstation中如果是【ef arc】是先不分這些東西，先畫出弧段，再間接用這些閉合的弧段和【build】或【clean】命令來生成polygon圖層。 所以按照這個(gè)原理，可以先build或clean一個(gè)polygon圖層【coverage里面的圖層，而非整個(gè)coverage】，然后再在這個(gè)圖層中用【ef poly】來添加多邊形，效果應(yīng)該一樣。類似的，【ef line】 不過下面還是按照先【ef arc】，把東西都畫到arc圖層中，再build或clean建立polygon。Arc: copy gz_a gz_a1備份一下剛才圈到的圖，便于以后單獨(dú)修改，不過build【建立拓?fù)洹渴墙⒌皆揷overage圖像的poly【polygon】或line圖層中，不影響剛?cè)Φ降腶rc圖層,不過后面會以gz_a1為基礎(chǔ)繼續(xù)操作Arc: build gz_a1建立拓?fù)洹灸J(rèn)為polygon多邊形】，之前剛裁剪到輪廓時(shí)已經(jīng)build【建立拓?fù)洹恳淮?，這時(shí)還要build是因?yàn)樾氯Φ降腶rc還未建立對應(yīng)的拓?fù)涞絧oly【polygon】圖層中Arcedit: clean gz_a1 gz_a1_clean1 50 20Clean【建立拓?fù)浜蛯?yīng)的屬性表到新的文件或輸入文件中】，第一個(gè)參數(shù)是輸入文件，第二個(gè)參數(shù)是輸出文件，第三個(gè)是剪斷懸掛弧的長度【大于它則不剪，小于才剪，包括閉合的弧段也在計(jì)算范圍內(nèi)】，第四個(gè)參數(shù)是各個(gè)體間的頂點(diǎn)的最小距離【小于它則合并，大于則不管】。單位都是對應(yīng)投影坐標(biāo)的單位，這里可以從【Arc：describe gz_a1】命令中看到其坐標(biāo)位置的單位是米【m】，而背景圖gz_OK.tif的分辨率只有30米，所以后面的兩個(gè)參數(shù)可能要設(shè)得大一點(diǎn)才能把所有懸掛弧剪掉。其實(shí)還有第五個(gè)參數(shù)，就是【polypolygon】或【line】，省略則為【polypolygon】更詳細(xì)的clean的各參數(shù)含義見后面附錄-下面的查錯(cuò)操作不需按順序，按需而用：-Arc: labelerrors gz_a1檢查每個(gè)多邊形【poly】是否有且只有一個(gè)lab點(diǎn)【因?yàn)閘abel點(diǎn)是標(biāo)志多邊形的屬性的，所以少于或大于1個(gè)lab點(diǎn)都是有錯(cuò)的】Arc: nodeerrors gz_a1 dangles檢查是否還有懸掛弧Arcedit: de node dangle設(shè)置要顯示的要素是懸掛弧Arcedit: de node 設(shè)置要顯示的要素是節(jié)點(diǎn)，要注意【Arcedit: de node dangle】Arcedit: nodecolor dangle 2設(shè)置dangle懸掛弧的顏色為2關(guān)于顏色的數(shù)字意思These are the default screen colors accessed by numbers 0 to 7 on most graphic terminals and plotters:Arcedit: sds 7sds【setdrawsymbol】，設(shè)置選中的要素的顏色Arcedit: dsds【drawselect】，把選中的要素的顏色畫出來-Arcedit: ef poly要編輯的要素是polygonArcedit: sel gz_a1# = 5選擇【gz_a1】文件中編號為5的多邊形，因?yàn)榍耙粋€(gè)命令是編輯polygon圖層 【Arcedit: ef poly】?！?】代表編號的標(biāo)識符。這樣就可以對選中的多邊形進(jìn)行刪除修改的等操作Arcedit: ef lab類似的，要編輯的要素是lab點(diǎn)- 至于線圖層的建拓?fù)?，類似，這里名為gz_l1。但注意的是，建線的拓?fù)鋾r(shí)，【build】不會打斷較差弧段，而【clean】會打斷。七、 修改懸掛弧和lab點(diǎn)重復(fù)上述步驟，不斷檢查懸掛弧錯(cuò)誤和lab點(diǎn)錯(cuò)誤，直到錯(cuò)誤均為0才算完成，最后輸出無錯(cuò)的面圖層gz_a1和線圖層gz_l1。八、 AML批處理實(shí)用 AML【Arc Macro Language】，這里AML是指ArcGIS批【/宏】處理語言。 一句話來說【/大致來說】，AML是一種可以將多個(gè)Workstation命令行中的命令保存下來的可以一次性執(zhí)行的文本文件，一般為【.aml】后綴，只要是文本格式就行。正是這樣，AML文件可以提高工作效率。更詳細(xì)關(guān)于AML的介紹見后面附錄。開始：Arc：&amlpath D:InterpretationamlFiles【&amlpath pathnames】設(shè)置aml文件的默認(rèn)工作目錄為【D:InterpretationamlFiles】，【&amlpath pathnames】命令可一次性設(shè)置多個(gè)目錄。如果要打開其他路徑的aml文件，則出了文件名之外，還需加上其所在路徑。Arc：&run openTif【&run】，也可縮寫為【&r】，是打開剛才設(shè)置的那些目錄些的以【.aml】結(jié)尾的aml文件，如果不是【.aml】結(jié)尾的文件，則不能省略后綴名。以下是【openTif.aml】文件的內(nèi)容【其中/*表示注釋，不執(zhí)行】：w D:Interpretationlcaedisp 9999 ec gz_ade arc labdrawimage gz_OK.tif composite 1 2 3 draw /*intersectarc all/*ef arc/*add 一個(gè)打開aml文件的命令就把一幅圖顯示出來了。只是這個(gè)設(shè)置和打開aml文件的使用就可以，用兩三條命令代替輸一堆命令了。除了打開圖像，其他一些設(shè)置如添加watch記錄文件、初始化環(huán)境等都可以用aml文件來提高操作效率。九、 ArcMap中成圖 打開Arcmap，加載上面得到的面圖層gz_a1和線圖層gz_l，最好是先轉(zhuǎn)成shp文件【方法可以是Export Data或者見下圖的to coverage/to shapefile】 給不同類型（TYPE）的多邊形和線賦上不同的顏色。 點(diǎn)擊layout，接著調(diào)整大小、加入題目、圖例、以tif或其他格式輸出圖像。最終得到文件，最好用【Arc：copy】命令，把兩個(gè)coverage文件復(fù)制到另外一個(gè)文件夾【直接右鍵復(fù)制會導(dǎo)致兩種情況中的一種：1雖然只有一個(gè)coverage文件，但公用的info文件夾會包含冗余信息；2沒把inbfo文件也復(fù)制過去，導(dǎo)致信息不全。所以最好還是copy命令】：-gz_a1【coverage格式的文件】gz_l1【coverage格式的文件】tif圖像【Arcmap成的圖】設(shè)置了相對路徑的Arcmap文件mxd，就是輸出tif格式圖像時(shí)的狀態(tài)【mxd文件一般不大】-十、 ArcMap中coverage文件與shp文件的轉(zhuǎn)換在ArcMap中不能直接編輯Coverage格式的文件，需要先轉(zhuǎn)換到shp格式才行，而其中的可選擇對應(yīng)關(guān)系：Coverage的label或pointshp中的pointCoverage的polygon或regionshp中的polygonCoverage的arcshp中的polyline同樣的，從shp格式要轉(zhuǎn)回coverage格式，可同時(shí)選中point圖層和polygon圖層，分別作為label和polygon【默認(rèn)的是region】輸出，輸出時(shí)會保留屬性表的所有信息。但有個(gè)缺點(diǎn)是，coverage格式轉(zhuǎn)到shp格式后，要拉動其中的相接的多邊形時(shí)，會出現(xiàn)空洞【即，本來是共用一條邊的，拉動了左邊，右邊卻依然還是原來位置】通過對比Workstation中的建立拓?fù)渑cArcMap、ArcCatalog建立拓?fù)涞牟煌?，可更容易理解建立拓?fù)涞母拍?。并可知道，對于建立拓?fù)鋪碚f，Workstation暫時(shí)還是比ArcMap9.3要好一點(diǎn)的，【W(wǎng)orkstation還有Arctools這個(gè)界面化的工具Arc：arctools】三、監(jiān)督分類【Supervised Classification】 監(jiān)督分類：就是認(rèn)為預(yù)先選好一些認(rèn)為是某種地物的區(qū)域【即“訓(xùn)練區(qū)”】，然后軟件按照某種分類方法進(jìn)行地物分類【最大似然法、波譜角法、最近距離法etc.】得到窗口如下： 上面，第6、7波段本來有信息的，但被用來作為訓(xùn)練區(qū)保存波段和分類結(jié)果保存波段了，如果想使最后兩個(gè)波段也參與分類分析，可以在開始分類前給圖形增加兩個(gè)波段【是真實(shí)的波段，不是顯示波段】。 還有，如果到這一步了，不繼續(xù)下去，關(guān)閉文件，再打開時(shí)會發(fā)現(xiàn)6、7波段已經(jīng)更改了，不再是原圖的6、7波段了，而且所選哪個(gè)波段作為顯示波段、哪個(gè)波段作為訓(xùn)練區(qū)波段和結(jié)果波段都被保留著，這說明到這步PCI會默認(rèn)把更改寫入到實(shí)際的圖像文件中了。所以注意備份好分類前的圖像文件。 下面，單擊【Trace & Close】后，要等一下，直到按鈕“陷下去”才能開始在主顯示窗口“圈圈”，與【Trace】不同的是【Trace】不會在鼠標(biāo)松開后自動閉合首尾點(diǎn)，一般用【Trace & Close】。 下面，【Undo】：撤銷上一步操作，但只能先后撤銷一步。下面，橡皮擦工具很好用，可以擦掉圈錯(cuò)的區(qū)域。按照各種類地物圈下去。序號主要地物類型圖像中 大概顏色分類后顯示顏色1闊葉林broadleaved forest鮮紅色，山中，成片分布鮮綠2針葉林coniferous forest暗紅色，山中，成片分布暗綠3草地grass淡紅色，可參照珠海校區(qū)后山的高爾夫球場的草地的顏色黃綠4灌木bush淡紅顏色，色雜，常與草地混在一起，山谷山間雜青5水田paddy珠三角較多，比較規(guī)則，【但基塘也很規(guī)則】墨綠6旱地dry land沒裸土那么白，色稍雜，暗黃7裸土exposed soil偏白，裸露山坡山腳，河灘黃色8居民地settlement淺綠色帶白，廣州城區(qū)很多淺白色9清潔水clean water深藍(lán)色，如圖像右上角藍(lán)色10渾濁水turbid water淺藍(lán)色，反射率較高藍(lán)白色11火燒跡地brule山上，灰黑色，清遠(yuǎn)北江河邊的山有些，但已與陰影搞混灰色Google Earth is nice now.僅供參考，且是R-4 G-3 B-2波段顯示下的顏色，不是的話無需參考圈完后，看下圖，Signaturestatistics特征統(tǒng)計(jì)?？梢钥疵恳活惖木?方差 和協(xié)方差矩陣。下圖，Signatureseparability 可分性看每種樣本與其他各樣本的可分性 (可分百分比)越大越能區(qū)分。下圖，Classify選擇算法 平行算法-最小距離法-最大似然法點(diǎn)擊classify分類。下圖,如果要重新打開圖像，查看分類后結(jié)果圖像，先要把波段選到保存的結(jié)果波段，這里是【7】，然后再在【Classify】中打開【Class Labeling】，可以看到分類后的結(jié)果圖。分類后大致結(jié)果【有較大錯(cuò)誤，需繼續(xù)修改】： 一次選擇訓(xùn)練區(qū)可能得到結(jié)果不太好【上圖的珠江口有些變成了植被，有問題，而背景色由于也被參與到分類了暫未找到去掉背景顏色的方法，所以顏色變了，但不影響圖像效果，但影響分類結(jié)果的統(tǒng)計(jì)】，需要選中有問題的訓(xùn)練區(qū)類型【只能選中整個(gè)類型，而貌似不能選中這類訓(xùn)練樣區(qū)的某一塊】，然后點(diǎn)【Clear Class（es）】按鈕，清除了這個(gè)分類，再找這個(gè)分類的訓(xùn)練區(qū)，再執(zhí)行分類算法。如此這般，直到符合要求為止。 通常為了避免有些地物被明顯分到其他分類中去，最好在整幅圖像有該地物的區(qū)域都“圈個(gè)圈”。 這里的“火燒跡地”接近于灰黑色，所以不少陰影也容易被分到火燒跡地中去了。第n次后：最大似然法：最小距離法：平行六面體法：無論用哪種方法，關(guān)鍵都是選好“訓(xùn)練樣區(qū)”和知道如何修改上次的訓(xùn)練樣區(qū)。分類完畢后最好截圖下來。把1-6波段這些圖層刪掉【刪掉前最好備份一個(gè)刪之前的pix文件，便于以后修改】。最終得到文件：-認(rèn)為分類最好的分類后的pix文件【只有一個(gè)圖層的，保存著分類結(jié)果】-#附錄【主要是Arcinfo Workstation的】附錄：ArcInfo Workstation相關(guān)概念復(fù)習(xí)概念理解：【摘自ArcGIS9.3產(chǎn)品白皮書】ArcGIS產(chǎn)品的基本架構(gòu)桌面GIS包括ArcReader、ArcView、ArcEditor、ArcInfo和ArcGIS擴(kuò)展模塊，功能方面ArcInfo ArcEditor ArcView ArcReader其中我們常用的ArcMap屬于ArcInfo的子軟件，而ArcInfo的Workstation相當(dāng)于命令行式的ArcMap，Arc/Info的原著格式包括coverage【這種矢量文件與GRID文件一樣是幾個(gè)文