J'ai une liste de géocroords (lat, long) et un shapefile avec différentes couches. Je veux être capable d'identifier à quelle couche appartient chaque coord.Conversion de la latitude du shapefile en vrai
Mais le shapefile (SHP) a poligons dans lequel la latitude et la longitude sont exprimées en nombre dans une plage bizarre, comme 120724,86008864 et 484497,34058312
Je sais que le fichier PRJ contient les informations sur la façon dont cette transformation a été fait, mais je ne semble pas obtenir comment. Voici:
PROJCS [ "RD_New", GEOGCS [ "GCS_Amersfoort", DONNÉE [ "D_Amersfoort", SPHEROID [ "Bessel_1841", 6377397.155,299.1528128]], PRIMEM [ "Greenwich", 0], UNITÉ [ » Degré ", 0.0174532925199432955]], PROJECTION [" Double_Stereographic "], PARAMETRE [" False_East ", 155000], PARAMETRE [" False_Northing ", 463000], PARAMETRE [" Central_Meridian ", 5.38763888888889], PARAMETRE [" Scale_Factor ", 0.9999079], PARAMETRE ["Latitude_Of_Origin", 52.15616055555555], UNIT ["Compteur", 1]]
La question concrète est de savoir comment transformer un point lat/long régulier en ceux du fichier de formes.
Travailler en Python, avec cette bibliothèque http://gdal.org/python/
Toute aide est appréciée.
Vous voulez obtenir les coordonnées? Ou analyser le fichier prj? –
Je peux analyser le fichier sans problème, mais les lat/long dans le fichier de formes sont dans une échelle ou une plage différente. Par exemple, une latitude de 120724.86008864 qui est manifestement erronée. J'ai des latitudes normales telles que 52.3605883. Donc je veux savoir quelle transformation je dois appliquer à mes coordonnées normales pour être comme celles du fichier. Ensuite, je serai en mesure de les croiser avec les poligones dans les couches. –