#1 Tue 06 December 2022 18:47
- Thais781
- Participant occasionnel
- Date d'inscription: 24 Oct 2022
- Messages: 25
Convertion geometry dans le bon repère
Bonjour,
Je galère un peut dans la transformation de mes données géométriques
J'ai un point sous le format suivant : POINT (-1.41 43.69)
Hors pour qu'il soit lisible (dans le bon repère) dans QGIS, il faudrait qu'il ait un format plutôt de ce type : POINT (343125.4 6296416.2)
J'ai essayé de jouer avec le CRS :
Code:
mon_df = pd.read_csv('mon_file.csv', encoding ='utf-8', dtype=str) mon_df['geometry'] = mon_df['geometry'].map(shapely.wkt.loads) mon_gdf = geopandas.GeoDataFrame(mon_df, geometry='geometry', crs="EPSG:2154").to_crs(2154)
Mais rien, je n'arrive pas à avoir le format désiré.
Si vous avez un idée, sur comment en python faire cette transformation, je suis preneuse ....
Thais
Hors ligne
#2 Tue 06 December 2022 19:47
- Thais781
- Participant occasionnel
- Date d'inscription: 24 Oct 2022
- Messages: 25
Re: Convertion geometry dans le bon repère
La nuit portant conseil, j'ai trouvé la solution.
Lors du chargement du fichier CSV et de la création du fichier geopandas associé, il faut spécifier le CRS utilisé dans le fichier que l'on charge.
Code:
mon_gdf = geopandas.GeoDataFrame(mon_df, crs="EPSG:4326", geometry=mon_geometry)
Ici le CRS du GPS
et ensuite on converti au format CRS souhaité (ici celui de la france)
Code:
mon_gdf = mon_gdf.to_crs(2154)
J'espère que cela sera utile à d'autres qui se poseraient la même question que moi ;-)
Thais
Dernière modification par Thais781 (Wed 07 December 2022 08:38)
Hors ligne