arcengine栅格数据使用总结

arcengine栅格数据使用总结 两个星期以来一直与栅格数据打交道，对AO的栅格部分应该有了一定的理解，下面是自己的一点体会，希望高手指教：－） 1、栅格数据的存储类型 栅格数据一般可以存储为ESRI GRID（由一系列文件组成），TIFF格式（包括一个TIF文件和一个AUX文件），IMAGINE Image格式 在AE中一般调用ISaveAs接口来保存栅格数据 2、栅格数据集和栅格编目的区别 一个栅格数据集由一个或者多个波段（RasterBand）的数据组成，一个波段就是一个数据矩阵。对于格网数据（DEM数据）和单波段的影像数据，表现为仅仅只有一个波段数据的栅格数据集，而对于多光谱影像数据则表现为具有多个波段的栅格数据集 栅格编目（RasterCatalog）用于显示某个研究区域内各种相邻的栅格数据，这些相邻的栅格数据没有经过拼接处理合成一副大的影像图 3、IRasterWorkspaceEx与IRasterWorkspace ,IRsterWorkspace2的区别 1).IRasteWorkspaceEx接口主要是用来读取GeoDatabase中的栅格数据集和栅格编目 2) . IRasterWorkspace ,IRsterWorkspace2主要是用来读取以文件格式存储在本地的栅格数据 4、加载栅格数据（以存储在本地的栅格数据文件为例） 1.直接用IRasterLayer接口打开一个栅格文件并加载到地图控件 IRasterLayer rasterLayer = new RasterLayerClass(); (fileName); 用IRasterDataset接口打开一个栅格数据集 IWorkspaceFactory workspaceFactory = new RasterWorkspaceFactory(); IWorkspace workspace; workspace = (inPath, 0); ; return; } IRasterWorkspace rastWork = (IRasterWorkspace)workspace; IRasterDataset rastDataset; rastDataset= (inName);; return; } 5、如何读取栅格数据的属性和遍历栅格数据 栅格数据的属性包括栅格大小，行数，列数，投影信息，栅格范围等等，见下面代码 （假设当前加载的栅格文件栅格值存储方式为：UShort类型） IRasterProps rasterProps = (IRasterProps)clipRaster; int dHeight = ;; ; oString()); } } 6、如何提取指定的范围的栅格数据 提取指定范围内的栅格数据通常用两种方法IRasterLayerExport(esriCarto), IExtractionOp, IExtractionOp2 ([url=]esriSpatialAnalyst[/url])，IRasterLayerExport接口提供的栅格数据提取功能有限，只能以矩形范围作为提取范围，而IExtractionOp接口提供了多边形，圆，属性，矩形等几种形式作为提取栅格数据. 1).IRasterLayerExport接口 IRasterLayerExport rLayerExport = new RasterLayerExportClass(); = rasterLayer; } Catch(Exception ex) { Throw new Argumention; } 2．IExtractionOp接口(调用此接口前，应该先检查空间许可) IExtractionOp extraction = new RasterExtractionOpClass(); try { IGeoDataset geoDataset = ((IGeoDataset)clipRaster, clipExtent, true); IRaster raster = geoDataset as IRaster; if (raster != null) { IWorkspaceFactory WF = new RasterWorkspaceFactoryClass(); IWorkspace rasterWorkspace = (_folder, 0); ISaveAs saveAs = (ISaveAs)raster; (“”, rasterWorkspace, "TIFF"); } } catch (Exception ex) { MessageBox..Show; } 7．栅格数据重采样 栅格数据的重采样主要基于三种方法：最邻近采样（NEAREST），双线性 ILINEAR）和三次卷积采样（CUBIC）。 （1）.最邻近采样：它用输入栅格数据中最临近栅格值作为输出值。因此，在重采 样后的输出栅格中的每个栅格值, 都是输入栅格数据中真实存在而未加任何改变的值。这种方法简单易用，计算量小，重采样的速度最快。 （2）.双线性采样：此重采样法取待采样点（x，y）点周围四个邻点，在y方向（或X方向）内插两次，再在x方向（或y方向）内插一次，得到（x，y）点的栅格值。 （3）.三次卷积采样：这是进一步提高内插精度的一种方法。它的基本思想是增加邻点来获 得最佳插值函数。取待计算点周围相邻的16个点，与双线性采样类似，可先在某一方向上内插，如先在x方向上，每四个值依次内插四次，再根据四次的计算结果在y方上内插，最终得到内插结果 代码示例：采用双线性采样 IRasterGeometryProc rasterGeometryProc = new RasterGeometryProcClass(); , newCellSize, clipRaster);

Tag标签: 栅格数据 ArcEngine 标签：栅格 重分类 分类：AE二次开发

public static IRasterLayer SetViewShedRenderer(IRaster pInRaster,string sField,string sPath) {

IRasterDescriptor pRD = new RasterDescriptorClass(); (pInRaster, new QueryFilterClass(), sField); IReclassOp pReclassOp = new RasterReclassOpClass(); IGeoDataset pGeodataset=pInRaster as IGeoDataset; IRasterAnalysisEnvironment pEnv = pReclassOp as IRasterAnalysisEnvironment; IWorkspaceFactory pWSF=new RasterWorkspaceFactoryClass(); IWorkspace pWS = (sPath, 0); = pWS; object objSnap = null; object objExtent = ; , ref objExtent, ref objSnap); = ; IRasterLayer pRLayer = new RasterLayerClass(); IRasterBandCollection pRsBandCol = pGeodataset as IRasterBandCollection; oString(); }

pFea = (); } = pTable; } else if (pLyr is IRasterLayer) { IRasterLayer pRlyr = pLyr as IRasterLayer; IRaster pRaster = ; IRasterProps pProp = pRaster as IRasterProps; = ; if == { IRasterBandCollection pBcol = pRaster as IRasterBandCollection; IRasterBand pBand = (0); ITable pRTable = ;

DataTable pTable = new DataTable(); for (int i = 0; i < i++)

ICursor pCursor= (null, false); IRow pRrow= (); while (pRrow != null) { DataRow pRow = (); for (int i =0 ;i{ pRow = (i).ToString () ; }

pRrow = (); } = pTable; }

}

IRasterWorkspace2 IRasterDataset CreateRasterDataset C# public IRasterDataset CreateFileRasterDataset(string directoryName, string fileName) { .......................................................................... ? ?? ?? ?? ?IWorkspace pWorkspace = 0); ? ?? ?? ?? ?IConversionOp pConversionOp = new RasterConversionOpClass(); ? ?? ?? ?? ?ISpatialReference pSpatialReference = ;

? ?? ?? ?? ?IGeoDataset pGeoDataset;