月度归档： 2016 年 6 月

关于C#中使用ArcGIS开发接口“ISaveAs2.SaveAs”无法写入数据的说明

在ArcGIS的二次开发中我们可以通过接口“ISaveAs2”的“SaveAs”方法将数据保存到本地。但是由于该方法将返回一个IDataset对象来绑定输出，因此在.NET开发时我们必须使用“Marshal.ReleaseComObject”方法来将该IDataset对象释放，否则数据将无法完整写入本地文件中。以下为官方说明：
Note, SaveAs will return a RasterDataset, to prevent from ISaveAs holding the output, .NET ReleaseCOMObject needs to be called to release the referene to the output raster dataset: System.Runtime.InteropServices.Marshal.ReleaseComObject(outRasterDS);

// 保存数据
IDataset dataset = saveAs2.SaveAs(fileName, workspace, "TIFF");
// 释放IDataset对象
Marshal.ReleaseComObject(dataset);

关于ArcGIS开发中读写图层文件的说明

ArcGIS的图层文件后缀名为“.lyr”，其保存了ArcGIS中地图图层的相关信息，如图层名称、图层数据、图层样式等。在ArcGIS的二次开发中我们可以通过接口“ILayerFile”实现对该文件的访问。

一、保存图层文件：

// 定义要保存图层文件的图层
ILayer layer;
// 定义要保存到的图层文件地址
string lyr = "C:/test.lyr";
// 创建图层文件类
ILayerFile layerFile = new LayerFileClass();
// 新建图层文件
layerFile.New(lyr);
// 写入图层数据
layerFile.ReplaceContents(layer);
// 保存图层文件
layerFile.Save();

二、读取图层文件：

// 定义要读取的图层文件地址
string lyr = "C:/test.lyr";
// 通过地图控件直接读取图层文件
axMapControl.AddLayerFromFile(lyr);
// 创建图层文件类
ILayerFile layerFile = new LayerFileClass();
// 打开图层文件
layerFile.Open(lyr);
// 通过图层文件类读取图层
axMapControl.AddLayer(layerFile.Layer);

关于ArcGIS开发中设置栅格分级专题的说明

在ArcGIS的二次开发中我们可以通过接口“IRasterClassifyColorRampRenderer”设置栅格图层的分级专题显示。以下介绍了两种常用的栅格分级方法：使用ArcGIS内置的分级方法和自定义间断点的分级方法。

一、使用内置分级方法：
ArcGIS内置了“自然间断点”、“相等间隔”、“几何间隔”、“标准差”、“分位数”等分级方法。

// 定义要分级显示的栅格图层
IRasterLayer rasterLayer;
// 定义分级方法名称
string classifyMethod;
// 获取栅格数据
IRaster raster = rasterLayer.Raster;
IRasterBandCollection rasterBandCollection = raster as IRasterBandCollection;
IRasterBand rasterBand = rasterBandCollection.Item(0);
// 计算栅格统计直方图
if (rasterBand.Histogram == null)
{
    rasterBand.ComputeStatsAndHist();
}
// 创建分级方法对象
IClassify classify;
if (classifyMethod == "自然间断点")
{
    classify = new NaturalBreaksClass();
}
else if (classifyMethod == "相等间隔")
{
    classify = new EqualIntervalClass();
}
else if (classifyMethod == "几何间隔")
{
    classify = new GeometricalIntervalClass();
}
else if (classifyMethod == "标准差")
{
    classify = new StandardDeviationClass();
    // 设置偏差数据
    IRasterStatistics rasterStatistics = rasterBand.Statistics;
    IDeviationInterval deviationInterval = classify as IDeviationInterval;
    deviationInterval.Mean = rasterStatistics.Mean;
    deviationInterval.StandardDev = rasterStatistics.StandardDeviation;
    deviationInterval.DeviationInterval = 1;
}
else if (classifyMethod == "分位数")
{
    classify = new QuantileClass();
}
else
{
    classify = new NaturalBreaksClass();
}
// 设置分级数
int numClasses = 3;
classify.Classify(ref numClasses);
// 创建分级渲染对象
IRasterClassifyColorRampRenderer rasterClassifyColorRampRenderer = new RasterClassifyColorRampRendererClass();
// 设置分级方式
IRasterClassifyUIProperties rasterClassifyUIProperties = rasterClassifyColorRampRenderer as IRasterClassifyUIProperties;
rasterClassifyUIProperties.ClassificationMethod = classify.ClassID;
// 设置分级数
rasterClassifyColorRampRenderer.ClassCount = 3;
// 设置分级样式1
Color color = Color.FromArgb(255, 0, 0);
IFillSymbol fillSymbol = new SimpleFillSymbolClass();
IRgbColor rgbColor = new RgbColorClass();
rgbColor.Red = color.R;
rgbColor.Green = color.G;
rgbColor.Blue = color.B;
rgbColor.Transparency = color.A;
fillSymbol.Color = rgbColor;
rasterClassifyColorRampRenderer.set_Symbol(0, fillSymbol as ISymbol);
// 设置分级样式2
color = Color.FromArgb(0, 255, 0);
fillSymbol = new SimpleFillSymbolClass();
rgbColor = new RgbColorClass();
rgbColor.Red = color.R;
rgbColor.Green = color.G;
rgbColor.Blue = color.B;
rgbColor.Transparency = color.A;
fillSymbol.Color = rgbColor;
rasterClassifyColorRampRenderer.set_Symbol(1, fillSymbol as ISymbol);
// 设置分级样式3
color = Color.FromArgb(0, 0, 255);
fillSymbol = new SimpleFillSymbolClass();
rgbColor = new RgbColorClass();
rgbColor.Red = color.R;
rgbColor.Green = color.G;
rgbColor.Blue = color.B;
rgbColor.Transparency = color.A;
fillSymbol.Color = rgbColor;
rasterClassifyColorRampRenderer.set_Symbol(2, fillSymbol as ISymbol);
// 设置栅格图层分级专题
IRasterRenderer rasterRender = rasterClassifyColorRampRenderer as IRasterRenderer;
rasterLayer.Renderer = rasterRender;
// 设置分级标签，必须在应用栅格专题后设置
rasterClassifyColorRampRenderer.set_Label(0, "0");
rasterClassifyColorRampRenderer.set_Label(1, "1");
rasterClassifyColorRampRenderer.set_Label(2, "2");

