Portál AbcLinuxu, 8. srpna 2025 22:30


Dotaz: Guplot - body nad krivkou nebo errorbar

masomlejn avatar 17.7.2008 22:18 masomlejn | skóre: 16
Guplot - body nad krivkou nebo errorbar
Přečteno: 172×
Odpovědět | Admin
Příloha:
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]
Nástroje: Začni sledovat (0) ?Zašle upozornění na váš email při vložení nového komentáře.

Odpovědi

masomlejn avatar 17.7.2008 22:20 masomlejn | skóre: 16
Rozbalit Rozbalit vše Re: Guplot - body nad krivkou nebo errorbar
Odpovědět | | Sbalit | Link | Blokovat | Admin
V nadpisu ma samozrejme byt "Gnuplot - body nad krivkou nebo errorbar", sorry za preklep
masomlejn avatar 25.7.2008 13:41 masomlejn | skóre: 16
Rozbalit Rozbalit vše Re: Gnuplot - body nad krivkou nebo errorbar
Odpovědět | | Sbalit | Link | Blokovat | Admin
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.

Založit nové vláknoNahoru

Tiskni Sdílej: Linkuj Jaggni to Vybrali.sme.sk Google Del.icio.us Facebook

ISSN 1214-1267, (c) 1999-2007 Stickfish s.r.o.