Vous utilisez le système de coordonnées personnalisé dans ArcGIS Desktop?

16

Je ne connais pas tellement les systèmes de coordonnées ...

Dans mon bureau, nous utilisons pour traiter des données spatiales provenant de sites archéologiques. Chaque site possède son propre système de coordonnées xyz (GCS). Trois axes cartésiens simples ortogonaux . Au cours des dernières années, nous avons géré ces données spatiales via un logiciel SIG (ArcGIS), sans utiliser de système de coordonnées spécifique (laissez-le simplement "indéfini")

Je voudrais savoir s'il existe un GCS conçu pour traiter de tels ensembles de données en utilisant un axe orthogonal cartésien simple, sans distorsions de grille du GCS typique. De plus, j'aimerais savoir si ce système est adapté à son utilisation dans une application de cartographie en ligne.

Soit dit en passant, nous gérons les environnements 2D (ArcMap) et 3D (ArcScene) et travaillons avec "mm" comme unité de base de longueur.

Si une telle chose n'existe pas, peut-être que quelqu'un sait comment la créer.

Maik
la source

Réponses:

17

En supposant que vos zones d'intérêt sont plutôt petites par rapport au globe, vous pouvez configurer une projection Mercator transversale personnalisée.

Vous devez connaître les coordonnées géographiques lat_0 et lon_0 de l'origine de votre CRS, et la direction des axes x et y:

+proj=tmerc +lat_0=51.4 +lon_0=7 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

Pour une utilisation avec Arcgis, le fichier .prj doit contenir:

PROJCS["Transverse_Mercator",GEOGCS["GCS_WGS_1984",DATUM["D_unknown",SPHEROID["WGS84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",51.4],PARAMETER["central_meridian",7],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]

S'ils ne sont pas parallèles aux directions Nord et Est, une rotation serait nécessaire en utilisant une projection de mercator oblique:

+proj=omerc +lat_0=51.4 +lonc=7 +alpha=-10 +k=1 +x_0=0 +y_0=0 +gamma=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

Le fichier .prj pour cette projection est:

PROJCS["Hotine_Oblique_Mercator_Azimuth_Center",GEOGCS["GCS_WGS_1984",DATUM["D_unknown",SPHEROID["WGS84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Hotine_Oblique_Mercator_Azimuth_Center"],PARAMETER["latitude_of_center",51.4],PARAMETER["longitude_of_center",7],PARAMETER["azimuth",-10],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]

Voici les deux ensemble (tmerc en bleu et omerc en rouge):

entrez la description de l'image ici

Les lat_0 et lon_0 (ou lonc) sont les coordonnées de votre origine CRS locale. Vous devez les explorer avec un GPS ou d'autres cartes. L'angle que vous devez calculer à partir des coordonnées de deux points connus sur votre axe X local.

La seule chose qui ne fonctionne pas comme prévu est le millimètre en unités. Cela devrait fonctionner, mais ne fonctionne pas dans QGIS actuel :-(

http://hub.qgis.org/issues/9414


Mise à jour

Les versions actuelles de QGIS peuvent désormais travailler avec des projections locales en utilisant mm comme unités.

AndreJ
la source