Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]


Groups > comp.graphics.apps.gnuplot > #3716

Re: Ephemerides data

From Karl Ratzsch <mail.kfr@gmx.net>
Newsgroups comp.graphics.apps.gnuplot
Subject Re: Ephemerides data
Date 2017-08-20 23:59 -0300
Organization solani.org
Message-ID <ondia2$rdd$1@solani.org> (permalink)
References <e132bb6e-fa89-4de4-b1a9-c12c13f352ad@googlegroups.com> <ondekk$p4g$1@solani.org>

Show all headers | View raw


Am 20.08.2017 um 22:56 schrieb Karl Ratzsch:
> Am 20.08.2017 um 10:36 schrieb robertvanpass@gmail.com:
>> Hi,
>>
>> I have difficulties to plot ephemerides data.
>> They are generated by SkytechX in the following format:
>>
>> Date;Time;Local R.A.;Local Dec.;
>> 8/20/2017;12:00:00 UTC+1.0;10h 35m 50.60s;+03° 47' 15.38";
>> 8/21/2017;12:00:00 UTC+1.0;10h 33m 13.63s;+04° 02' 23.70";
>> 8/22/2017;12:00:00 UTC+1.0;10h 30m 23.99s;+04° 20' 34.64";
>> ...
>>
>> The column "Local R.A." is in time format: Hours Minits Seconds,
>> should be converted to: Hours + Minits/60 + Seconds/3600.
>>
>> Column "Local Dec." is in angle format: Degrees Arcminits Arcseconds,
>> should be converted to Degrees + Arcminits/60 + Arcseconds/3600.
>>
>> So the first record should become:
>> 8/20/2017;12:00:00 UTC+1.0;10.587606;3.787606";
> 
> $dat << EOD
> 8/20/2017;12:00:00 UTC+1.0;10h 35m 50.60s;+03° 47' 15.38";
> 8/21/2017;12:00:00 UTC+1.0;10h 33m 13.63s;+04° 02' 23.70";
> 8/22/2017;12:00:00 UTC+1.0;10h 30m 23.99s;+04° 20' 34.64";
> EOD
> 
> set dataf sep "; \"
> set table
> 
> plot $dat us \
> 	(timecolumn(1,"%d/%m/%Y")):\
> 	(timecolumn(2,"%H:%M:%S")):\
> 	(timecolumn(4,"%Hh %Mm %Ss")/3600):\
> 	(($7+$8/60+$9/3600)) w table
> 
> gives
> 
>  1.53369e+009 43200 10.5974 3.78761
>  1.53636e+009 43200 10.5538 4.03992
>  1.53896e+009 43200 10.5067 4.34296
> 
> 	
> This ignores the timezone specifier, I'm not sure how to catch that.

Wasn't hard:

set dataf sep "; " # the backslash was nonsense
set table

plot $dat us \
(timecolumn(1,"%d/%m/%Y")):\
(timecolumn(2,"%H:%M:%S")+strcol(3)[4:8]*3600):\
(timecolumn(4,"%Hh %Mm %Ss")/3600):\
(($7+$8/60+$9/3600)) w table

Expects the format "UTC+x.y" or "UTC-x.y", anything else fails of
course.



Back to comp.graphics.apps.gnuplot | Previous | NextPrevious in thread | Next in thread | Find similar


Thread

Ephemerides data robertvanpass@gmail.com - 2017-08-20 06:36 -0700
  Re: Ephemerides data Karl Ratzsch <mail.kfr@gmx.net> - 2017-08-20 22:56 -0300
    Re: Ephemerides data Karl Ratzsch <mail.kfr@gmx.net> - 2017-08-20 23:59 -0300
      Re: Ephemerides data robertvanpass@gmail.com - 2017-08-21 15:01 -0700
        Re: Ephemerides data Karl Ratzsch <mail.kfr@gmx.net> - 2017-08-22 03:59 +0200
          Re: Ephemerides data robertvanpass@gmail.com - 2017-08-22 07:57 -0700
            Re: Ephemerides data Karl Ratzsch <mail.kfr@gmx.net> - 2017-08-22 19:09 +0200
              Re: Ephemerides data robertvanpass@gmail.com - 2017-08-22 14:09 -0700
                Re: Ephemerides data Karl Ratzsch <mail.kfr@gmx.net> - 2017-08-23 01:45 +0200
                Re: Ephemerides data robertvanpass@gmail.com - 2017-08-23 13:09 -0700

csiph-web