Gnuplot: One plot per file - bash

I'm trying to plot the 1st and 3rd columns of multiple files, where each file is supposed to be plotted to an own output.png.
My files have the following names:
VIB2--135--398.6241
VIB2--136--408.3192
VIB2--137--411.3725
...
The first number in the file name is an integer, which ranges from 135-162. The second number is just a decimal number and there is no regular spacing between the values.
Basically I want to do something like this
plot for [a=135:162] 'VIB2--'.a.'--*' u 1:3 w l
although this doesn't work, of course, since the ' * ' is just the placeholder I know from bash and I don't know, if there is something similar in gnuplot.
Furthermore, each of the files should be, as already said above, plotted to its own output.png, where the two numbers should be in the output name, e.g. VIB2--135--398.6241.png.
I tried to come up with a bash script for this, like (edited):
#!/bin/bash
for file in *
do
gnuplot < $file
set xtics 1
set xtics rotate
set terminal png size 1920,1080 enhanced
set output $file.png
plot "$file" u 1:3 w l
done
but I still get
gnuplot> 1 14 -0.05
^
line 0: invalid command
gnuplot> 2 14 0.01
^
line 0: invalid command
...
which are actually the numbers from my input file. So gnuplot thinks, that the numbers I want to plot are commands... ?? Also, when the end of the file is reached, I get the following error message
#PLOT 1
plot: an unrecognized command `0x20' was encountered in the input
plot: the input file `VIB2--162--496.0271' could not be parsed
I've seen a few questions similar to mine, but the solutions didn't really work for me and I cannot add a comment, since I do not have the reputation.
Please help me with this.

gnuplot < $file starts gnuplot and feeds it the content of $file as input. That means gnuplot will now try to execute the commands in the data file which doesn't work.
What you want is a "here document":
gnuplot <<EOF
set xtics 1
set xtics rotate
set terminal png size 1920,1080 enhanced
set output $file.png
plot "$file" u 1:3 w l
EOF
What this does is: The shell reads the text up to the line with solemn EOF, replaces all variables, puts that into a temporary file and then starts gnuplot feeding it the temporary file as input.
Be careful that the file names don't contain spaces, or set output $file.png will not work. To be safe, you should probably use set output "$file.png" but my gnuplot is a bit rusty.

Related

xmgrace: changing color of the XY graph in batch execution

I am using gracebat integrated into my bash script to make simple XY graph from 2D numerical data
gracebat input.xvg -hdevice PNG -fixed 800 600 -world 0 0.02 610 0.5 -printfile output.png -hardcopy
This simple execution without GUI produces the XY graph totally in black color!
would it be possible to quickly change to color of the line (I guess should be s0 color in the batch file) without using batch file for gracebat: e.g. via some option of gracebat or alternatively via editing of the input.xvg (via sed or any shell utility) specifying color or the line:
# s0 line color 2
Alternatively, if the grace has not good possibilities for graph customization, I would be grateful for simple example of the gnuplot batch execution.
What does xmgrace -help say? Please check the corresponding manuals and homepages and do a search on SO and the web before asking such questions. I don't now xmgrace and its documentation.
But a gnuplot command line would be using the command line option -e: (check help command line options).
gnuplot -e "set term pngcairo size 800,600; set output 'myFile.png'; plot 'myFile.dat' u 1:2 w l lc 'red' "

gnuplot format xK or xM

For the format of x-axis, current I am using the following command in pngcairo terminal:
set format x "%.sK"
Which recognizes numbers from 100K to 900K but when it gets to 1 million it prints "1K" instead of "1000K".
what is the command to automatically set the label to "xK" before 1 million and to "xM" after 1 million?
These kind of labels are controlled by gnuplot's own format specifiers (see doc for gprintf):
set format x '%.s%c'

Bash script for gnuplot

I have a c++ code that generates every time a variable number of ".txt" file.
I am using a bash script that loads the txt files and send them to a gnuplot script. Then I have to plot in the same figure always the same 2 columns of each of them (representing a altitude vs time plot).
My problem is that the x coordinate (1st column of txt files) does not finish always at the same value and, when the last value of x coordinate is not zero,I need to insert a cross in the gnuplot plot. Do you have an idea about how to do it?
This is how I am plotting the multiple files
plot for [object in objectnames] object.".txt" using 1:2 w l title object
where "objectnames" is a string with all the filenames sent by the BASH
Now, for plotting the CROSSES, I thought about two ways:
Send a array of coordinates to the gnuplot script and then set some label but GNUPLOT does not like arrays and I cannot create a variable number of label in this way.
The alternative is to read somehow the txt file inside the gnuplot script but I have no idea how to do it.
Thanks for the help
EDIT
BASH SCRIPT
export DATA_DIR=${1:-/home/../output}
export objectname
space=" "
OUT_FILES="$(find -L $DATA_DIR -name '*_Trajectory.txt')"
for file in $OUT_FILES; do
fig=$(basename "$file")
objectID=$(echo $fig| cut -d'_' -f 1)
# filenames concatenation
objectname=$objectname$objectID$space
done
export figNameAltitude="${file/${fig}/altitudePlot.png}"
gnuplot -e "folder='${DATA_DIR}';objectnames='${objectname}';figName='${figNameAltitude}';lineWidth='${linewidth}'" altitudePlot.gp
GNUPLOT SCRIPT ( I avoid to insert all the lines regarding the title, axis and so on)
plot for [object in objectnames] folder.'/'.object."_Trajectory.txt" using 1:2 w l lw lineWidth title object
This might help
set title "Energy vs. Time for Sample Data"
set xlabel "Time"
set ylabel "Energy"
plot "d1.dat" with lines
pause -1 "Hit Enter to continue"
execute the script file as **$**gnuplot filename
click here for more details

