John Snowデータセットをマップ座標に変換する方法
RのHistDataパッケージ(https://r-forge.r-project.org/R/?group_id=574)には、1854年のロンドンでのコレラの発生に関するJohn Snowのマップに関連するデータセットがあります。ウォルター・トブラーの監督の下で慎重にデジタル化された権威です。これらのデータセットの詳細については、John Mackenzieがhttp://www1.udel.edu/johnmack/frec480/cholera/cholera2.htmlで説明しています。 残念ながら、死亡、ポンプ、および道路の座標は、他のGISアプリケーションまたはRのマッピングソフトウェア(空間パッケージ、ggmapなど)に適したマップ座標ではなく、任意の座標系を使用します でhttp://freakonometrics.hypotheses.org/19213アーサー・シャルパンティエはからスノーデータのバージョンでggmapを使用しています http://www.rtwilson.com/downloads/SnowGIS_v2.zip。Cholera_Deaths.shpただし、このファイルには、私が記録した578人ではなく、489人の死亡のみがリストされていHistData::Snow.deathsます。 1つのアイデアは、平均と(x、y)座標の標準偏差との関係を見つけて線形に再スケーリングすることですが、おそらくもっと良い方法がありますか? ここに私が今まで試したものがあります > data(Snow.deaths, package="HistData") > D <- Snow.deaths[,2:3] > colMeans(D) x y 13.03312 11.69721 > var(D) x y x 3.8150987 0.3802654 y 0.3802654 2.7213828 Cholera_deathsファイルを読む > folder <- "C:/Dropbox/R/data/Snow/SnowGIS_v2/SnowGIS" > library(maptools) > deaths <- readShapePoints(file.path(folder, "Cholera_Deaths")) > head(deaths@coords) coords.x1 coords.x2 0 529308.7 181031.4 …