R googleway с данными переписи ABS - PullRequest
Новая
       10

R googleway с данными переписи ABS

0 голосов
/ 29 августа 2018

Я собираюсь проанализировать данные переписи в Австралии , используя googleway для составления тепловых карт.

Мой подход состоял в том, чтобы использовать подготовленные данные из googleway melbourne, который содержит столбец SA2_NAME, и объединить его с файлом формы ESRI из данных переписи после преобразования с помощью rgdal (код ниже). Проблема в том, что соединение по SA2_NAME не является уникальным для полилиний - некоторые области SA2 состоят из нескольких «подзон». Так что, похоже, это не очень хороший подход.

Лучшим подходом было бы преобразование данных формы ESRI sa2_shape ниже для получения полилиний в формате данных melbourne. Как это сделать?

Приведенный ниже код создает фрейм данных «мост» для использования при объединении данных melbourne из googleway с данными ABS, в которых ключевым полем является SA2_MAIN - как указано выше, проблема с этим подходом «хака» заключается в том, что что полилинии не являются уникальными SA2_NAME

library(tidyverse)
library(googleway)
library(rgdal)

shape_path <- "abs_data/sa2_esri_shapefile"
shape_file <- "SA2_2016_AUST"
sa2_shape <- readOGR(shape_path, shape_file)
sa2_df <- data.frame(sa2_shape$SA2_MAIN, sa2_shape$SA2_NAME)
names(sa2_df) <- c("SA2_MAIN", "SA2_NAME")
sa2_df <- sa2_df %>% semi_join(melbourne, by = "SA2_NAME")

1 Ответ

0 голосов
/ 29 августа 2018

Согласно комментарию SymbolixAU - используется sf для загрузки данных, и это работает, пока geometry не является пустым списком - см. Код ниже. 

library(tidyverse)
library(googleway)
library(sf)

shape_path <- "abs_data/sa2_esri_shapefile"
shape_file <- "SA2_2016_AUST"

shape_file_path <- paste0(shape_path, "/", shape_file, '.shp')
sa2_shape <- sf::st_read(shape_file_path)
sa2_shape <- sa2_shape %>% 
  filter(STATE_NAME == "Victoria",
         AREA_SQKM > 0)# This is important - otherwise google_map() will crash!
google_map() %>% 
  googleway::add_polygons(data = sa2_shape, 
                          polyline = "geometry", 
                          fill_colour = "SA2_NAME")

Victoria SA2 map

> sa2_shape %>% head()
Simple feature collection with 6 features and 6 fields
geometry type:  MULTIPOLYGON
dimension:      XY
bbox:           xmin: 143.6849 ymin: -37.68153 xmax: 143.951 ymax: -37.46847
epsg (SRID):    4283
proj4string:    +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs
   SA2_MAIN SA2_MAIN16         SA2_NAME STATE_CODE STATE_NAME AREA_SQKM                       geometry
1 201011001  201011001        Alfredton          2   Victoria   52.7111 MULTIPOLYGON (((143.7072 -3...
2 201011002  201011002         Ballarat          2   Victoria   12.3787 MULTIPOLYGON (((143.8675 -3...
3 201011003  201011003 Ballarat - North          2   Victoria   92.3577 MULTIPOLYGON (((143.853 -37...
4 201011004  201011004 Ballarat - South          2   Victoria   32.8541 MULTIPOLYGON (((143.8675 -3...
5 201011005  201011005        Buninyong          2   Victoria   51.5855 MULTIPOLYGON (((143.8533 -3...
6 201011006  201011006        Delacombe          2   Victoria   34.1608 MULTIPOLYGON (((143.7072 -3...
...