arcgis中對(duì)矢量和柵格數(shù)據(jù)進(jìn)行裁剪切割的方法

1、arcgis中對(duì)矢量和柵格數(shù)據(jù)進(jìn)行裁剪切割的方法引用2011-11-0713:21:46 字號(hào)：大中小 ArcGIS9:ArcToolbox-SpatailAnalystTools-Extraction :在ArcMap-SpatialAnalyst擴(kuò)展模塊工具條-RasterCalculator或者直接使用Mask 1、對(duì)矢量數(shù)據(jù)進(jìn)行裁減：Arctoolbox中，spatialtoolextractclip在InputFeatures中選擇被裁剪的圖層，在ClipFeatures中選擇裁剪形狀2、利用矩形對(duì)柵格數(shù)據(jù)裁剪：Arctoolbox中，datamanagementrasterclip

2、InputRaster中選擇被裁剪的柵格數(shù)據(jù)，設(shè)定好矩形四個(gè)頂點(diǎn)即可 3、在SpatialAnalystT001s中提供了多種對(duì)柵格數(shù)據(jù)的提取方法，Arctoolbox 中，SpatialAnalyst Toolsextraction包括提取值到點(diǎn)，根據(jù)屬性提取，用圓提取，用多邊形提取， 用掩模提取，用點(diǎn)提取，用矩形提取等，其中的用掩模提取功能可以讓我們通過不規(guī)則邊界來獲取需要的柵格數(shù)據(jù)。 可以先建一個(gè)polygon的圖層，再用這個(gè)圖層clip那個(gè)raster 加載空間分析模塊， option選項(xiàng)中設(shè)置mask圖層，范圍，cell尺寸, 用其柵格計(jì)算器計(jì)算: 圖層名.shape*1, 點(diǎn)擊ev

3、aluation, ok 如果不要求一定要多變形裁剪的話，建議使用clip datamanagementtoolbox-clip 這個(gè)工具可以保持原有像元信息，但是這里arcmap只提供了矩形切割。 這個(gè)約束應(yīng)該不會(huì)影響一般的應(yīng)用，否則的話建議使用圖像處理軟件中的mask 裁剪。 另外上面大家說得mask什么的方法，只能出來灰度圖，可能原先是設(shè)計(jì)成處理 grid的數(shù)據(jù)吧 用任意多邊形剪切柵格數(shù)據(jù)（矢量數(shù)據(jù)轉(zhuǎn)換為柵格數(shù)據(jù)） 在ArcCatlog下新建一個(gè)要素類（要素類型為：多邊形），命名為:在ArcMap中，加載柵格數(shù)據(jù)：例如、和 打開編輯器工具欄，開始編輯ClipPoly,根據(jù)要剪切的區(qū)域，繪

4、制一個(gè)任意形狀的多邊形。打開屬性表，修改多邊形的字段“ID”的值為1,保存修改，停止編輯。 打開空間分析工具欄 執(zhí)行命令：- 指定柵格大小：查詢要剪切的柵格圖層kunming的柵格大小，這里假設(shè)指定為1指定輸出柵格的名稱為路徑 執(zhí)行命令：- 構(gòu)造表達(dá)式：kunming*polyClip4-polyclip4,執(zhí)行柵格圖層：kunming 和用以剪切的柵格polyClip4之間的相乘運(yùn)算 兩種方法差不多，第一種可能更好些，第二種是自己發(fā)現(xiàn)的雖然不怎么好但也 分享下，希望有用呵呵 是否需要裁剪柵格圖象區(qū)域通過一個(gè)面狀的shapefile表達(dá)出來？ 如果可以，那么就很簡(jiǎn)單了。 在ArcMap中，調(diào)用

5、空間分析擴(kuò)展模塊，將你感興趣區(qū)的shapefile多邊形圖層 設(shè)置為掩膜，然后在柵格計(jì)算器中重新計(jì)算一下你的圖象，它就會(huì)沿掩膜裁出。 設(shè)置掩膜：空間分析工具條的下拉菜單option里面設(shè)置 用任意多邊形剪切柵格數(shù)據(jù)（矢量數(shù)據(jù)轉(zhuǎn)換為柵格數(shù)據(jù)） 在ArcCatlog下新建一個(gè)要素類（要素類型為：多邊形），命名為:在ArcMap中，加載柵格數(shù)據(jù)：例如、和 打開編輯器工具欄，開始編輯ClipPoly,根據(jù)要剪切的區(qū)域，繪制一個(gè)任意形狀的多邊形。打開屬性表，修改多邊形的字段“ID”的值為1,保存修改，停止編輯。 打開空間分析工具欄 執(zhí)行命令：- 指定柵格大?。翰樵円羟械臇鸥駡D層kunming的柵格大小

