.csv 파일의 3 개 열로 구성된 데이터 세트가 있습니다. 첫 번째 두 열은지도 좌표이며 세 번째 열은 해당지도 좌표에서 시추공에서 발견 된 아연의 백분율입니다. 거리에 따른 Zn 농도 변화를 보여주기 위해 등고선지도를 만들고 싶습니다. 모든 코드 예제는 행렬 형태로 데이터를 사용할 수 있지만 목록은 목록에 있습니다. 나는 밑에있는 그림을 그리는 여러 가지 다른 방법을 시도했다. 대부분의 메서드는 내 데이터의 레이아웃과 관련이있는 "개체 x를 찾을 수 없음"의 줄에 따라 오류 메시지를 표시합니다. 누구든지이 작업을 수행하는 방법을 알고 있습니까? 비슷한 데이터 세트를 아래에 추가했습니다. 사전에 도움을 주셔서 감사합니다. 홀리목록의 데이터를 사용하여 채워진 등고선 플롯 만들기
데이터 세트 :
Statsrep <- structure(list(X = c(156000L, 156010L, 156421L, 156450L, 156500L, 156700L, 158420L, 158646L, 158970L, 159050L, 159050L, 159130L, 159155L), Y = c(143630, 143980, 147260, 145000, 146000, 142800, 146700, 145207, 147170, 145200, 144800, 147815, 145890), Zn = c(2, 8, 4, 0, 3, 0, 2, 7, 12, 0, 4, 19, 0)), .Names = c("X", "Y", "Zn"), row.names = c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L), class = "data.frame")
Statsrep
코드 :
library(ggplot2)
Grade <- read.csv(file="filename.csv", header=TRUE, sep=",")
ggplot(Grade, aes(x$x="X", y$y="Y", z$z="Zn")) +
stat_contour()
library(lattice)
Grade <- read.csv(file="filename.csv", header=TRUE, sep=",")
levelplot(Grade ~x*y, data = Zn,
xlab = "Eastings", ylab = "Northings",
col.regions = terrain.colours)
Grade <- read.csv(file="filename.csv", header=TRUE, sep=",")
x$x <- X
y$y <- Y
z$z <- Zn
filled.contour(x$x, y$y, z$z, color = terrain.colours,
xlab = "Eastings", ylab = "Northings"),
plot.axes = {axis(1, seq(156000, 165000, by=1000)); axis(2, seq(142000, 150000, by=1000))},
key.title = title(main="Zn content\n(percent)"),
key.axes= axis(4, seq(0, 20, by = 2)))
아, 그래서 저는 잘못된 기능을 사용하고있었습니다. 이 방법이 효과적 일뿐 아니라 선들 사이에 색상을 채울 수 있습니까? –
시도해보세요.'g + stat_density2d (aes (채우기 = .. 밀도 ..), contour = FALSE, geom = "타일")' – csgillespie
고마워요! 두 개가 결합되면 그것은 내가 찾고있는 음모의 유형입니다. 그러나 밀도에 대한 열쇠는 아연의 가치를 제공하지 않으며, 그것은 e-08의 가치를 제공합니다. –