Re: [gdal-dev] gdal_calc.py: produce median raster

2014-11-25 Thread Marius Jigmond
That's because numpy.median expects one array and keyword arguments, not three arrays. It's unclear to me what gdal_calc.py does under the hood (I admit I did not read it) but the correct way to get the median of the three bands would be: --calc median([A,B,C], 0) The strange part is that--calc m

Re: [gdal-dev] Call for discussion on RFC 51: RasterIO() improvements : resampling and progress callback

2014-11-25 Thread Even Rouault
Mike, The gaussian resampling is the one already available for gdaladdo, that was added per http://trac.osgeo.org/gdal/ticket/2137 and http://trac.osgeo.org/gdal/changeset/15323 Basically : - for a downscaling of factor between 1 and 2 (and for upscaling as well), it will use a 3x3 kernel - for

Re: [gdal-dev] GDAL Driver for MongoDB

2014-11-25 Thread Even Rouault
Shuai, I may have the opportunity to integrate your driver in the future. I just wanted to make sure you had the right to contribute this work under the X/MIT license. From some headers like in https://github.com/mongogis/mongodb-gdal- driver/blob/master/mongo/ogr_mongo.h, I see Nanjing Universi

Re: [gdal-dev] Call for discussion on RFC 51: RasterIO() improvements : resampling and progress callback

2014-11-25 Thread Mike Toews
On 25 November 2014 at 05:34, Even Rouault wrote: >> Otherwise I have not much to say about this RFC. "Blurring" has two r's, > > Fixed Gaussian blur also has another complexity, since it's output dimensions depend on the fftconvolve mode (using scipy.signal docs http://docs.scipy.org/doc/scipy-0

Re: [gdal-dev] gdal_calc.py: produce median raster

2014-11-25 Thread Simen Langseth
Dear Dr Chris Yesson Thanks for your quick response. I got the following error: ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()` Perhaps the expression for calc seems to be revised... On Tue, Nov 25, 2014 at 4:59 PM, Chris Yesson wrote:

Re: [gdal-dev] gdal_calc.py: produce median raster

2014-11-25 Thread Chris Yesson
Dear Simen, You can refer to each band explicitly using the --A_band option. Try something like this gdal_calc.py -A infile --A_band 1 -B infile --B_band 2 -C infile --C_band 3 --outfile outfile --calc median(A,B,C) - Chris *Dr Chris Yesson* *Institute of Zoology, Zoological Society of London*