[gdal-dev] Cascaded VRTs

2015-08-11 Thread Fabian Schindler
Hi list, Premise: I’m tinkering to add support for datasets from multiple sources (like different subdatasets) in EOxServer (Python, MapScript, …) Versions: GDAL 1.9.2, MapServer 6.2.2 The sample dataset I try visualize is a MODIS/HDF 4 file with 7+ subdatasets (see [1] for a description for

Re: [gdal-dev] Rasdaman VRT generation script

2014-08-26 Thread Fabian Schindler
Even, I was not aware that creating a VRT was that easy... shame on me. My script would have handled connection parameters and created a query string, but I guess there is not enough added value. Thanks for the info Fabian On 2014-08-25 16:06, Even Rouault wrote: Selon Fabian Schindler

[gdal-dev] Rasdaman VRT generation script

2014-08-25 Thread Fabian Schindler
Hi folks, I created a small script that generates VRT files that act as wrappers to add georeferencing to rasdaman arrays (which do not have this feature). I thought I could add this somewhere under /swig/python/samples/, or is this too specific for a certain driver? Regards, Fabian

Re: [gdal-dev] LAS support

2013-11-28 Thread Fabian Schindler
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Rui, IMHO point clouds are special form of vector data. But still, an OGR driver would not be a viable solution. Did you have a look at PDAL (http://www.pointcloud.org/)? It is similar to GDAL and seems to have a built-in driver for LAS files.

[gdal-dev] Patch: Improvements to the rasdaman driver

2013-11-18 Thread Fabian Schindler
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Devs, I implemented a couple of improvements of the GDAL rasdaman driver. Ticket with a detailed description of the changes and patch are available here: http://trac.osgeo.org/gdal/ticket/5298 Of course, I appreciate any feedback and I hope the

Re: [gdal-dev] GeoTIFF with LZW compression Overviews error

2013-02-14 Thread Fabian Schindler
Frank, I'm assuming you mean gdaladdo is run on y.tif, the output of gdalwarp, right? Yes, sorry, that's what I meant. I tried the steps with a manufactured input file of the same size and did not encounter the problem on my system. The most likely important difference is that I am not using

[gdal-dev] GeoTIFF with LZW compression Overviews error

2013-02-13 Thread Fabian Schindler
Dear devs, I'm experiencing strange errors when I try to add overviews to specific GeoTIFF files with internal tiling and LZW compression. This error shows only on some (about 3-5%) of the input files. This is how I process the file: $gdal_translate -gcp 0 0 19.09 49.22 -gcp 100 0 20.1

[gdal-dev] Reworking the rasdaman driver

2013-02-12 Thread Fabian Schindler
Dear devs, I've been tasked to improve or rewrite the current rasdaman driver implementation as it currently suffers a number of drawbacks: * it is quite slow (compared to other drivers and other rasdaman connection techiques) * it does not support creation of new rasters within the

Re: [gdal-dev] GCP reprojection

2012-12-21 Thread Fabian Schindler
Hi devs, Sorry for the noise, I figured it out. Seems like it worked all along with order 1, but I was just to blind to see. Regards, Fabian On 12/18/2012 04:12 PM, Fabian Schindler wrote: Sorry for the traffic, just forgot to say that with the command line utilities it works as expected

[gdal-dev] GCP reprojection

2012-12-18 Thread Fabian Schindler
Hey devs, I've spent a full day to find a solution to this, but unfortunately was not able to solve it, so you are my last hope :) I'm trying to reproject an image with GCPs to a CRS specified by WKT. I used `gdal.AutoCreateWarpedVRT` to retrieve the output image size + geotransform in

Re: [gdal-dev] GCP reprojection

2012-12-18 Thread Fabian Schindler
gdalwarp -order 1 tmp.tif output.tif Thanks, Fabian On 12/18/2012 04:06 PM, Fabian Schindler wrote: Hey devs, I've spent a full day to find a solution to this, but unfortunately was not able to solve it, so you are my last hope :) I'm trying to reproject an image with GCPs to a CRS specified

Re: [gdal-dev] TIFF overviews issue

2012-12-03 Thread Fabian Schindler
Perfect, your workaround works for me. As I'm mostly dealing with small fairly small images, the performance decrease may not be a concern for me. I'll stick with this solution until I can definitely say that the target system has libtiff 4.0. Thanks Even! On 12/01/2012 11:37 AM, Even

[gdal-dev] TIFF overviews issue

2012-11-30 Thread Fabian Schindler
Dear devs, I'm experiencing an strange issue concerning overviews when I try to add overviews via Python. Essentially I want to create overviews for Greyscale images encoded in RGB, but on one overview level I always get colored pixels in a specific rectangle area of the image. I created

Re: [gdal-dev] Error handlers in Python bindings

2012-07-12 Thread Fabian Schindler
Thanks for the answer, but unfortunately this does not work on my machine: from osgeo import gdal def f(err_class, err_code, msg): ... print err_class, err_code, msg ... gdal.PushErrorHandler(f) Traceback (most recent call last): File stdin, line 1, inmodule File

Re: [gdal-dev] Error handlers in Python bindings

2012-07-12 Thread Fabian Schindler
Thanks Even, I feared that this was the case. But the `UseException` function looks promising. Regards, Fabian On 07/12/2012 10:59 AM, Even Rouault wrote: Selon lpinnerluke.pin...@environment.gov.au: The python function needs to accept three arguments, error class, error code and error

[gdal-dev] Error handlers in Python bindings

2012-07-11 Thread Fabian Schindler
Hey folks, I'd like to know how I could set a custom error handling function in the GDAL Python bindings? Currently I was able to turn off error messages like so: gdal.PushErrorHandler(CPLQuietErrorHandler) What I want is something like this: def f(error): ...logger.error(error) ...

Re: [gdal-dev] [NEW] Multi-threaded warping

2012-06-12 Thread Fabian Schindler
Etienne, The main advantage of OpenMP is that you can compile the code with (-fopenmp on gcc) or without OpenMP support but the output *should* stay the same. So you can produce code that is highly optimized for multithreaded environments (with little overhead) and still have backwards

[gdal-dev] Bug in rasdaman dataset

2011-11-11 Thread Fabian Schindler
Dear devs, I'd like to draw attention to a bug in the rasdaman dataset (to be more specific: its GetGeoTransform method). I created a new ticket [1] and also added a resolving patch. I hope it will be applied to the trunk soon. Thanks, Fabian [1] http://trac.osgeo.org/gdal/ticket/4341