[OSGeoJapan-discuss] 【R：教えて下さい】ポイントをラスタ変換すると中心がずれる

[shimada yasu]

すみません、どなたかご教授下さい。

地図表示をしたく、Rで格子状のポイントデータをラスタ変換すると
ラスタの中心にポイントが来ないのですが、、、これは

（1）以下のスクリプトが間違っている？
（2）そもそも経度の距離は緯度に依存するから、平均的なラスタ解像度とは、当然一致しない？
（3）そのほか・・・

ご助言を頂けると幸いです。

RとGISも中途半端な理解なので、根本的に間違っているかも知れませんが、
数日考えてもどうにも解決できず投稿した次第です。

----ここから----
library(maptools)
library(rgdal)
library(raster)

rm(list=ls(all=TRUE))

#ダミーデータ（データフレーム）作成
x <- rep(seq(121,140,1),10)
y <- numeric(0);for ( i in 21:40 ){ y <- c(y, rep(i,20))}
z <- c(1:(length(x)*length(y)))

dd <- data.frame(x=x, y=y, z=z)

#座標を指定
coordinates(dd) <- ~x+y
#projectionを指定
proj4string(dd) <- CRS("+init=epsg:4326")

#領域設定
ext <- extent(dd)
r <- raster(ext, ncol=20, nrow=20)

#ラスタ変換
dd.map <- rasterize(dd, r, field=dd$z)

#左下隅を拡大（ラスタの中心に＋がない・・・）
plot(dd.map, xlim=c(120,125), ylim=c(21,25))
par(new=T)
plot(dd, xlim=c(120,125), ylim=c(21,25))
------ここまで-----

-- 
島田泰夫/一般財団法人日本気象協会/環境・エネルギー事業部/環境影響評価室/〒170-6055豊島区東池袋3-1-1サンシャイン60-55F
/TEL:03-5958-8160 FAX:03-5958-8157 /shimada ＠ jwa.or.jp
/携帯:080-8018-1567
-------------- next part --------------
HTMLの添付ファイルを保管しました...
URL: <http://lists.osgeo.org/pipermail/osgeojapan-discuss/attachments/20171015/15c0bccd/attachment.html>