空間データの読み込み

R で地理空間データを可視化する

Charlotte Wickham

Assistant Professor at Oregon State University

ニューヨーク郡の中央値所得

  • センサストラクトは、ほぼ同程度の人口の区域
  • 空間オブジェクト:
    • センサストラクトのポリゴン
    • より大きな近隣区域のポリゴン
    • 水域のポリゴン

nyc_map.png

R で地理空間データを可視化する

手順

  • 近隣区域と水域を表すシェープファイルを読み込む
  • 2 つの異なる座標系を整合
  • データフレームのデータを SpatialPolygonsDataFrame に結合
  • 公開用に地図を仕上げる

nyc_map.png

R で地理空間データを可視化する

シェープファイルの読み込み

  • ベクタデータ: 点・線・ポリゴンで表すデータ
  • 最も一般的な形式はシェープファイル
R で地理空間データを可視化する

シェープファイルの読み込み

# rgdal::readOGR() はベクタ形式を読み込む
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")

ドライバ: ESRI Shapefile の OGR データソース 
ソース: "water", レイヤー: "water-areas"
20 フィーチャ
フィールドは 5 個
R で地理空間データを可視化する

結果の確認

summary(water)

SpatialPolygonsDataFrame クラスのオブジェクト
座標:
        min       max
x -74.04731 -73.90866
y  40.68419  40.88207
投影あり: FALSE 
...

plot(water)

nyc_lines.png

R で地理空間データを可視化する

ラスタファイルの読み込み

library(rgdal)  # rgdal::readGDAL() はラスタ形式を sp オブジェクトに読み込む
library(raster) # raster::raster() はラスタ形式を raster オブジェクトに読み込む

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")
R で地理空間データを可視化する

結果の確認

total_pop

class       : RasterLayer 
dimensions  : 3120, 7080, 22089600  (nrow, ncol, ncell)
resolution  : 0.008333333, 0.008333333  (x, y)
extent      : -125, -66, 24, 50  (xmin, xmax, ymin, ymax)
coord. 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 
names       : uspop300 
values      : 0, 65535  (min, max)
R で地理空間データを可視化する

Vamos praticar!

R で地理空間データを可視化する

Preparing Video For Download...