Re: [gdal-dev] Gdal_grid and lake bathymetry

2019-11-21 Thread Jan Hartmann
These are important options. It would be nice if they could be implemented in gdal_grid itself. Not sure how difficult that would be, though. Jan On 11/21/2019 5:21 PM, Rahkonen Jukka (MML) wrote: Hi Joaquim, GMT surface probably has all the features that I need. I followed the tutorial and

Re: [gdal-dev] ogr2ogr -tps with more than 1000 control points

2019-08-13 Thread Jan Hartmann
m Auftrag von Jan Hartmann Gesendet: Freitag, 9. August 2019 21:18 An: Even Rouault Cc: gdal-dev@lists.osgeo.org Betreff: Re: [gdal-dev] ogr2ogr -tps with more than 1000 control points Yes, I already noticed that. But a run with thousand gcp's is still quite fast. I am looking for a smart wa

Re: [gdal-dev] ogr2ogr -tps with more than 1000 control points

2019-08-09 Thread Jan Hartmann
ing common points on the borders of adjacent parts, and adding them to the control points. If you have any ideas, let me know. Jan On 8/9/2019 8:19 PM, Even Rouault wrote: On vendredi 9 août 2019 20:08:31 CEST Jan Hartmann wrote: Thanks Even. GDAL was built with Armadillo, so I'm going to

Re: [gdal-dev] ogr2ogr -tps with more than 1000 control points

2019-08-09 Thread Jan Hartmann
ault wrote: On vendredi 9 août 2019 13:31:57 CEST Jan Hartmann wrote: Hi all, I'm converting large historical vector maps with ogr2ogr, using -gcp and -tps to rubbersheet the polygons with control points. Works fine up to about a thousand control points, each in the form of "-gcp fromx f

[gdal-dev] ogr2ogr -tps with more than 1000 control points

2019-08-09 Thread Jan Hartmann
. Is this a bug? BTW, I am using gdal 2.1.0 Regards, Jan Hartmann Amsterdam ___ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] FW: Don't we have any ideas for GSoC 2017?

2017-03-01 Thread Jan Hartmann
erdam ,together with Edzer Pebesma, of the University of Muenster, the leading specialist of spatial programming in R. Regards, Jan Hartmann http://www.uni-muenster.de/Geoinformatics/en/institute/staff/index.php/119/Edzer_Pebesma http://www.uva.nl/over-de-uva/organisatie/medewerkers/conten

Re: [gdal-dev] gdal_grid: number of data points to use

2014-01-18 Thread Jan Hartmann
dedicated interpolation programs like Surfer, and I always found it very useful. Funding would be possible. Jan On 01/18/2014 01:50 PM, Even Rouault wrote: Le samedi 18 janvier 2014 12:48:23, Jan Hartmann a écrit : gdal_grid has the option to set a minimum and maximum number of data points to

[gdal-dev] gdal_grid: number of data points to use

2014-01-18 Thread Jan Hartmann
gdal_grid has the option to set a minimum and maximum number of data points to use for interpolation. How are these data points selected from the total number of points? Are they always the nearest points to the grid-point that is interpolated? Jan _

Re: [gdal-dev] gdalwarp -tps

2013-08-29 Thread Jan Hartmann
The -tps option can only be used with control points added to an unreferenced raster. It's done this way: gdal_translate \ -gcp scanx1 scany1 worldx1 wordy1 \ -gcp scanx2 scany2 worldx1 worldy1 \ ... etc \ in.tif out.tif gdalwarp -tps out.tif world.tif A georeferenced raster

Re: [gdal-dev] Geotransform - rotation and translation

2013-08-27 Thread Jan Hartmann
The only way would be then to concatenate the transformations: first rotate and then translate half a pixel width. The computation can be simplified by matrix multiplication: http://www.mathplanet.com/education/geometry/transformations/transformation-using-matrices Jan That won't work, since i

Re: [gdal-dev] concatenated VRT file

2013-08-07 Thread Jan Hartmann
information for each individual file. Perhaps it isn't possible; in that case I'll hard-copy the whole bunch. Jan On 6-8-2013 18:18, Eli Adam wrote: On Tue, Aug 6, 2013 at 2:21 AM, Jan Hartmann <mailto:j.l.h.hartm...@uva.nl>> wrote: Not sure if this is possible, but perhaps s

[gdal-dev] concatenated VRT file

2013-08-06 Thread Jan Hartmann
Not sure if this is possible, but perhaps someone here knows a trick. I use a VRT file as an index to a large set of detailed, tiled scans of the Netherlands, built with gdalbuildvrt. I recently discovered how to transform this map on the fly, so that it exactly aligns to older maps in differ

Re: [gdal-dev] GeomTransform in VRT-file

