[gdal-dev] Gdal_grid with a mask layer?

2016-10-06 Thread Rahkonen Jukka (MML)
Hi, I have depth data from a lake as scattered points and I would like to convert them into raster DEM. Unfortunately this lake is not a rectangular, north oriented one but it has a curvy shoreline and lots of islands. In my first trials the resulting DEM leaks away from the water to the dry g

[gdal-dev] GPKG: Drop table but how to drop the rest?

2016-09-29 Thread Rahkonen Jukka (MML)
Hi, I wonder if it could be possible to do something like "DROP CASCADING" in GeoPackage and delete the spatial index by the same. It will for sure be unusable once the table is dropped. I feel that SQLite does not support that because cascading drop seems to require that parent and child table

Re: [gdal-dev] Preserve Source Raster Resolution When Merging

2016-08-25 Thread Rahkonen Jukka (MML)
Sam Wilson wrote: > Hi, > I'm attempting to merge (or maybe the proper term is mosaic) overlapping > raster datasets of varying resolutions. I'd like to have the higher > resolution datasets override the lower resolution datasets where available. > I've been able to get close using various me

Re: [gdal-dev] Issues with MapInfo Conversion to GeoJson using ogr2ogr

2016-08-16 Thread Rahkonen Jukka (MML)
Suman Humagain wrote: > Hi folks, > I am converting the Mapinfo files to GeoJson format using the ogr2ogr tool. > The things are working great and its converting fine but i am getting > different polylines than i have drawn in the mapinfo. Some extra polylines > are created in the intersections.

Re: [gdal-dev] Motions: Promote 1.11.5 RC 1, 2.0.3 RC2 and and 2.1.1 RC2 for release

2016-07-12 Thread Rahkonen Jukka (MML)
+1 for all three. -Jukka Rahkonen- Lähettäjä: Even Rouault<mailto:even.roua...@spatialys.com> Lähetetty: ‎12.‎7.‎2016 12:15 Vastaanottaja: Rahkonen Jukka (MML)<mailto:jukka.rahko...@maanmittauslaitos.fi> Aihe: Fwd: [gdal-dev] Motions: Promote 1.11.5

Re: [gdal-dev] GDAL 2.0.2 & gpkg woes

2016-05-02 Thread Rahkonen Jukka (MML)
Alan Stewart wrote: >Unfortunately, makes no difference, with my application code or ogr2ogr. I >also tested the new libspatialite-4.2.0a DLL >against my old GDAL build. Works >fine. Hi, Perhaps you should try to GPKG driver with GDAL that is built without Spatialite. As far as I know GPKG d

Re: [gdal-dev] How to add more than one column to a shapefile using ogrinfo

2016-05-02 Thread Rahkonen Jukka (MML)
Andrea wrote: >Hi Jukka, >Rahkonen Jukka (MML) wrote >> Do you know some database that accepts SQL "ALTER TABLE foo ADD COLUMN >> bar1, bar2"? At least SQLite and PostgreSQL don't. That kind of usage is >> not mentioned in a bunch of docum

Re: [gdal-dev] How to add more than one column to a shapefile using ogrinfo

2016-05-01 Thread Rahkonen Jukka (MML)
Andrea Borruso wrote: >thank you very much. I had read that page, but stupidly I did not think >that "ALTER TABLE tablename ADD [COLUMN] columnname columntype" meant >strictly that I could only use it in this way. >My wrong intepretion was that I could use "alter table add" SQL command >and its v

Re: [gdal-dev] How to append to all layers with ogr2ogr?

2016-02-01 Thread Rahkonen Jukka (MML)
More intormation: Reason may be that GML files do not have same layers but the second one has one layer more: ogrinfo /vsizip/K2344R_mtk.zip ERROR 1: Random access not supported for /vsizip Had to open data source read-only. INFO: Open of `/vsizip/K2344R_mtk.zip' using driver `GML' success

