Mekansal verileri içe aktarma

R ile Coğrafi Verileri Görselleştirme

Charlotte Wickham

Assistant Professor at Oregon State University

New York County'de medyan gelirler

Nüfus sayımı bölgeleri, kabaca aynı nüfusa sahip alanlardır
Mekansal nesneler:
- Nüfus sayımı bölgesi poligonları
- Daha büyük mahalle poligonları
- Su alanı poligonları

Yöntem

Mahalleleri ve su yollarını tanımlayan shapefile'ları içe aktarın
İki farklı koordinat sistemini eşleyin
Bir veri çerçevesinden veriyi SpatialPolygonsDataFrame içine birleştirin
Haritayı yayıma hazır hale getirin

Bir shapefile'ı içe aktarma

Vektör veri: nokta, çizgi, poligon ile tanımlanan veriler
En yaygın biçim: shapefile

Bir shapefile'ı içe aktarma

# rgdal::readOGR() vektör biçimlerini okur
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")

Sürücü: ESRI Shapefile ile OGR veri kaynağı 
Kaynak: "water", katman: "water-areas"
20 öznitelik var
5 alanı var

Sonucu kontrol etme

summary(water)

SpatialPolygonsDataFrame sınıfından nesne
Koordinatlar:
        min       max
x -74.04731 -73.90866
y  40.68419  40.88207
Projeksiyonlu: FALSE 
...

plot(water)

Bir raster dosyayı içe aktarma

library(rgdal)  # rgdal::readGDAL() raster biçimlerini sp nesnelerine okur
library(raster) # raster::raster() raster biçimlerini raster nesnelerine okur

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")

Sonucu kontrol etme

total_pop

sınıf        : RasterLayer 
boyutlar     : 3120, 7080, 22089600  (nrow, ncol, ncell)
çözünürlük   : 0.008333333, 0.008333333  (x, y)
kapsam       : -125, -66, 24, 50  (xmin, xmax, ymin, ymax)
koord. ref.  : +proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0 
veri kaynağı : /Users/wickhamc/Documents/Projects/courses-visualizing-geospatial-
              data-in-r/data/census_grids/usgrid_data_2000/uspop300.tif 
adlar        : uspop300 
değerler     : 0, 65535  (min, max)

Haydi pratik yapalım!

R ile Coğrafi Verileri Görselleştirme