Блог

Мне интересно, есть ли какой — либо способ построить 2D-изображение.plot() под драпировкой.plot (на пустом пространстве x~y) — или на любых графических объектах в этом пустом пространстве (например, точки, из которых оценивается TPS). По сути, объединение двух фигур, но на драпировке.перспектива сюжета. Пример кода ниже:

 #create data
a = rnorm(100, mean = 0, sd = 1)
b = rnorm(100, mean = 0, sd = 1)
x = cbind(a,b)
matrix = as.matrix(x, ncol=2)
c = rnorm(100, mean = 0.1, sd = 0.05)
c = c(c)

library(fields)
#fit thin plate spline
fit = Tps(matrix, c)
#generate surface
out = predictSurface(fit)
#2D plot of tps
image.plot(out, 
       xlim = c(-5,5),
       ylim = c(-5,5))
#3D plot of tps
drape.plot(out,zlab = "", xlab = "", ylab = "",
       zlim = c(-0.25,0.25), 
       horiz = T, box = T,
       xlim = c(-5,5),
       ylim = c(-5,5))