GDAL ComputeRasterMinMax .aux.xml

这个函数的第一个参数bApproxOK,TRUE表示精确统计,速度慢,图像的所有像元都遍历一边,FALSE表示粗略统计,速度快,但是不一定准确。函数还有个返回值,是CPLErr类型,如果返回CE_None表示计算成功,其他的可以参考GDAL的相关说明。

第二个参数pdfMinMax,是用来存储统计出来的最小值和最大值。这个数组一般是一个double [2],第0个表示最小值,第1个表示最大值。

切记,调用了函数GDALRasterBand::ComputeRasterMinMax ()之后,GDAL不会把最大值和最小值写入.aux.xml中。如果想把最大值和最小值写入.aux.xml中,需要手动调用:
    CPLsnprintf( szValue, sizeof(szValue), "%.14g", dfMin );
    SetMetadataItem( "STATISTICS_MINIMUM", szValue );


    CPLsnprintf( szValue, sizeof(szValue), "%.14g", dfMax );
    SetMetadataItem( "STATISTICS_MAXIMUM", szValue );

通常情况下此函数和GDALRasterBand::GetMinimum()、GDALRasterBand::GetMaximum()配套使用。

只有直接调用函数GDALRasterBand::GetMinimum()和GDALRasterBand::GetMaximum() 获取不成功时才需要调用ComputeRasterMinMax函数。

并可以将ComputeRasterMinMax的计算的结果存入.aux.xml中:下次就不需要再统计一遍了 ,如果有.aux.xml文件,且.aux.xml中有最大最小值,直接调用函数GDALRasterBand::GetMinimum()和GDALRasterBand::GetMaximum()就能获取到准确的结果。

代码逻辑如下:

      int bGotMin, bGotMax;
      double adfMinMax[2];
      adfMinMax[0] = poBand->GetMinimum( &bGotMin );
      adfMinMax[1] = poBand->GetMaximum( &bGotMax );
      if( ! (bGotMin && bGotMax) )
      {
         GDALComputeRasterMinMax(poBand, FALSE, adfMinMax);
         CPLsnprintf( szValue, sizeof(szValue), "%.14g", dfMin );
         SetMetadataItem( "STATISTICS_MINIMUM", szValue );
         CPLsnprintf( szValue, sizeof(szValue), "%.14g", dfMax );
         SetMetadataItem( "STATISTICS_MAXIMUM", szValue );
      }



阅读更多
个人分类: GDAL
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