ich habe da mal wieder eine Frage. Auch Dank der Hilfe, welche mir im Forum zuteil wurde komme ich als Anfänger sehr gut voran. Meine heutige Frage gilt dem Thema 3D surface plots mittels plotly:
Ich habe mit den Testdaten, welche den reellen Werte im Gelände nicht ganz unähnlich sind, versucht die gegen die Tiefe in einem 3D Surface zusammenzufassen siehe https://plot.ly/r/3d-surface-plots/. Erstens ist mir bei den gezeigten Beispielen (Nr. 1 Basic 3D Surface Plot) nicht ganz klar welche Daten hierfür für die Berechnung verwendet werden.
Dann habe ich versucht mit meinen eigenen Daten die Beispiele 2 (2D Kernel Density Estimation) und 3 (Multiple Surfaces)nachzubauen. Vor allem letzteren ist für mich von enormes Interesse, erlaubt diese Art von Beispielen die räumliche Darstellung von horizontalen Schichtungen. Also im meinem Fall der Durchlässigkeitsbeiwert (sog. Kf Wert m/s) oder auch andere Daten. Diese Art von Darstellung wäre für die Auswertung solcher Daten eine enorme Hilfe!
Leider funktioniert das ganze nicht so wie in dem Beipiel angeführt. Beim Beispiel 2 (2D Kernel Density Estimation auf der Webseite),rechnet R das Gitter aus, gibt aber nur die Tiefe in der Legende an. Das Gitter ist leer.
Im Fall des Beispiels Nr. 3 (Multiple Surfaces), wird der Plot zwar dargestellt, nur die Darstellung sieht etwas schräg aus.
Meine Fragen wären:
1. Weiß jemand welche Daten hinter dem Plot in Beipiel 1 (Basic 3D Surface Plot) und in welcher Form diese für den Plot aufbereitet wurden.
2. Für den Nachbau des Beipiels Nr. 2 (2D Kernel Density Estimation) habe ich folgende eigene Daten verwendet:
Code: Alles auswählen
x <-runif(301, min=0.000001, max=0.00001)
y <-runif(301, min=0.000001, max=0.00001)
z <- seq(0,30,0.1)
z<- as.matrix(z)
p <- list( x = x , y = y, z = z)
p<- as.matrix(p)
p <- plot_ly(x = p$x, y = p$y, z = z) %>% add_surface()
p
Testdaten evtl. in einem falschen Format vor? Was könnte der Grund für das leere Gitter sein? Außerdem nützt mir
in der Legende die Höhe nichts. Ich möchte eher das die sog. Kf- Wert - x,y in der Legende angezeigt werden.
3. Beim Nachbau des Beispiels Nr. 3 (Multiple Surfaces):
Code: Alles auswählen
# erste 1m
z <-c(
c(8.753440e-06, 3.087675e-06, 0.0),
c(2.354936e-06, 4.172422e-06, 0.1),
c(4.809585e-06, 8.902687e-06, 0.2),
c(7.300328e-06, 1.498156e-06, 0.3),
c(2.918452e-06, 1.915362e-06, 0.4),
c(5.005054e-06, 2.203559e-06, 0.5),
c(4.474008e-06, 5.454518e-06, 0.6),
c(7.727593e-06, 5.699172e-06, 0.7),
c(1.679446e-06, 8.038964e-06, 0.8),
c(2.969207e-06, 1.764419e-06, 0.9),
c(8.288681e-06, 2.674241e-06, 1.0)
)
#plotting
dim(z) <- c(11,3)
z2 <- z + 1
z3 <- z - 1
p <- plot_ly(showscale = FALSE) %>%
add_surface(z = ~z) %>%
add_surface(z = ~z2, opacity = 0.98) %>%
add_surface(z = ~z3, opacity = 0.98)
p
Testdaten liegen könnte.
Zum Abschluss: Mein Ziel ist es Bohrlochdaten (2-4 Bohrungen) in der Tiefe 3 dimensional darzustellen. Ist die
Herangehensweise meinerseits mit Plotly dafür geeignet oder liegt hier ein grundsätzliches Missverständnis vor?
Vielen Dank für Eure Hilfe und die Zeit, welche Ihr auch dafür nehmt.
PS: Trotz mehrfacher Wiederholung gelang es mir nicht den Code mit dem zugehörigen Button einzufügen. Was vorher kein Problem war.