Ruimtelijke data inlezen

Georuimtelijke data visualiseren in R

Charlotte Wickham

Assistant Professor at Oregon State University

Mediaaninkomens in New York County

Censusgebieden hebben ongeveer evenveel inwoners
Ruimtelijke objecten:
- Censustract-polygons
- Grotere buurtpolygons
- Waterpolygonen

Stappenplan

Shapefiles van buurten en waterwegen inlezen
Twee verschillende coördinatensystemen afstemmen
Data uit een data frame mergen in een SpatialPolygonsDataFrame
Een kaart polijsten tot publicatieklaar

Een shapefile inlezen

Vectordata: data als punten, lijnen, polygons
Shapefile is het meest gangbare formaat

Een shapefile inlezen

# rgdal::readOGR() leest vectorformaten in
library(rgdal) 
library(sp)

dir()

"water"

dir("water")

"water-areas.dbf" "water-areas.prj" 
"water-areas.shp" "water-areas.shx"

water <- readOGR("water", "water-areas")

OGR-databron met driver: ESRI Shapefile 
Bron: "water", laag: "water-areas"
met 20 features
Bevat 5 velden

Resultaat controleren

summary(water)

Object van klasse SpatialPolygonsDataFrame
Coördinaten:
        min       max
x -74.04731 -73.90866
y  40.68419  40.88207
Geprojecteerd: FALSE 
...

plot(water)

Een rasterbestand inlezen

library(rgdal)  # rgdal::readGDAL() leest rasterformaten in als sp-objecten
library(raster) # raster::raster() leest rasterformaten in als rasterobjecten

dir()

"usgrid_data_2000"   "usgrid_data_2000_1"

dir("usgrid_data_2000")

"metadata"     "usarea00.tif"  "usba00.tif" "usfb00.tif" "usgrid-2000-variables.xls" "usp2500.tif" 
"uspop300.tif" "uspov00.tif"  "uspvp00.tif"

total_pop <- raster("usgrid_data_2000/uspop300.tif")

Resultaat controleren

total_pop

class       : RasterLayer 
afmetingen  : 3120, 7080, 22089600  (nrow, ncol, ncell)
resolutie   : 0.008333333, 0.008333333  (x, y)
extent      : -125, -66, 24, 50  (xmin, xmax, ymin, ymax)
coörd. ref. : +proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0 
data source : /Users/wickhamc/Documents/Projects/courses-visualizing-geospatial-
              data-in-r/data/census_grids/usgrid_data_2000/uspop300.tif 
namen       : uspop300 
waarden     : 0, 65535  (min, max)

Laten we oefenen!

Georuimtelijke data visualiseren in R