proj4string ne semble pas définir complètement les informations de projection ESRI

Question

J'ai un fichier de formes que j'ai importé dans R en utilisant readOGR du paquetage rgdal. J'y travaille un peu, comme ajouter des informations d'attribut, etc., puis je l'exporte à nouveau sous forme de fichier de formes ESRI, avec un nouveau nom. Cependant, lorsque j'apporte à la fois l'original et le nouveau fichier de formes dans ArcGIS, il me dit que le CRS ne correspond pas.

Ainsi, en notant que tous les paramètres de projection restent les mêmes, mais les projections et coordonner les noms de système sont différents, et la donnée nom est tombé, mes questionas sont:

1. est le deuxième CRS la même que le premier?

2. Si oui, pourquoi les noms ont-ils changé et pourquoi ArcGIS ne les reconnaît plus comme identiques?

3. Sinon, comment cela a-t-il été changé?

4. La proj4string peut-elle être modifiée pour être plus spécifique, et si oui, pourquoi readOGR ne l'a-t-elle pas déjà fait pour conserver toutes les informations?

je peux utiliser les nouveaux fichiers de formes très bien, mais il serait bon de savoir que le CRS est identique à l'original. Et, bien sûr, je pourrais le définir de nouveau dans ArcGIS, mais une partie de la motivation pour travailler dans R est d'éviter le besoin de pointer et de cliquer pour de nombreux fichiers.

J'apprécie toute idée ou éclaircissement.

Voici les informations de projection d'origine à partir ArcGIS:

Projected Coordinate System: NAD_1983_HARN_Transverse_Mercator 
Projection: Transverse_Mercator 
False_Easting: 520000.00000000 
False_Northing: -4480000.00000000 
Central_Meridian: -90.00000000 
Scale_Factor: 0.99960000 
Latitude_Of_Origin: 0.00000000 
Linear Unit: Meter 
Geographic Coordinate System: GCS_North_American_1983_HARN 
Datum: D_North_American_1983_HARN 
Prime Meridian:  Greenwich 
Angular Unit: Degree 

Voici le proj4string de R, qui est également d'accord avec le proj4string donné pour cette projection à www.spatialreference.org pour epsg: 3071 et aussi pour SR-ORG: 7396.

+proj=tmerc +lat_0=0 +lon_0=-90 +k=0.9996 +x_0=520000 +y_0=-4480000 +ellps=GRS80 +units=m +no_defs 

Lorsque j'utilise writeOGR pour exporter le SpatialPolygonsDataFrame avec le proj4string ci-dessus, puis le ramener dans ArcGIS, les informations de projection est donnée comme suit, et ne sera plus reconnu comme l'original.

Projected Coordinate System: Transverse_Mercator 
Projection: Transverse_Mercator 
false_easting: 520000.00000000 
false_northing: -4480000.00000000 
central_meridian: -90.00000000 
scale_factor: 0.99960000 
latitude_of_origin: 0.00000000 
Linear Unit: Meter 
Geographic Coordinate System: GCS_GRS 1980(IUGG, 1980) 
Datum: D_unknown 
Prime Meridian:  Greenwich 
Angular Unit: Degree 

Answer 1

Peut-être pas une réponse définitive, mais je posté cette question sur la liste R-sig-Geo service, et a obtenu quelques solutions possibles de travail autour. Pour l'instant, j'ai simplement utilisé un script R pour écraser le fichier .prj avec une copie de l'original, et cela semble fonctionner correctement. Il a également été suggéré d'utiliser un paquet appelé arcgisbinding pour faire le pont entre ArcGIS et R (et peut-être une solution similaire serait-elle disponible pour QGIS?). Je n'ai pas vérifié la solution d'arcgisbinding, mais plus d'informations peuvent être trouvées à l'article de blog here et dans la documentation de paquet here.