spatstat

    0熱度

    2回答

    要了解spatstat中的一些繪圖,我試圖使用分割函數繪製數據庫粘膜以創建一個多元繪圖。我得到了兩個不同的粘膜小區,但默認情況下它們都具有相同的符號系統。 library(spatstat) plot(split(mucosa)) 現在兩個次要情節,我想自定義每個插曲在顏色,尺寸方面,形狀等我怎樣才能做到這一點?請指教。 我試圖通過PCH和CEX和col兩個元素的矢量(因爲我知道這是兩種類

    0熱度

    2回答

    幫助density`功能「時間迭代沒有超過iterMax;增加DTT或增加iterMax」中的R爲density功能 如何增加iterMax? 我正在使用庫spatstat。 這是我的腳本: D660 < - 密度(取消標記(control23_network.ppp),660) 線性網絡上我的觀點模式有: 1969年點 線性網絡20126頂點和21363線 圍窗口:矩形= [-87.63141,

    0熱度

    1回答

    我有大約8K分段的每個20米的網絡。隨機泊松過程在這個網絡上創建一個實現。接下來,我想添加另一個泊松過程,但第二個過程的點不應該落在已經包含來自第一個過程的點的段上。 Q1>很大程度上,我古玩知道,如果每線段的點可以被定義。我的理解是這是不可能的,因爲泊松過程的性質。但是可能有一個可選的參數來限制每個段的點數?我知道可以限制整個linnet對象上的點數,但是我想知道這是否可以用於每個linnet段

    0熱度

    1回答

    我正在閱讀spatstat教科書中的文本,並嘗試使用ppm學習模型擬合。 我創建了carteisan座標作爲協變量的模型。然後我想看看在模型上只有一個協變量的影響, model1 = ppm(chicago_ppp ~ x+y) plot(effectfun(model1, covname = "x")) ,但我得到的錯誤 Error in effectfun(model1, covname =

    0熱度

    1回答

    我有一個點模式,有標記。它是使用connected.ppp函數創建的,也是在spatstat包中。我想只繪製標記等於2的點。 我認爲要做到這一點,需要使用which.marks參數,但輸出仍然是包含所有標記的圖。這些標記肯定存在,因爲當我設置use.marks = FALSE時,標記不再使用。 plot.ppp(testconn, use.marks = TRUE, which.marks = 1

    0熱度

    1回答

    我學習做一些地塊在spatstat的情節音階符號(圓圈)。我創建了一種只有一種泊松點過程,並附加了一些標記 - 一個數值。現在,當我繪製時,情節默認情況下是圓形的。這些圈子是相對比例的,但這些圈子在我的網絡上太大了。我想問是否有一些縮放因子來控制在點實現中繪製的這些圓的大小。 所以,我從提取點圖像處理點,並使用他們的X,Y,我能夠與繪圖功能的CEX參數來控制。但是這並不能控制圖例的大小。所以我試圖

    0熱度

    1回答

    我想子集一個linnet對象來獲取一個linnet對象列表。在幫助文件中,我發現這可以通過owin對象來實現。但我得到的錯誤,因爲語法不清晰,我也不清楚是在調用子集中使用snip參數。 我想基於來自標記來用戶屬性拆分紅雀或基於窗口分割紅雀對象(每個窗口是從shape文件多邊形) 可能有人請建議如何子集一個紅雀根據一個因子/窗口獲得一個linnet對象列表: 在線段上使用用戶屬性/因子/標記。 標記

    0熱度

    1回答

    我想從一個函數中創建一個im對象,該函數根據相對於固定線的距離來分配一個值,即,如果該點離開如果這條線靠近該線,則該線的值接近零並且相反。我創建了以下功能: library(spatstat) W <- owin(xrange = c(-0.5, 4.5), yrange = c(-0.5, 3.5)) covValue <- function(x, y) { # definin

    0熱度

    1回答

    我想用包Spatstat創建協變量,如果點在多邊形之外,我有一個shp文件作爲多邊形以便爲協變量賦值。我的問題是在座標投影上存在問題,因爲當我完成這個過程時,我有一個奇怪的數字,請有人幫助我。我的代碼: #Creating covariate Z: W<- owin(xrange=c(767768.3,768773.6),yrange=c(517335.8,518044.8)) xp<-c(

    0熱度

    1回答

    我有以下DF: name x y A -47 134 B -11 311 C 49 100 D -40 138 E -33 233 F 30 134 我使用x &ÿCOORDS以產生沃羅努瓦剖分由封閉矩形: library(spatstat) library(deldir) rectangle <- owin(c(