<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">すみません、どなたかご教授下さい。</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">地図表示をしたく、Rで格子状のポイントデータをラスタ変換すると</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">ラスタの中心にポイントが来ないのですが、、、これは</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">（1）以下のスクリプトが間違っている？</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">（2）そもそも経度の距離は緯度に依存するから、平均的なラスタ解像度とは、当然一致しない？</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">（3）そのほか・・・</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">ご助言を頂けると幸いです。</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">RとGISも中途半端な理解なので、根本的に間違っているかも知れませんが、</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">数日考えてもどうにも解決できず投稿した次第です。</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">----ここから----</div><div class="gmail_default"><div class="gmail_default"><font face="verdana, sans-serif">library(maptools)</font></div><div class="gmail_default"><font face="verdana, sans-serif">library(rgdal)</font></div><div class="gmail_default"><font face="verdana, sans-serif">library(raster)</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif">rm(list=ls(all=TRUE))</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif">#ダミーデータ（データフレーム）作成</font></div><div class="gmail_default"><font face="verdana, sans-serif">x <- rep(seq(121,140,1),10)</font></div><div class="gmail_default"><font face="verdana, sans-serif">y <- numeric(0);for ( i in 21:40 ){ y <- c(y, rep(i,20))}</font></div><div class="gmail_default"><font face="verdana, sans-serif">z <- c(1:(length(x)*length(y)))</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif">dd <- data.frame(x=x, y=y, z=z)</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif">#座標を指定</font></div><div class="gmail_default"><font face="verdana, sans-serif">coordinates(dd) <- ~x+y</font></div><div class="gmail_default"><font face="verdana, sans-serif">#projectionを指定</font></div><div class="gmail_default"><font face="verdana, sans-serif">proj4string(dd) <- CRS("+init=epsg:4326")</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif">#領域設定</font></div><div class="gmail_default"><font face="verdana, sans-serif">ext <- extent(dd)</font></div><div class="gmail_default"><font face="verdana, sans-serif">r <- raster(ext, ncol=20, nrow=20)</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif">#ラスタ変換</font></div><div class="gmail_default"><font face="verdana, sans-serif">dd.map <- rasterize(dd, r, field=dd$z)</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif">#左下隅を拡大（ラスタの中心に＋がない・・・）</font></div><div class="gmail_default"><font face="verdana, sans-serif">plot(dd.map, xlim=c(120,125), ylim=c(21,25))</font></div><div class="gmail_default"><font face="verdana, sans-serif">par(new=T)</font></div><div class="gmail_default"><font face="verdana, sans-serif">plot(dd, xlim=c(120,125), ylim=c(21,25))</font></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">------ここまで-----</div></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div>-- <br><div class="gmail_signature"><div dir="ltr">島田泰夫/一般財団法人日本気象協会/環境・エネルギー事業部/環境影響評価室/〒170-6055豊島区東池袋3-1-1サンシャイン60-55F<br>/TEL:03-5958-8160 FAX:03-5958-8157 /<a href="mailto:shimada@jwa.or.jp" style="color:rgb(17,85,204)" target="_blank">shimada@jwa.or.jp</a><br>/携帯:080-8018-1567<br></div></div>
</div>