arcengine栅格数据使用总结 下载本文

内容发布更新时间 : 2024/4/25 10:04:12星期一 下面是文章的全部内容请认真阅读。

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(); rasterLayer.CreateFromFilePath(fileName); // fileName指存本地的栅格文件路径 axMapControl1.AddLayer(rasterLayer, 0); 2. 用IRasterDataset接口打开一个栅格数据集 IWorkspaceFactory workspaceFactory = new RasterWorkspaceFactory(); IWorkspace workspace; workspace = workspaceFactory.OpenFromFile(inPath, 0); //inPath栅格数据存储路径 if (workspace == null) { Console.WriteLine(\return; } IRasterWorkspace rastWork = (IRasterWorkspace)workspace; IRasterDataset rastDataset; rastDataset= rastWork.OpenRasterDataset(inName);//inName栅格文件名 if (rastDataset == null) { Console.WriteLine(\return; } 5、如何读取栅格数据的属性和遍历栅格数据 栅格数据的属性包括栅格大小,行数,列数,投影信息,栅格范围等等,见下面代码 (假设当前加载的栅格文件栅格值存储方式为:UShort类型) IRasterProps rasterProps = (IRasterProps)clipRaster; int dHeight = rasterProps.Height;//当前栅格数据集的行数 int dWidth = rasterProps.Width; //当前栅格数据集的列数 double dX = rasterProps.MeanCellSize().X; //栅格的宽度 double dY = rasterProps.MeanCellSize().Y; //栅格的高度 IEnvelope extent=rasterProps.Extent; //当前栅格数据集的范围 rstPixelType pixelType=rasterProps.PixelType; //当前栅格像素类型 IPnt pntSize = new PntClass(); pntSize.SetCoords(dX, dY); IPixelBlock pixelBlock = clipRaster.CreatePixelBlock(pntSize); IPnt pnt = new PntClass(); for (int i = 0; i < dHeight; i++) for (int j = 0; j < dWidth; j++) { pnt.SetCoords(i, j); clipRaster.Read(pnt, pixelBlock); if (pixelBlock != null) { object obj = pixelBlock.GetVal(0, 0, 0); MessageBox.Show( Convert.ToUInt32(obj).ToString()); } } 6、如何提取指定的范围的栅格数据 提取指定范围内的栅格数据通常用两种方法IRasterLayerExport(esriCarto), IExtractionOp, IExtractionOp2 ([url=]esriSpatialAnalyst[/url]),IRasterLayerExport接口提供的栅格数据提取功能有限,只能以矩形范围作为提取范围,而IExtractionOp接口提供了多边形,圆,属性,矩形等几种形式作为提取栅格数据. 1).IRasterLayerExport接口 IRasterLayerExport rLayerExport = new RasterLayerExportClass(); rLayerExport.RasterLayer = rasterLayer;// rasterLayer指当前加载的栅格图层 rLayerExport.Extent = clipExtent;//clipExtent指提取栅格数据的范围 if (proSpatialRef != null) rLayerExport.SpatialReference = proSpatialRef;// proSpatialRef当前栅格数据的投影信息 IWorkspaceFactory pWF = new RasterWorkspaceFactoryClass(); try { IWorkspace pRasterWorkspace = pWF.OpenFromFile(_folder, 0);// _folder指栅格文件保存路径 IRasterDataset outGeoDataset = rLayerExport.Export(pRasterWorkspace, code, strRasterType); //调用ISaveAs接口将导出的数据集保存 …………………….. } Catch(Exception ex) { Throw new Argumention(ex.Message); } 2.IExtractionOp接口(调用此接口前,应该先检查空间许可) IExtractionOp extraction = new RasterExtractionOpClass(); try { IGeoDataset geoDataset = extraction.Rectangle((IGeoDataset)clipRaster, clipExtent, true); IRaster raster = geoDataset as IRaster; if (raster != null) { IWorkspaceFactory WF = new RasterWorkspaceFactoryClass(); IWorkspace rasterWorkspace = WF.OpenFromFile(_folder, 0); ISaveAs saveAs = (ISaveAs)raster; saveAs.SaveAs(“Result.tif”, rasterWorkspace, \ } } catch (Exception ex) { MessageBox..Show(Ex.message); } 7.栅格数据重采样 栅格数据的重采样主要基于三种方法:最邻近采样(NEAREST),双线性 ILINEAR)和三次卷积采样(CUBIC)。