如何處理大型連雀網目標

Question

我正在嘗試使用整個城市網絡進行特定分析，我知道這個分析非常龐大。我也將其設置爲稀疏網絡。

library(maptools) 
library(rgdal) 
StreetsUTM=readShapeSpatial("cityIN_UTM") 

#plot(StreetsUTM) 
library(spatstat) 
SS_StreetsUTM =as.psp(StreetsUTM) 
SS_linnetUTM = as.linnet(SS_StreetsUTM, sparse=TRUE) 

> SS_linnetUTM 
Linear network with 321631 vertices and 341610 lines 
Enclosing window: rectangle = [422130.9, 456359.7] x [4610458, 
4652536] units 
> SS_linnetUTM$sparse 
[1] TRUE 

我有以下問題：

1. 花了15-20分鐘打造psp對象
2. 花了近5個小時打造linnet對象我想每次
3. 到分析它爲一個點模式或信封，R崩潰

我明白我應該儘量減少網絡的大小，但是：

1. 我想知道是否有一個聰明的方法來克服這個問題。重新調整幫助？
2. 我怎麼能把它放在更多的處理能力？
3. 我也很好奇知道spatstat是否可以使用parallel包
4. 最後，spatstat的網絡大小有什麼限制。

ř崩潰當我使用從Spatstat書說明書

ř崩潰：

1. KN <- linearK(spiders, correction="none");在我的網絡上（linnet）當然是

2. envelope(spiders, linearK, correction="none", nsim=39);我的網絡上

我不認爲RAM的問題，我有16GB內存和2.5GHz的雙核i5處理器的SSD的機器上。

請有人指導我。

Answer 1

請更具體地說明您使用的命令。

你有沒有建立使用as.linnet.psp一個psp對象linnet對象（在這種情況下，網絡的連接一定要猜測，這可能需要很長的時間），還是你有關於網絡的連通性信息你傳遞給linnet()命令？

究竟是什麼命令來「分析點模式或信封」導致崩潰，什麼樣的崩潰？

spatstat中的線性網絡代碼是仍在開發中的研究代碼。更快的K函數算法將很快發佈。

Answer 2

我只能用Simplify Geometries工具中的Douglas-Peucker算法簡化QGIS中的網絡來解決此問題。所以這是shapefile中線性網絡幾何的輕微折中。