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

[shimada yasu]

水谷様

ありがとうございました。

> ​この部分の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
-------------- next part --------------
HTMLの添付ファイルを保管しました...
URL: <http://lists.osgeo.org/pipermail/osgeojapan-discuss/attachments/20171015/a9f347ec/attachment-0001.html>
-------------- next part --------------
テキスト形式以外の添付ファイルを保管しました...
ファイル名: clipbord_R_res19.png
型:         image/png
サイズ:     18432 バイト
説明:       無し
URL:        <http://lists.osgeo.org/pipermail/osgeojapan-discuss/attachments/20171015/a9f347ec/attachment-0001.png>