ArcGIS中修改柵格圖的像元值

1、在處理圖象數(shù)據(jù)時(shí)，我們經(jīng)常會(huì)碰到要求修改柵格圖象象元值的問(wèn)題，比如說(shuō)DEM圖的部分?jǐn)?shù)據(jù)錯(cuò)誤，我們要進(jìn)行修改；再比如說(shuō)柵格圖象中有些與周圍均勻色彩不一致的錯(cuò)誤斑塊要更正等等，那我們?nèi)绾蝸?lái)處理這一類問(wèn)題呢？現(xiàn)我以一DEM柵格圖（名字為eldodem）為例，現(xiàn)在我要修改它的部分象元值，總結(jié)出以下三種方法，大家可以參考一下。第一種方法不大實(shí)用，但可借鑒，第二三種方法針對(duì)的條件不一致，大家可以在具體情況下進(jìn)行選擇。 一、直接運(yùn)用轉(zhuǎn)換，思路簡(jiǎn)單，易操作，但實(shí)用性1、柵格ascii文件柵格這種方法是先將柵格圖用工具 直接轉(zhuǎn)成ascii文件，然后在ascii文本文件中直接修改需要修改的象元的值，修改好后又用工

2、具 轉(zhuǎn)換成柵格圖。這種方法可行，但是不實(shí)用，因?yàn)槲覀円阉鞯街付ǖ南笤孟癫皇悄敲春?jiǎn)單。那么有沒有別的比較好的方法呢？二、在柵格計(jì)算器中操作，方法靈活，可操作性強(qiáng)，實(shí)用性強(qiáng)準(zhǔn)備工作先要在option中設(shè)置保留的柵格范圍，通常情況下默認(rèn)的為相交后的部分，這里我們要保留整個(gè)DEM，所以要改為以下設(shè)置： 2、通過(guò)點(diǎn)的位置修改點(diǎn)象元值21* 問(wèn)題一 、要是我想修改圖象中的指定的行列的點(diǎn)的象元值（比如說(shuō)把第100行，200列的點(diǎn)的值修改為0），那該如何操作呢？這時(shí)候我們可以在柵格計(jì)算器中輸入以下公式：con($rowmap = 99 & $colmap = 199),0,eldodem.img)

3、(柵格圖的編號(hào)是從0行0列開始的)22* 問(wèn)題二、那若是要修改指定行列范圍內(nèi)的柵格的象元值，比如說(shuō)把第101行，251列到第401行，301列的部分的象元值改為0，又改如何進(jìn)行呢？con($rowmap < 400 & $rowmap > 100 & $colmap < 300 & $colmap > 250),0,eldodem.img) 結(jié)果如下 23 * 問(wèn)題三、這些行和列是給定的，但是實(shí)際情況往往不是這樣的，那若看到一些要修改的點(diǎn)，我如何通過(guò)以上的方法來(lái)修改它的屬性值呢，這時(shí)要解決的首要問(wèn)題是如何獲知這個(gè)點(diǎn)的行列值，我們可以通過(guò)以下方法來(lái)

4、解決：（1）準(zhǔn)備工作，設(shè)置spatial analyst中的option，設(shè)置所需范圍，象素以及捕捉 這個(gè)確定后再次打開設(shè)定捕捉，使生成的柵格和原柵格象元重合（2）在柵格計(jì)算器中分別做出對(duì)應(yīng)的行列柵格圖行柵格圖的生成：在柵格計(jì)算器中輸入以下公式 r=$rowmapc=$colmap（在這里r-行，c-列）（3）然后可以用identify工具， 對(duì)要修改的點(diǎn)的行列號(hào)進(jìn)行查詢 如圖，所要修改的點(diǎn)位于第166行，129列，若要修改此點(diǎn)象元值為6，可以在柵格計(jì)算器中輸入以下公式：con($rowmap = 165 & $colmap = 128),6,eldodem.img)3、通過(guò)點(diǎn)面文件來(lái)

5、修改象元值*問(wèn)題、現(xiàn)在我們有一幅柵格圖，由于某種需要要修改其中的部分點(diǎn)單元格的象元值或者是修改部分斑塊的象元值，這時(shí)候該如何實(shí)現(xiàn)呢？本人覺得以下方法可行。思路主要是（1）先建一點(diǎn)shp或面shp圖層，（2）然后在點(diǎn)或面圖層上數(shù)字化需要修改的點(diǎn)或面，（3）接著將點(diǎn)或面轉(zhuǎn)換成柵格圖，（4）在柵格計(jì)算器中進(jìn)行替換運(yùn)算?，F(xiàn)我仍以eldodem為例，假設(shè)其中有部分象元的值是錯(cuò)誤的，要對(duì)其進(jìn)行修改，（同時(shí)進(jìn)行點(diǎn)面轉(zhuǎn)換操作）具體步驟如下：（1）新建一點(diǎn)文件point_11,然后在上面數(shù)字化出需要修改的點(diǎn)。新建一面文件eldodem，數(shù)字化出需要修改的面。 （2）將point_11和eldodem在featu

6、re to raster工具中分別轉(zhuǎn)換成柵格文件poin1和和eldo1.（3）在柵格計(jì)算器中進(jìn)行運(yùn)算假設(shè)現(xiàn)在要把個(gè)別點(diǎn)象元值修改為0，我們?cè)跂鸥裼?jì)算器中填入以下的公式：con(IsNull(poin1),eldodem.img,0) 然后點(diǎn)右下角的Evaluate運(yùn)行因?yàn)辄c(diǎn)象元值很小，放大可以看到我們修改后的微小的黑點(diǎn)對(duì)于面狀象元值的修改，相似的，輸入以下公式：con(IsNull(eldo1),eldodem.img,0) 然后運(yùn)行，得到的圖如下 （大家很可能得到的圖與此圖有異，要在屬性的symbol里設(shè)置一下，使其為拉伸樣式） 可以看到，修改過(guò)的部分變?yōu)楹谏?、修改一定經(jīng)緯度范圍內(nèi)的柵格象元值比如將經(jīng)度大于105度30分，緯度