Question

我是R和GIS菜鸟，非常感谢您对我的问题的帮助。我需要在R中创建一个热图，并且具有一个数据集，该数据集具有华盛顿特区都市区的某些属性，经度和纬度列。我使用ggmap导入华盛顿特区的地图，然后在其上添加了我的数据集的图。地图看起来很棒。现在我需要创建此绘图的shapefile并以ArcGIS格式导出它。我环顾四周，但我找不到教程或帮助页面的确切内容。这是我使用的代码：

library(ggmap)
library(ggplot2)

myLocation <- "Washington DC"
#Getting a map of Washington DC metro area from google API
myMap <- get_map(location=myLocation, source="google", maptype="terrain", crop="FALSE")
ggmap(myMap) #maps the Washington DC map
zoom=5

#This plot the heatmap
ggmap(myMap) + 
  geom_density2d(data = mydata, 
                 aes(x = longitude, y = latitude), size = 0.3) + 
  stat_density2d(data = mydata, 
                 aes(x = longitude, y = latitude, fill = ..level.., 
                     alpha = ..level..), size = 0.01, bins = 16, geom = "polygon") + 
  scale_fill_gradient(low = "green", high = "red") + 
  scale_alpha(range = c(0, 0.3), guide = FALSE)

这个命令对我有用。我现在迷失了如何基于DC地图和我绘制的点以及如何导出它来创建shapefile。我探索了图书馆的shapefile，却找不到我想要的东西。

这就是我的数据集“mydata”看起来像

ID_Var         latitude      longitude
223573820000001 38.8529910  -77.3640670
435153420000004 39.0964730  -76.5175040
579885949000002 38.9926000  -76.8358700
218543557000001 39.3958180  -76.4634520
578842012000004 38.9113450  -76.7900130
214727140000003 39.4023900  -76.4168250
214726688000001 39.4455920  -76.6112660
214760285000002 39.2898340  -76.6021910
214738051000001 39.2529070  -77.0197160
214727140000002 39.4023900  -76.4168250

使用long和lat在R中创建shapefile并将其导出

0 个答案: