New Reply

[Gnuplot] density mapping help needed

 
Share Thread
Sep11-12, 11:22 PM   #1
 

[Gnuplot] density mapping help needed


Hi all,

I used below bash command to plot my data (.dat file)

Code:
#!/bin/bash

(
echo 'set term jpeg'
echo 'set style data lines'
echo 'set yrange [0:200]'
echo 'set xrange[0:200]' 
echo 'set pm3d map'
echo 'set palette defined (-2 "yellow", 0 "green", 2 "red")'


for f in "$@"
do
#	echo "Processing $f" >&2
	pic_name="${f%.dat}.jpg"
	echo "set out '$pic_name'"
	echo "splot '$f'"
done) | gnuplot
but I got this warning:

Warning: Single isoline (scan) is not enough for a pm3d plot.
Hint: Missing blank lines in the data file? See 'help pm3d' and FAQ.
As I know, to map using gnuplot for 3d density mapping, I should have output with blank lines to separate different x values, like this

Code:
-180.0 -180.0 0.0
-180.0 -179.0 0.2
[...]
-180.0  180.0 0.0

-179.0 -180.0 0.4
-179.0 -179.0 0.2
[...]
-179.0  180.0 0.0

[...]
but, if I have a few thousand lines, with 200 .dat files, what can I do to expedite my job? Anyone can help...thank you in advance.
PhysOrg.com science news on PhysOrg.com

>> City-life changes blackbird personalities, study shows
>> Origins of 'The Hoff' crab revealed (w/ Video)
>> Older males make better fathers: Mature male beetles work harder, care less about female infidelity
Sep12-12, 12:52 AM   #2
 
Hey s_hy and welcome to the forums.

One idea I have for you is to run a script (or write some code to create an executable) to just remove every instance of a double or more than double CR-LF (carriage return line-feed) and just set it to a single CR-LF.
Sep12-12, 12:59 AM   #3
 
chiro,

how can I do that? gnuplot is new to me..still not familiar with the command
Sep12-12, 01:10 AM   #4
 

[Gnuplot] density mapping help needed


It's not a gnuplot specific suggestion: it's a general one that can be implemented on your language (scripting/non-script/whatever) and platform of your choice.

All you want to do is basically read in every line of the file that is relevant and whenever you encounter a file with multiple CR-LF (carriage return-line feed) values just remove all of them but one.

If you want to know what these values are use a hex-editor and look at the raw data and see if it's just a one or two byte code and you'll see what the values are.

The simplest way to do this without even worring about the CR-LF values is to just have a loop that reads in one line at a time and then if the line is blank forget about it and if not write it to another file.

When the loop terminates you'll have another file without any extra spaces in the data and you can use that.

It's basically just an open statement, one loop that reads until end of file, read line, check if empty or with only spaces and if so don't add it to output but if not write that line to the output file, loop terminates, close input file, flush output to existing file close handle, program exits and you have a file that has no gaps.

If you have any platform that allows you to read in one line at a time and write one line at a time you're looking at something like 10-20 lines for the entire thing and will do its magic pretty much on the spot.
New Reply

Similar discussions for: [Gnuplot] density mapping help needed
Thread Forum Replies
Show that an analytic mapping is an open mapping Calculus & Beyond Homework 6
plot some graphs to use in LaTeX Math & Science Software 5
An open mapping is not necessarily a closed mapping in functional analysis Calculus 4
gnuplot Computing & Technology 3
gnuplot Computing & Technology 0