Как сохранить поле SHP для использования в качестве ковариаты в lppm?
У меня есть записи о мертвых животных вдоль дорог (полилинии SHP), и мне нужно прогнозировать эти записи по ряду ковариат, включая тип дороги (поле SHP). Мне нужно преобразовать шейп-файл в объект linnet, и я хотел бы использовать поле "тип дороги" как ковариату в lppm. Предусмотренные решения: 'linfun': я не могу понять, как адаптировать эту функцию к этому случаю; 'linim' (растр): это приведет к неприятностям на перекрестке с несколькими типами дорог.
Вот пример сети дорог: https://wetransfer.com/downloads/c581143a04569b9dadfaac0543d58f5020180522161621/ec5223
и код для чтения:
library(spatstat)
# to read the SHP including the linear network and the field "type":
library(maptools)
library(sp)
library(rgdal)
oldfolder2<-getwd()
setwd(system.file("shapes", package="maptools"))
L<-readShapeSpatial("Road-sample.shp") # maybe need to write the complete file path
setwd(oldfolder2)
plot(L)
Мне не удалось адаптировать функцию 'linfun'; я должен извлечь координаты X, Y? В этом случае параметры "x,y" будут такими же, как "seg, tp".