Я хотел бы построить несколько плоскостей, каждое из которых является неравенством. После того, как я намету все самолеты, я хотел бы объединить их и покрасить область внутри этих линий. Изображение рисует много трехмерных линий и окрашивает область внутри - это то, что я пытаюсь сделать.
Мои данные выглядят следующим образом:
df <- structure(list(z = c(0, 0.06518, 0.08429, -0.01659, 0, 0.06808,
0.12383, -1, -0.01662, 0.28782, 0, -0.09539, 0.04255, 0.09539,
-0.13361, -0.28782, -0.14468, -0.19239, 0.10642), x = c(1, 0.02197,
0.03503, -0.02494, 0, 0.04138, 0.17992, 0, -0.02482, 0.1122,
0, 0.01511, 0.0011, -0.01511, -0.06699, -0.1122, -0.06876, 0.12078,
0.10201), y = c(0, 0.08735, 0.09927, 0.03876, -1, 0.22114, -0.00152,
0, 0.03811, -0.07335, 0, -0.03025, 0.07681, 0.03025, -0.23922,
0.07335, -0.25362, -0.09879, 0.05804), value = c(5801L, 135L,
162L, 109L, 4250L, 655L, 983L, 4500L, 108L, 1594L, 4400L, 540L,
147L, 323L, 899L, 1023L, 938L, 1627L, 327L)), .Names = c("z",
"x", "y", "value"), class = "data.frame", row.names = c(NA, -19L
))
Каждая строка представляет собой уравнение вида: z + x + y < value
. x
- горизонтальное значение, y
- вертикальное значение, а z
- глубина. z может быть решено быть: -x - y + value > z.
Пределы системы координат:
x <- z <- seq(-6000, 6000, by = 1)
y <- seq(-4000, 4000, by = 1)
Итак, из каждой строки я хотел бы нарисовать самолет. Затем я хотел бы объединить все эти плоскости и заполнить значения внутри строк. Результат должен выглядеть как многосторонний неравный кубик. Или уродливый бриллиант.
Я играл с функциями rgl
и persp
, но я не уверен, с чего начать. Я открыт для других рекомендаций по программному обеспечению.
Опираясь на один из примеров из persp3d
:
x <- seq(-6000, 6000, by = 1)
z <- x
y <- seq(-4000, 4000, by = 1)
f <- function(x, y) <- { r <- -x - y + value > z } # stuck here, can you handle an inequality here?
z <- outer(x, y, f)
open3d()
bg3d("white")
material3d(col = "black")
persp3d(x, y, z, col = "lightblue",
xlab = "X", ylab = "Y", zlab = "z")
Я понимаю, что это довольно большие пределы. Если это помогает уменьшить их, не стесняйтесь, или увеличьте sequence(..., by = )
.