Wieder mal: Gnuplot-Frage (Kugelkoordinaten etc.)

[Tobias P. (hubertus)]

Hallo Kollegen

ich möchte wieder mal was spezielles mit Gnuplot darstellen lassen.
Und zwar habe ich das 3D Richtdiagramm einer Antenne gemessen. Siehe 
Anhang, desired.png ist das 'soll'-Bild. Leider kann die Software, mit 
der das Bild erstellt wurde, nur erbärmliches bmp und jpg (!) 
exportieren, d.h. ich muss aus den Messdaten selber versuchen eine Figur 
zu erstellen da ich ein Vektorbild will.

Die Messdaten der Antenne habe ich hier deponiert:
http://files.kooltek.net/antennapattern/patch_starlab_meas.txt
Sie sind in Kugelkoordinaten (1. Spalte Frequenz, 2. Spalte Phi, 3. 
Spalte Theta, 4. Spalte Gewinn).

Mit diesem Gnuplot-Code erzeuge ich dann meine Figur:

set mapping spherical
set pm3d depthorder
#set pm3d scansbackward
#set pm3d implicit
#unset hidden3d
set ticslevel 0
unset colorbox
set size square
unset tics
unset border

splot 'patch_starlab_meas.txt' u 2:($3-pi/2):($4+30):4 w pm3d notitle

daraus resultiert pattern_1.png, was schon mal nicht so schlecht 
aussieht! wenn ich die Figure aber nun drehe (die "Vorderseite" sieht ja 
nicht so spannend aus, ich möchte die Rückseite sehen!) dann wirds 
komisch. Siehe pattern_2.png. Offenbar ist diese Oberfläche nicht 
geschlossen und hat an den zwei rot markierten Stellen Lücken, was etwas 
sonderbar aussieht. Kann man da Gnuplot irgendwie sagen, es soll diese 
Löcher einfach schliessen? Ich habe das schon versucht, indem ich den 
Messdaten den letzten Datensatz für Phi = pi eingefügt habe, aber das 
hat nichts gebracht.

Bin mal gespannt ob ihr ein paar Ideen habt ;-)
Die Frage ist auch noch, ob man die Oberfläche allenfalls mit einem 
Gitternetz ausstatten könnte, damit man die Krümmung besser sieht.

[Max M. (jens2001)]

Kenn mich mit GnuPlot nicht aus, hab aber mal in die Doku gesehen.

Imho passen deine Datensätze nicht zu der Art wie GP sie erwartet.

Wenn ich die Doku richtig verstanden habe verwendet GP die Buchstaben 
Phi und Theta anders als Wikipedia und ich.

Die Werte für den Polarwinkel werden vom Equator aus angegeben. 
Wertebereich +Pi/2 bis -Pi/2

GP erwartet als ersten Wert den Azimutwinkel (+Pi bis -Pi)

[Tobias P. (hubertus)]

Hallo Max
korrekt, das ist so. Und darum subtrahiere ich bei meinem Theta ja noch 
pi/2, um auf den richtigen Wertebereich zu kommen.

[Max M. (jens2001)]

Tobias P. schrieb:
> subtrahiere ich bei meinem Theta ja noch
> pi/2, um auf den richtigen Wertebereich zu kommen.

Bleibt ImhO immer noch falsch.

Aus der GP-Doku
"For a spherical coordinate system, the data occupy two or three columns 
(or using entries). The rst two
are interpreted as the azimuthal and polar angles theta and phi (or 
"longitude" and "latitude"), in the units
specied by set angles. The radius r is taken from the third column if 
there is one, or is set to unity if
there is no third column. The mapping is:
x = r * cos(theta) * cos(phi)
y = r * sin(theta) * cos(phi)
z = r * sin(phi)
Note that this is a "geographic" spherical system, rather than a "polar" 
one (that is, phi is measured from
the equator, rather than the pole)."