Nous utilisons des cookies pour vous garantir la meilleure expérience sur notre site. Si vous continuez à utiliser ce dernier, nous considèrerons que vous acceptez l'utilisation des cookies. J'ai compris ! ou En savoir plus !.
banniere

Le portail francophone de la géomatique


Toujours pas inscrit ? Mot de passe oublié ?
Nom d'utilisateur    Mot de passe              Toujours pas inscrit ?   Mot de passe oublié ?

Annonce

Suite à un problème technique intervenu entre le 22 et le 23 mars, nous avons du procéder dans la soirée du 25 mars, à la restauration de la base de données du 24 mars (matinée).

En clair, nous avons perdu vos contributions et inscriptions du dimanche 24 et du lundi 25 mars.
Nous vous prions de nous excuser.

#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

 

Pied de page des forums

Powered by FluxBB