Gast
#6063998
Liebe Gemeinde,
bislang habe ich ortsaufgelöste Messdaten immer mit Gnuplot als Heatmap
geplottet.
Das Datenformat sieht folgendermaßen aus:
[code]
[x] [y] [z]
0 0 8
1 0 5
2 0 3
0 1 4
1 1 5
2 1 4
0 2 8
1 2 7
2 2 6
[\code]
[x] und [y] sind dabei die Ortskoordinaten und [z] der Messwert.
In Gnuplot kann man ja mit set cbrange[a:b] den Datenbereich für die
Palette angeben und dann einfach mit
[code]
set view map
splot 'data.dat' using 1:2:3 with image palette
[\code]
die Heatmap erstellen.
Nun möchte ich das ganze vollständig in der Sprache "R" implementieren,
da ich dort auch noch weitere Prozessierung durchführe.
Leider scheitere ich kläglich bei der Verwendung der image-Funktion:
[code]
image(x, y, z, zlim, xlim, ylim,
col = hcl.colors(12, "YlOrRd", rev = TRUE),
add = FALSE, xaxs = "i", yaxs = "i", xlab, ylab,
breaks, oldstyle = FALSE, useRaster, …)
[\code]
Ich erzeuge aus der eingelesenen Tabelle die Variablen x,y und z:
[code]
ip <- read.table("1.dat")
x <- ip[,1]
y <- ip[,2]
z <- ip[,3]
[\code]
Der Aufruf der image-Funktion ist dann wie folgt:
[code]
image(x, y, z, 170, 1200, 1200,col = hcl.colors(12, "YlOrRd", rev =
TRUE),add = FALSE, xaxs = "i", yaxs = "i",breaks, oldstyle = FALSE)
[\code]
Die Fehlerausgabe lautet:
[code]
Fehler in image.default(x, y, z, 170, 1200, 1200, col = hcl.colors(12,
:
ansteigende 'x' und 'y' Werte erwartet
[\code]
Ich kann meine Daten aber nicht stetig ansteigen lassen, oder nicht?
Hat jemand von euch ein probates Mittel, um in R eine vernünftige
Heatmap zu erstellen, die mit meinen Daten arbeiten kann?
Beste Grüße
Jiri
