<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">この部分の</span>ncolとnrowは19になると思います。</div><div>print(r)とするとresolutionで解像度が確認できます。</div><div>この場合、ポイントはラスタの中心ではなく、ラスタの隅となります。</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: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">地図表示をしたく、<wbr>Rで格子状のポイントデータをラスタ変換すると</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）以下のスクリプトが間違っている？</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">（3）そのほか・・・</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">RとGISも中途半端な理解なので、<wbr>根本的に間違っているかも知れませんが、</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"><div class="gmail_default"><font face="verdana, sans-serif">library(maptools)</font></div><div class="gmail_default"><font face="verdana, sans-serif">library(rgdal)</font></div><div class="gmail_default"><font face="verdana, sans-serif">library(raster)</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif">rm(list=ls(all=TRUE))</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif">#ダミーデータ（データフレーム）作成</font></div><div class="gmail_default"><font face="verdana, sans-serif">x <- rep(seq(121,140,1),10)</font></div><div class="gmail_default"><font face="verdana, sans-serif">y <- numeric(0);for ( i in 21:40 ){ y <- c(y, rep(i,20))}</font></div><div class="gmail_default"><font face="verdana, sans-serif">z <- c(1:(length(x)*length(y)))</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif">dd <- data.frame(x=x, y=y, z=z)</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif">#座標を指定</font></div><div class="gmail_default"><font face="verdana, sans-serif">coordinates(dd) <- ~x+y</font></div><div class="gmail_default"><font face="verdana, sans-serif">#projectionを指定</font></div><div class="gmail_default"><font face="verdana, sans-serif">proj4string(dd) <- CRS("+init=epsg:4326")</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif">#領域設定</font></div><div class="gmail_default"><font face="verdana, sans-serif">ext <- extent(dd)</font></div><div class="gmail_default"><font face="verdana, sans-serif">r <- raster(ext, ncol=20, nrow=20)</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif">#ラスタ変換</font></div><div class="gmail_default"><font face="verdana, sans-serif">dd.map <- rasterize(dd, r, field=dd$z)</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="verdana, sans-serif">#左下隅を拡大（ラスタの中心に＋がない・・・）</font></div><div class="gmail_default"><font face="verdana, sans-serif">plot(dd.map, xlim=c(120,125), ylim=c(21,25))</font></div><div class="gmail_default"><font face="verdana, sans-serif">par(new=T)</font></div><div class="gmail_default"><font face="verdana, sans-serif">plot(dd, xlim=c(120,125), ylim=c(21,25))</font></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">------ここまで-----</div></div><span class="HOEnZb"><font color="#888888"><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br></div>-- <br><div class="m_2560801062928531410gmail_signature"><div dir="ltr">島田泰夫/一般財団法人日本気象協会/環境・エネルギー事業部/<wbr>環境影響評価室/〒170-6055豊島区東池袋3-1-<wbr>1サンシャイン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></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>