如何将John Snow数据集转换为地图坐标
在R的HistData数据包(https://r-forge.r-project.org/R/?group_id=574)中,我具有与约翰·斯诺(John Snow)1854年伦敦霍乱疫情地图有关的数据集。我相信它们在Walter Tobler的监督下经过精心数字化处理后,具有权威性。John Mackenzie在http://www1.udel.edu/johnmack/frec480/cholera/cholera2.html上描述了有关这些数据集的一些详细信息。 不幸的是,死亡,水泵和街道的坐标使用任意坐标系,而不是适用于R中其他GIS应用程序或地图软件(空间包,ggmap等)的地图坐标。 在http://freakonometrics.hypotheses.org/19213中, Arthur Charpentier将ggmap与来自http://www.rtwilson.com/downloads/SnowGIS_v2.zip的John Snow数据版本一起使用 。Cholera_Deaths.shp但是,该文件仅列出489人死亡,而不是我在中记录的578人HistData::Snow.deaths。 一种想法是找到均值与(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 …