Проблема с Географическим CRS при выполнении KernelUD (пакет adehabitatHR) - PullRequest
Новая
       34

Проблема с Географическим CRS при выполнении KernelUD (пакет adehabitatHR)

0 голосов
/ 11 июня 2019

Я пытаюсь создать ядро ​​с пакетом AdehabitatHR, используя пространственные данные, и, похоже, у меня проблема с системами координат (CRS).

Я преобразовал свой набор данных в SpatialPoint DataFrame, указал координаты (долготу и широту) и установил CRS ("+ proj = longlat + datum = WGS84").

Когда я пытаюсь выполнить kernelUD (из пакета adehabitatHR), я получаю ошибку "Географический CRS, предоставленный для несоответствующих данных".

Это мой набор данных: 

> str(data)
'data.frame':   4537 obs. of  4 variables:
 $ ID_geoYear: chr  "6098001_1" "6098001_1" "6098001_1" "6098001_1" ...
 $ dtime     : POSIXct, format:  ...
 $ lat       : num  -4.71 -2.89 -0.24 -1.9 3.18 3.97 -1.74 -0.14 0.65 
 -1.59 ...
 $ long      : num  -31.3 -32.3 -32.7 -33 -31.8 ...

> head(data)
       ID_geoYear        dtime        lat   long
45062  6098001_1 2007-10-13 13:52:00 -4.71 -31.29
45063  6098001_1 2007-10-14 01:55:00 -2.89 -32.33

Я указываю столбцы широты и долготы: 

>coord <- data[,c("long", "lat")]

Я создаю SpatialPointDataFrame: 

>df <- SpatialPointsDataFrame(coord, data, coords.nrs = numeric(0))

Я установил CRS моего SpatialPointDataFrame: 

> proj4string(df) <- CRS("+proj=longlat +datum=WGS84")

Кажется, у моих столбцов долготы и широты есть нормальные значения и нет NA: 

> df@bbox 
        min   max
 long -44.55  1.45
 lat  -33.30 37.81

> unique(is.na(df@data$lat))
[1] FALSE
> unique(is.na(df@data$long))
[1] FALSE

Когда я пытаюсь сделатьЯдро, я получаю сообщение об ошибке: 

> KUD <- adehabitatHR::kernelUD (df, h = "href")
Warning in adehabitatHR::kernelUD(df, h = "href") :
xy should contain only one column (the id of the animals) 
id ignored
Error in `proj4string<-`(`*tmp*`, value = CRS(pfs1)) : 
Geographical CRS given to non-conformant data: -106.217881356 
110.727881356

Я читал, что обычно вы получаете эту проблему, когда у вас странные позиции долготы или широты, или если вы переключаете долготу и широту, но я изучил свои данные иэто не так.Итак, я не знаю, в чем проблема !!

Заранее спасибо!

Марта.

...