значения широты и долготы за пределами границ в преобразованном файле формы с помощью ggplot - PullRequest
0 голосов
/ 30 сентября 2019

Я пытаюсь составить карту округов Пеннильвании, используя данные из PA Geospatial Data Clearinghouse . Я читаю в файле формы и затем конвертирую его в фрейм данных примерно так:

setwd(proj_path)
dsn <- "Data/PaCounty2019_07/PaCounty2019_07.shp"
map <- readOGR(dsn)
map_df <- tidy(map)

Затем я рисую его так:

ggplot() +
geom_path(data = map_df, aes(x = long, y = lat, group = group)) +
labs(title = "ggplot map of pa counties")

Это работает, но масштаб широта / долготаочень странноШирота идет от -100 000 до 150 000, а долгота идет от -2e + 05 до 23 + 05. Если я попытаюсь нанести на карту какие-либо точки, похоже, что все они имеют значение [0,0], потому что значения широты и долготы для точек, конечно, все <180. <a href="https://i.stack.imgur.com/tL5Ln.png" rel="nofollow noreferrer"> Изображение карты PA со странным масштабом здесь

Почему масштаб настолько искажен? Как я могу это исправить?

1 Ответ

0 голосов
/ 30 сентября 2019

Вам необходимо преобразовать шейп-файл в нормальную широту / долготу проекции. Сделайте это до определения map_df ...

map <- spTransform(map, CRS("+proj=longlat +ellps=WGS84 +datum=WGS84"))

enter image description here

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...