Проблема с системой координат при экспорте полигонов из растра dBBMM - PullRequest
Новая
       42

Проблема с системой координат при экспорте полигонов из растра dBBMM

0 голосов
/ 14 декабря 2018

Я пытаюсь вычислить распределение использования по данным отслеживания животных, создавая динамическую модель движения Броуновского моста с функцией brownian.bridge.dyn, используя объект перемещения. Вот скрипт для созданияпереместите объект и затем растр с помощью brownian.bridge.dyn 

PR <- move(x=bird$Long, y=bird$Lat, time=as.POSIXct(muleys$PosTime, format="%Y%m%d %H:%M:%S"), proj=CRS("+proj=longlat +ellps=WGS84"),data=bird, animal=bird$ID)
PR_t <- spTransform(x=PR, CRSobj="+proj=aeqd +ellps=WGS84", center=TRUE)
PR_dbbmm <- brownian.bridge.dyn(object=PR_t, location.error=22, window.size=35, ext=0.4, margin=7, dimSize=30, time.step=5)

Теперь мне нужно создать указанные контуры и экспортировать их в виде многоугольников в формате .shp для визуализации в QGIS. Для этого я нашел следующий скрипт: 

dbbmm.px <- as(PR_dbbmm, "SpatialPixelsDataFrame")
image(dbbmm.px)
dbbmm.ud <- new("estUD",dbbmm.px)
dbbmm.ud@vol = FALSE
dbbmm.ud@h$meth = "dBBMM"

shp50 <- getverticeshr(dbbmm.ud, percent=50, standardize=TRUE)
class(shp50)
map.ps50 <- SpatialPolygons2PolySet(shp50)
diss.map.50 <- as.PolySet(map.ps50, projection = "LL", zone = "35N")
diss.map.p50 <- PolySet2SpatialPolygons(diss.map.50, close_polys = TRUE)
diss.map.p50p <- spTransform(x=diss.map.p50, CRSobj="+ellps=WGS84 +proj=longlat +datum=WGS84 +no_defs", center=TRUE)
diss.map.p50 <- diss.map.p50p
data50 <- data.frame(PID = 1)
diss.map.p50 <- SpatialPolygonsDataFrame(diss.map.p50, data = data50)
writeOGR(diss.map.p50, dsn = ".", layer="PR_dBBMM_contour50", driver = "ESRI Shapefile")

Я бы хотел, чтобы полигоны были в проекции WGS84 .Сценарий, кажется, работает и создает полигоны.Когда я загружаю файл .shp в QGIS, кажется, что полигоны находятся в точном месте.Я загружаю исходные точки данных, используемые для расчета dBBMM в проекции WGS84, а также некоторые .shp с границами страны. Кажется, что все слои перекрывают НО, когда вы видите многоугольники, созданные с помощью описанного выше сценария, точки похожи на одну точку позади, поэтому вам нужно увеличивать масштаб, пока не увидите нормальное распределение точек на карте играницы страны. Похоже, что полигон намного выше, чем другие слои, поэтому вам нужно значительно уменьшить масштаб, чтобы увидеть их правильно.Это выглядит очень странно, и я понятия не имею, как это исправить и что не так с используемым мной сценарием.

Любая помощь будет принята с благодарностью.Обратите внимание, что я новичок в R

...