r - Troubles to calibrate and to filter occurrences with flexsdm package - Stack Overflow

admin管理员组

文章数量:1346043

I am using flexsdm to create Species Distribution Models but there are some problems I am facing.

When I am creating the calibration area using flexsdm it returns a table 0 x 0

> # Define calibration area using a buffer method
> calibration_area <- calib_area(
+   data = occurrences_bin,
+   x = "longitude",
+   y = "latitude",
+   method = c("buffer", width = 500),  # Try a larger width
+   crs = crs(selected_raster_stack)
+ )
> summary(calibration_area)
< table of extent 0 x 0 >

I have already checked if there is any NA value in the 'occurrences_bin' and 'selected_raster_stack'.

> any(is.na(occurrences_bin$longitude))
[1] FALSE
> any(is.na(occurrences_bin$latitude))
[1] FALSE
> occurrences_bin <- na.omit(occurrences_bin)
> any(is.na(values(selected_raster_stack)))
[1] FALSE
> selected_raster_stack <- trim(selected_raster_stack)

#I have already checked if there crs for 'occurrences_bin' and 'selected_raster_stack'. 

# Load species occurrence data
occurrences <- fread("Databases/Occurences Folder/CleannedOcc_Classified_Nearctic.txt", 
                     stringsAsFactors = TRUE, encoding = "UTF-8")

# Convert occurrences to spatial object
wgs84<-"+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0" 
occurrences_sf <- st_as_sf(occurrences, coords = c("longitude", "latitude"), crs = wgs84)
# Convert species column to binary presence/absence format
# Use the sf object to create a binary presence/absence table
occurrences_bin <- dcast(occurrences_sf, 
                         st_coordinates(occurrences_sf)[,1] + st_coordinates(occurrences_sf)[,2] ~ species,
                         value.var = "species", fun.aggregate = length)

# Convert to data.table, if not already
setDT(occurrences_bin)

> # Recheck raster stack
> crs(occurrences_sf) ==  crs(selected_raster_stack)
[1] TRUE

# Geographical filtering using Moran's I method
occ_filtered <- occfilt_geo(
  data = occurrences_bin,
  x = "longitude",
  y = "latitude",
  env_layer = selected_raster_stack,
  method = c('moran'),
  prj = crs(selected_raster_stack)
)

Error in if (any(mor <= method[2])) {:
Missing value where TRUE/FALSE is required
Calls: occfilt_geo -> occfilt_geo_0
Execution interrupted

本文标签： rTroubles to calibrate and to filter occurrences with flexsdm packageStack Overflow