Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

[Gnuplot] density mapping help needed

  1. Sep 11, 2012 #1
    Hi all,

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

    Code (Text):

    #!/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:

    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 (Text):

    -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.
     
    Last edited: Sep 11, 2012
  2. jcsd
  3. Sep 12, 2012 #2

    chiro

    User Avatar
    Science Advisor

    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.
     
  4. Sep 12, 2012 #3
    chiro,

    how can I do that? gnuplot is new to me..still not familiar with the command
     
  5. Sep 12, 2012 #4

    chiro

    User Avatar
    Science Advisor

    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.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: [Gnuplot] density mapping help needed
  1. Need help in Fortran (Replies: 4)

Loading...