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 );
}