我总是得到错误
# polygon that to be window
neighborhoods <- st_read("neighborhoods/neighborhoods.shp")
# convert CRS to planar projection as recommended by (https://stackoverflow.com/questions/59597078/use-sf-polygon-object-as-window-in-spatstat)
neighborhoods_3857 <- st_transform(neighborhoods, crs = 3857)
# point that to be PPP spatstat
trees <- st_read("trees/trees.shp")
# convert to planar projection
trees_3857 <- st_transform(trees, crs = 3857)
问题是，"trees_3857"没有以EPSG3857坐标表示的数据帧列，因此"trees_3857"的功能列没有符合 EPSG 3857 的 x 和 y 列
q <- ppp(x=?, y=?, win=neighborhoods_3857)
我做了什么但错误
z <- as.ppp(trees_3857, win=neighborhoods_3857)
Error in as.ppp.sf(trees_3857, win = neighborhoods_3857): unused argument (win = neighborhoods_3857)
Traceback:
您可以从数据营地免费获取数据。
https://assets.datacamp.com/production/repositories/738/datasets/96a72364e69d872645038b3a6dc7c0dbcb1114d6/neighborhoods.zip
https://assets.datacamp.com/production/repositories/738/datasets/08a3684dc4d538d59ba051a64a834166883ab5d1/trees.zip
尽管您希望将数据转换为spatstat包中的类"ppp"的对象，但错误消息指示问题源于函数as.ppp.sf，该函数是sf包的一部分。
错误消息显示unused argument: win，这意味着函数无法识别或接受参数win。
只是为了使其更具挑战性，功能as.ppp.sf没有记录，也不可见......通过键入sf:::as.ppp.sf我们可以看到函数体并找出该函数只有一个参数，因此它不接受任何窗口信息。
这不是通用函数as.ppp设计为在spatstat包中工作的方式，所以如果你在spatstat中寻找help(as.ppp)，你感到困惑也就不足为奇了。
最好的解决方案是(如 Ege Rubak 建议的那样(转换点坐标，然后添加窗口信息：
z <- as.ppp(trees_3857)
Window(z) <- as.owin(neighborhoods_3857)
转换as.ppp和as.owin将使用sf中的代码执行，因此我不能保证它们会起作用。赋值Window(z) <-将在spatstat代码中执行，请参见help("Window<-.ppp")。