r - 如何在 persp 中只有其他边界

Question

我想做以下事情

set.seed(1)
x <- seq(-10, 10, length= 600)
y <- x
f <- function(x, y) { r <- sqrt(x^2+y^2); 10 * sin(r)/r }
z <- outer(x, y, f)
persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = "lightblue")

但是网格太薄，所以边框占主导地位(您可以使用 border = NA 确认它确实是边框而不是照明。它变成蓝色):



解决这个问题的一种方法当然是使用不太精细的网格(例如，如果您将 length = 600 更改为 length= 50，它看起来非常令人愉快，实际上是 ?persp 中的示例)。但我想要与这个精细网格完全相同的形状和平滑度。我只是不想绘制所有的边框，例如可能只有它们的 1/5(或者我认为我可以自定义的一半)。

Answer 1

绘制平滑形状然后在其顶部绘制网格的一个问题是，您可以通过形状看到另一侧的网格。为了解决这个问题，您可以首先在白色对象的顶部绘制类(class)网格，这意味着您看不到网格的背面，将结果保存到文件中。

x <- seq(-10, 10, length=50)
y <- x
z <- outer(x, y, f)
png("top.png")
print(persp(x, y, z, theta = 30, phi = 30, expand = 0.5, border="black", col="white"))
dev.off()

然后，您可以绘制平滑的图像，然后是所有白色完全透明的网格。

x <- seq(-10, 10, length= 600)
y <- x
f <- function(x, y) { r <- sqrt(x^2+y^2); 10 * sin(r)/r }
z <- outer(x, y, f)
png("bottom.png")
print(persp(x, y, z, theta = 30, phi = 30, expand = 0.5, border="lightblue", col="lightblue"))
dev.off()

par(oma=c(0, 0, 0, 0), mar=c(0, 0, 0, 0))
library(png)
top.img <- readPNG("top.png")
top.img[,,4][top.img[,,1] + top.img[,,2] + top.img[,,3] > 2] <- 0

plot.new()
rasterImage(bottom.img, 0, 0, 1, 1)
rasterImage(top.img, 0, 0, 1, 1)