Как получить баллы (лат. И лонг.), Которые будут отображаться на карте шейп-файлов?

Question

Я пытаюсь нанести точки с координатами возникновения насилия на карту шейп-файла Нью-Йорка. Я применил скрипт Google, чтобы превратить исходные адреса улиц набора данных в координаты широты и долготы. Я экспортировал его .csv в RStudio. Очистите данные немного больше, удалив ненужные столбцы и значения NA + изменив Lat. столбец к цифре c.

Кажется, я все сделал правильно, пока не наложил точки на карту шейп-файла. Когда я запускаю следующий код, он возвращает координаты отдельно от карты (см. Изображения в приложении). То есть они не наслоены друг на друга, поэтому я в конечном итоге смогу использовать их для создания картограммы. Кроме того, изображение точки ниже, кажется, не показывает все координаты широты / долготы в наборе данных. В общей сложности около 500 инцидентов с предоставленными данными о координатах разбросаны по всему штату Нью-Йорк. Я не очень уверен в том, что показано, но это, вероятно, топи c для другого вопроса.

library(data.table)
library(sp)
library(rgdal)
library(ggplot2)


df_2 <- Gun_Violence_Clean_3 %>%
  select(-`Incident ID`, -Operations, -`City 2`, -Combine, -`Lat & Long`) %>%
  na.omit()
df_2

df_2$Lat <- as.numeric(df_2$Lat)

coordinates(df_2) = c("Lat","Long")
crs.geo1 = CRS("+proj=longlat")  
proj4string(df_2) = crs.geo1  

plot(df_2, pch = 20, col = "steelblue")


New_York = readOGR(dsn = "./NYS Boundaries", layer = "new-york-state-city-and-town- 
boundaries")
plot(New_York)
points(df_2, pch = 20, col = "orange")

Воспроизводимые данные:

structure(list(`Incident ID` = c(1664753, 1664770, 1664768, 1664751, 
1664723, 1664721), `Incident Date` = c("23-Apr-20", "22-Apr-20", 
"22-Apr-20", "22-Apr-20", "22-Apr-20", "22-Apr-20"), State = c("New 
York", 
"New York", "New York", "New York", "New York", "New York"), 
`City Or County` = c("Buffalo", "Schenectady", "Schenectady", 
"Albany", "Brooklyn", "Corona (Queens)"), Address = c("50 block of 
Langmeyer Ave", 
"1009 McClyman St", "1013 McClyman St", "200 block of Second Ave", 
"255 Havemeyer St", "225-37 Murdock Ave"), `# Killed` = c(1, 
0, 0, 0, 0, 0), `# Injured` = c(0, 0, 1, 1, 0, 1), Operations = 
c("N/A", 
"N/A", "N/A", "N/A", "N/A", "N/A"), `City 2` = c("Buffalo", 
"Schenectady", "Schenectady", "Albany", "Brooklyn", "Corona (Queens)"
), Combine = c("50 block of Langmeyer Ave Buffalo", "1009 McClyman St 
Schenectady", 
"1013 McClyman St Schenectady", "200 block of Second Ave Albany", 
"255 Havemeyer St Brooklyn", "225-37 Murdock Ave Corona (Queens)"
), `Lat & Long` = c("42.92484, -78.815534", "#ERROR!", 
"42.80176729999999, -73.9331919", 
"42.6390962, -73.77026289999999", "40.7079479, -73.95942509999999", 
"40.703342, -73.731005"), Lat = c("42.92484000000", "#ERROR!", 
"42.80176730000", "42.63909620000", "40.70794790000", "40.70334200000"
), Long = c(-78.815534, NA, -73.9331919, -73.7702629, -73.9594251, 
-73.731005)), row.names = c(NA, -6L), class = c("tbl_df", 
"tbl", "data.frame"))

Очки Карта Нью-Йорка

Снимок экрана набора данных

Где я получил шейп-файл штата Нью-Йорк

Большое спасибо за вашу помощь.

Answer 1

Использование этого ответа в качестве ориентира: значения широты и долготы за пределами границ в преобразованном файле формы с использованием ggplot

Вот минимальное минимальное решение:

library ( rgdal) библиотека (sp)

df_2 <- structure(list(`Incident Date` = c("23-Apr-20", "22-Apr-20", 
"22-Apr-20", "22-Apr-20", "22-Apr-20"), `City Or County` = c("Buffalo", 
"Schenectady", "Albany", "Brooklyn", "Corona (Queens)"), `# Killed` = c(1, 
0, 0, 0, 0), `# Injured` = c(0, 1, 1, 0, 1), Lat = c(42.92484, 
42.8017673, 42.6390962, 40.7079479, 40.703342), Long = c(-78.815534, 
-73.9331919, -73.7702629, -73.9594251, -73.731005)), row.names = c(NA, 
-5L), na.action = structure(c(`2` = 2L), class = "omit"), class = c("tbl_df", 
"tbl", "data.frame"))

#read shape file
New_York = readOGR(dsn = "./NYS Boundaries", layer = "Counties")

#Transform the shape file coordinates to Lat/Longitude
NY <-spTransform(New_York, CRS("+proj=longlat +ellps=WGS84 +datum=WGS84"))
plot(NY, xlab="long")
axis(1)  #Check the coordinates
axis(2)
#plot the points
points(x=df_2$Long, y=df_2$Lat, pch = 20, col = "orange")