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