［R］SpatialLinesDataFrameオブジェクトから座標データを簡単に抜き出す
ggplot2パッケージに含まれる関数fortifyを使うと、SpatialLinesDataFrameオブジェクトから座標データを簡単なデータフレーム形式に変換してくれる。以下に、線が20含まれるシェープファイルmap.shpを使用した例を示す。
> library(maptools)
> library(ggplot2)
> shp <- readShapeLines("map.shp")
> shp
（大量にデータが表示される）
> length(shp) # 線の数
[1] 20
> df <- fortify(shp) # データフレーム形式に変換
> is.data.frame(df)
[1] TRUE
> head(df, 3)
long lat order piece group id
1 ・・・
2 ・・・
3 ・・・
> tail(df, 3)
long lat order piece group id
○○ ・・・
○○ ・・・
○○ ・・・
含まれる線の識別は列idを見ればよい。上記の例でいうと、列idは1から始まり線が変わることにidが1ずつ増えて最終行のidは20になる。
« SpatialLinesDataFrameオブジェクトから座標データを抜き出す。 | トップページ | 1回の繰り返し計算（ループ）で分散を求める »
「R（ggplot2）」カテゴリの記事
- ［R］ggplot2パッケージによる図を、簡単に並べて表示する(2014.12.16)
- ［R］ggplot2による図の枠線内の描画範囲（軸の範囲）をxlim、ylimで指定をした値きっかりにしたい(2014.12.09)
- ［ggplot2］点を描画する(2014.04.04)
- ［R］SpatialLinesDataFrameオブジェクトから座標データを簡単に抜き出す(2014.01.25)
- ［ggplot2］複数のデータを重ねて図に描画する(2014.01.17)
「R（数値計算）」カテゴリの記事
- ［R］複数のパラメーターと定数を持つ関数の値が最大・最小となるパラメーターを求める(2026.03.30)
- ［R］重み付き線形最小二乗法を行う(2026.03.20)
- ［R］重み付き非線形最小二乗法を行う(2026.03.06)
- ［R］計算機イプシロンを求める(2025.10.29)
- ［R］複数の引数を持つ関数の値の最小値（最大値）を求める(2024.12.07)
「R（maptools）」カテゴリの記事
- ［R］SpatialLinesDataFrameオブジェクトから座標データを簡単に抜き出す(2014.01.25)
- SpatialLinesDataFrameオブジェクトから座標データを抜き出す。(2014.01.24)
- エラーメッセージ「Not polygon shapes」(2014.01.24)
« SpatialLinesDataFrameオブジェクトから座標データを抜き出す。 | トップページ | 1回の繰り返し計算（ループ）で分散を求める »
コメント