微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

如何在我的地区创建此特定地图

如何解决如何在我的地区创建此特定地图

我使用以下代码来创建特定地图。不幸的是,由于多次尝试,我无法制作所需的地图。因此,我要求帮助我创建这张地图。 The required map is attached herethe following code gives the map is here 。需要进行哪些更改才能制作所需的地图。请注意,我在附加代码中没有错误,我想修改代码以提供所需的地图。

library(tidyverse)
 library(raster)
    
#obtain map from GADM.org
pak <- getData('GADM',country='Pakistan',level=1)

map <- ggplot(pak) + geom_sf()

#Object for ggplot
#Need a name for joins
pak$name = pak$NAME_1
pak_ggplot <- map_data(pak)
#map_data transform name in region

#I create a random variable for color in map
zone <- data.frame(region = unique(pak$NAME_1))
zone$other_name = runif(nrow(zone),1,10)

#Join ggplot object with variable for color
pak_ggplot <- left_join(pak_ggplot,zone)

dist <- c('Astore','Bunji','Gupis','Chilas','Gilgit','Skardu')
lat <- c(35.357,35.6431,36.2274,35.4222,35.8819,35.3247)
long <- c(74.8624,74.6342,73.4421,74.0946,74.4643,75.551)

points <- cbind.data.frame(dist,lat,long)


mapa <- ggplot(data = pak_ggplot,aes(x = long,y = lat,group = group)) +
  geom_polygon(aes(fill = other_name)) +
  scale_fill_continuous(low="white",high="red")  +
  labs(title = "Pakistan Map")

mapa + geom_point(data=points,aes(x=long,y=lat,group=1)) +
  geom_text(data=points,label=dist,group=1),size=2)

pak2 <- pak_ggplot %>% filter(region=='northern Areas')

mapa2 <- ggplot(data = pak2,group = 1)) +
  geom_polygon() +
  scale_fill_continuous(low="white",high="red")  +
  labs(title = "Pakistan Map")

mapa2 + geom_point(data=points,size=2)

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。