Проблемы проецирования диапазона ядра ядра UTM tpeqd - PullRequest
Новая
       77

Проблемы проецирования диапазона ядра ядра UTM tpeqd

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

Я работаю, чтобы вычислить домашний диапазон леопардов, используя пакеты adehabitatHR и ks.Я успешно получил пространственный многоугольник и экспортировал его с помощью writeOGR, чтобы получить шейп-файл, но когда я пытаюсь просмотреть его в Google Планета Земля Про (я тупой пользователь Mac и не имею доступа к Arc), и он проецируется в зону UTM38, даже когда проекция, которую я дал, составляет 40N.Мои исходные точки данных GPS проецируются правильно, как и другой шейп-файл, предоставленный мне моими сотрудниками, который имеет проекцию: 

"+proj=tpeqd +lat_1=37.5464870070476 +lon_1=58.7026354113759 +lat_2=37.387284424216 +lon_2=58.8159213873777 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"

Я незнаком с этим форматом tqeqd, но он работает в ГИС.Вот результаты того, что я пробовал: spTransform to 40N = проекты в зоне 37 на севере Москвы;spTransform to 40S = проекты в зоне 38 к югу от Мадагаскара;скопировать аргумент CRS в формате tpeqd = проекты в Северной Японии.

Это может быть проблемой в моем домашнем расчете дальности, но я не знаю, с чего начать, или почему преобразование в 40N спроецирует это так далеко.Я приведу свой код преобразования ниже и буду рад даже простому руководству по прогнозам и аргументам CRS.Я также могу добавить любой соответствующий код по запросу.

## Create Shapefile of KDE shape
Kaveh_isovol@proj4string <- CRS("+init=epsg:32740")
Kaveh_95vol@proj4string <- CRS("+init=epsg:32740")
writeOGR(Kaveh_isovol,dsn="output", layer="Kaveh_core",driver="ESRI 
 Shapefile",overwrite=TRUE)
writeOGR(Kaveh_95vol,dsn="output", layer="Kaveh_95",driver="ESRI 
 Shapefile", overwrite=TRUE)
...