2013-07-08 Thread Jan Hartmann
, Jan Hartmann пишет: I'm struggling with the order of the Geotransform section in a vrt file. I would like to rotate an image, and computed the six parameters for the affine transformation. I can transform points exactly, using the following formula: x2 = a.x1 + b.y1 + c y2 = d.x1 + e.y1 + f

[gdal-dev] GeomTransform in VRT-file

2013-07-04 Thread Jan Hartmann
I'm struggling with the order of the Geotransform section in a vrt file. I would like to rotate an image, and computed the six parameters for the affine transformation. I can transform points exactly, using the following formula: x2 = a.x1 + b.y1 + c y2 = d.x1 + e.y1 + f Can anyone tell me wh

Re: [gdal-dev] NTV2 Grid Shift files and gdalwarp

2013-06-25 Thread Jan Hartmann
Hi Andre, As far as I know, the transformation from one coordinate system to another goes as follow: - convert source coordinate to latlon - compute grid shift, either with +towgs84 parameters or by a grid shift raster - convert resulting latlon to target coordinate Ntv2 is a Canadian form

[gdal-dev] NTV2 Grid Shift files and gdalwarp

2013-06-25 Thread Jan Hartmann
Hi all, I've been experimenting with grid-shift files in ntv2 format to rubbersheet old maps. I created a working ntv2-file as described in http://lists.maptools.org/pipermail/proj/2013-January/006539.html and added an epsg definition for PROJ and PostGIS. The thing works with cs2cs and with

Re: [gdal-dev] Reducing the size of a GeoTIFF file with significant blackfill

2013-06-17 Thread Jan Hartmann
Hi Peter, You can get a compress-ratio to about 10% of the original size with the GDAL option -co COMPRESS=JKEG -co PHOTOMETRIC=YCBCR. Perhaps even more for you sort of data. Of course that compression is lossy, so for archival purposes it is useless. For viewers, you won't note the differenc

Re: [gdal-dev] Feature request: gauss and other interpolations in gdalwarp

2013-04-10 Thread Jan Hartmann
Hi folks, just a general question: I opened this thread asking to implement a few gdaladdo filters in gdalwarp, and am happy to see that happen now. Would it be possible to add a few new filters to gdalwarp/gdaladdo? I'm thinking about the gauss filter for gdalwarp, and the unsharp mask filter

Re: [gdal-dev] Feature request: gauss and other interpolations in gdalwarp

2013-03-28 Thread Jan Hartmann
right up there with being able to access the warp API from Python. Is this at all possible? Should I submit a feature request on Trac, or what? Just let me know! Jack. -- mathuin at gmail dot com On Mon, Feb 18, 2013 at 8:17 AM, Jan Hartmann <mailto:j.l.h.hartm...@uva.nl>> wrote:

[gdal-dev] Feature request: gauss and other interpolations in gdalwarp

2013-02-18 Thread Jan Hartmann
Hi devs, Would it be possible to add gauss and other interpolations to gdalwarp? At the moment I georeference large scans to 2000*2000 tiles at the most detailed scale, and then create 2000*2000 tiles at resolutions of 2, 4 6 etc times the original scale, using gdaladdo and gauss or other i

[gdal-dev] Feature request: gauss and other interpolations in gdalwarp

2013-02-18 Thread Jan Hartmann
Hi devs, Would it be possible to add the gauss and other interpolations to gdalwarp? At the moment I georeference large scans to 2000*2000 tiles at the most detailed scale, and then create 2000*2000 tiles at resolutions of 2, 4 6 etc times the original scale, using gdaladdo and gauss or othe

[gdal-dev] Creating a nadgrid shift file

