Блог

У меня есть данные о широте / долготе для разрешений на строительство в городе Чикаго, и я надеюсь перейти по ним к почтовым индексам. У меня есть список всех почтовых индексов Чикаго, которые я могу использовать, но я не уверен, как решить эту проблему.

У меня есть две переменные, chi_zips_map которая содержит геопространственные данные через tigris пакет, и mapping_data_19 которая содержит мои широтные / длинные данные.

Я смог успешно сопоставить как почтовые индексы, так и данные по широте / длине с помощью листовки, используя следующий код:

 chi_zips_map %>%
  leaflet() %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addPolygons(weight = '0.5', color = 'purple') %>%
  addCircleMarkers(data = mapping_data_19, 
                   lng = ~longitude, 
                   lat = ~latitude, 
                   radius = 2, 
                   color = ~permitpal(permit_type))
  

и мы получили следующую карту, которая великолепна:

Карта points.zip

Возможно ли для меня добавить столбец почтового индекса в mapping_data_19 фрейм данных на основе процессов, которые я выполнил до сих пор?