arcgis中对矢量和栅格数据进行裁剪切割的方法

arcgis中对矢量和栅格数据进行裁剪切割的方法[引用2011-11-0713:21:46]字号：大中小ArcGIS9:ArcToolbox--SpatailAnalystTools--ExtractionArcGIS8.X:在ArcMap--SpatialAnalyst扩展模块工具条--〉RasterCalculator或者直接使用Mask1、对矢量数据进行裁减：Arctoolbox中，spatialtoolextractclip在InputFeatures中选择被裁剪的图层，在ClipFeatures中选择裁剪形状2、利用矩形对栅格数据裁剪:Arctoolbox中，datamanagementrasterclipInputRaster中选择被裁剪的栅格数据，设定好矩形四个顶点即可3、在SpatialAnalystTools中提供了多种对栅格数据的提取方法，Arctoolbox中，SpatialAnalystToolsextraction包括提取值到点，根据属性提取，用圆提取，用多边形提取，用掩模提取，用点提取，用矩形提取等，其中的用掩模提取功能可以让我们通过不规则边界来获取需要的栅格数据。可以先建一个polygon的图层，再用这个图层clip那个raster加载空间分析模块，option选项中设置mask图层，范围，cell尺寸，用其栅格计算器计算：图层名.shape*1,点击evaluation，ok如果不要求一定要多变形裁剪的话，建议使用clipdatamanagementtoolbox---clip这个工具可以保持原有像元信息，但是这里arcmap只提供了矩形切割。这个约束应该不会影响一般的应用，否则的话建议使用图像处理软件中的mask裁剪。另外上面大家说得mask什么的方法，只能出来灰度图，可能原先是设计成处理grid的数据吧(2)用任意多边形剪切栅格数据(矢量数据转换为栅格数据)2.1在ArcCatlog下新建一个要素类（要素类型为:多边形），命名为：ClipPoly.shp2.2在ArcMap中，加载栅格数据：例如kunming.img、和ClipPoly.shp2.3打开编辑器工具栏，开始编辑ClipPoly，根据要剪切的区域，绘制一个任意形状的多边形。打开属性表，修改多边形的字段“ID”的值为1，保存修改，停止编辑。2.4打开空间分析工具栏执行命令：空间分析－转换--要素到栅格指定栅格大小：查询要剪切的栅格图层kunming的栅格大小，这里假设指定为1指定输出栅格的名称为路径2.5执行命令:空间分析-栅格计算器2.6构造表达式：[kunming]*[polyClip4-polyclip4]，执行栅格图层：kunming和用以剪切的栅格polyClip4之间的相乘运算两种方法差不多，第一种可能更好些，第二种是自己发现的虽然不怎么好但也分享下，希望有用呵呵(1)是否需要裁剪栅格图象区域通过一个面状的shapefile表达出来？如果可以，那么就很简单了。在ArcMap中，调用空间分析扩展模块，将你感兴趣区的shapefile多边形图层设置为掩膜，然后在栅格计算器中重新计算一下你的图象，它就会沿掩膜裁出。设置掩膜：空间分析工具条的下拉菜单option里面设置(2)用任意多边形剪切栅格数据(矢量数据转换为栅格数据)2.1在ArcCatlog下新建一个要素类（要素类型为:多边形），命名为：ClipPoly.shp2.2在ArcMap中，加载栅格数据：例如kunming.img、和ClipPoly.shp2.3打开编辑器工具栏，开始编辑ClipPoly，根据要剪切的区域，绘制一个任意形状的多边形。打开属性表，修改多边形的字段“ID”的值为1，保存修改，停止编辑。2.4打开空间分析工具栏执行命令：空间分析－转换--要素到栅格指定栅格大小：查询要剪切的栅格图层kunming的栅格大小，这里假设指定为1指定输出栅格的名称为路径2.5执行命令:空间分析-栅格计算器2.6构造表达式：[kunming]*[polyClip4-polyclip4]，执行栅格图层：kunming和用以剪切的栅格polyClip4之间的相乘运算文章二、使用ArcGIS的RasterCaclutaor剪切栅格文件如果我们需要把一个栅格文件（可以是tif、grid等arcgis所能识别的）的一部分特定感兴趣区截取出来，而且是任意形状的那种，就可以使用arcgis的RasterCaclutaor功能实现。利用RasterCaclutaor功能我们也可以把某个栅格图像与某一个矢量文件的相重叠部分单独切出来。两者的步骤基本一致，都是需要把矢量文件先转成栅格格式，然后利用进行象元间的加减乘除操作。步骤：1．将矢量文件和栅格文件同时添加到arcmap当中。2．先将矢量文件转为栅格文件”polygonRaster-polygonraster”Spatialanalysis工具条下的convert,FeaturestoRaster3.Spatialanalysis工具条打开RasterCaclutaor，在计算框内输入：[Topo.sid]-[polygonRaster-polygonraster]。Tips：这里采用的是两图层相剪操作，是因为polygonRaster–polygonrasterd栅格文件的Value为0，这样两图层象素间Value值对应相减，得到的是Topo.sid的Value值。如果说polygonRaster–polygonrasterd的value值为1.为了获得Topo.sid的value值则应该选择两图层象素值相乘操作。有时候这样的图层相剪操作会不成功，提示“Error:Failedtoevaluatethecalculatorexpression”。目前我还不太清楚具体原因是什么，在ESRI支持中心找到一篇文章，提到了可能的原因（语法错误、像元大小不一致、大小写不一致等）（=knowledgebase.techarticles.articleShow&d=18406）TheerrormessageFailedtoevaluatethecalculatorequationappearsafterclickingEvaluateontheRasterCalculatordialogbox.CauseThereareseveralcausesforthisproblem:•Youhaveasyntaxerrorinthecalculatorexpression.•Theextentandcellsizeofyourinputrasterisdifferentthantheextentsetforyouranalysisresults.•Theexpressionwillnotbeevaluatediftherearecarriagereturnsbeforeoraftertheexpression.•TheRasterCalculatoriscasesensitive.•Therearetoomanyopengrids.Thereisalimitof50openrastersofthegridformatatanyonetime.SolutionorWorkaroundSomeoftheseissueshavebeenaddressedinArcGISVersion8.1.2.Thecorrespondingworkaroundsfortheabovecausesare:Checkyoursyntax.TypeMapAlgebraintheindexoftheonlineHelpsystemforinformationonMapAlgebraanditssyntax.MakesuretheextentandcellsizeofyouranalysisresultissetappropriatelyintheExtentandcellsizetabsoftheOptionsdialogbox.Theextentofyouranalysisresultsshouldbesetthesameas,orencompass,theextentofthelayeryouareusingintheRasterCalculator.Donotusecarriagereturnsinyourexpression.Allowtheexpressiontoword-wrapifitislengthy.ThisproblemhasbeenresolvedinArcGIS8.1.2.Makesurethelayernametypedintotheexpressionboxcase-matchesthenameofthelayerintheLayerslist.Forexample:Layername:ElevationExpressiontyped:slope([elevation])Thisexpressionwillnotwork.Thefirstletterofelevationshouldbecapitalized.ThisproblemhasbeenresolvedinArcGIS8.1.2.Removeanyun-neededgridsfromtheArcMapTableofContents.Thetotalnumberofopengrids,includingtheonebeingcreatedbytheSpatialAnalystfunction,cannotexceed50.Thislimitdoesnotapplytonon-gridrasters.AworkaroundistoconvertgridstorastersinArcCatalogwiththeExportRasterToDifferentFormatoption