[gdal-dev] Numpy and ogr

2014-11-26 Thread Elliot Hallmark
Hello, I was wondering if the python api allows for vectorizing functions over arrays of geometry object (ogr objects?) through numpy. Like can I call overlaps on two (N,1) geomtry arrays and get an array of bools? Is it using python loops to iterate through, or is it applying it as a ufunc? Tha

Re: [gdal-dev] Multi-band JPEG2000 support?

2014-11-26 Thread Even Rouault
Le mercredi 26 novembre 2014 22:08:43, jplacy2 a écrit : > Thanks Even. Is it true that a recompile from source is necessary to > support other drivers, namely the JasPer JP2 option? I haven't been able > to find any windows binaries that make that driver available. But perhaps > that option only s

Re: [gdal-dev] Multi-band JPEG2000 support?

2014-11-26 Thread jplacy2
Thanks Even. Is it true that a recompile from source is necessary to support other drivers, namely the JasPer JP2 option? I haven't been able to find any windows binaries that make that driver available. But perhaps that option only supports three bands as well? I haven't been able to find a def

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

2014-11-26 Thread Pieter Kempeneers
Simen, In case you can not find a solution with gdal_calc.py, you could try to use pkfilter from pktools (http://pktools.nongnu.org). Once installed (currently under Linux only), you can produce your median raster as follows: pkfilter -i infile -o outfile -dz 1 -f median The option -dz 1 indicat

Re: [gdal-dev] Multi-band JPEG2000 support?

2014-11-26 Thread Even Rouault
Le mercredi 26 novembre 2014 17:31:34, jplacy2 a écrit : > Hi Folks- > > I'm attempting to generate some 4-band JP2 files using the GDAL 1.11 Win64 > binary acquired from gisinternals.com. > > I quickly ran into a problem... JP2 driver doesn't like four bands: > > gdal_translate --debug on -of J

Re: [gdal-dev] How to get additional information about JPEG 2000 using GDAL + ECW/JPG2 Intergraph SDK

2014-11-26 Thread Even Rouault
Le mercredi 26 novembre 2014 15:58:21, Tomer a écrit : > I want to be able to display additional properties of JPEG 2000 that are > normally not available in programs like GDALINFO (with the SDK integrated). > > Those properties are for example: > PROGRESSION=LRCP/RLCP/RPCL: > CODESTREAM_ONLY=YES/

[gdal-dev] Multi-band JPEG2000 support?

2014-11-26 Thread jplacy2
Hi Folks- I'm attempting to generate some 4-band JP2 files using the GDAL 1.11 Win64 binary acquired from gisinternals.com. I quickly ran into a problem... JP2 driver doesn't like four bands: gdal_translate --debug on -of JP2OpenJPEG test.tif test.jp2 GDAL: GDALOpen(test.tif, this=0061CB

[gdal-dev] How to get additional information about JPEG 2000 using GDAL + ECW/JPG2 Intergraph SDK

2014-11-26 Thread Tomer
I want to be able to display additional properties of JPEG 2000 that are normally not available in programs like GDALINFO (with the SDK integrated). Those properties are for example: PROGRESSION=LRCP/RLCP/RPCL: CODESTREAM_ONLY=YES/NO: LEVELS=n: LAYERS=n: PRECINCT_WIDTH=n: PRECINCT_HEIGHT=n: TI

Re: [gdal-dev] curved polygon -> polygon

2014-11-26 Thread Even Rouault
Le mercredi 26 novembre 2014 12:17:38, Martin Landa a écrit : > Hi, > > I wonder why curved polygons are not approximated to polygons somehow > automatically. Eg. conversion such data fails > > $ ogr2ogr -f PostgreSQL PG:dbname=ruian_test > /vsicurl/http://vdp.cuzk.cz/vymenny_format/soucasna/2014

[gdal-dev] curved polygon -> polygon

2014-11-26 Thread Martin Landa
Hi, I wonder why curved polygons are not approximated to polygons somehow automatically. Eg. conversion such data fails $ ogr2ogr -f PostgreSQL PG:dbname=ruian_test /vsicurl/http://vdp.cuzk.cz/vymenny_format/soucasna/20141031_OB_564567_UKSH.xml.gz Parcely -overwrite Warning 1: Geometry to be inse

Re: [gdal-dev] Motion: Adopt RFC 50: OGR field subtypes

2014-11-26 Thread Even Rouault
Hi, > So: > > Motion : I move to adopt RFC 50: OGR field subtypes > > http://trac.osgeo.org/gdal/wiki/rfc50_ogr_field_subtype > > Starting with my +1 I declare this motion passed with support from JukkaR, TamasS, FrankW and myself. Even -- Spatialys - Geospatial professional services http:

[gdal-dev] Warping with +proj=ob_tran

2014-11-26 Thread Knut-Frode Dagestad
Hi, I have some NetCDF files tagged with this proj4 string: proj4=+proj=ob_tran +o_proj=longlat +lon_0=-40 +o_lat_p=22 +R=6.371e+06 +no_defs As ob_tran seems not supported, is there any way I can warp this onto a regular (unrotated) lonlat-grid? There is a long, 3 yr old discussion on th

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

2014-11-26 Thread Mike Toews
Hi Simen, Try this: gdal_calc.py -A rgb.tif --outfile ouput.tif --calc "median(A, axis=0)" The "axis=0" parameter is to perform a median calculation along the first dimension, which is the band dimension. Lean more here: http://docs.scipy.org/doc/numpy/reference/generated/numpy.median.html -Mik