Can't transform coordinates, source layer has no

[mhugo]

When trying to import https://github.com/BRGM/gml_application_schema_toolbox/blob/master/samples/BRGM_raw_database_observation_waterml2_output.xml

Maybe it lacks a way to set the source SRS ...

[sgrellet]

you already have the SRS

49.9506859240222 2.3597603712661

[mhugo]

Exact. There is something else GDAL does not like apparently ...

[mhugo]

@rouault apparently if 'srcSRS' is not set, the conversion keeps giving me errors and no geometries are created in the sqlite file. When I set srcSRS to 4326, it works fine.

Example script to reproduce: https://github.com/BRGM/gml_application_schema_toolbox/blob/master/tests/test_load_in_qgis.py#L36

Does that mean the input SRS must be supplied ?

[rouault]

If you run with gdal.VectorTranlate() with reproject = True and a dstSRS, then the source layers with geometries must have a SRS set, or it must be explictly set with srcSRS

Here the issue is that GDAL doesn't recognize the GML srs as it is set on the pos element (which is valid, but unusual) and not on the Point

Looking if I can come up with a GDAL fix

[rouault]

Fixed in GDAL per https://trac.osgeo.org/gdal/ticket/6962

[mhugo]

Thanks !

[rouault]

Wondering if the combo in the plugin shouldn't have an item that would completely disable reprojection ("Source SRS") ?

[mhugo]

I've added a 'Source SRS" button in https://github.com/BRGM/gml_application_schema_toolbox/commit/70bbd0fe61da6d18731e8d2e22b6f7dac0c1b1f5