Re: [gdal-dev] How to append to all layers with ogr2ogr?

2016-02-01 Thread Rahkonen Jukka (MML)
Hi, Did you try to add another GML file into existing GeoPackege? First run OK: ogr2ogr -f gpkg -update -append mtk.gpkg /vsizip/K2344R_mtk.zip Secong run gives an error: ogr2ogr -f gpkg -update -append mtk.gpkg /vsizip/K2433L_mtk.zip ERROR 1: Layer Kunta already exists, and -append not specifie

Re: [gdal-dev] GeoPDF content inverted as GeoTIF

2015-12-14 Thread Rahkonen Jukka (MML)
Hi, There is something odd in the PDF. It is without coordinate system and perhaps because of that the top-left coordinates are (0,0) and north coordinate is increasing from top to bottom. Probably the "Invalid length for GPTS object" error does not have any effect on this. gdalinfo loopmap.

Re: [gdal-dev] Motion: Commit access to David Adler

2015-12-03 Thread Rahkonen Jukka (MML)
+1 -Jukka Rahkonen- Mateusz Loskot wrote: Motion: Committer Status for David Adler --- ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] Question: Creating overviews

2015-12-01 Thread Rahkonen Jukka (MML)
Hi, That is normal. Overviews are kind of support files for the main image and geotransform can be stolen from there. With internal overviews in TIFF that's also the only possibility. If you need standalone downsampled images create them with gdal_translate or gdalwarp. -Jukka Rahkonen- Gane

Re: [gdal-dev] strange grid pattern when reprojecting

2015-11-30 Thread Rahkonen Jukka (MML)
Hi, Sorry for being unclear, I meant a sample of source data that could be used for reproducing the issue, not a sample of output with a grid artifact. -Jukka Rahkonen- > Hi, > I've tried adding -et 0 and the grid feature is still there. Here is an > example: https://www.dropbox.com/s/teywuev

Re: [gdal-dev] strange grid pattern when reprojecting

2015-11-27 Thread Rahkonen Jukka (MML)
Hi, Gdalwarp is the tool to use and usually it works well but your case is not the most common one "lunar south polar DEMs in stereographic projection to gnomonic projection." Best if you could put a sample image available somewhere but perhaps but there are some documents about advanced trick

Re: [gdal-dev] Proposed Block iterator API

2015-11-25 Thread Rahkonen Jukka (MML)
Hi, I was wondering if this idea could resolve also this old ticket https://trac.osgeo.org/gdal/ticket/1233 -Jukka Rahkonen- Even Rouault wrote: Hi, I was wondering if there wasn't a need for a block iterator API for code that wants to process raster efficiently and iterates over blocks f

Re: [gdal-dev] Motion: Adopt RFC 60: Improved round-tripping in OGR

2015-11-05 Thread Rahkonen Jukka (MML)
+1 -Jukka Rahkonen- >Hi, >Since no remarks have been done on the latest proposal, I move to adopt RFC >60 : Improved round-tripping in OGR >https://trac.osgeo.org/gdal/wiki/rfc60_improved_roundtripping_in_ogr >Starting with my +1, >Even __

Re: [gdal-dev] Bad design in -clipsrs?

2015-10-25 Thread Rahkonen Jukka (MML)
bbox of a feature - no cross - for sure no need to clip. -Jukka- Lähettäjä: Even Rouault<mailto:even.roua...@spatialys.com> Lähetetty: ‎25.‎10.‎2015 12:02 Vastaanottaja: gdal-dev@lists.osgeo.org<mailto:gdal-dev@lists.osgeo.org> Kopio: Rahkone

Re: [gdal-dev] Gdaltransform hangs sometimes after being idle

2015-09-05 Thread Rahkonen Jukka (MML)
Hi, I can reproduce now. GDAL is 64-bit 2.1-dev from http://download.gisinternals.com/sdk/downloads/release-1600-x64-gdal-mapserver.zip 1. gdaltransform -s_srs epsg:4326 -t_srs epsg:3067 30 65 641428.429705435 7211811.30748109 0 --wait 5 minutes 31 70 652590.036387019 7770880.82231984 0 --now Ct

