Re: [gdal-dev] Raster statistics

2017-08-03 Thread Chris Waigl
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 ). Chris -- Christine (Chris) Wa

Re: [gdal-dev] Question in orthographic projection usage with PROJ4

2017-01-31 Thread Chris Waigl
> On Jan 31, 2017, at 7:46 AM, > niels.vonfestenberg-packisch@rohde-schwarz.com wrote: > > Hi all, > I would like to transform local cartesian coordinates into longitude and > latitude. According to all the documentation I could find, I tried doing this > (in C++, specifying a center fo

Re: [gdal-dev] gdalwarp - dataset shifted several metres compared to arcgis

2017-01-26 Thread Chris Waigl
, but that seems not so likely. Still, also EPSG:32232, which uses WGS 72. If you have rally old data it might be that. In any event, I’d advise inspecting the output of gdalinfo for both your target and your 2006 version. Chris -- Christine (Chris) Waigl - cwa...@alaska.edu - +1-907-474