通过以上方法设置“IClassify”和“IRasterClassifyColorRampRenderer”接口后，我们也可以继续设置需要分级的栅格数据，然后手动获取内置分级方法计算出来的间断点值：

// 设置需要分级的栅格数据
IRasterRenderer rasterRender = rasterClassifyColorRampRenderer as IRasterRenderer;
rasterRender.Raster = raster;
rasterRender.Update();
// 获取对应索引的间断点值
double value = rasterClassifyColorRampRenderer.get_Break(0);

获取了内置分级方法计算出来的间断点值后，我们可以再次使用自定义间断点方法去进行栅格分级。

二、自定义间断点方法：

// 定义要分级显示的栅格图层
IRasterLayer rasterLayer;
// 定义间断点值集合
IList<double> breakPoints;
// 定义分段颜色集合
IList<Color> colors;
// 定义分级标签
IList<string> labels;
// 创建分级渲染对象
IRasterClassifyColorRampRenderer rasterClassifyColorRampRenderer = new RasterClassifyColorRampRendererClass();
IRasterRenderer rasterRender = rasterClassifyColorRampRenderer as IRasterRenderer;
// 设置分级数
rasterClassifyColorRampRenderer.ClassCount = colors.Count;
// 设置分级间断点
for (int i = 0, length = breakPoints.Count; i < length; i++)
{
    rasterClassifyColorRampRenderer.set_Break(i, breakPoints[i]);
}
// 设置各分级样式
for (int i = 0, length = colors.Count; i < length; i++)
{
    Color color = colors[i];
    IFillSymbol fillSymbol = new SimpleFillSymbolClass();
    IRgbColor rgbColor = new RgbColorClass();
    rgbColor.Red = color.R;
    rgbColor.Green = color.G;
    rgbColor.Blue = color.B;
    rgbColor.Transparency = color.A;
    fillSymbol.Color = rgbColor;
    rasterClassifyColorRampRenderer.set_Symbol(i, fillSymbol as ISymbol);
}
// 设置栅格图层分级专题
rasterLayer.Renderer = rasterRender;
// 设置分级标签，必须在应用栅格分级专题后设置
for (int i = 0, length = labels.Count; i < length; i++)
{
    rasterClassifyColorRampRenderer.set_Label(i, labels[i]);
}

关于Visual Studio中解决“无法嵌入互操作类型……，请改用适用的接口”的说明

今天在使用Visual Studio进行COM组件二次开发时，出现了“无法嵌入互操作类型……，请改用适用的接口”的错误：
vs-embedded-interoperability-type-1
解决方法很简单，选中项目中引入的COM组件库，鼠标右键，选择属性，把“嵌入互操作类型”设置为False即可：
vs-embedded-interoperability-type-2

关于ArcGIS开发中栅格代数计算的说明

在ArcGIS的二次开发中我们可以通过接口“IMapAlgebraOp”的实现类“RasterMapAlgebraOpClass”进行栅格数据之间的代数运算，即实现栅格计算器的功能。

// 定义要计算的栅格数据
IRaster raster1;
IRaster raster2;
// 定义结果栅格保存文件
string file = "C:/result.tif";
// 获取结果栅格文件所在文件夹
string folder = Path.GetDirectoryName(file);
// 获取结果栅格文件名
string fileName = Path.GetFileName(file);
// 定义结果工作空间
IWorkspaceFactory workspaceFactory = new RasterWorkspaceFactoryClass();
IWorkspace workspace = workspaceFactory.OpenFromFile(folder, 0);
// 定义地图代数运算操作类
IMapAlgebraOp mapAlgebraOp = new RasterMapAlgebraOpClass();
// 设置地图代数运算结果工作空间
IRasterAnalysisEnvironment rasterAnalysisEnvironment = mapAlgebraOp as IRasterAnalysisEnvironment;
rasterAnalysisEnvironment.OutWorkspace = workspace;
// 绑定栅格数据到地图代数运算操作类中
IGeoDataset geoDataset1 = raster1 as IGeoDataset;
IGeoDataset geoDataset2 = raster2 as IGeoDataset;
mapAlgebraOp.BindRaster(geoDataset1, "raster1");
mapAlgebraOp.BindRaster(geoDataset2, "raster2");
// 定义计算函数，栅格数据名称使用中括号包围
string calc = "[raster1] + [raster2]";
// 开始栅格计算获取计算结果栅格
IRaster raster = mapAlgebraOp.Execute(calc) as IRaster;
// 保存结果栅格为TIFF文件
ISaveAs2 saveAs2 = raster as ISaveAs2;
IDataset dataset = saveAs2.SaveAs(fileName, workspace, "TIFF");
// 释放保存结果数据，完成TIFF文件的写入
Marshal.ReleaseComObject(dataset);
return raster;
 ©  LuYutao