Valores de NA en mi trama después de cambiar la resolución, la extensión y el origen

Question

Estoy trabajando con dos rásteres que difieren en su origen, la extensión y la resolución. Tengo una batimetría de trama, con una resolución muy alta (x=0.0008333333, y=0.0008333333) y una gran extensión espacial. También tengo un mar de temperatura de la superficie de ráster, que tiene mucho más gruesa de la resolución (x=0.04166667, y=0.04166667). Tanto los rásteres tienen la misma proyección (longlat, datum=WGS84). Tanto los rásteres también tiene 'NA' valores que corresponden a masas de tierra.

Los datos de la batimetría:

https://www.ngdc.noaa.gov/dem/squareCellGrid/download/4956

Sst de datos: 25 de junio-2 de julio (2016) de 8 días de composite, recortada a -131.8333, -127.6667, 50.125, 54.54166 (xmin, xmax, ymin, ymax)

https://oceancolor.gsfc.nasa.gov/cgi/l3?per=8D&prd=NSST_sst.nc&sen=A&res=9km&num=24&ctg=Standard&date=25Jun2016

Me han manipulado la batimetría de trama (con R) para que coincida con la medida, el origen, y la resolución de la temperatura superficial del mar de trama.

bathy2<-projectRaster(bathymetry, sst, method="bilinear")

Hice un marco de datos para que yo pudiera comparar la temperatura superficial del mar de los valores y de la batimetría de los valores en cada píxel. Entonces me di cuenta de que tengo un problema. Algunos de los píxeles que tiene 'NA' los valores de la batimetría, pero los valores numéricos para la SST. ¿Cómo puede ser esto, dado que las masas de tierra debe estar exactamente en el mismo lugar? El original de la batimetría de trama tiene una resolución mucho mayor que la SST de trama, por lo que su esquema de las masas de tierra debe ser precisa.

Es este un resultado de la utilización de projectRaster(), o un problema con los datos de la misma? De todos modos puedo solucionarlo?

Answer 1

La batimetría de trama que se nos da como ejemplo no es en todo el mundo, tal SST de la trama. Por lo tanto, si intenta proyecto de la trama de SST con projectRaster() obtendrá la mayoría de NA valores en remuestreada de datos, debido a que la función de remuestreo a la nueva resolución/medida/CRS. Primero, los cultivos de SST raster y después de eso, el proyecto de la batimetría de raster:

library(raster)

bathy <- raster("~/Downloads/British_Columbia_DEM_5076/british_columbia_3sec.asc")

crs(bathy) <- "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs "

SST <- raster("~/Downloads/A20163372016344.L3m_8D_NSST_sst_9km.nc")

SST2 <- crop(SST, bathy)

plot(SST2, legend=F)

bathy2 <- projectRaster(bathy, SST2, method = "bilinear")

plot(bathy2, col=grey(1:100/100))
plot(SST2, add=T, alpha=0.6,legend=F)

---

NA de datos

Desde la descripción del producto:

El digital de costa utilizados en el desarrollo de la Columbia Británica DEM fue generado por la fusión de vector de las costas de la NOAA Electrónico Cartas de navegación (ENCs) y de Recursos Naturales de Canadá (NRCAN), a continuación, editado basado en ESRI mundial de la capa de imágenes. El final digital de la costa fue convertida a formato xyz, con elevación de cero y punto espaciado a 10 metros de distancia. El digital de la costa también fue convertido a un polígono y en última instancia, una trama para el enmascaramiento de la topografía y la eliminación de datos interpolados.

df <- as.data.frame(stack(bathy2, SST2))

dim(df[complete.cases(df),])[1]
[1] 8070

length(df$british_columbia_3sec[!is.na(df$british_columbia_3sec)])
[1] 8710

length(df$Sea.Surface.Temperature[!is.na(df$Sea.Surface.Temperature)])
[1] 8264

test <- overlay(is.na(bathy2),is.na(SST2),fun=sum)

test[test!=2] <- NA ;test[test==2] <- 1

plot(test, col="yellow",legend=F)
plot(bathy2, col=scales::alpha("red",0.5),legend=F,add=T)
plot(SST2, col=scales::alpha("blue",0.5),legend=F,add=T)

Así que usted puede esperar de valores válidos en las costas o la parte interior del terreno. Los NA valores proviene de la falta de coincidencia entre los valores válidos de la SST en la costa. Si quieres completar los valores de SST, necesita combinar la batimetría con una superficie de DEM.