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

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

关于ArcGIS开发中计算两个几何对象之间距离的说明

在ArcGIS的二次开发中我们可以通过接口“IProximityOperator”的“ReturnDistance”方法来计算一个几何对象与另一个几何对象之间的距离。该接口被BezierCurve、CircularArc、EllipticArc、Envelope、GeoEllipse、GeoPolygon、GeoPolyline、MultiPatch、Multipoint、Point、Polygon、Polyline等几何对象实现。

// 定义用于计算的几何对象
IPoint point1 = new PointClass();
point1.X = 120;
point1.Y = 30;
// 定义计算距离的几何对象
IPoint point2 = new PointClass();
point2.X = 110;
point2.Y = 35;
// 获取IProximityOperator接口
IProximityOperator proximityOperator = point1 as IProximityOperator;
// 计算并返回几何对象point1到point2的距离
return proximityOperator.ReturnDistance(point2);

关于C#中读取ArcGIS本地切片缓存的说明

在ArcGIS中有两种存储在本地文件系统中的切片缓存格式:松散型(Exploded)、紧凑型(Compact)。

一、切片参数:
两种切片缓存格式的顶层文件及文件夹逻辑是一致的:
arcgis-tile-cache-1
在“_alllayers”文件夹中存储所有切片缓存,并以不同级别分别存放对应的图片数据:
arcgis-tile-cache-2
在“conf.xml”文件中记录切片的投影类型、级别参数、存储格式等信息:
arcgis-tile-cache-3
SpatialReference:投影参数
TileOrigin:原点坐标
TileCols:切片宽度
TileRows:切片高度
LODInfos:各级别参数
LevelID:级别号
Scale:比例尺
Resolution:像素分辨率
TileImageInfo:图片格式
CacheStorageInfo:缓存方式

在“conf.cdi”文件中记录切片缓存的范围:
arcgis-tile-cache-4
XMin:最小X坐标
YMin:最小Y坐标
XMax:最大X坐标
YMax:最大Y坐标
SpatialReference:投影参数

二、松散型:
每幅切片按照其所在行列单独存放在对应的图片文件中:
arcgis-tile-cache-5
读取代码如下:

// 切片缓存文件夹
string path = "C:/tiles";
// 切片所在级别
int level = 0;
// 切片所在行号
int row = 0;
// 切片所在列号
int col = 0;
// 切片图像格式
string format = "png";
// 解析切片图像文件路径
string l = string.Format("L{0:D2}", level);
string r = string.Format("R{0:X8}", row);
string c = string.Format("C{0:X8}", col);
string img = string.Format("{0}/_alllayers/{1}/{2}/{3}.{4}", path, l, r, c, format);
// 获取切片图像数据
return Image.FromFile(img);

三、紧凑型:
分块存储在切片数据文件bundle和索引文件bundlx中:
arcgis-tile-cache-6
每个bundle文件最多存储16384幅切片,每个bundlx文件大小固定为81952字节。bundlx文件中前后分别有16个字节与切片索引无关,中间剩余的81920字节以每5个字节的频率重复,构成了对bundle文件的索引信息。
读取代码如下:

// 切片缓存文件夹
string path = "C:/tiles";
// 切片所在级别
int level = 0;
// 切片所在行号
int row = 0;
// 切片所在列号
int col = 0;
// 解析切片所在区块的行列号
int rowGroup = 128 * (row / 128);
int colGroup = 128 * (col / 128);
// 解析切片数据文件和索引文件路径
string l = string.Format("L{0:D2}", level);
string r = string.Format("R{0:X5}", rowGroup);
string c = string.Format("C{0:X5}", colGroup);
path = string.Format("{0}/_alllayers/{1}/{2}{3}", path, l, r, c);
FileStream bundlx = new FileStream(path + ".bundlx", FileMode.Open);
FileStream bundle = new FileStream(path + ".bundle", FileMode.Open);
// 定位切片索引位置
long offset = 128 * (col - colGroup) + (row - rowGroup);
bundlx.Seek(16 + 5 * offset, SeekOrigin.Begin);
// 读取切片数据位置
byte[] buffer = new byte[5];
bundlx.Read(buffer, 0, buffer.Length);
// 定位切片数据位置
offset = 0;
for (int i = 0, length = buffer.Length; i < length; i++)
{
    offset += (buffer[i] & 0xff) * (long)Math.Pow(256, i);
}
bundle.Seek(offset, SeekOrigin.Begin);
// 读取切片数据长度
buffer = new byte[4];
bundle.Read(buffer, 0, buffer.Length);
offset = 0;
for (int i = 0, length = buffer.Length; i < length; i++)
{
    offset += (buffer[i] & 0xff) * (long)Math.Pow(256, i);
}
// 读取切片图像数据
byte[] tiles = new byte[offset];
bundle.Read(tiles, 0, tiles.Length);
bundlx.Close();
bundle.Close();
MemoryStream memoryStream = new MemoryStream(tiles);
return Image.FromStream(memoryStream);