2013-01-17 Thread Jan Hartmann
(Sorry for crossposting, don't know if it is a gdal or proj problem ) This has been discussed on this list some time ago. I have been able to create a working nadgrid file with gdal, but it doesn't return the right coordinates. I did it as follows: - create a 4 band, float32 geotif, each pixe

Re: [gdal-dev] does gdal support multiple simultaneous writers to raster

2013-01-14 Thread Jan Hartmann
ssues with write consistency but of a different nature. The easiest approach is to assign different spatial regions to each worker server and run these linearly. Best regards, On Jan 12, 2013 8:05 AM, "Jan Hartmann" <mailto:j.l.h.hartm...@uva.nl>> wrote: I would be interes

Re: [gdal-dev] does gdal support multiple simultaneous writers to raster

2013-01-12 Thread Jan Hartmann
This is about the same acceleartion (2-3 times) as I got on jobs running for a few days. My impression is that distributed tile processing would give much more dramatic results. On 01/12/2013 03:57 PM, Even Rouault wrote: Le samedi 12 janvier 2013 15:08:55, Jan Hartmann a écrit : You

Re: [gdal-dev] does gdal support multiple simultaneous writers to raster

2013-01-12 Thread Jan Hartmann
I would be interested in an implementation. I'm preparing a proposal to georeference the complete cadastral map of the Netherlands in 1832 at a 10 cm/pixel scale with Cloud facilities. Gdalwarp is the central piece of software, and distributed processing capabilities would be very important. Co

Re: [gdal-dev] does gdal support multiple simultaneous writers to raster

2013-01-12 Thread Jan Hartmann
You probably know this, but there is an option to let gdalwarp use more cores: -wo NUM_THREADS=ALL_CPUS. It gives some improvement, but not really staggering. Splitting up operations over individual tiles would really fasten up things. Even if I use only one VM, I can define 32 cores, and it wo

Re: [gdal-dev] does gdal support multiple simultaneous writers to raster

2013-01-11 Thread Jan Hartmann
I would like to make a further suggestion: what about an option for gdalwarp to create *external* tiles, instead of one big output file? That's what I am doing now: just create one big georeferenced raster, and split it up in tiles, mostly 2000*2000. Accessing those tiles with a tileindex is re

Re: [gdal-dev] does gdal support multiple simultaneous writers to raster

2013-01-11 Thread Jan Hartmann
I'v tried this in a Cloud environment with multiple VMs writing to the same shared file location. I didn't use tiles, so the results where negative. It would be a great asset for Cloud environments, however. Is it in principle possible to build such a facility into GDAL, and if yes, what kind o

Re: [gdal-dev] Compression artifacts

2012-12-27 Thread Jan Hartmann
e 4GB limit you might need to add "-co bigtiff=if_safer". Good luck, Trent On Thu, Dec 27, 2012 at 5:38 AM, Jan Hartmann <mailto:j.l.h.hartm...@uva.nl>> wrote: Can anyone tell me wath is happening here? I have a large three-band raster (3.8G) look

Re: [gdal-dev] Compression artifacts

2012-12-27 Thread Jan Hartmann
On 12/27/2012 03:11 PM, Even Rouault wrote: Le jeudi 27 décembre 2012 15:11:23, Jan Hartmann a écrit : On 12/27/2012 02:41 PM, Even Rouault wrote: I've just pushed a fix since building overviews with "--config COMPRESS_OVERVIEW JPEG -- config PHOTOMETRIC_OVERVIEW YCBCR" when

Re: [gdal-dev] Compression artifacts

2012-12-27 Thread Jan Hartmann
On 12/27/2012 02:41 PM, Even Rouault wrote: I've just pushed a fix since building overviews with "--config COMPRESS_OVERVIEW JPEG -- config PHOTOMETRIC_OVERVIEW YCBCR" when there's a .msk didn't work. OK Even, thanks. This means that I should use gdal-svn? Not necessarily. I think you can wo

Re: [gdal-dev] Compression artifacts

2012-12-27 Thread Jan Hartmann
On 12/27/2012 02:27 PM, Even Rouault wrote: The artifacts are due to the lossy characteristics of JPEG compression that are particularly visible on borders between nodata pixel and valid pixels. An option would be to use nearblack to create a .msk mask file : nearblack -setmask -o big_with_mask

Re: [gdal-dev] Compression artifacts

2012-12-27 Thread Jan Hartmann
artifacts can be prevented? On 12/27/2012 01:44 PM, Dmitry Baryshnikov wrote: Hi Jan, You can cut result raster with the geometry (shape file) using gdalwarp with cutline option (http://gdal.org/gdalwarp.html). Best regards, Dmitry 27.12.2012 16:38, Jan Hartmann ?: Can anyone tell me wath

[gdal-dev] Compression artifacts

2012-12-27 Thread Jan Hartmann
Can anyone tell me wath is happening here? I have a large three-band raster (3.8G) looking like this: http://mapserver.sara.nl/map2a.png I compress it (from information from this list a few days ago) with: gdal_translate -co COMPRESS=JPEG -co PHOTOMETRIC=YCBCR big.tif small.tif gdaladdo --conf

Re: [gdal-dev] Q: on gdalwarp of ecw file

2012-12-18 Thread Jan Hartmann
Thanks Jukka, I always wanted to learn Finnish :-). Almost as difficult as Dutch. Cheers, Jan On 12/18/2012 01:16 PM, Jukka Rahkonen wrote: Here is a fresh table comparing resampling methods and cgi vs. fast-cgi with Mapserver and with raster map data. http://latuviitta.org/Apps4Finland-2012

Re: [gdal-dev] Q: on gdalwarp of ecw file

2012-12-18 Thread Jan Hartmann
wrote: I've found setting PHOTOMETRIC=YCBCR to allow much higher compression when using jpeg compression in gtifs. Mike -- Michael Smith US Army Corps Remote Sensing GIS/Center From: Jan Hartmann mailto:j.l.h.hartm...@uva.nl>> Date: Tuesday, December 18, 2012 5:47 AM To: Juk

Re: [gdal-dev] Q: on gdalwarp of ecw file

2012-12-18 Thread Jan Hartmann
Jan On 12/18/2012 12:33 PM, Jukka Rahkonen wrote: Jan Hartmann uva.nl> writes: Hi Even, are there any benchmarks to compare (uncompressed) gtif with the three formats above? My production maps are always tiled to 2000*2000 pixels, all zoomlevels precomputed. V

Re: [gdal-dev] Q: on gdalwarp of ecw file

2012-12-18 Thread Jan Hartmann
On 12/17/2012 09:14 PM, Even Rouault wrote: The performance vs size when zooming out is impressive with ECW, but it does not look like I have ability to create ECW. Are there better formats then GTif in the list below? JPEG2000 or MrSID are somewhat similar to ECW. Hi Even, are there any bench

Re: [gdal-dev] gdalwarp performance with many GCPs

2012-12-12 Thread Jan Hartmann
On 12/12/2012 01:54 PM, Knut-Frode Dagestad wrote: The GCPs are completely regular, every 40 pixel in each direction. Warping quality is very good (only slow), so the GCPs are probably also ok. I would appreciate very much if you could download this file and test the command below: http://d

Re: [gdal-dev] gdalwarp performance with many GCPs

2012-12-12 Thread Jan Hartmann
and promising. For the mentioned file it takes about 2 minutes with "-et 5" (low accuracy), and 12 minutes without this switch. Without -tps it takes less than 0.5 seconds. The machine is quite decent with 16 GB memory and SSD disk. Best regards from Knut-Frode On 12. des. 2012

Re: [gdal-dev] gdalwarp performance with many GCPs

2012-12-12 Thread Jan Hartmann
Hi Knut, What do you mean by "very slow"? I regularly use gdalwarp -tps on much larger rasters with a few thousand gcp-s, and never noticed unacceptable delays. Do you have very little physical memory, or many parallel processes running? Jan On 12/12/2012 01:12 PM, Knut-Frode Dagestad wrote

Re: [gdal-dev] EPSG 8.0 Upgrade

2012-12-06 Thread Jan Hartmann
Hi Frank, where do I find which datums have new datum shift solutions? Jan On 12/06/2012 01:18 AM, Frank Warmerdam wrote: Folks, At the request of Howard Butler, I have upgraded PROJ.4, libgeotiff, and GDAL to use the EPSG 8.0 database in development "trunk". This was accomplished, as usual b

[gdal-dev] Python programs on multi-core machines

2012-08-13 Thread Jan Hartmann
I'm working on multi-core VMss in a Cloud environment, that access their data on a central dataserver via NFS. Parallellizing jobs for different map sheets gives huge accelerations for C-programs like gdaladdo, but there seems to be a problem with Python-based programs like rgb2pct.py. Consider

[gdal-dev] gdaltransform with gcps from file

2012-06-23 Thread Jan Hartmann
Silly question, perhaps. I use gdaltransform with gcps given at the command line. According to the documentation, gcps can also be put into the parameter to the program. What is the format of that file? Most of the time I get an error message like "file fomat not supported". With four numbers

Re: [gdal-dev] Warping without shear

2011-12-17 Thread Jan Hartmann
;t be too difficult to adapt the -order 1 flag so as to leave out the shearing. Jan On 12/16/2011 08:21 PM, Joaquim Luis wrote: On 16-12-2011 18:50, Jan Hartmann wrote: On 16-12-2011 18:25, Joaquim Luis wrote: Jan, What you want is to apply a transformation called (at least on the M

Re: [gdal-dev] Warping without shear

2011-12-16 Thread Jan Hartmann
On 16-12-2011 18:25, Joaquim Luis wrote: Jan, What you want is to apply a transformation called (at least on the Matlab parlance) "Linear conformal" which an afine transform but without the the shear terms. Mirone has it implemented and if your image has the GCPS embedded the operation is rea

[gdal-dev] Warping without shear

2011-12-16 Thread Jan Hartmann
Hi all, Is it possible to do a first order GCP-warp *without* shearing? I have a rectangular scan with very imprecise control points, and I want it georeferenced without any distortion. With a first order gdalwarp I get back a map in the form of a parallellogram. Is there a way to get back th

Re: [gdal-dev] adding georeferencing information

2011-09-22 Thread Jan Hartmann
2011 12:37, Jan Hartmann <mailto:j.l.h.hartm...@uva.nl>> wrote: Found the answer to this one: use the "listgeo" and "geotifcp" utilities from the Geotiff library (http://trac.osgeo.org/geotiff/). Did you not receive the email I sent to the list on 6 Septembe

Re: [gdal-dev] adding georeferencing information

2011-09-22 Thread Jan Hartmann
Thanks Jean-Claude, I'll try it out as soon as I manage to compile GDAL trunk. It's still broken at the moment ... Jan On 09/22/2011 02:26 PM, Jean-Claude Repetto wrote: Le 22/09/2011 13:37, Jan Hartmann a écrit : Found the answer to this one: use the "listgeo" and "

[gdal-dev] Using the alpha channel to mask parts of an image

2011-09-22 Thread Jan Hartmann
Hi, I have a set of unreferenced RGB images with black borders I want to cut off. I do this by adding a fourth alpha band to the image when warping it, and filling this with 255 within the data rectangle, and 0 outside, using gdal_rasterize ( I have the exact georeferenced boundaries of the

Re: [gdal-dev] adding georeferencing information

2011-09-22 Thread Jan Hartmann
ss there are binary versions for other distributions. Alternatively, you could use fwtools or osgeo4w. Jan On 09/05/2011 04:46 PM, Jan Hartmann wrote: Hi all, Is there a quick way to add back georeferencing information to a geotif file that has been processed by Imagemagick? I have lots of

Re: [gdal-dev] mapfile config options - scan directory

2011-09-20 Thread Jan Hartmann
Hi all, What is the relationship between the "PER_DATASET_ALPHA" mask of the RGB bands in a raster dataset, and the Alpha Channel? I created a four-band tiff file with gdalwarp from an RGB original. On creation, the fourth band has all values set to 255, and the whole picture is visible. Howe

Re: [gdal-dev] Numerical instability with thin plate spline transform

2011-09-17 Thread Jan Hartmann
x3Dataset]+0x30): undefined reference to `GDALDataset::CloseDependentDatasets()' gdaldem.o:(.data.rel.ro._ZTV22GDALColorReliefDataset[vtable for GDALColorReliefDataset]+0x30): undefined reference to `GDALDataset::CloseDependentDatasets()' Jan On Thu, Sep 8, 2011 at 5:23 AM, Jan H

[gdal-dev] Numerical instability with thin plate spline transform

2011-09-08 Thread Jan Hartmann
Not sure whether this can be considered a bug, so I give it for what it is worth. I'm doing thin plate spline transformation from one set of projected coordinates to another. Both sets have values between -60 and 60 (meters). A typical set of gcps looks like: -gcp 62402 -74383 18191

[gdal-dev] adding georeferencing information

2011-09-05 Thread Jan Hartmann
Hi all, Is there a quick way to add back georeferencing information to a geotif file that has been processed by Imagemagick? I have lots of georeferenced images that I want to sharpen up a bit, but ImageMagick throws the georeference headers away. Is there a way to transport them from the ori

[gdal-dev] Re: Warping an image with both a geotransform and control points

2011-07-07 Thread Jan Hartmann
Hi all, This has come up before, but I'm still struggling for an answer. I want to rubbersheet an already projected image with control points, so as to remove small errors.On a suggestion by Frank (http://lists.osgeo.org/pipermail/gdal-dev/2009-February/019775.html) I produced a VRT file of t

Re: [gdal-dev] Re: gdal_rasterize 1.8.0 options

2011-02-07 Thread Jan Hartmann
On 2011-02-07 15:45, Frank Warmerdam wrote: On 11-02-07 06:25 AM, Jan Hartmann wrote: Is there a place to put al the links on the stere-sterea subject together? In my experience, discussions about projection parameters tend to get fragmented, and even if a solution has been found, old

Re: [gdal-dev] Re: gdal_rasterize 1.8.0 options

2011-02-07 Thread Jan Hartmann
As far as the Dutch projection is concerned (EPSG:28992), Peter's summary is correct: it first projects the ellipsoid to the sphere and then projects the sphere stereographically to the plane. I don't know if that is the way PROJ4 computes the transformation, but I can garantee that using "ster

Re: [gdal-dev] Re: gdal_rasterize 1.8.0 options

2011-02-07 Thread Jan Hartmann
For the Dutch case: in 2005 the PROJ-definition for the Dutch national system (which uses "Double Stereographic") was changed from "stere" to "sterea". See the thread starting at: http://lists.osgeo.org/pipermail/gdal-dev/2005-November/007045.html. The thread refers to a bug ticket at http://b

Re: [gdal-dev] gdalwarp with gauss resampling

2011-01-13 Thread Jan Hartmann
Oh, you mean gdalbuildvrt? I say it in your previous posting. I did try it out and it solves my problem. No need at all for a separate gauss-flag to gdalwarp. The whole idea of chaining pipelined operations this way is brilliant! Thanks, Jan On 01/13/11 22:18, Even Rouault wrote: is i

Re: [gdal-dev] gdalwarp with gauss resampling

2011-01-13 Thread Jan Hartmann
On 01/12/11 20:17, Frank Warmerdam wrote: Yes, you can dump overviews to a standalone file using the "dumpoverviews" commandline application. I don't think it is built by default so you might need to go into gdal/apps and make it explicitly. Yes, this works after compiling: I get a stand

[gdal-dev] gdalwarp with gauss resampling

2011-01-12 Thread Jan Hartmann
Is it possible to use gauss resampling with gdalwarp? The current resampling options (cubic, bilinear etc) work very well as long as the pixels in input and output rasters are about the same size, but not very well when downsampling to a much coarser scale. Gdaladdo has gauss resampling, but on

Re: Fwd: Re: [gdal-dev] FWTools and GDAL 1.7.0

2011-01-06 Thread Jan Hartmann
Three months ago I put an announcement on the MapServer-Users list about a shell script I wrote to compile the whole MapServer-GDAL-PostGIS-PLR suite, with and without PHP/Python, in a separate user directory using only two commands. Separate installations can coexist in separate directory tree

Re: [gdal-dev] Frank Eats Crow - or Fixing PixelIsPoint Interpretation

2010-11-16 Thread Jan Hartmann
On 11/16/10 02:40, Frank Warmerdam wrote: Folks, After 15 years of denial, and bullheadedness it is now time for me to admit my interpretation of PixelIsPoint in GeoTIFF files is wrong. To that end I have prepared a brief RFC discussing how I intend to fix it. I'd appreciate comment before

Re: [gdal-dev] Re: Obtain more information with Ogrinfo

2010-05-10 Thread Jan Hartmann
The manual page says: *-al*: List all features of all layers (used instead of having to give layer names as arguments). This is not a really clear answer for Antonio's question: when I read this the first time, I thought this flag was meant to display information about more layers, not

Re: [gdal-dev] Re: Thin Plate Spline

2010-05-07 Thread Jan Hartmann
precision of sub-pixel. Attached a ZIP with the source image and C source code used to generate it. Command line used to warp : gdalwarp -tps src_grid.tif dst_grid.tif Best regards, Even Le Tuesday 04 May 2010 15:14:34 Jan Hartmann, vous avez écrit : Thanks, Jukka, that's exactly wh

Re: [gdal-dev] Re: Thin Plate Spline

2010-05-04 Thread Jan Hartmann
Thanks, Jukka, that's exactly what I did. I used QGIS to test the coordinates in both images. Jan On 05/04/10 15:08, Jukka Rahkonen wrote: Jan Hartmann uva.nl> writes: you can see a screenshot of the original image (right) and the georeferenced one (left, rotated almost 90

Re: [gdal-dev] Thin Plate Spline

2010-05-04 Thread Jan Hartmann
t happens in between? Jan On 05/01/10 00:21, Even Rouault wrote: It should work... Perhaps you can share how to reproduce what you see Le Saturday 01 May 2010 00:17:06 Jan Hartmann, vous avez écrit : Hi, I'm rubber-sheeting old, deformed maps using the thin plate spline option of gdalwarp

Re: [gdal-dev] Thin Plate Spline

2010-05-04 Thread Jan Hartmann
h segment. In that case, I ended up transforming each segment and then stitching the segments together - it worked remarkably well. Hope that helps, best wishes, Peter Jan Hartmann wrote: Hi, I'm rubber-sheeting old, deformed maps using the thin plate spline option of gdalwarp (-tps)

[gdal-dev] Thin Plate Spline

2010-04-30 Thread Jan Hartmann
Hi, I'm rubber-sheeting old, deformed maps using the thin plate spline option of gdalwarp (-tps). My impression was that the control points would be transformed to exactly their georeferenced locations, but in the georeferenced image the control poimts don't quite align with the gcp-coordinat

[gdal-dev] transform PostGIS geometries with control points, like gdalwarp -gcp

2010-04-19 Thread Jan Hartmann
(Sorry for cross-posting, this is about both Gdal and PostGIS) Hi all, Would it be possible / make sense to add a procedure to PostGIS to transform geometries based on control points? We now have the base transforms: scale rotate and translate, and IMHO it would not be very complicated to add

Re: [gdal-dev] Converting contour lines (Shape files) into DEM

2010-03-30 Thread Jan Hartmann
Frank, is there a description of the way the smoothing is done, or should I look it up in the code? BTW: there is a description of the algorithm used in GRASS to create contour lines: https://trac.osgeo.org/grass/browser/grass/trunk/raster/r.contour/README This looks much more adapted to the

[gdal-dev] GDAL WMS driver

2010-03-21 Thread Jan Hartmann
On 2010-03-19 23:00, Tamas Szekeres wrote: You might want to try to use the GDAL WMS/TMS driver , as something like: Very interesting. How about using gdalwarp to georeference WMS images over the Internet? Any hints how to optimise this, as against using

[gdal-dev] gdal --with-php

2010-03-18 Thread Jan Hartmann
Hi all, What is the status of the php-bindings for gdal (./configure --with-php)? I found a message from Frank (2006) stating that it is not well supported (http://lists.osgeo.org/pipermail/gdal-dev/2006-March/008258.html), but that there are php-bindings for OGR, developed by DMSolutions (h

Re: [gdal-dev] Motion: Commit access for

2010-02-26 Thread Jan Hartmann
On 26-2-2010 18:08, Frank Warmerdam wrote: Frank Warmerdam wrote: Motion: Commit access on the GDAL/OGR project is extended to Brian Case. --- Brian Case (Winkey in irc/trac) has been using and contributing to GDAL/OGR for some time. He has provided a number of patches, including those in t

Re: [gdal-dev] Re: Add GCP's to a file

2010-02-19 Thread Jan Hartmann
On 19-2-2010 11:46, Manuel Massing wrote: Hi, Gdal_translation has on option [-gcp pixel line easting northing [elevation]]* this should also work using the VRT output format, which has the advantage that the raster data doesn't need to be written to an intermediate file (VRT "wr

Re: [gdal-dev] Color-palette file for use in rgb2pct

2010-02-01 Thread Jan Hartmann
Hi Frank, could this be added to the documentation? Jan On 1-2-2010 12:37, Frank Warmerdam wrote: Luuk Schaminee wrote: Hi, I am trying to convert a 24 bit RGB TIFF file to a 8 bit TIFF file. Without the –pct option this works perfectly but I want to control the created colors with a palett

Re: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-18 Thread Jan Hartmann
On 17-1-2010 22:02, Mateusz Loskot wrote: Does that mean that I can use ogrinfo on a gzipped archive, like gdalinfo (http://trac.osgeo.org/gdal/wiki/UserDocs/ReadInZip)? Yes, it does but...as long as OGR driver performs filesystem operations using VSI*L API. The problem is that only fe

Re: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-16 Thread Jan Hartmann
Thanks Frank, I'll try and let you know. BTW: Perhaps your idea of GML-output for OGR would solve the Metadata-problem for OGR. Everyone can add whatever data they want to the GML-file, as long as OGR knows what parts to retrieve to reconstruct the vector map. Jan On 16-Jan-10 21:53, Frank

Re: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-16 Thread Jan Hartmann
On 16-1-2010 16:03, Mateusz Loskot wrote: Jan Hartmann wrote: Yes, that is clear, thanks. I see that at the moment only raster files are supported. Would it make sense to do this for vector formats too? The VSI layer is available to all parts of GDAL and OGR. If you scan source

Re: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-16 Thread Jan Hartmann
Yes, that is clear, thanks. I see that at the moment only raster files are supported. Would it make sense to do this for vector formats too? I am thinking of a dump from a large PostGIS database I had to upgrade from a 32 to a 64 bits server. I didn't like the pgdump format, as I got in all sor

Re: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-16 Thread Jan Hartmann
Oh, that makes a difference indeed! When I read Frank's long answer, this was the only point I didn't like. How does this mechanism work? Jan On 16-1-2010 12:33, Even Rouault wrote: I think Frank meant With Even's work, it is *now* possible for many drivers to to transparently access compres

Re: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-16 Thread Jan Hartmann
On 16-1-2010 3:01, Frank Warmerdam wrote: Instead, if this is a goal of archiving, I'd suggest archiving the original data (in a possibly arcane format), and a copy in a more accessable format likely to still be usable decades later. That's what we are doing now. It's OK for practical purpose

Re: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-15 Thread Jan Hartmann
endor neutral; and it be possible to retrieve all information in the future. Peter Jan Hartmann wrote: I was thinking along the same lines, but more in the direction of OGR as an archival standard. I have been working with archives and stored maps all my life and am now busy with lots of di

Re: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-15 Thread Jan Hartmann
I was thinking along the same lines, but more in the direction of OGR as an archival standard. I have been working with archives and stored maps all my life and am now busy with lots of digital historical maps. For raster maps the best storing format is Geotiff, from which all other formats ca

Re: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-14 Thread Jan Hartmann
On 13-1-2010 21:19, Mateusz Loskot wrote: IMHO, it's misunderstanding to consider OGR fully featured data model and I/O engine to read, write, process and analyse spatial vector data, especially if performance is a critical factor. IMHO, there are too many compromises in OGR. OK, that is

Re: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-13 Thread Jan Hartmann
On 13-1-2010 15:49, Ari Jolma wrote: Jan Hartmann wrote: Just curious, would it make sense / be possible to implement indexing in OGR, something like a generalized version of Mapserver's shptree, the "quadtree-based spatial index for a shapefiles"? http://mapserve

Re: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-13 Thread Jan Hartmann
On 13-1-2010 2:33, Mateusz Loskot wrote: OGR does not provide any spatial indexing layer common to various vector datasets. For many simple formats it performs the brute-force selection. Just curious, would it make sense / be possible to implement indexing in OGR, something like a genera

Re: [gdal-dev] (article) social organization of the open source R project

2010-01-06 Thread Jan Hartmann
I would everyone advise to read Frederick P. Brooks: The Mythical Man-Month. Essays on Software Engineering Brooks was manager of IBM's OS/360, the classical mammoth operating system of the nineteen-sixties. He wrote this book in 1975 (it is still in print) about what can go wrong in designing

Re: [gdal-dev] Calibrating colors of a series of maps

2009-12-26 Thread Jan Hartmann
ing here. Jan On 26-12-2009 17:51, Frank Warmerdam wrote: Jan Hartmann wrote: Hi, I put the following question on the GRASS user list. It isn't possible with GRASS (or the programs mentioned by Glynn), but if someone on this list can help me further, I would be very grateful. I have been

[gdal-dev] Calibrating colors of a series of maps

2009-12-26 Thread Jan Hartmann
s lots of scanned old maps in the near future. Jan Jan Hartmann wrote: > I am working with series of old topographical maps, ranging from 10 to > 600 map sheets per map. Each sheet is discolored a bit from age. How can > I process a complete series in such a way that there are no dif

Re: [gdal-dev] Remove Old Generation Python bindings for 1.7 release

2009-12-15 Thread Jan Hartmann
I'm not a developer, but I have had the same problems in my builds (I must compile everything on a Linux-cluster as non-root in non-standard locations, with up-to-date versions of Python, that kind of thing). Removing the old bindings would remove one further source of headaches, for me at le

Re: [gdal-dev] OGR DXF Driver, /me blogs

2009-12-05 Thread Jan Hartmann
"If only one in ten cities in the world made serious use of free gis software and provided enough financial support for one core developer it would have a huge boosting effect". Frank, do you realise in your Canadian basement how many cities there are in the world? What are you trying to b

[gdal-dev] Warping an already georeferenced image with control points

2009-02-18 Thread Jan Hartmann
of conflicting georeferencing information. And to me it doesn't matter anyway, as that "GCP-ed" image is just temporary. After the gdalwarp step I get an image with more traditional/well-supported georeferencing, I can delete it. Should gdalinfo (and possibly other programs) use

Re: [gdal-dev] Re: [Fwd: Re: Warping an already georeferenced image with control points

2009-02-15 Thread Jan Hartmann
to me it doesn't matter anyway, as that "GCP-ed" image is just temporary. After the gdalwarp step I get an image with more traditional/well-supported georeferencing, I can delete it. Should gdalinfo (and possibly other programs) use the GCPs to calculate the image extents? Mayb

Re: [gdal-dev] Re: [Fwd: Re: Warping an already georeferenced image with control points

2009-02-13 Thread Jan Hartmann
that adds the capability of adding the 4 corners as 4 additional GCPs with the "-add_corners_as_gcps" option. You might give it at try (it applies cleanly on GDAL 1.6.0 and trunk). I'm not sure yet if it is valuable enough to include it in baseline gdal_translate.cpp. Best r

Re: [Fwd: Re: [gdal-dev] Warping an already georeferenced image with control points

2009-02-11 Thread Jan Hartmann
e (or not), perhaps with adequate funding. Everything is better than talking to a blind wall. Sorry again Frank, Jan Jan Hartmann wrote: No, this doesn't seem to work. Although the coordinate system *is* retained when I add -a_srs to the gcps, the boundaries of the input file are discarded: th

Re: [gdal-dev] Re: [Fwd: Re: Warping an already georeferenced image with control points

2009-02-09 Thread Jan Hartmann
slate. If gdal_translate could keep this information in the target file, even when gcps have been added, everything would be fine (and my impression is it should). Jan Rahkonen Jukka wrote: Jan Hartmann uva.nl> writes: Jukka Rahkonen wrote: >> Hi, >> >> I am not su

Re: [gdal-dev] Re: [Fwd: Re: Warping an already georeferenced image with control points

2009-02-09 Thread Jan Hartmann
Jukka Rahkonen wrote: Jan Hartmann uva.nl> writes: No, this doesn't seem to work. Although the coordinate system *is* retained when I add -a_srs to the gcps, the boundaries of the input file are discarded: they are reset to the pixel dimensions of the file. Is this a bug? Ja

[Fwd: Re: [gdal-dev] Warping an already georeferenced image with control points

2009-02-07 Thread Jan Hartmann
No, this doesn't seem to work. Although the coordinate system *is* retained when I add -a_srs to the gcps, the boundaries of the input file are discarded: they are reset to the pixel dimensions of the file. Is this a bug? Jan Frank Warmerdam wrote: Jan Hartmann wrote: Hi, I am tryi

  1   2   >