[gdal-dev] Raster statistics
bontepaarden at gmail.com
Thu Aug 3 11:51:22 PDT 2017
Thanks Chris for your reply.
I forgot to mention I'm not using GDAL with Python.
I use it with C++ and/or C#.
*Paul Meems *
Release manager, configuration manager
and forum moderator of MapWindow GIS.
Owner of MapWindow.nl - Support for
Dutch speaking users.
*The MapWindow GIS project has moved to GitHub
Download the latest MapWinGIS mapping engine.
Download the latest MapWindow 5 open source desktop application.
2017-08-03 20:05 GMT+02:00 Chris Waigl <cwaigl at alaska.edu>:
> I would not use gdal for this particular task. I presume you have the band
> data in a 2D numpy array. Then I’d get the 80th percentile for example with
> np.percentile() and use a boolean expression to generate a mask for the
> array (droneraster > perc80value ).
> Christine (Chris) Waigl - cwaigl at alaska.edu - +1-907-474-5483
> <(907)%20474-5483> - Skype: cwaigl_work
> Geophysical Institute, UAF, 903 Koyukuk Drive, Fairbanks, AK 99775-7320,
> On Aug 3, 2017, at 5:43 AM, Paul Meems <bontepaarden at gmail.com> wrote:
> I have a drone raster file which I want to use for some calculation.
> Before the calculation, I need to loose some extreme values.
> I want to do something like a percentile calculation where you get all
> values, order them and loose the top 10%.
> For this, I need to get all values first which can be slow when using a
> large file.
> I looked at the statistics (band.GetStatistics) but that doesn't work well.
> I thought I could use 2 times the standard deviation added to the mean to
> get roughly 97%.
> But with these statistics:
> Mean + 2*std is larger than the max.
> So I moved to the histogram. It is also very fast, but I'm not sure how to
> use it.
> I have this:
> 256 buckets from 1.53322 to 33.248:
> 410 77 66 66 65 58 56 45 42 87 57 72 61 65 68 70 73 82 93 ...
> Does this mean, bucket 1 = 410 that I have 410 pixels of value 1.53322 and
> the second bucket means I have 77 pixels between 1.53322 and 1.657? 1.657 =
> 1.53322 + ((33.248 - 1.53322)/256)
> Is this a good approach? Or can/should I use a different one.
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the gdal-dev