Re: [gdal-dev] Multidimensional raster support in GDAL

2017-11-21 Thread Ari Jolma
Even Rouault kirjoitti 21.11.2017 klo 18:51: > Somehow the user must be able to write > > dataset_1 = GDAL.Open("WCS:URL?coverage=coverage_a") > dataset_2 = GDAL.Open("WCS:URL?coverage=coverage_a") > > So that for example dataset_1 is band c from coverage_a and dataset_2 is > band d from t

Re: [gdal-dev] Multidimensional raster support in GDAL

2017-11-21 Thread Even Rouault
> > And the following GetCoverage request succeeds and returns 3 bands > > > > SERVICE=WCS&VERSION=1.0.0&REQUEST=GetCoverage&COVERAGE=multi&FORMAT=GEOT IF > > F_8&BBOX=15,48,16,49&bands=9,5,1&CRS=EPSG:4326&WIDTH=5&HEIGHT=5 > So it's a misconfiguration at 194.66.252.155 Yes, probably an issue wit

Re: [gdal-dev] Multidimensional raster support in GDAL

2017-11-21 Thread Ari Jolma
Even Rouault kirjoitti 21.11.2017 klo 16:06: > and different time steps as bands. My reading of the code and doc of the current driver is that time steps are treated as subdatasets, and not bands. Yes. My mistake. > I'm not sure if there are many people > using it since I believe the fu

Re: [gdal-dev] Multidimensional raster support in GDAL

2017-11-21 Thread Even Rouault
On mardi 21 novembre 2017 12:23:31 CET Ari Jolma wrote: > I'd like to get back to this topic a bit since I learned something from > the standards that I believe was not discussed in the earlier discussion > > https://lists.osgeo.org/pipermail/gdal-dev/2017-October/047464.html > > In WCS 1 the ran

Re: [gdal-dev] Extracting keys & values from OSM Planet

2017-11-21 Thread Nikos Alexandris
* Even Rouault [2017-11-20 15:59:51 +0100]: On lundi 20 novembre 2017 15:48:20 CET Nikos Alexandris wrote: * Even Rouault [2017-11-07 22:37:02 +0100]: >> 2. Is OGR handling well the conversion from .osm to ESRI Shapefiles? > >Yes, but within the limits ot the shapefile, and particularly .dbf

Re: [gdal-dev] Vertical datum conversion for SRTM data

2017-11-21 Thread Michael Smith
Egm96 grid shift file available from http://download.osgeo.org/proj/vdatum/egm96_15/egm96_15.gtx Mike -- Michael Smith Remote Sensing/GIS Center US Army Corps of Engineers From: gdal-dev on behalf of Even Rouault Date: Tuesday, November 21, 2017 at 6:18 AM To: Subject: Re: [

Re: [gdal-dev] Vertical datum conversion for SRTM data

2017-11-21 Thread Even Rouault
On mardi 21 novembre 2017 10:25:42 CET Nikos Alexandris wrote: > Dear list, can GDAL handle the following exercise? > > Of interest is the Mean Sea Level, based on NASA’s SRTM. I am new in this > topic. Is the following equation correct? > > SRTM(MSL) = SRTM(WGS84) - MSL(WGS84) [*] > > > Is it

[gdal-dev] Multidimensional raster support in GDAL

2017-11-21 Thread Ari Jolma
I'd like to get back to this topic a bit since I learned something from the standards that I believe was not discussed in the earlier discussion https://lists.osgeo.org/pipermail/gdal-dev/2017-October/047464.html In WCS 1 the range of the data - remember the main conceptual division to the dom

Re: [gdal-dev] Merge many raster using maximum value where overlaps

2017-11-21 Thread Rutger
Sorry, i had the code embedded as raw text. It appears to be left out of my previous post, even though it worked in the preview. Here it is: max Python Regards, Rutger -- Sent from: http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html

[gdal-dev] Vertical datum conversion for SRTM data

2017-11-21 Thread Nikos Alexandris
Dear list, can GDAL handle the following exercise? Of interest is the Mean Sea Level, based on NASA’s SRTM. I am new in this topic. Is the following equation correct? SRTM(MSL) = SRTM(WGS84) - MSL(WGS84) [*] Is it “correct” to recommend the https://hydrosheds.cr.usgs.gov/ data set for Mean Sea

Re: [gdal-dev] Merge many raster using maximum value where overlaps

2017-11-21 Thread Rutger
You could also consider using the VRT pixel functions (since 2.2). I'm not sure about the performance when using the Python pixel functions (or perhaps with Numba), i haven't done any benchmarks. It might be a 'nice to have' if some statistics like max,min,std are added to the default functions if

Re: [gdal-dev] Extracting keys & values from OSM Planet

2017-11-21 Thread Nikos Alexandris
* Even Rouault [2017-11-20 15:59:51 +0100]: On lundi 20 novembre 2017 15:48:20 CET Nikos Alexandris wrote: * Even Rouault [2017-11-07 22:37:02 +0100]: >> 2. Is OGR handling well the conversion from .osm to ESRI Shapefiles? > >Yes, but within the limits ot the shapefile, and particularly .dbf