R语言绘制等值线和等高线

数据读取

data1 <- read.table(‘forams-Plio.txt’,sep = ‘\t’,header = TRUE)
head(data1)

str查看数据结构后，知道数据是242个，那么说明高是242个高，那咱们看242能被2和121进行相乘。涉及的矩阵为2成121。那么我想设计成30乘以8的矩阵，就可以删掉2个元素，暂时只是删掉一个。如下所示。

#设置高程值

str(data1)
zz=data1[-1,]
str(zz)
z1=zz$NumSp
z=matrix(z1, 10, 10)
str(z)
str(x)
#x坐标，每个网格为10米分辨率，方向由南向北
x <- 10*(1:nrow(z))
#y坐标，每个网格为10米分辨率，方向由东向西
y <- 10*(1:ncol(z))

二维可视化：栅格+等高线

par(mar=rep(0.5,4))
image(x, y, z, col=terrain.colors(100), axes=F)
contour(x, y, z, levels=seq(from=min(z), to=max(z), by=10),axes=T, add=T)
library(png)
logo <- readPNG(‘E:/logo.png’)

layout(z,heights=x,widths=y)