RとPostGIS

RとPostGIS

　ここでは、Open GIS Simple Features Specification for SQL に準拠したPostGIS を　R 内より利用して、GISの空間検索、空間解析および空間分析の方法を説明していくことにする。

PostGIS とは†

PostgreSQL 用のオープンソースの空間データベースエンジン。カナダの Refractions 社が開発。PostgreSQL 上で大容量の地理空間ベクトルデータの蓄積・検索・解析・地図投影変換が可能。WKT を格納するフィールドのデフォルト名は"the_geom"。

WKT(Well Known Text)†

PostGIS は国際標準の地理情報格納形式であるOpenGIS シンプルフィーチャをベースにしたWKTというテキスト形式で図形情報を格納している*1。

シェープファイルの変換†

PostGISのユーティリティツール shp2pgsql を利用してSQL 文を生成する。

shp2pgsql シェープファイル名 テーブル名 データベース名 > SQLファイル 

使用例*2

createdb columbusdb #データベース作成
createland plpgsql columbusdb #PL/PGSQL の利用を可能にする
psql -d columbusdb -f <PostgreSQLのインストール先のディレクトリ>/share/contrib/postgis.sql #PostGIS の利用を可能にする
shp2pgsql columbus columbus columbusdb > columbus.sql #SQL文生成
psql -d columbusdb -f columbus.sql                    #テーブル作成 

Rとの接続†

rgdal†

読み込み†

library(rgdal)
meuse = readOGR("PG:dbname=postgis", "meuse2")
meuse[1:10,] # first ten rows
proj4string(meuse)

WKT(Well Known Text)†

　rgeos の利用

書込み†

RODBC + rgeos --- [R-sig-Geo] Connection between R and PostGISより引用†

ここではRODBC を利用して、PostgreSQL に接続する。ただし、この方法は 処理できるサイズの小さなものに限られる。

library(rodbc) #パッケージのロード
con <- odbcConnect("gisdb")
 ## query
 sql <- "SELECT ST_AsText(the_geom) from gis_data.gadm2 where name_2='Berlin'"
 ## execute query
d <- readWKT(sqlQuery(con, sql, stringsAsFactors=F)[1,1])

関連リンク†

• feature request: rpostgresql does not support wkb/wkt

RPostgreSQL†

Rdbi/RdbiPgSQL†

なかまさんのWin32 バイナリがあるが、PostgreSQL のバージョンが違うので利用できない？ ビルドも面倒なので今はパス。

凸包†

convexhull

ユニオン†

GeomUnion または memGeomUnion

ディゾルブ†

ユニオンの応用

オーバーレイ解析†

ポリゴン†

• intersection
• symdifference
• difference

オンデマンド・トポロジー†

α版でサポート予定

道路ネットワーク解析†

空間ウェイトマトリックス†

オンデマンド・トポロジーを利用

空間データマイニング†

リンク†

• PL/R & PostGIS
  • Boston GIS: Tutorial Booklet 上記の PDF 版
• Open Source Online Geospatial Interfaces to Multiscale Indicators WebGIS での利用
• Connecting PostGIS with R ODBC と　OGR の利用例
• An example of using PostGIS from R ODBC & QGIS
  • countries share a common boundary
• PostGIS-R - How to use R functions from PostgreSQL
  • How to use R functions from PostgreSQL deldir

• Quick Intro to R and PL/R - Part 1 - Postgres OnLine Journal

参考書†

• PostGIS in Action PLR が取り上げられる予定。2011年2月発行。//www.amazon.co.jp/Postgis-Action-Regina/dp/1935182269]] でも入手可能。Early Access 版は全章が完了。ソースコードも公開。翻訳の話もあるような。
  • Chapter 10 rgdal, sp, spatial dataframe への PostGIS データの読み込み、R による作図など。
  • 第２版 R 関連は？

最新の20件