ArcGISEngine实现栅格数据裁剪
ArcGIS空间分析扩展模块提供了对栅格数据进行裁剪抽取操作的功能,在ArcGISEngine中可以调用相关的GP工具来对栅格数据操作。
- 调用空间分析ArcToolbox中Extract By Polygon工具
privatevoid ExtractByPolygon()
{
IVariantArray parameters =
newVarArrayClass();
string sDocPath =
@”C:\filerastergdb.gdb\raster”;
//地图文档路径
string sDocPath1 =
@”C:\filerasteroutput.gdb\outraster”;
parameters.Add(sDocPath);
string ss1 =
“‘545160.9444639276.889′;’545216.242 4639266.57’;’545227.355 4639308.375′”;;
parameters.Add(ss1);
parameters.Add(sDocPath1);
parameters.Add(“INSIDE”);
IGeoProcessor2 gp =
newGeoProcessorClass();
gp.OverwriteOutput =
true;
IGeoProcessorResult pResult =
null;
try
{
pResult =gp.Execute(“ExtractByPolygon”, parameters,
null);
IWorkspaceFactory pFWSF =
new ESRI.ArcGIS.DataSourcesGDB.FileGDBWorkspaceFactoryClass();
ESRI.ArcGIS.Geodatabase.IWorkspace pWS = pFWSF.OpenFromFile(@”C:\filerasteroutput.gdb”,
0) ;
ESRI.ArcGIS.Geodatabase.IRasterWorkspaceEx pRasterWs = pWS
asESRI.ArcGIS.Geodatabase.IRasterWorkspaceEx;
IRasterDataset pRasterDataset =pRasterWs.OpenRasterDataset(“outraster”);
IRasterLayer pRasterLayer =
newRasterLayerClass();
pRasterLayer.CreateFromDataset(pRasterDataset);
axMapControl1.AddLayer(pRasterLayer);
string mess =
null;
for (int k = 0; k < gp.MessageCount; k++)
{
mess+= gp.GetMessage(k) + “\n”;
}
MessageBox.Show(mess);
}
catch (Exception err)
{
string msg = err.Message;
string mess =
null;
for (int k = 0; k < gp.MessageCount; k++)
{
mess+= gp.GetMessage(k) + “\n”;
}
MessageBox.Show(mess);
}
}
- 调用IExtractionOp接口对栅格数据裁剪
publicvoid ExtractionRasterData(IGeoDataset pGeoDataset, ESRI.ArcGIS.Geometry.IPolygon
pPolygone, bool binside)
{
ESRI.ArcGIS.SpatialAnalyst.IExtractionOp pExtractionOp =
new ESRI.ArcGIS.SpatialAnalyst.RasterExtractionOpClass();
ESRI.ArcGIS.Geodatabase.IGeoDataset pGeoOutput = pExtractionOp.Polygon(pGeoDataset, pPolygone,binside);
}
转载自:https://blog.csdn.net/SYDBC/article/details/17120085