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

上传人:b****2 文档编号:24168907 上传时间:2023-05-25 格式:DOCX 页数:6 大小:19.74KB
下载 相关 举报
arcgis中对矢量和栅格数据进行裁剪切割的方法.docx_第1页
第1页 / 共6页
arcgis中对矢量和栅格数据进行裁剪切割的方法.docx_第2页
第2页 / 共6页
arcgis中对矢量和栅格数据进行裁剪切割的方法.docx_第3页
第3页 / 共6页
arcgis中对矢量和栅格数据进行裁剪切割的方法.docx_第4页
第4页 / 共6页
arcgis中对矢量和栅格数据进行裁剪切割的方法.docx_第5页
第5页 / 共6页
点击查看更多>>
下载资源
资源描述

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

《arcgis中对矢量和栅格数据进行裁剪切割的方法.docx》由会员分享,可在线阅读,更多相关《arcgis中对矢量和栅格数据进行裁剪切割的方法.docx(6页珍藏版)》请在冰豆网上搜索。

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

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

arcgis中对矢量和栅格数据进行裁剪切割的方法[引用2011-11-0713:

21:

46]

 

 

字号:

大中小

 

ArcGIS9:

  ArcToolbox-->SpatailAnalystTools-->Extraction

:

在ArcMap-->SpatialAnalyst扩展模块工具条--〉Raster Calculator或者直接使用Mask

 

1、对矢量数据进行裁减:

Arctoolbox中,spatialtool>extract>clip

在InputFeatures中选择被裁剪的图层,在ClipFeatures中选择裁剪形状

2、利用矩形对栅格数据裁剪:

Arctoolbox中,datamanagement>raster>clip

InputRaster中选择被裁剪的栅格数据,设定好矩形四个顶点即可

3、在SpatialAnalystTools中提供了多种对栅格数据的提取方法,Arctoolbox中,SpatialAnalyst

Tools>extraction包括提取值到点,根据属性提取,用圆提取,用多边形提取,用掩模提取,用点提取,用矩形提取等,其中的用掩模提取功能可以让我们通过不规则边界来获取需要的栅格数据。

 

可以先建一个polygon的图层,再用这个图层clip那个raster

 

加载空间分析模块,

option选项中设置mask图层,范围,cell尺寸,

用其栅格计算器计算:

 

图层名.shape*1,

点击evaluation,

ok

 

如果不要求一定要多变形裁剪的话,建议使用clip

datamanagementtoolbox---clip

这个工具可以保持原有像元信息,但是这里arcmap只提供了矩形切割。

这个约束应该不会影响一般的应用,否则的话建议使用图像处理软件中的mask裁剪。

另外上面大家说得mask什么的方法,只能出来灰度图,可能原先是设计成处理grid的数据吧

 

 

 

 

 

(2)

用任意多边形剪切栅格数据(矢量数据转换为栅格数据)

在ArcCatlog下新建一个要素类(要素类型为:

多边形),命名为:

在ArcMap中,加载栅格数据:

例如、和

打开编辑器工具栏,开始编辑ClipPoly,根据要剪切的区域,绘制一个任意形状的多边形。

打开属性表,修改多边形的字段“ID”的值为1,保存修改,停止编辑。

打开空间分析工具栏  

执行命令:

<空间分析>-<转换>--<要素到栅格>

指定栅格大小:

查询要剪切的栅格图层kunming的栅格大小,这里假设指定为1

指定输出栅格的名称为路径

执行命令:

<空间分析>-<栅格计算器>

构造表达式:

[kunming]*[polyClip4-polyclip4],执行栅格图层:

kunming和用以剪切的栅格polyClip4之间的相乘运算

两种方法差不多,第一种可能更好些,第二种是自己发现的虽然不怎么好但也分享下,希望有用 呵呵

 

(1)

是否需要裁剪栅格图象区域通过一个面状的shapefile表达出来

如果可以,那么就很简单了。

在ArcMap中,调用空间分析扩展模块,将你感兴趣区的shapefile多边形图层设置为掩膜,然后在栅格计算器中重新计算一下你的图象,它就会沿掩膜裁出。

 

设置掩膜:

空间分析工具条的下拉菜单>option里面设置

 

(2)

用任意多边形剪切栅格数据(矢量数据转换为栅格数据)

在ArcCatlog下新建一个要素类(要素类型为:

多边形),命名为:

在ArcMap中,加载栅格数据:

