ArcGIS Engine栅格数据使用总结

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('Could not open the workspace.');
return;
}
IRasterWorkspace rastWork = (IRasterWorkspace)workspace;
IRasterDataset rastDataset;
rastDataset= rastWork.OpenRasterDataset(inName);//inName栅格文件名
if (rastDataset == null)
{
Console.WriteLine('Could not open the raster dataset.');
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接