ホーム  •  ニュース  •  フォーラム  •  アカウント情報  •  サイト内検索  •  新規登録
 ログイン
ユーザ名:

パスワード:


パスワード紛失

新規登録
 カウントダウンブロック
カウントダウンイベントはありません
 メニュー
メイン
   地球科学のための R
     地図に行政区分を入れる		投稿するにはまず登録を

スレッド表示 | 新しいものから 前のトピック | 次のトピック | 下へ
投稿者 スレッド
石戸谷
投稿日時: 2004-11-25 17:39
研究室メンバ
登録日: 2004-7-30
居住地:
投稿: 13
地図に行政区分を入れる
適当なライブラリを入れた後、例えば、

source("~/lib/R/filledcontour.R")
library(akima)
library(maps)
library(maptools)

map <- read.shape("~/lib/R/SHP_MAP/japan_jdg.shp")

x <- seq(139,140,0.1)
y <- 35.2 + runif(11)
z <- 1:11

filled.contour(interp(x,y,z))

plot(map, fg="transparent", add=T)

とする。ただし、japan_jdg.shp は、http://www.esrij.com/gis_data/ から取得する。filled.contour() は、地図データを載せられるようにパッチを当てたものである。
森
投稿日時: 2004-11-29 19:51
登録日: 2004-7-29
居住地: 地球
投稿: 303
Re: 地図に行政区分を入れる
石戸谷さんのために、特定の行政区分だけ描くようにしたいと思っています。その方が速く描けるからです。

しかし、どうもうまくいきません。Windows であれば、このデータでうまく行くかもしれません。しかし、行政区分の名前が全て日本語になっているため、Linux 上の R ではうまくいきません。

そこで、国土地理院のデータの地球地図日本からbnda.zip をダウンロードし、同様に shape 形式のファイルを取り出しこれを使おうと思いました。

しかし、依然として、たとえば、TOKYO だけのデータを描く方法がわかりません。

mapselect <- function(map=map,PREF="TOKYO"){
  len <- length(map$att$NAM[map$att.data$NAM==PREF])
  for( i in 1:len ){
    plot(map$Shapes[map$att.data$NAM==PREF][[i]]$vert,type='l',add=T)
  }
}

を用意しておいて、
 library(maptools)
 data <- read.shape("bnda.shp")

として、
mapselect(data,"KANAGAWA")

とでもすればうまくいくような気がしますが、add=T が有効にききません。どうすればいいのか…。

いずれにしろ、日本地図全体をいちいち描くようにすれば行政区分を入れることは可能です。これを読んでわからなければコンタクトをとって下さい。
森
投稿日時: 2005-12-18 20:54
登録日: 2004-7-29
居住地: 地球
投稿: 303
Re: 地図に行政区分を入れる
各都道府県を選んで地図に描く方法は別にあるようです。自分では試していません。次のリンクを参照して試してみて下さい。(特にM君)

青木先生(群馬大)による都道府県の地図を扱うためのページ
スレッド表示 | 新しいものから 前のトピック | 次のトピック | トップ

投稿するにはまず登録を
 

WWW を検索 meteorology.jp を検索

Powered by XOOPS 2.0 © 2001-2006 The XOOPS Project, Maitained by A. Mori
FI Theme :: XOOPS 2 Theme by ImageSquare :: Costomized by matchan and A.Mori