例如、和

打开编辑器工具栏,开始编辑ClipPoly,根据要剪切的区域,绘制一个任意形状的多边形。

打开属性表,修改多边形的字段“ID”的值为1,保存修改,停止编辑。

打开空间分析工具栏  

执行命令:

<空间分析>-<转换>--<要素到栅格>

指定栅格大小:

查询要剪切的栅格图层kunming的栅格大小,这里假设指定为1

指定输出栅格的名称为路径

执行命令:

<空间分析>-<栅格计算器>

构造表达式:

[kunming]*[polyClip4-polyclip4],执行栅格图层:

kunming和用以剪切的栅格polyClip4之间的相乘运算

 

文章二、

 使用ArcGIS的RasterCaclutaor剪切栅格文件

如果我们需要把一个栅格文件(可以是tif、grid等arcgis所能识别的)的一部分特定感兴趣区截取出来,而且是任意形状的那种,就可以使用arcgis的RasterCaclutaor功能实现。

利用RasterCaclutaor功能我们也可以把某个栅格图像与某一个矢量文件的相重叠部分单独切出来。

两者的步骤基本一致,都是需要把矢量文件先转成栅格格式,然后利用进行象元间的加减乘除操作。

步骤:

1. 将矢量文件和栅格文件同时添加到arcmap当中。

2. 先将矢量文件转为栅格文件”polygonRaster-polygonraster”

Spatialanalysis工具条下的convert,FeaturestoRaster

3.Spatialanalysis工具条打开RasterCaclutaor,在计算框内输入:

[]-[polygonRaster-polygonraster]。

 

Tips:

这里采用的是两图层相剪操作,是因为polygonRaster–polygonrasterd栅格文件的Value为0,这样两图层象素间Value值对应相减,得到的是的Value值。

如果说polygonRaster–polygonrasterd的value值为1.为了获得的value值则应该选择两图层象素值相乘操作。

有时候这样的图层相剪操作会不成功,提示“Error:

 Failedtoevaluatethecalculatorexpression”。

目前我还不太清楚具体原因是什么,在ESRI支持中心找到一篇文章,提到了可能的原因(语法错误、像元大小不一致、大小写不一致等)

(&d=18406)

 

Theerrormessage"Failedtoevaluatethecalculatorequation"appearsafterclickingEvaluateontheRasterCalculatordialogbox.

Cause

Thereareseveralcausesforthisproblem:

•Youhaveasyntaxerrorinthecalculatorexpression.

•Theextentandcellsizeofyourinputrasterisdifferentthantheextentsetforyouranalysisresults.

•Theexpressionwillnotbeevaluatediftherearecarriagereturnsbeforeoraftertheexpression.

•TheRasterCalculatoriscasesensitive.

•Therearetoomanyopengrids.Thereisalimitof50openrastersofthegridformatatanyonetime.

SolutionorWorkaround

SomeoftheseissueshavebeenaddressedinArcGISVersion

Thecorrespondingworkaroundsfortheabovecausesare:

Checkyoursyntax.TypeMapAlgebraintheindexoftheonlineHelpsystemforinformationonMapAlgebraanditssyntax.

MakesuretheextentandcellsizeofyouranalysisresultissetappropriatelyintheExtentandcellsizetabsoftheOptionsdialogbox.Theextentofyouranalysisresultsshouldbesetthesameas,orencompass,theextentofthelayeryouareusingintheRasterCalculator.

Donotusecarriagereturnsinyourexpression.Allowtheexpressiontoword-wrapifitislengthy.

ThisproblemhasbeenresolvedinArcGIS

Makesurethelayernametypedintotheexpressionboxcase-matchesthenameofthelayerintheLayerslist.Forexample:

Layername:

Elevation

Expressiontyped:

slope([elevation])

Thisexpressionwillnotwork.Thefirstletterof"elevation"shouldbecapitalized.

ThisproblemhasbeenresolvedinArcGIS

Removeanyun-neededgridsfromtheArcMapTableofContents.Thetotalnumberofopengrids,includingtheonebeingcreatedbytheSpatialAnalystfunction,cannotexceed50.

Thislimitdoesnotapplytonon-gridrasters.AworkaroundistoconvertgridstorastersinArcCatalogwiththeExport>RasterToDifferentFormatoption

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 工作范文 > 行政公文

copyright@ 2008-2022 冰豆网网站版权所有

经营许可证编号:鄂ICP备2022015515号-1