Management of longitude/latitude with GIS data
Not fully sure this is a real incident or if this is related to the way I manage the spatial reference in my projects... By the way, I usually define my projection with ImportFromEPSG. I didn't had any problems when I build mesh with UTM projection but I had to build a model in long/lat so I indicate ImportFromEPSG(4326) and then the problem occurred. I tried to use GIS file preliminary converted in 4326, or just indicate an output_SRS in 4326 but the mesh was not at the right position. There was always a switch between X and Y values... I noticed there is a particular mention in field.py in the init of class Raster where there was a specific mention of _tools.osr.OAMS_TRADITIONAL_GIS_ORDER.
After adding domain_srs.SetAxisMappingStrategy(seamsh._tools.osr.OAMS_TRADITIONAL_GIS_ORDER) in my mesh script everything went well. It's possible to add this each time but maybe this could be directly implemented in geometry.py
Regards