Re: [gdal-dev] SRS of the Elastic Search geometries?

2015-09-01 Thread Rahkonen Jukka (MML)
Even Rouault wrote: >> >> Probably that's it, the explicit SRS. I started with data that have native >> SRS epsg:3067 and OpenJUMP JML format which has no means for holding SRS. >> Conversion into ES gives always an empty layer, just the schema gets >> inserted but no geometries, even is I use -s_

Re: [gdal-dev] Question

2015-08-26 Thread Rahkonen Jukka (MML)
Hi, I made my test with GDAL 2.0 binaries from gisinternals by placing your osmconf.ini into \bin\gdal-data. Any change you do for the selected tags in osmconf.ini will have on effect on " ogrinfo hamburg-latest.osm.pbf points -so ". Don't do anything before you see that your osmconf.ini is use

Re: [gdal-dev] Gdal_translate considers --config means -co

2015-06-18 Thread Rahkonen Jukka (MML)
Even Rouault wrote > Le jeudi 18 juin 2015 16:35:37, Jukka Rahkonen a écrit : >> Hi, >> >> Spotted this from gis.stackexchange >> http://gis.stackexchange.com/questions/151153/is-gdal-translate-ignori >> ng-th >> e-cpl-tmpdir-switch >> >> I repeated the test on Windows with GDAL 2.0 and it seems

Re: [gdal-dev] Relating to GeoLocation tool of Geospatial PDF Document

2015-04-08 Thread Rahkonen Jukka (MML)
Hi, The OGC_BP option is actually created by TerraGo and I am surprised that it does not work for you with the TerraGo Toolbar because that is the native way for adding georeferencing info into GeoPDF (tm) files. In some of my early tests I found problems with interoperability: Only ISO: OK f

Re: [gdal-dev] New error messages when writing Spatialite

2015-03-24 Thread Rahkonen Jukka (MML)
Even Rouault wrote: > Le mardi 24 mars 2015 14:14:25, Jukka Rahkonen a écrit : >> Hi, >> >> I noticed that ogr2ogr from GDAL 2.0-dev is nowadays printing error >> messages which I have not seen before: >> >> "CreateSpatialIndex() error: either "points"."GEOMETRY" isn't a >> Geometry column or

Re: [gdal-dev] Shapefile driver and support for Boolean attributes

2015-03-18 Thread Rahkonen Jukka (MML)
Hi, I am testing new Boolean data type in OpenJUMP which is saved as Logical when data are exported into shapefiles. Kosmo GIS supports Logical already and Geotools and gvSIG as well. I said Logical is well defined but that is not totally true. Some implementations seem to use space to mean not

Re: [gdal-dev] Ogrtindex: schema consistency check too strict?

2015-03-17 Thread Rahkonen Jukka (MML)
- Lähettäjä: Even Rouault [mailto:even.roua...@spatialys.com] Lähetetty: 17. maaliskuuta 2015 17:52 Vastaanottaja: Rahkonen Jukka (MML) Kopio: gdal-dev@lists.osgeo.org Aihe: Re: [gdal-dev] Ogrtindex: schema consistency check too strict? Jukka, Accepting both Integer and Integer64 as the same could

Re: [gdal-dev] Affine parameters in the gdal datamodel doc

2015-01-09 Thread Rahkonen Jukka (MML)
Even Rouault wrote: >> Hi, >> >> I am trying to understand the geotransform parameters by reading >> http://www.gdal.org/gdal_datamodel.html but I am confused. In the >> Affine GeoTransform section there are formulas >> >> Xgeo = GT(0) + Xpixel*GT(1) + Yline*GT(2) >> Ygeo = GT(3) + Xpi

<    1   2   3   4