Tengo un problema que muchos usuarios de ArcGIS, QGIS y R de la República Checa y Eslovaquia pueden encontrar. En ArcGIS, estoy trabajando EPSG: 5514, Sistema de coordenadas proyectadas, definido como:
S-JTSK_Krovak_East_North
WKID: 5514 Authority: EPSG
Projection: Krovak
False_Easting: 0.0
False_Northing: 0.0
Pseudo_Standard_Parallel_1: 78.5
Scale_Factor: 0.9999
Azimuth: 30.28813975277778
Longitude_Of_Center: 24.83333333333333
Latitude_Of_Center: 49.5
X_Scale: -1.0
Y_Scale: 1.0
XY_Plane_Rotation: 90.0
Linear Unit: Meter (1.0)
Geographic Coordinate System: GCS_S_JTSK
Angular Unit: Degree (0.0174532925199433)
Prime Meridian: Greenwich (0.0)
Datum: D_S_JTSK
Spheroid: Bessel_1841
Semimajor Axis: 6377397.155
Semiminor Axis: 6356078.962818189
Inverse Flattening: 299.1528128Cuando cargo mi archivo en R, mi proj4string(data_in_5514) espectáculos:
"+proj=krovak +lat_0=49.5 +lon_0=24.83333333333333 +k=0.9999 +x_0=0 +y_0=0 +ellps=bessel +units=m +no_defs"Tengo la misma información cuando exporto mis archivos desde R usando writeRaster() función de raster paquete. Sin embargo, cuando quiero abrir mis datos modificados en R de nuevo en ArcGIS, mi ArcGIS se bloquea, y si tengo suerte, tengo un mensaje de advertencia :
Advertencia, ¡extensión incoherente!
que significa ( http://support.esri.com/em/knowledgebase/techarticles/detail/37732 ):
Este mensaje de error se devuelve si los datos tienen una definición de proyección incorrecta.
En concreto, los datos están en un sistema de coordenadas proyectadas con unidades de pies, metros u otras unidades lineales, pero la proyección se ha definido con un Sistema de Coordenadas Geográficas (SCG) que tiene unidades en grados decimales.
Por lo tanto, es evidente que la proyección de mis datos ha cambiado, pero ¿cuándo y cómo?
¿Cómo puedo transferir correctamente los datos entre R y ArcGIS utilizando Krovak EPSG 5514? También he encontrado, que este PCS también aparece como EPSG: 102067 bajo ArcGIS, sin embargo, R no reconoce este código.
EDITAR :
datos de la muestra: http://ulozto.sk/xNHrBv2D/r3-zip
como salida estaba utilizando varios formatos. Pude opet Gtiff una vez, .img neidther una vez, .grd parece no ser reconocido directamente por ArcGIS.
datos de exportación:
writeRaster(aaa, filename="r3.img", format="HFA", overwrite=TRUE)
# export data
writeRaster(aaa, filename="r3.tif", format="GTiff", overwrite=TRUE)
# export data
writeRaster(aaa, filename="r3", overwrite=TRUE) # creates a .grd
