У меня есть файл .tif из ArcGIS, представляющий собой карту местоположений с сеткой и значениями глубины z для средней точки каждой сетки. Я пытаюсь импортировать этот растровый файл в R, используя библиотеку «растровых данных», и вычислять наклон в каждом месте, используя библиотеку «SDMTools». Однако у меня возникли проблемы с проекцией, так как она не в широте / долготе, и я не могу понять, как получить доступ к фактическим местоположениям с соответствующими глубинами и уклонами.
Пока что я сделал следующее: '' '
library(raster)
grid <- raster(here::here("data/arcGIS/GridDepthCPUE1.tif"))
plot(grid) #map looks right, but x/y axis not in lat/lon
grid
class : RasterLayer
dimensions : 272, 252, 68544 (nrow, ncol, ncell)
resolution : 2800, 2800 (x, y)
extent : 753206.3, 1458806, 535587.9, 1297188 (xmin, xmax, ymin, ymax)
crs : +proj=aea +lat_1=55 +lat_2=65 +lat_0=50 +lon_0=-154 +x_0=0 +y_0=0 +datum=NAD83
+units=m +no_defs +ellps=GRS80 +towgs84=0,0,0
source : /Users/laurenwild/Desktop/UAF/Thesis/R-Markdown-GitHub/UAF_PhD_Ch3_Pm_Movement/data/ArcGIS/GridDepthCPUE1.tif
names : GridDepthCPUE1
values : 40, 3426 (min, max)
attributes :
ID Rowid COUNT
from: 40 0 7
to : 3426 2536 1
' ''
Тогда, если я пойду дальше и запущу функцию slope () следующим образом:'' '
library(SDMTools)
grid.slope <- slope(grid)
plot(grid.slope)
class : RasterLayer
dimensions : 272, 252, 68544 (nrow, ncol, ncell)
resolution : 2800, 2800 (x, y)
extent : 753206.3, 1458806, 535587.9, 1297188 (xmin, xmax, ymin, ymax)
crs : +proj=aea +lat_1=55 +lat_2=65 +lat_0=50 +lon_0=-154 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0
source : memory
names : GridDepthCPUE1
values : 0, 0.3361841 (min, max)
attributes :
ID Rowid COUNT
from: 40 0 7
to : 3426 2536 1
Как видите, значения "экстента" имеют минимальное / максимальное значения: 753206,3, 1458806, 535587,9, 1297188 (xmin, xmax, ymin, ymax). Я хотел бы, чтобы они были в лат / лоне, я думаю! Я хочу, чтобы crs = 4326.
Кроме того, моя конечная цель - извлечь местоположения с их значениями глубины и уклона для использования в линейных моделях.
Мне сложно понять, как получить доступ к необработанным данным из объекта grid.slope с помощью операторов @ и $.
Есть ли способ сделать это? Любая помощь будет потрясающей!