你好所有潜在的帮手,
我从SpatialPolygonDataFrame包中获得了一个tigris对象,我希望在创建ppp对象时使用它作为多边形窗口。以下是我尝试过的:
# Needed packages
library(spatstat)
library(sf)
# Download geospatial data for Lee county in Alabama (home of the great Auburn University by the way!)
county <- tigris::county_subdivisions(state = "Alabama
我试着用克罗斯做空间点模式分析。我希望从我的数据集“鸟”中比较物种A和物种B。感兴趣的变量是species_name。它只有两个级别- species A和species B。
下面是我尝试将species_name设置为多类型对象的代码。
marks(birds) <- factor("species_name") #where birds is a ppp object
Kcross(birds,i = "species A", j = "species B")
但是,我遇到了错误消息:No points have mark i = s
我总是犯错误
# 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(neighborhoo
我想使用我的空间数据集(sp_ds)作为ppp (spatstat Point Pattern object)来进行完整的空间测试,我尝试这样做: # Packages
library(spatstat)
library(sf)
library(sp)
library(raster) ## maybe a cause of pointDistance() function, not used yet
# Open spatial data set in GitHub
sp_ds<-read.csv("https://raw.githubusercontent.com/Lepr
我想写一个文件搜索代码,我不知道我正在搜索的目录是否有子目录,我想检查它,这样我就不会得到这样的错误: error 267目录名称无效:'C:/Path/ to / directory‘。我写了一段类似这样的代码,如果它找到了文件,它就会中断并停止程序,如果没有,它就会下一层,依此类推。
filename = raw_input('> ')
path = 'C:/Path/Of/Directory/You/Want/To/Search/In'
fldr = os.listdir(path)
for f in fldr:
p = path
我正在尝试从CRC spatstat书(第17章)中了解一个小的linnet对象的linearK估计,当我使用linearK函数时,spatstat抛出了一个错误。我已经在下面r代码的注释中记录了这个过程。错误如下所示。
Error in seq.default(from = 0, to = right, length.out = npos + 1L) : 'to' cannot be NA, NaN or infinite
我不知道如何解决这个问题。我正在遵循这个过程:
# I have data of points for each data of the week
# d
我正在编写一个C代码来打开txt文件,并在上面读取两行,然后打印它工作了1018次的值,然后给出“分段错误”,我尝试过刷新缓冲区,但它不起作用。
while(running) {
i = 0;
if ((fptr = fopen("pwm.txt", "r")) == NULL) {
printf("Error! File cannot be opened.");
// Program exits if the file pointer returns NULL.
exit(1);