Dotaz: Guplot - body nad krivkou nebo errorbar

[masomlejn]

Příloha:

• gama-nokey.png (31424 bytů)

Bry vecer, snazim se pri vykreslovani grafu v gnuplotu (4.2) dostat body nad (do popredi) krivky tak aby ji prekryvaly. Jde o to ze mam treba body tvorene "prazdnyma" ctvereckama, a ty jsou protnuty errorbarem a krivkou (pokud na ni lezi). Viz obrazek. Nevypada to moc dobre, kdyz jsou tam ty ctvercky, kundicky, a pentagony preskrtnute errorbarem. Bohuzel nemuzu pouzit jen "vyplnene" symboly pro body...je tam moc dat, ktere je treba rozlisit. Nevite, prosim vas, jak je dostat nad tu krivku? Diky, Milan.

Obrazek je generovan nasledujicimi dvema skripty. Upraveni vzheldu obrazku (viewgama.gp):

set terminal wxt enhanced
set timestamp "generated at %H:%M:%S on date %d-%m-%Y by `whoami`"
###########################################################
# set common properties for all output terminals          #
###########################################################
set style line 1 lt 1 lw 2
set style line 2 lt -1 lw 2
set border 15 ls 2
set bar 1.75
set size 1,1
set pointsize 2
Tabs=273.15
Tmax=385.15
Tmin=273.15
Tincr=20
fxaxis(x)=1/x

# xtics
t1=Tmin
t2=t1+Tincr
t3=t2+Tincr
t4=t3+Tincr
t5=t4+Tincr
t6=t5+Tincr
t7=t6+Tincr

# minor xtics
tm1=Tmin + Tincr/2
tm2=t1+Tincr/2
tm3=t2+Tincr/2
tm4=t3+Tincr/2
tm5=t4+Tincr/2
tm6=t5+Tincr/2

###########################################################
# set primary graph  (x1axis in 1/T)                      #
###########################################################
set grid
set xrange [fxaxis(Tmax):fxaxis(Tmin)]
set xtics 0.0026,0.0002,0.00365 nomirror
set mxtics 2
set xlabel "{/=22 1000.T ^{ -1} / K ^{ -1}}"
set format x "%.1s"
set mytics 2
set ylabel "{/=22 ln {/Symbol=22 \g@^{/=22 \245}_{1}}}"
set key outside noreverse
call 'params.dat'
############################################################
# set secondary graph  (x2axis in C^{o})                   #
############################################################
set x2tics ('0' fxaxis(t1),"20" fxaxis(t2),"40" \
fxaxis(t3),"60" fxaxis(t4),"80" fxaxis(t5),"100"\
fxaxis(t6))
set x2label "{/=22 t / ^{o}C}"
replot
############################################################
# write modified (pretty  output to *.eps file          #
############################################################
set notimestamp
set terminal postscript enhanced monochrome "Times-Romain" 22
set output 'gama.eps'
replot
set output 'gama-nokey.eps'
set key default
set key off
replot
pause -1

A pro samotne vykresleni params.dat(tento skrip je volan predchozim skriptem viewgama.gp):

plot "work.tmp" index 0:0 u ($$5/$$1):2:3 notitle w yerrorbars 1,\
 "work.tmp" index 0:0 u ($$5/$$1):2 title "92HWA/OLS" w p ps 1.75 pt 4

replot "work.tmp" index 1:1 u ($$5/$$1):2:3 notitle w yerrorbars 1, \
 "work.tmp" index 1:1 u ($$5/$$1):2 title "99BEN/DOH" w p ps 1.75 pt 6

replot "work.tmp" index 2:2 u ($$5/$$1):2:3 notitle w yerrorbars 1, \
 "work.tmp" index 2:2 u ($$5/$$1):2 title "16DAV" w p ps 1.75 pt 8

replot "work.tmp" index 3:3 u ($$5/$$1):2:3 notitle w yerrorbars 1, \
 "work.tmp" index 3:3 u ($$5/$$1):2 title "31GRO/SAY" w p ps 1.75 pt 10

replot "work.tmp" index 4:4 u ($$5/$$1):2:3 notitle w yerrorbars 1, \
 "work.tmp" index 4:4 u ($$5/$$1):2 title "41VER/WER" w p ps 1.75 pt 12

replot "work.tmp" index 5:5 u ($$5/$$1):2:3 notitle w yerrorbars 1, \
 "work.tmp" index 5:5 u ($$5/$$1):2 title "80BAN/YAL" w p ps 1.75 pt 14

replot "work.tmp" index 6:6 u ($$5/$$1):2:3 notitle w yerrorbars 1, \
 "work.tmp" index 6:6 u ($$5/$$1):2 title "03BER/DOH" w p ps 1.75 pt 5

replot "work.tmp" index 7:7 u ($$5/$$1):2:3 notitle w yerrorbars 1, \
 "work.tmp" index 7:7 u ($$5/$$1):2 title "04BER/DOH" w p ps 1.75 pt 7

replot (36.391525069038+(-28.248913125533)/(1/(x*298.15))+(-29.889623983045)*log(1/(x*298.15))+(0)*(1/(x*298.15)))\
title "FIT with 3 params" with lines ls 1
set pointsize 2
set xrange [0.0026:0.0036]

Odpovědi

[masomlejn]

V nadpisu ma samozrejme byt "Gnuplot - body nad krivkou nebo errorbar", sorry za preklep

[masomlejn]

Tak si odpovim Nasel jsem to po tydnu hledani. V pointtype (zkaratka pt) je treba zvolit cislo z intervalu 70 az 75. Pak jsou symboly po exportu do eps (set term postscript enhanced monochrome) vyplnene, a tim padem nepreskrtnute errorbarem.