标签： C#

C#开发指南

关于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]);
}

关于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;

关于SuperMap Objects .NET中计算沿线某个距离点方向的说明

我们可以通过查找指定距离点所在的前后两个线节点来计算线上某个距离点的方向：

// 定义要计算的线
GeoLine geoLine = new GeoLine();
// 定义沿线距离
double distance = 100;
// 定义要计算的剩余距离
double tempDistance = distance;
// 遍历线上所有点
for (int i = 0; i < geoLine.PartCount; i++)
{
    Point2Ds point2Ds = geoLine[i];
    for (int j = 1; j < point2Ds.Count; j++)
    {
        // 计算两个节点的距离并修正剩下需要沿线的距离
        Point2D point2D1 = point2Ds[j - 1];
        Point2D point2D2 = point2Ds[j];
        tempDistance -= Geometrist.Distance(new GeoPoint(point2D1), new GeoPoint(point2D2));
        // 判断是否到达指定的距离
        if (tempDistance <= 0)
        {
            // 根据当前所在前后节点计算并返回方向
            double dx = point2D2.X - point2D1.X;
            double dy = point2D2.Y - point2D1.Y;
            // 东北方向
            if (dx > 0 && dy > 0)
            {
                return Math.Atan(dx / dy) * 180d / Math.PI;
            }
            // 东南方向
            else if (dx > 0 && dy < 0)
            {
                return 180d - (Math.Atan(dx / -dy) * 180d / Math.PI);
            }
            // 西南方向
            else if (dx < 0 && dy < 0)
            {
                return 180d + (Math.Atan(-dx / -dy) * 180d / Math.PI);
            }
            // 西北方向
            else if (dx < 0 && dy > 0)
            {
                return 360d - (Math.Atan(-dx / dy) * 180d / Math.PI);
            }
            // 正北方向
            else if (dx == 0 && dy > 0)
            {
                return 0;
            }
            // 正南方向
            else if (dx == 0 && dy < 0)
            {
                return 180;
            }
            // 正东方向
            else if (dx > 0 && dy == 0)
            {
                return 90;
            }
            // 正西方向
            else if (dx < 0 && dy == 0)
            {
                return 270;
            }
            // 其他
            else
            {
                return 0;
            }
        }
    }
}
// 没有线则返回默认方向
return 0;

关于SuperMap Objects .NET中点串拟合飞行轨迹的说明

在三维GIS开发中，场景飞行是一种常用的展示手段，平滑的飞行轨迹能使三维展示更加流畅。下面介绍在SuperMap Objects .NET中如何将一系列的坐标点串拟合成平滑的飞行轨迹线：

// 定义轨迹点串，包括X、Y、Z三维坐标
double[][] coordinates = new double[99][];
// 定义相关点集
Point2Ds point2Ds = new Point2Ds();
Point2Ds heightItemPoint2Ds = new Point2Ds();
Point2Ds heightPoint2Ds = new Point2Ds();
// 定义曲线拟合段数
int cardinalSegment = 20;
// 遍历所有轨迹点
foreach (double[] coordinate in coordinates)
{
    // 根据X、Y坐标创建高度曲线
    Point2D point2D = new Point2D(coordinate[0], coordinate[1]);
    point2Ds.Add(point2D);
    heightItemPoint2Ds.Add(point2D);
    GeoCardinal heightItemGeoCardinal = new GeoCardinal(heightItemPoint2Ds);
    // 将高度曲线转换为几何线对象
    GeoLine heightItemGeoLine = heightItemGeoCardinal.ConvertToLine(cardinalSegment);
    // 将高度曲线长度和Z坐标定义为轨迹曲线的坐标
    point2D = new Point2D(heightItemGeoLine.Length, coordinate[2]);
    heightPoint2Ds.Add(point2D);
}
// 获取水平轨迹曲线
GeoCardinal geoCardinal = new GeoCardinal(point2Ds);
GeoLine geoLine = geoCardinal.ConvertToLine(cardinalSegment);
// 获取垂直轨迹曲线
GeoCardinal heightGeoCardinal = new GeoCardinal(heightPoint2Ds);
GeoLine heightGeoLine = heightGeoCardinal.ConvertToLine(cardinalSegment);
// 获取轨迹点集
Point2Ds allPoint2Ds = geoLine[0];
Point2Ds allHeightPoint2Ds = heightGeoLine[0];
// 拟合轨迹曲线
Point3Ds allPoint3Ds = new Point3Ds();
for (Int32 i = 0; i < allHeightPoint2Ds.Count; i++)
{
    allPoint3Ds.Add(new Point3D(allPoint2Ds[i].X, allPoint2Ds[i].Y, allHeightPoint2Ds[i].Y));
}
GeoLine3D geoLine3D = new GeoLine3D(allPoint3Ds);
// 返回轨迹曲线
return geoLine3D;
 ©  LuYutao