[R-es] Problema utilizando la función dist2isobath (marmap package) con un sistema de coordenadas proyectadas

[Isa García Barón]

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!