Alex/Gnuplot

Aus Wiki CCC Göttingen
Wechseln zu: Navigation, Suche

... einfach etwas Beispielcode zu Gnuplot, so als Gedankenstütze

<highlightSyntax language="gnuplot"> call="`echo $CALL`"

  1. !rm -rf list.txt
  2. !scp -C alexander@db0smg.de:wspr_data/list.txt .
  3. !cat list.txt | sort | uniq > list1.txt

cmd=sprintf("%s%s%s",'egrep "^[0123456789-]* [0123456789:]* ',call,'.*" list1.txt | head -n1 | perl qth_cor.pl | head -n1')

  1. print cmd

loc=system(cmd)

  1. print loc

cmd=sprintf("%s%s%s%s%s",'egrep "^[0123456789-]* [0123456789:]* ',call,'.*" list1.txt | perl qth_cor.pl | perl conv_nt.pl "',loc,'" > l.txt')

  1. print cmd

system(cmd)


cmd=sprintf("%s%s%s",'cat list1.txt | perl qth_cor.pl | perl conv_nt.pl "', loc,'" > ll.txt')

  1. print cmd

system(cmd)

cmd=sprintf("%s%s%s",'cat world_50m.txt | perl conv_nt.pl "',loc,'" >aw.txt')

  1. print cmd

system(cmd)

cmd=sprintf("%s%s%s",'cp -- aw.txt "aw_',loc,'.txt"')

  1. print cmd

system(cmd)


cmd=sprintf("%s%s%s",'egrep "^[0123456789-]* [0123456789:]* ',call,'.*" list1.txt | sort -n -k12 | cut -d " " -f 12|tail -n1|sed -e "s/\n//"')

  1. print cmd

max_d=system(cmd)

cmd=sprintf("%s%s%s",'egrep "^[0123456789-]* [0123456789:]* ',call,'.*" list1.txt | sort -n -k12 | cut -d " " -f 9 | tail -n1|sed -e "s/\n//"')

  1. print cmd

power=system(cmd)


min_t="`head -n1 list1.txt | cut -d " " -f 1,2`" max_t="`tail -n1 list1.txt | cut -d " " -f 1,2`"

set title sprintf ("Hearing: %s with: %sW on 30meter\n %s max distance\nStart %s - End %s",call, power,max_d,min_t,max_t)

set size sq set polar set rrange [0:(max_d+500)] set angles deg set grid polar set terminal png size 1000,1000

plot 'aw.txt' u ((90-$3)):4 w l t "",\

       'll.txt' u  ((90-$3)):4 w points pt 4 t "",\
       'l.txt' u  ((90-$3)):4 w l t "",\
       'l.txt'  u  ((90-$3)):4:5 with labels t ""

</highlightSyntax>


<highlightSyntax language="gnuplot">

cmd=sprintf("%s%s%s",'egrep "^[0123456789-]* [0123456789:]* ',call,'.*" list1.txt | head -n1 | perl qth_cor.pl | head -n1')

  1. print cmd

loc=system(cmd)

  1. print loc

cmd=sprintf("%s%s%s%s%s",'egrep "^[0123456789-]* [0123456789:]* ',call,'.*" list1.txt | perl qth_cor.pl | perl conv_nt.pl "',loc,'" > l.txt')

  1. print cmd

system(cmd)


  1. hier sind noch alle stationen bei ... auch welche die nur senden

cmd=sprintf("%s%s%s",'cat list1.txt | perl qth_cor.pl | perl conv_nt.pl "', loc,'" > ll.txt')

  1. print cmd

system(cmd)

cmd=sprintf("%s%s%s",'cat world_50m.txt | perl conv_nt.pl "',loc,'" >aw.txt')

  1. print cmd

system(cmd)

cmd=sprintf("%s%s%s",'cp -- aw.txt "aw_',loc,'.txt"')

  1. print cmd

system(cmd)


cmd=sprintf("%s%s%s",'egrep "^[0123456789-]* [0123456789:]* ',call,'.*" list1.txt | sort -n -k12 | cut -d " " -f 12|tail -n1|sed -e "s/\n//"')

  1. print cmd

max_d=system(cmd)

cmd=sprintf("%s%s%s",'egrep "^[0123456789-]* [0123456789:]* ',call,'.*" list1.txt | sort -n -k12 | cut -d " " -f 9 | tail -n1|sed -e "s/\n//"')

  1. print cmd

power=system(cmd)


min_t="`head -n1 list1.txt | cut -d " " -f 1,2`" max_t="`tail -n1 list1.txt | cut -d " " -f 1,2`"

set title sprintf ("Hearing: %s with: %sW on 30meter\n %s max distance\nStart %s - End %s",call, power,max_d,min_t,max_t)

set size sq set polar set rrange [0:(max_d+500)] set angles deg set grid polar set terminal png size 1000,1000

plot 'aw.txt' u ((90-$3)):4 w l t "",\

       'll.txt' u  ((90-$3)):4 w points pt 4 t "",\
       'l.txt' u  ((90-$3)):4 w l t "",\
       'l.txt'  u  ((90-$3)):4:5 with labels t "

</highlightSyntax>