Как сохранить поле 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".

0 ответов

Другие вопросы по тегам