Parametric gnuplot Scripts

When I try to write a shell script such as
gnuplot
cd '/home/cagirici/test'
plot "case1.test" with linespoints
and run it in sudo, the script stops after running gnuplot.
In each .test file, I have 12 columns. I won't plot them all at once but choose what to plot. For instance, if I choose to plot the success percentage, it should plot 1:2 and 1:8. But if I choose error percentage, it should plot 1:4 and 1:10 (it is basically i:(i+6).)
I need to write a scripts that takes two inputs: i) file name ii) what to plot.
Notice that I also want to choose the line type, point type and output shell. (if there is a configuration to set specific line and point types as default values, please tell me this as well)
Can I write this script in Ubuntu? Or should I use a different method.
My script would look something like
case "$2" in
("success percentage") plot "$1" using 1:2, plot "$1" using 1:8
("error percentage") plot "$1" using 1:4, plot "$1" using 1:10
You can make this work by creating a Gnuplot script with variables and then passing values to these variables at run time from command line or your bash script.
Lets call the Gnuplot script Myscript.gp:
set terminal pngcairo
set output 'graph.png'
plot file u x:y w linespoints
Here file and x:y can be varried when you invoke Gnuplot from your bash script. A possible call to Myscript.gp is as follows:
gnuplot -e "x=1;y=2;file = \"mydata.dat\"" Myscript.gp
In the snippet of your bash script that you provide this may look like:
case "$2" in
("success percentage") gnuplot -e "x=1; y=2; file = \"$1\"" Myscript.gp
Below I'm also providing a dummy mydata.dat file which you can use to try things.
1.54 23.66 43.66
1.75 26.25 46.25
1.92 30.20 40.20
2.08 34.46 44.46
2.44 42.08 42.08
2.78 46.81 46.81
3.03 51.10 41.10
3.70 52.99 42.99
4.17 56.15 46.15
4.76 59.34 49.34
If you are calculating the column numbers (x and y) in your bash script then you can also use variables inside the plot command.
gnuplot -e "x=$x;y=$y;file = \"mydata.dat\"" Myscript.gp
Below I show the output once using:
gnuplot -e "x=1;y=2;file = \"mydata.dat\"" Myscript.gp
and once
gnuplot -e "x=1;y=3;file = \"mydata.dat\"" Myscript.gp
I had faced similar problem 2 years back, then my senior team members decided to write a tool for plotting graphs using a simple csv and plotting specifications and tool will plot graphs for you,
Initially it was designed for sar data on linux so named RSAR but in general it is a csv plotter.)
if required you can contact me nachiket.kate90#gmail.com

Gnuplot tic frequency

For a very large data-set, how can gnuplot be used to only put tic marks/labels on the x axis for just the first and last data point?
With gnuplot 4.6 and up, you can use the commands
stats 'data.dat'
set xtics \
(sprintf("%.2g",STATS_min_x) STATS_min_x, \
sprintf("%.2g",STATS_max_x) STATS_max_x)
plot 'data.dat'
With other versions of gnuplot, you can use this similar sequence of commands:
# this setting makes sure we don't make an output right away
set terminal unknown
plot 'data.dat'
set xtics \
(sprintf("%.2g",GPVAL_DATA_X_MIN) GPVAL_DATA_X_MIN, \
sprintf("%.2g",GPVAL_DATA_X_MAX) GPVAL_DATA_X_MAX)
set terminal <actual terminal>
replot
The set xtics command takes comma-separated pairs of strings with data values, all inside parentheses.
(Here I assume you want the minimum and maximum, not the first and last, data points.)
For more info you can run these at the gnuplot command line:
help set format
help set stats
show variables all
I'm adding another answer because my other one really answered a different question, but may be useful. To answer the question,
How to mark the x tics for the first and last data points?
here is my method:
#!/usr/bin/env gnuplot
set terminal png
set output 'test.png'
# define a function to get the first and last values.
# this assumes the file has not changed since first running 'stats',
# (and contains at least one data point)
# `plot/stats 'data.dat' using (firstlast($x))`
# should be interchangable with
# `plot/stats 'data.dat' using x
# that is, the resulting variables should be the same
firstlast(x) = ($0==0) ? (first=$1, last=$1, $1) : (last=$1, $1)
# run stats to find the first and last values
# just on x data column
stats 'data.dat' u (firstlast($1)) nooutput
# set the x tics to the first and last x points
set xtics \
(sprintf("%.2g (first)", first) first, \
sprintf("%.2g (last)", last) last)
print first
print last
plot 'data.dat'
I used this example data file:
data.dat
1 1
2 2
3 3
4 2
0 3
and got this output:

Resources