ArcGIS Engine 中对栅格数据的波段信息统计 (转)

先打开栅格文件所在的工作空间(文件),然后获取其所有的波段,访问每一个波段
有时候波段中已经有直方图或统计信息,有时候没有这些信息,可以使用ComputeStatsAndHist()函数对其进行计算
(数据量较大时,可能耗时较长)

IWorkspaceFactory workspaceFactory = new RasterWorkspaceFactory();
IWorkspace workspace;
workspace = workspaceFactory.OpenFromFile(filePath, 0);

IRasterWorkspace rastWork = (IRasterWorkspace)workspace;
IRasterDataset rastDataset = rastWork.OpenRasterDataset(pathBox.Text.Substring(pathBox.Text.LastIndexOf("\\") + 1));
IRasterDataset2 rd2 = rastDataset as IRasterDataset2;
IRaster raster = rd2.CreateFullRaster();
IRasterBandCollection rbc = (IRasterBandCollection)raster;
for (int i = 0; i < rbc.Count; i++)
{
IRasterBand rb = rbc.Item(i);
bool tmpBool ;
rb.HasStatistics(out tmpBool);
if(!tmpBool)
rb.ComputeStatsAndHist();
IRasterHistogram rh = rb.Histogram;
IRasterStatistics rs = rb.Statistics;
}

from: http://blog.csdn.net/rrrrssss00/article/details/6585211