6、，這里假設(shè)指定為1 指定輸出柵格的名稱為路徑 執(zhí)行命令：- 構(gòu)造表達(dá)式：kunming*polyClip4-polyclip4,執(zhí)行柵格圖層：kunming 和用以剪切的柵格polyClip4之間的相乘運(yùn)算 ，、一. 又早一、 使用ArcGIS的RasterCaclutaor剪切柵格文件 如果我們需要把一個(gè)柵格文件（可以是tif、grid等arcgis所能識(shí)別的）的一 部分特定感興趣區(qū)截取出來，而且是任意形狀的那種，就可以使用arcgis的 RasterCaclutaor功能實(shí)現(xiàn)。利用RasterCaclutaor功能我們也可以把某個(gè)柵 格圖像與某一個(gè)矢量文件的相重疊部分單獨(dú)切出來。兩者的步驟

7、基本一致，都是 需要把矢量文件先轉(zhuǎn)成柵格格式，然后利用進(jìn)行象元間的加減乘除操作。 步驟： 1. 將矢量文件和柵格文件同時(shí)添加到arcmap當(dāng)中。 2. 先將矢量文件轉(zhuǎn)為柵格文件polygonRaster-polygonraster” Spatialanalysis工具條下的convert,FeaturestoRaster 3. Spatialanalysis工具條打開RasterCaclutaor,在計(jì)算框內(nèi)輸入： -polygonRaster-polygonraster。 Tips:這里采用的是兩圖層相剪操作，是因?yàn)閜olygonRaster-polygonrasterd 柵格文件的Valu

8、e為0,這樣兩圖層象素間Value值對(duì)應(yīng)相減，得到的是的Value值。 如果說polygonRaster-polygonrasterd的value值為1.為了獲得的value值則應(yīng)該選擇兩圖層象素值相乘操作。 有時(shí)候這樣的圖層相剪操作會(huì)不成功， 提示“Error:Failedtoevaluatethecalculatorexpression0目前我還不太清楚具體原因是什么，在ESRI支持中心找到一篇文章，提到了可能的原因（語法錯(cuò)誤、像元大小不一致、大小寫不一致等） （&d=18409 TheerrormessageFailedtoevaluatethecalculatorequationapp

9、ears afterclickingEvaluateontheRasterCalculatordialogbox.Cause Thereareseveralcausesforthisproblem: ?Youhaveasyntaxerrorinthecalculatorexpression. ?Theextentandcellsizeofyourinputrasterisdifferentthantheextentsetforyouranalysisresults. ?Theexpressionwillnotbeevaluatediftherearecarriagereturnsbeforeo

10、raftertheexpression. ?TheRasterCalculatoriscasesensitive. ?Therearetoomanyopengrids.Thereisalimitof50openrastersofthegridformatatanyonetime. SolutionorWorkaround SomeoftheseissueshavebeenaddressedinArcGISVersion Thecorrespondingworkaroundsfortheabovecausesare: Checkyoursyntax.TypeMapAlgebraintheinde

11、xoftheonlineHelpsystemforinformationonMapAlgebraanditssyntax. MakesuretheextentandcellsizeofyouranalysisresultissetappropriatelyintheExtentandcellsizetabsoftheOptionsdialogbox.Theextentofyouranalysisresultsshouldbesetthesameas,orencompass,theextentofthelayeryouareusingintheRasterCalculator.Donotusec

12、arriagereturnsinyourexpression.Allowtheexpressiontoword-wrapifitislengthy. ThisproblemhasbeenresolvedinArcGIS Makesurethelayernametypedintotheexpressionboxcase-matchesthenameofthelayerintheLayerslist.Forexample: Layername:Elevation Expressiontyped:slope(elevation) Thisexpressionwillnotwork.Thefirstletterofelevationshouldbecapitalized. ThisproblemhasbeenresolvedinArcGIS Removeanyun-neededgridsfromtheArcMapTableofContents.Thetotalnumberofopengrids,inclu