Isa García Barón
2018-Sep-28 14:00 UTC
[R-es] Problema utilizando la función dist2isobath (marmap package) con un sistema de coordenadas proyectadas
Tengo un grid creado en un sistema de coordenadas proyectado (WGS84 zona 30N) y estoy intentando extraer la distancia más cercana a la costa desde de todas las celdas del grid, asi como a otras isóbatas utilizando la función "dist2isobath" del paquete marmap. Para ello reproyecté los datos del ETOPO1 (batimetría oceánica) a la misma proyección que los datos de mi grid y después he intentado utilizar la función "dist2isobath", pero me da un error. Aparentemente sólo funciona con coordenadas geográficos (ver el error a continuación). ¿Sabeis cómo puedo resolverlo o si existe otra función con la que trabajar con coordenadas proyectadas? Aquí os paso un ejemplo de los datos y el script, asi como el error que obtengo: head(grid) # ejemplo de mis datos del grid Lon Lat1 -124195.7 49866522 -120195.7 49866523 -116195.7 49866524 -112195.7 49866525 -108195.7 49866526 -104195.7 4986652 summary(etopo1) # ETOPO1 ya reproyectado y convertido a la clase "bathy" # Bathymetric data of class 'bathy', with 1329 rows and 709 columns# Latitudinal range: 4306614.28 to 5618264.28 (4306614.28 N to 5618264.28 N)# Longitudinal range: -550748.25 to 1203531.75 (550748.25 W to 1203171.75 E)# Cell size: 79258.1 minute(s) # Depth statistics:# Min. 1st Qu. Median Mean 3rd Qu. Max. NA's# -5586.23 -4004.18 -99.56 -1392.81 244.61 2928.77 141413 # First 3 columns and rows of the bathymetric matrix: # 4306614.27700283 4308466.88999718 4310319.50299153 # -550748.253167697 NA NA NA # -549427.259191793 NA NA NA # -548106.265215889 NA NA NA ## Y aquí utilizao la función para obtener las distancias: DistCoast <- dist2isobath(etopo1, grid[,1:2], isobath=0) Error in .pointsToMatrix(p) : longitude < -360 Gracias por la ayuda! -- *-----------------------------------------------------------------------------------------------------* *Isabel García Barón* Email: xanamb en gmail.com PhD Student at AZTI Foundation - AZTI Fundazioa Marine Ecosystems Functioning Herrera Kaia, Portualdea z/g 20110 - Pasaia (Gipuzkoa) [[alternative HTML version deleted]]