Hi Jukka,

jratike80 wrote
> You can have only one SRID on a layer so you must either split your data
> by
> the epsg field or transform each record to same SRID. Otherwise your
> request
> was almost complete.

I know, but my problem was the source!


jratike80 wrote
> ogrinfo -dialect sqlite -sql "select ST_GeomFromText('POINT('|| x || ' '
> ||
> y || ')',CAST(epsg AS integer)) from coords" coords.csv

And thank you to your example I have found the solution that was right for
me:

ogrinfo -dialect sqlite -sql "select ST_TRANSFORM(ST_GeomFromText('POINT('||
x || ' ' || y || ')',CAST(epsg AS integer)),4326) from coords" coords.csv 

It gives me 

OGRFeature(SELECT):0
  POINT (15.0 36.1447180988163)

OGRFeature(SELECT):1
  POINT (9.0 36.1447180988163)

Than I have unified all, without splitting my source file (as you have
written I have trasformed each record to same SRID).

You and ogr are great :)




--
Sent from: http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html
_______________________________________________
gdal-dev mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Reply via email to