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

Isa García Barón x@n@mb @ending from gm@il@com
Vie Sep 28 16:00:01 CEST 2018


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 using gmail.com
PhD Student at AZTI Foundation -  AZTI Fundazioa
Marine Ecosystems Functioning
Herrera Kaia, Portualdea z/g 20110 - Pasaia (Gipuzkoa)

	[[alternative HTML version deleted]]



Más información sobre la lista de distribución R-help-es