<div dir="ltr">水谷です。<div><br></div><div>この場合のraster関数は、第1引数を四隅の範囲として、</div><div>その範囲をncol、nrowの数で分割した解像度のラスタを作成します。</div><div>rasterize関数は、raster関数で作成した空のラスタにポイントの値を</div><div>入れ込むというイメージになります。</div><div>以下のコードではいかがでしょうか？</div><div><br></div><div><div>e<-extent(c(120.5, 140.5, 20.5, 40.5))</div><div>r <- raster(e, ncol=20, nrow=20)</div><div>dd.map <- rasterize(dd, r, field=dd$z)</div><div>plot(dd.map, xlim=c(120,125), ylim=c(21,25))</div><div>plot(dd,add=T)</div></div><div><br></div><div>水谷</div></div><div class="gmail_extra"><br><div class="gmail_quote">2017年10月15日 17:56 shimada yasu <span dir="ltr"><<a href="mailto:yas.shimada35@gmail.com" target="_blank">yas.shimada35@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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">日曜日にも係わらずありがとうございます。</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">（1）ポイントは、<wbr>当初ダミーで設定したとおりの場所に表示されまし（添付）</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">（2）ただ、点はラスタの中心でなく、四隅に位置する・・ので、<wbr>ラスタがどの点を代表しているのかが分からない。。。</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">　私が欲しかったのは、下記のような「<wbr>ラスタのど真ん中に点がある」図なのですが。。。（<wbr>考え方が間違っているかな。。）</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">|              |</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><div class="gmail_default">|      +      |</div><div><div class="gmail_default">|              |</div></div><div><div class="gmail_default">------------     <br></div></div></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"><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_extra"><div class="gmail_quote">2017年10月15日 17:18 Mizutani Takayuki <span dir="ltr"><<a href="mailto:mizutani.takayuki@gmail.com" target="_blank">mizutani.takayuki@gmail.com</a>></span>:<div><div class="h5"><br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><span style="font-size:14px">水谷です。</span><div style="font-size:14px"><br></div><div style="font-size:14px">ポイントのプロット部分を以下のようにadd=Tで重ねると大丈<wbr>夫でした。</div><div style="font-size:14px">（par(new=T)で重ねると上手くいかないのかもしれませ<wbr>ん。わかりません）</div><div style="font-size:14px"><br></div><div style="font-size:14px"><span class="m_-7304054549269228701gmail-m_7637094393477911161m_2430877528391582202gmail-im">r <- raster(ext, ncol=19, nrow=19)</span><span class="m_-7304054549269228701gmail-m_7637094393477911161m_2430877528391582202gmail-im">dd.map <- rasterize(dd, r, field=dd$z)</span><div>plot(dd.map, xlim=c(120,125), ylim=c(21,25))</div><div>plot(dd,add=T)</div></div></div><div class="gmail_extra"><br><div class="gmail_quote">2017年10月15日 16:52 shimada yasu <span dir="ltr"><<a href="mailto:yas.shimada35@gmail.com" target="_blank">yas.shimada35@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><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">ありがとうございました。</div><div class="gmail_default"><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">> ​この部分のncolとnrowは19になると思います。</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">> print(r)とするとresolutionで解像度が確認で<wbr>きます。</div><div style="font-family:verdana,sans-serif;font-size:small"><br></div><div style="font-family:verdana,sans-serif;font-size:small">早速、やってみました。</div><div style="font-family:verdana,sans-serif;font-size:small">解像度は1となりましたが、plotしたところ、中心にはなりせ<wbr>んでした・・・</div><div style="font-family:verdana,sans-serif;font-size:small"><br></div><div style="font-family:verdana,sans-serif;font-size:small"><div>-----<br></div><div style="font-family:arial,sans-serif"></div><div>r <- raster(ext, ncol=19, nrow=19)</div><div>-----<br></div></div><div><div><font face="verdana, sans-serif">class       : RasterLayer </font></div><div><font face="verdana, sans-serif">dimensions  : 19, 19, 361  (nrow, ncol, ncell)</font></div><div><font face="verdana, sans-serif">resolution  : 1, 1  (x, y)</font></div><div><font face="verdana, sans-serif">extent      : 121, 140, 21, 40  (xmin, xmax, ymin, ymax)</font></div><div><font face="verdana, sans-serif">coord. ref. : NA </font></div></div><div style="font-family:verdana,sans-serif;font-size:small">-----</div></div><div class="gmail_extra"><br><div class="gmail_quote">2017年10月15日 15:30 Mizutani Takayuki <span dir="ltr"><<a href="mailto:mizutani.takayuki@gmail.com" target="_blank">mizutani.takayuki@gmail.com</a>></span>:<div><div class="m_-7304054549269228701gmail-m_7637094393477911161m_2430877528391582202h5"><br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">島田さん<div><br></div><div>水谷です。</div><div><span style="font-family:verdana,sans-serif;font-size:14px">r <- raster(ext, ncol=20, nrow=20)</span><br></div><div><span style="font-family:verdana,sans-serif;font-size:14px"><br></span></div><div><span style="font-family:verdana,sans-serif;font-size:14px"><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small;display:inline">​​</div>この部分の</span>ncolとnrowは19になると思います。</div><div>print(r)とするとresolutionで解像度が確認で<wbr>きます。</div><div>この場合、ポイントはラスタの中心ではなく、ラスタの隅となりま<wbr>す。</div><div><br></div><div>--</div><div>水谷貴行</div></div><div class="gmail_extra"><br><div class="gmail_quote">2017年10月15日 12:03 shimada yasu <span dir="ltr"><<a href="mailto:yas.shimada35@gmail.com" target="_blank">yas.shimada35@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div style="font-family:verdana,sans-serif;font-size:small">すみません、どなたかご教授下さい。</div><div style="font-family:verdana,sans-serif;font-size:small"><br></div><div style="font-family:verdana,sans-serif;font-size:small">地図表示をしたく、Rで格子状のポイントデータをラスタ変換する<wbr>と</div><div style="font-family:verdana,sans-serif;font-size:small">ラスタの中心にポイントが来ないのですが、、、これは</div><div style="font-family:verdana,sans-serif;font-size:small"><br></div><div style="font-family:verdana,sans-serif;font-size:small">（1）以下のスクリプトが間違っている？</div><div style="font-family:verdana,sans-serif;font-size:small">（2）そもそも経度の距離は緯度に依存するから、平均的なラスタ<wbr>解像度とは、当然一致しない？</div><div style="font-family:verdana,sans-serif;font-size:small">（3）そのほか・・・</div><div style="font-family:verdana,sans-serif;font-size:small"><br></div><div style="font-family:verdana,sans-serif;font-size:small">ご助言を頂けると幸いです。</div><div style="font-family:verdana,sans-serif;font-size:small"><br></div><div style="font-family:verdana,sans-serif;font-size:small">RとGISも中途半端な理解なので、根本的に間違っているかも知<wbr>れませんが、</div><div style="font-family:verdana,sans-serif;font-size:small">数日考えてもどうにも解決できず投稿した次第です。</div><div style="font-family:verdana,sans-serif;font-size:small"><br></div><div style="font-family:verdana,sans-serif;font-size:small">----ここから----</div><div><div><font face="verdana, sans-serif">library(maptools)</font></div><div><font face="verdana, sans-serif">library(rgdal)</font></div><div><font face="verdana, sans-serif">library(raster)</font></div><div><font face="verdana, sans-serif"><br></font></div><div><font face="verdana, sans-serif">rm(list=ls(all=TRUE))</font></div><div><font face="verdana, sans-serif"><br></font></div><div><font face="verdana, sans-serif">#ダミーデータ（データフレーム）作成</font></div><div><font face="verdana, sans-serif">x <- rep(seq(121,140,1),10)</font></div><div><font face="verdana, sans-serif">y <- numeric(0);for ( i in 21:40 ){ y <- c(y, rep(i,20))}</font></div><div><font face="verdana, sans-serif">z <- c(1:(length(x)*length(y)))</font></div><div><font face="verdana, sans-serif"><br></font></div><div><font face="verdana, sans-serif">dd <- data.frame(x=x, y=y, z=z)</font></div><div><font face="verdana, sans-serif"><br></font></div><div><font face="verdana, sans-serif">#座標を指定</font></div><div><font face="verdana, sans-serif">coordinates(dd) <- ~x+y</font></div><div><font face="verdana, sans-serif">#projectionを指定</font></div><div><font face="verdana, sans-serif">proj4string(dd) <- CRS("+init=epsg:4326")</font></div><div><font face="verdana, sans-serif"><br></font></div><div><font face="verdana, sans-serif">#領域設定</font></div><div><font face="verdana, sans-serif">ext <- extent(dd)</font></div><div><font face="verdana, sans-serif">r <- raster(ext, ncol=20, nrow=20)</font></div><div><font face="verdana, sans-serif"><br></font></div><div><font face="verdana, sans-serif">#ラスタ変換</font></div><div><font face="verdana, sans-serif">dd.map <- rasterize(dd, r, field=dd$z)</font></div><div><font face="verdana, sans-serif"><br></font></div><div><font face="verdana, sans-serif">#左下隅を拡大（ラスタの中心に＋がない・・・）</font></div><div><font face="verdana, sans-serif">plot(dd.map, xlim=c(120,125), ylim=c(21,25))</font></div><div><font face="verdana, sans-serif">par(new=T)</font></div><div><font face="verdana, sans-serif">plot(dd, xlim=c(120,125), ylim=c(21,25))</font></div><div style="font-family:verdana,sans-serif;font-size:small">------ここまで-----</div></div><span class="m_-7304054549269228701gmail-m_7637094393477911161m_2430877528391582202m_-1592478299379076876gmail-HOEnZb"><font color="#888888"><span class="m_-7304054549269228701gmail-m_7637094393477911161m_2430877528391582202m_-1592478299379076876gmail-m_-7817052739310242743HOEnZb"><font color="#888888"><div style="font-family:verdana,sans-serif;font-size:small"><br></div>-- <br><div class="m_-7304054549269228701gmail-m_7637094393477911161m_2430877528391582202m_-1592478299379076876gmail-m_-7817052739310242743m_2560801062928531410gmail_signature"><div dir="ltr">島田泰夫/一般財団法人日本気象協会/環境・エネルギー事業部/<wbr>環境影響評価室/〒170-6055<a href="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" target="_blank">豊島区東池袋3-1-</a><a href="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" target="_blank">1</a>サン<wbr>シャイン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>
</font></span></font></span></div><span class="m_-7304054549269228701gmail-m_7637094393477911161m_2430877528391582202m_-1592478299379076876gmail-HOEnZb"><font color="#888888">
<br>______________________________<wbr>_________________<br>
OSGeoJapan-discuss mailing list<br>
<a href="mailto:OSGeoJapan-discuss@lists.osgeo.org" target="_blank">OSGeoJapan-discuss@lists.osgeo<wbr>.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailma<wbr>n/listinfo/osgeojapan-discuss</a><br></font></span></blockquote></div><br></div>
</blockquote></div></div></div><div><div class="m_-7304054549269228701gmail-m_7637094393477911161m_2430877528391582202h5"><br><br clear="all"><span class="m_-7304054549269228701gmail-m_7637094393477911161HOEnZb"><font color="#888888"><div><br></div>-- <br><div class="m_-7304054549269228701gmail-m_7637094393477911161m_2430877528391582202m_-1592478299379076876gmail_signature"><div dir="ltr">島田泰夫/一般財団法人日本気象協会/環境・エネルギー事業部/<wbr>環境影響評価室/〒170-6055<a href="https://maps.google.com/?q=%E8%B1%8A%E5%B3%B6%E5%8C%BA%E6%9D%B1%E6%B1%A0%E8%A2%8B3-1&entry=gmail&source=g" target="_blank">豊島区東池袋3-1</a>-1サン<wbr>シャイン60-55F<br>/TEL:03-5958-8160 FAX:03-5958-8<a href="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" target="_blank">1</a>57 /<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>
</font></span></div></div></div></div><span class="m_-7304054549269228701gmail-m_7637094393477911161HOEnZb"><font color="#888888">
<br>______________________________<wbr>_________________<br>
OSGeoJapan-discuss mailing list<br>
<a href="mailto:OSGeoJapan-discuss@lists.osgeo.org" target="_blank">OSGeoJapan-discuss@lists.osgeo<wbr>.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailma<wbr>n/listinfo/osgeojapan-discuss</a><br></font></span></blockquote></div><br></div>
</blockquote></div></div></div><div><div class="h5"><br><br clear="all"><div><br></div>-- <br><div class="m_-7304054549269228701gmail-m_7637094393477911161gmail_signature"><div dir="ltr">島田泰夫/一般財団法人日本気象協会/環境・エネルギー事業部/<wbr>環境影響評価室/〒170-6055豊島区東池袋3-1-1サン<wbr>シャイン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></div></div></div>
<br>______________________________<wbr>_________________<br>
OSGeoJapan-discuss mailing list<br>
<a href="mailto:OSGeoJapan-discuss@lists.osgeo.org">OSGeoJapan-discuss@lists.<wbr>osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss" rel="noreferrer" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/osgeojapan-<wbr>discuss</a><br></blockquote></div><br></div>