我试图用来levelplot绘制一个简单的数字高程模型(DEM).
这是我的代码:
r1 = raster("ned10dem.tif")
e = extent(460000,480000,4555000,4567500)
rr1 = crop(r1,e)
p = levelplot(rr1, scales=list(x=list(at=seq(450000,480000,4000))),
margin=F, cuts=200,
col.regions = terrain.colors(350,alpha=1),
colorkey=list(space="bottom"),
xlab="Easting(m)", ylab="Northing(m)")
plot(p)
Run Code Online (Sandbox Code Playgroud)
情节最终看起来像这样:

我无法弄清楚的是如何增加colorkey和x轴之间的空间,使得colorkey不会覆盖x轴标签.
我正在使用 R markdown 并使用 Rstudio 来针织。
我有以下 R markdown 文件:
---
title: "Untitled"
author: ""
date: ""
output: html_document
---
```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```
```{r}
#Load libraries
library(knitr)
library(ggplot2)
library(reshape2)
library(chron)
library(gridExtra)
#Get current directory
directory = getwd()
setwd(directory)
readname = function()#prompt user for the name and number of the site
{
n <- readline(prompt="Site name and number: ")
return(n)
}
prjName=readname()
print(prjName)
prjName2="TEST NAME"
prjName2
Run Code Online (Sandbox Code Playgroud)
该代码要求用户提供站点名称和编号。我可以在 R markdown 中运行该块,并通过控制台获取用户输入并打印出来。它也作为字符串存储在全局环境中。然而,当我编织文档时,我得到以下结果:
看来,在编织文档时,我无法打印用户输入的字符串......
任何对此的帮助将不胜感激。