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

[Mizutani Takayuki]

島田さん

水谷です。
r <- raster(ext, ncol=20, nrow=20)

この部分のncolとnrowは19になると思います。
print(r)とするとresolutionで解像度が確認できます。
この場合、ポイントはラスタの中心ではなく、ラスタの隅となります。

--
水谷貴行

2017年10月15日 12:03 shimada yasu <yas.shimada35 ＠ gmail.com>:

> すみません、どなたかご教授下さい。
>
> 地図表示をしたく、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
>
> _______________________________________________
> OSGeoJapan-discuss mailing list
> OSGeoJapan-discuss ＠ lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>
-------------- next part --------------
HTMLの添付ファイルを保管しました...
URL: <http://lists.osgeo.org/pipermail/osgeojapan-discuss/attachments/20171015/410d9762/attachment-0001.html>