遙感地學(xué)分析-圖像解譯(PCI幾何糾正 ArcInfoWorkStation目視解譯 PCI監(jiān)督分類)【Beta10】.doc_第1頁
遙感地學(xué)分析-圖像解譯(PCI幾何糾正 ArcInfoWorkStation目視解譯 PCI監(jiān)督分類)【Beta10】.doc_第2頁
遙感地學(xué)分析-圖像解譯(PCI幾何糾正 ArcInfoWorkStation目視解譯 PCI監(jiān)督分類)【Beta10】.doc_第3頁
遙感地學(xué)分析-圖像解譯(PCI幾何糾正 ArcInfoWorkStation目視解譯 PCI監(jiān)督分類)【Beta10】.doc_第4頁
遙感地學(xué)分析-圖像解譯(PCI幾何糾正 ArcInfoWorkStation目視解譯 PCI監(jiān)督分類)【Beta10】.doc_第5頁
免費(fèi)預(yù)覽已結(jié)束,剩余36頁可下載查看

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡介

幾何校正圖像解譯圖像監(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)楸恍U膱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)榇中U?,圖像會按照對應(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è)文

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

最新文檔

評論

0/150

提交評論