Re: [gdal-dev] python-gdal ComputStatistics return just two values
Lucena, Ivan
ivan.lucena at pmldnet.com
Fri Mar 27 15:11:42 EDT 2009
Actually...
GetStatistics dows the same as ComputeRasterStatistcs()
That is confusing...
> -------Original Message-------
> From: Lucena, Ivan <ivan.lucena at pmldnet.com>
> Subject: Re: [gdal-dev] python-gdal ComputStatistics return just two values
> Sent: Mar 27 '09 19:05
>
> Even,
>
> I see. You are right. With the absence of ComputeRasterStatistics() I got confused and I made the wrong
> assumption. BTW, why there is no wrapper for ComputeRasterStatistics()? That would work then:
>
> >>> bd1.SetStatistics(bd2.ComputeRasterStatistics())
>
> Not a big deal. The work around is easy.
>
> Thank you very much,
>
> Ivan
>
>
> > -------Original Message-------
> > From: Even Rouault <even.rouault at mines-paris.org>
> > Subject: Re: [gdal-dev] python-gdal ComputStatistics return just two values
> > Sent: Mar 27 '09 18:45
> >
> > Hi Ivan,
> >
> > I don't think there's a bug.
> >
> > The C API for ComputeBandStats is :
> >
> > CPLErr CPL_DLL CPL_STDCALL
> > GDALComputeBandStats( GDALRasterBandH hBand, int nSampleStep,
> > double *pdfMean, double *pdfStdDev,
> > GDALProgressFunc pfnProgress,
> > void *pProgressData );
> >
> > So, the python API is consistant with it.
> >
> > So it is expected that "bd1.SetStatistics(bd2.ComputeBandStats())" will fail
> > as it is provided with 2 values instead of the expected 4.
> >
> > Best regards,
> > Even
> >
> > Le Friday 27 March 2009 19:02:14 Lucena, Ivan, vous avez écrit :
> > > Hi there,
> > >
> > > I believe that ComputeBandStats() should return 4 values min, max, mean and
> > > stddev) [1] but the Python API is returning just mean and stddev.
> > >
> > > python:
> > > >>> bd.ComputeBandStats()
> > >
> > > (9.0071713001936118, 11.350094331448236)
> > >
> > > >>> bd.ComputeRasterMinMax()
> > >
> > > (-1.7899999618530273, 31.149999618530273)
> > >
> > > gdalinfo:
> > >
> > > Band 1 Block=360x1 Type=Float32, ColorInterp=Gray
> > > Min=-1.790 Max=31.150
> > > Minimum=-1.790, Maximum=31.150, Mean=9.007, StdDev=11.350
> > > Unit Type: unspecified
> > > Metadata:
> > > STATISTICS_MINIMUM=-1.789999961853
> > > STATISTICS_MAXIMUM=31.14999961853
> > > STATISTICS_MEAN=9.0071713001936
> > > STATISTICS_STDDEV=11.350094331448
> > >
> > > Is that a bug?
> > >
> > > See what happens in that case:
> > > >>> bd1.SetStatistics(bd2.ComputeBandStats())
> > >
> > > Traceback (most recent call last):
> > > File "<interactive input>", line 1, in <module>
> > > File "E:\Python25\Lib\site-packages\osgeo\gdal.py", line 729, in
> > > SetStatistics return _gdal.Band_SetStatistics(*args)
> > > TypeError: Band_SetStatistics() takes exactly 5 arguments (2 given)
> > >
> > > My best regards,
> > >
> > > Ivan
> > >
> > > [1] http://www.gdal.org/gdal_8h.html#a93b6b6ea6e71017ce25524e1a9ef1e3
> > >
> > > _______________________________________________
> > > gdal-dev mailing list
> > > gdal-dev at lists.osgeo.org
> > > http://lists.osgeo.org/mailman/listinfo/gdal-dev
> >
> >
> >
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
>
More information about the gdal-dev
mailing list