Блог

В настоящее время я работаю над анализом локальных перемещений животных с использованием динамических моделей движения броуновского моста, рассчитанных с помощью пакета move (версия 4.0.4). В прошлом код R, с которым я работал, функционировал без какого-либо предупреждения. Недавно, когда я проверяю проекцию объекта перемещения, а затем преобразую проекцию в AEQD, я получаю следующее предупреждение при использовании:

proj4string(move.ind) [1] " proj=longlat datum=WGS84 no_defs" Warning message: In proj4string(move.ind) : CRS object has comment, which is lost in output

Однако, когда я создаю растр, чтобы перепроектировать результаты в последующее использование:

 > x.ind = raster(xmn=min(ind$NewEastingUTM), xmx=max(ind$NewEastingUTM), ymn=min(ind$NewNorthingUTM), ymx=max(ind$NewNorthingUTM), crs = CRS("  proj=aeqd  datum=WGS84  lon_0=-79.29058  lat_0=25.73195"))
 

> proj4string(x.ind) [1] " proj=aeqd lat_0=25.73195 lon_0=-79.29058 x_0=0 y_0=0 datum=WGS84 units=m no_defs"

это предупреждение не появляется.

Я предполагаю, что предупреждение связано с моими версиями sp (1.4-4) и rgdal (1.5-27).

Теперь мой вопрос заключается в том, можно ли «игнорировать» это предупреждение или есть решение о том, как правильно работать с проекциями в свете обновленных пакетов sp и rgdal?

Для воспроизводимости я добавляю заголовок и str моего df ниже.

              Datetime station elasmo  location      Lat
1 2016-05-01 12:17:00  124083 e11454 new barge 25.74273
2 2016-05-01 12:19:00  124083 e11454 new barge 25.74273
3 2016-05-01 12:21:00  124083 e11454 new barge 25.74273
4 2016-05-01 12:33:00  124083 e11454 new barge 25.74273
5 2016-05-01 12:59:00  124083 e11454 new barge 25.74273
6 2016-05-01 13:07:00  124083 e11454 new barge 25.74273
        Lon NewEastingUTM NewNorthingUTM RE
1 -79.30235      670280.8        2848289  0
2 -79.30235      670280.8        2848289  0
3 -79.30235      670280.8        2848289  0
4 -79.30235      670280.8        2848289  0
5 -79.30235      670280.8        2848289  0
6 -79.30235      670280.8        2848289  0

'data.frame':   3543 obs. of  9 variables:
 $ Datetime      : POSIXct, format: "2016-05-01 12:17:00" ...
 $ station       : chr  "124083" "124083" "124083" "124083" ...
 $ elasmo        : chr  "e11454" "e11454" "e11454" "e11454" ...
 $ location      : chr  "new barge" "new barge" "new barge" "new barge" ...
 $ Lat           : num  25.7 25.7 25.7 25.7 25.7 ...
 $ Lon           : num  -79.3 -79.3 -79.3 -79.3 -79.3 ...
 $ NewEastingUTM : num  670281 670281 670281 670281 670281 ...
 $ NewNorthingUTM: num  2848289 2848289 2848289 2848289 2848289 ...
 $ RE            : num  0 0 0 0 0 0 0 0 0 0 ...```