gpt4 book ai didi

r - 在 R map 中包含比例和坐标

转载 作者:行者123 更新时间:2023-12-04 11:07:59 24 4
gpt4 key购买 nike

我正在使用以下代码在 R 中绘制 map :

library(maps)
library(mapdata)

map('worldHires',c('UK'),
xlim=c(-10,-1),
ylim=c(56,59.5))
points(-5.65,56.7233,col=2,pch=18)

导致

enter image description here

如何在 map 上添加比例尺以及在 map 周围显示纬度和经度的方框?类似于:

enter image description here

最佳答案

maps 包还有一个函数 map.axesmap.scales

library(maps)
library(mapdata)

map('worldHires',c('UK'), xlim=c(-10,-1),
ylim=c(56,59.5))

points(-5.65,56.7233,col=2,pch=18)
map.axes()
map.scale()

给你

enter image description here

编辑1

如果您查看 map.axes 的代码,您会发现它只是调用了 axisbox。一种方法是手动执行此操作并调整标签。

library(maps)
library(mapdata)

m <- map('worldHires',c('UK'), xlim=c(-10,-1),
ylim=c(56,59.5), mar=c(4, 4, 4, 4))

points(-5.65,56.7233,col=2,pch=18)

xat <- pretty(m$range[1:2])
xlab <- paste0(xat, " E")

yat <- pretty(m$range[3:4])
ylab <- paste0(yat, " N")

box()
axis(1, at=xat, labels=xlab)
axis(2, las=TRUE, at=yat, labels=ylab)
axis(3, at=xat, labels=xlab)
axis(4, las=TRUE, at=yat, labels=ylab)
map.scale(y=59)

这将为您提供具有调整标签的相同图片 enter image description here

编辑2:

或者按照@Pascal 的建议使用raster::degreeLabels* 函数。

library(maps)
library(mapdata)
library(raster)

m <- map('worldHires',c('UK'), xlim=c(-10,-1),
ylim=c(56,59.5), mar=c(4, 4, 4, 4))

points(-5.65,56.7233,col=2,pch=18)

xat <- pretty(m$range[1:2])
xlab <- parse(text=degreeLabelsEW(xat))

yat <- pretty(m$range[3:4])
ylab <- parse(text=degreeLabelsNS(yat))


box()
axis(1, at=xat, labels=xlab)
axis(2, las=TRUE, at=yat, labels=ylab)
axis(3, at=xat, labels=xlab)
axis(4, las=TRUE, at=yat, labels=ylab)

map.scale(y=59)

enter image description here

关于r - 在 R map 中包含比例和坐标,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/26114962/

24 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com