Как преобразовать координаты в распознаваемое измерение для карт листовок?

Я новичок в функциях сопоставления в R и связанных пакетах, поэтому любая помощь приветствуется. Вот что я делаю: я использую листовку на панели инструментов, которую создал с помощью shiny. Я пытаюсь сопоставить имеющиеся у меня данные, но буклет не распознает их, потому что мои широта и долгота измеряются в футах, а не в десятичных градусах (пример: x = 33,867022, y = -112,1409750 ).

Я знаю, что у меня есть данные широты и долготы в футах (пример: x = 632086,70, y = 1043211,35), потому что я легко смог отобразить это в ArcMap, просто используя функцию «отобразить xy». Затем я смог преобразовать формат измерения с помощью функции «вычислить геометрию» для преобразования из футов в десятичные градусы. Я пытаюсь автоматизировать свои карты и предпочел бы НЕ использовать ArcMap для этого преобразования, поскольку я хотел бы остаться на одной платформе (R).

Мои возможные варианты (как мне кажется) и поэтому вопросы следующие:

1. Преобразуйте мои данные в то, что буклет может распознать. Кто-нибудь знает код, шаги и / или пакет, который мне нужно использовать для преобразования измерения моих широтных и длинных данных из футов в десятичные градусы?

2. Сообщить буклету, чтобы он считывал мои данные в другом измерении. Есть ли способ (код, упаковка и т. д.), чтобы буклет читал мои баллы в футах, а не в десятичных градусах?

Или какие-либо другие предложения???


r leaflet gis shiny coordinate-transformation
person SeaSun    schedule 27.03.2017    source источник
comment
Какой объект вы пытаетесь построить, это объект data.frame или spatial*? Если это пространственный объект, вы можете использовать функцию sp::spTransform() - пример здесь   -  person SymbolixAU    schedule 28.03.2017
comment
Итак, теперь ваш вопрос звучит так: «Как мне преобразовать data.frame координат в широту/долготу». На StackOverflow есть несколько примеров того, как это сделать. Если вы все еще застряли, перепишите свой вопрос, приведя пример ваших данных и предпринятых вами попыток.   -  person SymbolixAU    schedule 28.03.2017
comment
Однако вы можете пойти по пути варианта 2 в зависимости от того, как вы получаете данные в R?   -  person SymbolixAU    schedule 28.03.2017
comment
Какая у карты CRS (система отсчета координат)?? Измерение в футах слишком широкое, и вам нужно знать о CRS, чтобы проецировать координаты в EPSG: 4326 (широта-долгота равнопрямоугольные). Узнав об этом, взгляните на proj4 для R.   -  person IvanSanchez    schedule 28.03.2017

Ответы (1)


arrow_upward
0
arrow_downward

Спасибо за всю предоставленную помощь... Мой коллега понял это, мы использовали приведенный ниже код.

proj4string(d) <- CRS("+proj=tmerc +lat_0=31+lon_0=-111.9166666666667+k=0.9999 +x_0=213360 +y_0=0 +ellps=GRS80 +datum=NAD83+to_meter=0.3048006096012192+no_defs")
CRS.new <- CRS("+init=epsg:4326") # WGS 84
d.ch102650 <- spTransform(d, CRS.new)

@IvanSanchez комментарий CRS был полезен.

person SeaSun    schedule 29.03.2017