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

[Mizutani Takayuki]

水谷です。

ポイントのプロット部分を以下のようにadd=Tで重ねると大丈夫でした。
（par(new=T)で重ねると上手くいかないのかもしれません。わかりません）

r <- raster(ext, ncol=19, nrow=19)dd.map <- rasterize(dd, r, field=dd$z)
plot(dd.map, xlim=c(120,125), ylim=c(21,25))
plot(dd,add=T)

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

> 水谷様
>
> ありがとうございました。
>
> > ​この部分のncolとnrowは19になると思います。
> > print(r)とするとresolutionで解像度が確認できます。
>
> 早速、やってみました。
> 解像度は1となりましたが、plotしたところ、中心にはなりせんでした・・・
>
> -----
> r <- raster(ext, ncol=19, nrow=19)
> -----
> class       : RasterLayer
> dimensions  : 19, 19, 361  (nrow, ncol, ncell)
> resolution  : 1, 1  (x, y)
> extent      : 121, 140, 21, 40  (xmin, xmax, ymin, ymax)
> coord. ref. : NA
> -----
>
> 2017年10月15日 15:30 Mizutani Takayuki <mizutani.takayuki ＠ gmail.com>:
>
> 島田さん
>>
>> 水谷です。
>> 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-
>>> <https://maps.google.com/?q=%E8%B1%8A%E5%B3%B6%E5%8C%BA%E6%9D%B1%E6%B1%A0%E8%A2%8B3-1-1&entry=gmail&source=g>
>>> 1
>>> <https://maps.google.com/?q=%E8%B1%8A%E5%B3%B6%E5%8C%BA%E6%9D%B1%E6%B1%A0%E8%A2%8B3-1-1&entry=gmail&source=g>
>>> サンシャイン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
>>>
>>
>>
>
>
> --
> 島田泰夫/一般財団法人日本気象協会/環境・エネルギー事業部/環境影響評価室/〒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/1a809c20/attachment.html>