Используя shapefiles vignette ofpretstat, мне удалось создать список объектов psp, чтобы каждая функция shapefile была элементом в этом списке psp с атрибутами в виде меток. Конкретно, я использовал раздел 3.2.4 в виньетка следующимКак преобразовать список объектов psp в linnet
> x<- readShapeSpatial("includes/roadshapefile.shp")
> out <- lapply([email protected], function(z) { lapply([email protected], as.psp) })
> dat <- [email protected]
> for(i in seq(nrow(dat)))
+ out[[i]] <- lapply(out[[i]], "marks<-", value=dat[i, , drop=FALSE])
> roadlist <- do.call("c", out)
Теперь я хотел бы создать целую сеть linnet из этого списка УЧС объектов. Есть ли прямая функция или мне нужно прокручивать каждый из объектов psp, чтобы сделать каждый из них linnet, а затем каким-то образом объединить его с полной сетью? as.linnet не будет принимать этот список psp в качестве допустимого ввода, потому что его класс является простым списком.
Я чувствую, что мне не хватает небольшого шага. Может ли кто-нибудь мне помочь.
Как только у меня есть linnet, я хочу создать шаблон точки с помощью linfun и получить доступ к этим атрибутам (меток), чтобы сделать неоднородную интенсивность точек.
Благодаря
Adrian Baddeley недавно отправил код Roger Bivand с кодом 'linnet', который поместил его в версию разработки« maptools », доступную в R-forge: http://r-forge.r-project.org/R/?group_id = 943, возможно, вы можете использовать это, чтобы легче импортировать linnet. Я не уверен, как вы можете сохранить атрибуты и поместить их в 'linfun' или' linim'. Мне нужно будет изучить это. –
Спасибо. Я прочитаю вашу ссылку. Не могли бы вы, возможно, рассказать мне, как я могу генерировать неоднородную интенсивность в сети, используя spartstat? Я хочу принести по крайней мере 2 атрибута из шейп-файла. Или может быть только один атрибут до того, как столбцы будут отброшены? – BKS
Спасибо за ссылку. Я нашел новую функцию для создания linnet из объекта SpatialLines, но, как вы сказали, я потерял атрибуты :( – BKS