package App::SourcePlot;
=head1 NAME
App::SourcePlot - implements Souce Plot application
=head1 SYNOPSIS
use App::SourcePlot;
App::SourcePlot::run_sourceplot_gui();
=head1 DESCRIPTION
This module contains the implementation of the Source Plot application,
which can be launched using the L<sourceplot> command.
Source Plot is a simple astronomical source plotter designed to
display a plot of astronomical sources on adjustable axes.
=cut
use strict;
#use warnings;
our $VERSION = '1.26';
use Config::IniFiles;
use Tk;
use Tk::Balloon;
use Tk::FileSelect;
use App::SourcePlot::Plotter::Tk;
use DateTime;
use DateTime::Format::Strptime;
use App::SourcePlot::Source;
use File::HomeDir;
use File::ShareDir qw/dist_file/;
use File::Spec;
use Tk::AstroCatalog;
use Astro::PAL;
use Astro::Telescope;
use Math::Trig;
use Astro::Coords::Planet 0.05;
#global variables that will be used....
my $locateBug = 0;
my $MW;
my $TIME;
my @SOURCE_LIST = ();
my @planets = map {ucfirst($_)} Astro::Coords::Planet::planets();
my $CATALOG_OPEN;
my $EDIT_OPEN;
my $LAST_COMMAND;
my @UNDO_LIST;
my $undoBut;
my $cBut;
my $TimeLap = 30000; # time between white dot updates
my $dotSizeX = 4;
my $dotSizeY = 4;
my @axes = ('Time', 'Elevation', 'Air Mass', 'Azimuth', 'Parallactic Angle');
my %defaults = (
TEL => 'JCMT',
XAXIS => 'Time',
YAXIS => 'Elevation',
TIME => '1:30:00'
);
my $plotter;
my $TEL;
my $X_AXIS;
my $Y_AXIS;
my $DATE;
my $telObject;
my ($minX,$minY);
my ($maxX,$maxY);
my $NUM_POINTS = 97;
my $BUSY = 0;
my $TIMER;
my @COLOR_LIST = qw/#ffaaaa #00ff00 #ff55ff #ffff00 #00ffff #ff00ff #ffffff #ff5555 #55ff55 #55ffff #ffff55/;
my $COLOR_INDEX = 0;
my $EditWin = {};
my $CatWin = {};
#setup dimensions
my $xplot = 590;
my $xborder = 50;
my $yplot = 550;
my $yborder = 30;
my $RESPONSE;
my $H_LIGHT = undef;
my $H_WIDTH = 3;
my $optBut;
my $dateBut;
my $eBut;
my $defaults = undef;
my $balloon;
=head1 METHODS
=over 4
=item B<run_sourceplot_gui>
Initializes the Source Plot GUI application and enters the Tk main loop.
=cut
sub run_sourceplot_gui {
# setting up default values for Source Plot Options
my $defaults_file = File::Spec->catfile(File::HomeDir->my_home(),
'.splotcfg');
if (-e $defaults_file) {
# Set the "fallback" section to allow reading of defaults
# files from previous versions of SourcePlot.
$defaults = new Config::IniFiles(-file => $defaults_file,
-fallback => 'Options');
}
unless (defined $defaults) {
print STDERR $_, "\n" foreach @Config::IniFiles::errors;
$defaults = new Config::IniFiles();
$defaults->SetFileName($defaults_file);
}
$TEL = $defaults->val('Options', 'TEL', $defaults{'TEL'});
$X_AXIS = $defaults->val('Options', 'XAXIS', $defaults{'XAXIS'});
$Y_AXIS = $defaults->val('Options', 'YAXIS', $defaults{'YAXIS'});
$TIME = $defaults->val('Options', 'TIME', $defaults{'TIME'});
$telObject = new Astro::Telescope($TEL);
##### global windows
$MW = new MainWindow();
$MW->positionfrom('user');
$MW->geometry('+0+90');
$MW->title('Source Plot');
$MW->resizable(0, 0);
$MW->iconname('Source Plot');
$MW->update;
$EditWin->{'Changeable'} = 1;
$CatWin->{'Changeable'} = 0;
#create the balloon used for help
$balloon = $MW->Balloon();
# set the date to the current date
my ($ss, $mm, $hh, $md, $mo, $yr, $wd, $yd, $isdst) = gmtime(time);
$mo++; # this catches the month up to the current date
$mo = '0'.$mo if length($mo) < 2;
$md = '0'.$md if length($md) < 2;
$yr += 1900;
$DATE = "$yr\/$mo\/$md";
my $canFrame = $MW->Frame(
-takefocus => 1,
-width => $xplot,
);
# Create a canvas and calculate the world to pixel ratio
$plotter = new App::SourcePlot::Plotter::Tk($xplot, $yplot, 1, 1, $canFrame);
$plotter->setBackground('black');
$plotter->worldCenter($xborder, $yplot-$yborder);
$plotter->usingWorld(1);
# Set up the X Y coordinate feedback
#my ($XPOS, $YPOS);
#$plotter->monitorXY( \$XPOS, \$YPOS );
my $buttonFrame = $MW->Frame(
-takefocus => 1,
-width =>$xplot
);
my $exitBut = $buttonFrame->Button(
-text => 'Exit',
-width => 8,
-command => sub{destroy $MW;}
)->pack(-side=>'right');
$balloon->attach($exitBut,
-balloonmsg => "Press to Exit program");
my $printBut = $buttonFrame->Button(
-text => 'Print',
-width => 8,
-command => sub{
my $choice = 'printer';
my $Top = $MW->Toplevel();
$Top->title('Source Plot Printer Options');
$Top->resizable(0,0);
my $radf = $Top->Frame()->pack(-padx=>10, -pady=>10, -side=>'top');
my $tf = $Top->Frame(-relief=>'groove', -borderwidth =>2)->pack(-padx=>10, -side=>'top');
# create the radio button for the file option
my $radbut = $radf->Radiobutton( -text => "To File",
-value => 'file',
-variable => \$choice );
$radbut->grid(-column=>0, -row=>0, -padx => 10);
# create the print radio button
my $radbut2 = $radf->Radiobutton( -text => "To Printer",
-value => 'printer',
-variable => \$choice );
$radbut2->grid(-column=>1, -row=>0, -padx => 10);
# create the file entry
$tf->Label (
-text => "Filename",
-justify => 'right'
)->grid(-sticky => 'e', -column=>2, -row=>0);
my $fileEnt = $tf->Entry(-relief=>'sunken',
-width=>20)->grid(-column=>3, -row=>0, -padx =>10, -pady=>5);
# create the printer options box
$tf->Label (
-text => "Printer Command",
-justify => 'right'
)->grid(-column=>2, -row=>1, -sticky=>'e');
my $printerEnt = 'lp';
$tf->Entry(-relief=>'sunken',
-textvariable=>\$printerEnt,
-width=>20)->grid(-column=>3, -row=>1, -padx =>10, -pady=>5);
# create the print button
my $bf = $Top->Frame()->pack(-side=>'bottom', -pady=>10, -padx=>10);
$bf->Button(
-text => 'Print',
-width => 8,
-command => sub {
tagOnOff('white', 1, 'dark green');
tagOnOff('black', 1, 'white');
foreach my $source (@SOURCE_LIST) {
tagOnOff($source->name(), 1, 'dark green');
}
if ($choice eq 'printer') {
$plotter->printCanvas($choice, $printerEnt, $MW);
} else {
$plotter->printCanvas($choice, $fileEnt->get, $MW);
}
tagOnOff('white', 1, 'white');
tagOnOff('black', 1, 'black');
foreach my $source (@SOURCE_LIST) {
tagOnOff($source->name(), 1, $source->color());
}
destroy $Top;
}
)->pack(-side=>'right', -padx => 10);
my $optBut = $bf->Button(
-text => 'Cancel',
-width => 8,
-command => sub {destroy $Top;}
)->pack(-side=>'right', -padx =>10);
$Top->update;
$Top->grab;
})->pack(-side=>'right');
$balloon->attach($printBut,
-balloonmsg => "Press to print the graph");
$optBut = $buttonFrame->Button(
-text => 'Options',
-width => 8,
-command => \&changeOpt
)->pack(-side=>'left');
$balloon->attach($optBut,
-balloonmsg => "Press to change the x and y axes");
$dateBut = $buttonFrame->Button(
-text => 'UT Date',
-width => 8,
-command => \&changeDate
)->pack(-side=>'left');
$balloon->attach($dateBut,
-balloonmsg => "Press to change the UT date");
print "made it to just before the menu button, planet\n" if $locateBug;
$eBut = $buttonFrame->Button(
-text => 'Edit',
-width => 8,
-command => sub {
if (!$EDIT_OPEN) {
$EDIT_OPEN = 1;
&editSource ();
}
}
)->pack(-side=>'left');
$balloon->attach($eBut,
-balloonmsg => "Press to open the Edit Window");
$cBut = $buttonFrame->Button(
-text => 'Catalog',
-width => 8,
-command => sub {
if (!$CATALOG_OPEN) {
open_catalog();
}
}
)->pack(-side=>'left');
$balloon->attach($cBut,
-balloonmsg => "Press to open the Catalog Window");
$canFrame->grid(qw/-sticky nsew/);
$buttonFrame->grid(qw/-sticky nsew/);
print "made it to just before the windows come up\n" if $locateBug;
$MW->update;
$eBut->configure(-state=>'disabled');
&editSource ();
open_catalog();
$EDIT_OPEN = 1;
plot();
calcTime();
print "made it to just before the main loop\n" if $locateBug;
MainLoop;
}
=item B<open_catalog>
Creates an Astro::Catalog object by reading the JCMT catalog which is
distributed with this module. This catalog is then passed to the
constructor of Tk::AstroCatalog to open a catalog window.
Also sets the C<$CATALOG_OPEN> variable and disables the catalog button.
=cut
sub open_catalog {
my $astrocat = new Astro::Catalog(Format => 'JCMT',
File => dist_file('App-SourcePlot', 'jcmt.cat'),
ReadOpt => {incplanets => 0});
my $catalog = new Tk::AstroCatalog($MW, -addCmd => \&addCommand,
-onDestroy => \&reset, -upDate => \&update_status,
-catalog => $astrocat);
$catalog->fillWithSourceList('full');
$CATALOG_OPEN = 1;
$cBut->configure(-state=> 'disabled');
}
=item B<reset>
Resets status relating to the catalog window.
Clears the C<$CATALOG_OPEN> variable and enables the catalog button.
=cut
sub reset
{
$CATALOG_OPEN = 0;
$cBut->configure(-state=>'normal');
}
=item B<addCommand>
This is the subroutine which is provided to C<Tk::AstroCatalog> to be called
when a source should be added to the display.
It recieves an C<Astro::Coords> object and uses it to construct an enclosing
C<App::SourcePlot::Source> object.
=cut
sub addCommand
{
my $selected = shift;
if(@{$selected} != ())
{
my $source;
$LAST_COMMAND = 'Add';
@UNDO_LIST = ();
$undoBut->configure(-state => 'normal') if $EDIT_OPEN;
foreach my $coords (@$selected)
{
$source = new App::SourcePlot::Source($coords);
if(!isWithin($source, @SOURCE_LIST))
{
my $s = $source->copy();
$s->color(getColor());
push (@SOURCE_LIST, $s);
push (@UNDO_LIST, $source->copy());
}
}
if ($EDIT_OPEN) {
fillWithSourceList ($EditWin, 'full');
}
plot();
}
}
=item B<waitForResponse>
Waits until the C<$RESPONSE> variable has been changed.
=cut
sub waitForResponse {
$RESPONSE = 'NOTHING';
while ($RESPONSE == 'NOTHING') {
$MW->update;
}
return 0 if $RESPONSE == -1;
return $RESPONSE;
}
=item B<isWithin>
Checks whether a source is already part of a list. The source is compared
to sources from the list by invoking the C<summary> method on their
C<Astro::Coords> objects.
next if isWithin($source, @list);
=cut
sub isWithin {
my $element = shift;
my @array = @_;
my $len = @array;
foreach (@array) {
if ($element->coords()->summary() eq $_->coords()->summary()) {
return 1;
}
}
return 0;
}
=item B<remove>
Removes entries from the list which match the given source. The comparison
is performed in the same way as C<isWithin>.
remove($source, \@list)
=cut
sub remove {
my $element = shift;
my $array = shift;
my $len = @$array;
my @temp;
my $flag = 0;
for (my $index = 0; $index < $len; $index++) {
if ($element->coords()->summary() eq $$array[$index]->coords()->summary()) {
$flag = -1;
} else {
$temp[$index+$flag] = $$array[$index];
}
}
@$array = @temp;
}
=item B<update_status>
Invokes the update method of the main window.
=cut
sub update_status {
$MW->update;
}
=item B<changeDate>
Changes the date to a new date.
=cut
sub changeDate {
$dateBut->configure(-state => 'disabled');
my @months = qw/January February March April May June July August
September October November December/;
my $name;
my $Top = $MW->Toplevel;
$Top->title('Source Plot UT Date');
$Top->resizable(0,0);
my $topFrame = $Top->Frame(-relief=>'groove', -borderwidth =>2)->pack(-padx=>10, -pady=>10);
# create the day entry
$topFrame->Label (
-text => "Day:"
)->grid(-column=>0, -row=>0);
my $dayEnt = $topFrame->Entry(-relief=>'sunken',
-width=>10)->grid(-column=>1, -row=>0, -padx =>10, -pady=>5);
# create the month menu button
$topFrame->Label (
-text => "Month:"
)->grid(-column=>0, -row=>1, -padx =>5, -pady=>5);
my $strp = new DateTime::Format::Strptime(pattern => '%Y/%m/%d',
time_zone => 'UTC', on_error => 'croak');
my $dt = $strp->parse_datetime($DATE);
my $monthEnt = $months[$dt->month() - 1];
my $mb = $topFrame->Menubutton( -text=>$monthEnt, -relief => 'raised', -width => 10);
foreach $name (@months) {
$mb->command(-label =>$name, -command=> sub{
$mb->configure( -text => $name );
$monthEnt = $name;
});
}
$mb->grid(-column=>1, -row=>1, -padx =>10, -pady=>5, -sticky=>'w');
# create the year entry
$topFrame->Label (
-text => "Year:"
)->grid(-column=>0, -row=>2, -padx =>5, -pady=>5);
my $yearEnt = $topFrame->Entry(-relief=>'sunken',
-width=>10)->grid(-column=>1, -row=>2, -padx =>10, -pady=>5);
$yearEnt->bind('<KeyPress-Return>' =>sub {
my $strp = new DateTime::Format::Strptime(
pattern => '%Y %B %d', on_error => 'croak');
my $dt = $strp->parse_datetime($monthEnt . ' ' .
$dayEnt->get() . ' ' . $yearEnt->get());
$DATE = $dt->strftime('%Y/%m/%d');
destroy $Top;
});
# create the update subroutine
my $complete = sub {
my $strp = new DateTime::Format::Strptime(
pattern => '%B %d %Y', on_error => 'croak');
my $dt = $strp->parse_datetime($monthEnt . ' ' .
$dayEnt->get() . ' ' . $yearEnt->get());
$DATE = $dt->strftime('%Y/%m/%d');
foreach my $source (@SOURCE_LIST) {
$source->erasePoints();
}
plot();
};
# create the apply button
my $F = $Top->Frame->pack();
my $buttonF = $F->Frame->pack(-side => 'left', -padx => 5, -pady => 10);
my $okBut = $buttonF->Button(
-text => 'Apply',
-command => $complete
)->pack(-side=>'left');
$okBut->bind('<KeyPress-Return>' => $complete);
# create the accept button
$buttonF = $F->Frame->pack(-side => 'right', -padx => 5, -pady => 10);
my $okBut = $buttonF->Button(
-text => 'Accept',
-command => sub {&$complete; destroy $Top;}
)->pack(-side=>'right');
$okBut->bind('<KeyPress-Return>' => sub {&$complete; destroy $Top;});
# create the cancel button
my $canBut = $buttonF->Button(
-text => 'Cancel',
-command => sub{destroy $Top;}
)->pack(-side=>'right');
$canBut->bind('<KeyPress-Return>' =>sub {destroy $Top;});
# Closing the window should reset $dateBut.
$Top->bind('<Destroy>', sub {
my $widget = shift;
return unless $widget == $Top;
$dateBut->configure(-state => 'normal');
});
$dayEnt->insert(0, $dt->day());
$yearEnt->insert(0, $dt->year());
$MW->update;
}
=item B<changeOpt>
Displays a window allowing the options to be changed.
=cut
sub changeOpt {
my $name;
my $telEnt = $TEL;
my ($tb, $tb2);
$optBut->configure(-state => 'disabled');
my $Top = $MW->Toplevel;
$Top->title('Source Plot Options');
$Top->resizable(0,0);
my $topFrame = $Top->Frame(-relief=>'groove', -borderwidth =>2)->pack(-padx=>10, -pady=>10);
# place the telescope menuButton
$topFrame->Label (
-text => "Telescope:"
)->grid(-column=>0, -row=>0);
my $f = $topFrame->Frame()->
grid(-column=>1, -row=>0, -padx =>10, -pady=>5, -sticky=>'w');
my $f2 = $topFrame->Frame()->
grid(-column=>1, -row=>1, -padx =>10, -pady=>5, -sticky=>'w');
$tb = $f->Menubutton( -text=>$telEnt, -relief => 'raised',
-width => 15)->pack(-side=>'left');
$tb->cascade(-label => 'A - C', -underline => 0);
$tb->cascade(-label => 'D - F', -underline => 0);
$tb->cascade(-label => 'G - I', -underline => 0);
$tb->cascade(-label => 'J - L', -underline => 0);
$tb->cascade(-label => 'M - O', -underline => 0);
$tb->cascade(-label => 'P - R', -underline => 0);
$tb->cascade(-label => 'S - U', -underline => 0);
$tb->cascade(-label => 'V - X', -underline => 0);
$tb->cascade(-label => 'Y - Z', -underline => 0);
my $cm = $tb->cget(-menu);
my $ac = $cm->Menu;
my $df = $cm->Menu;
my $gi = $cm->Menu;
my $jl = $cm->Menu;
my $mo = $cm->Menu;
my $pr = $cm->Menu;
my $su = $cm->Menu;
my $vx = $cm->Menu;
my $yz = $cm->Menu;
$tb->entryconfigure('A - C', -menu => $ac);
$tb->entryconfigure('D - F', -menu => $df);
$tb->entryconfigure('G - I', -menu => $gi);
$tb->entryconfigure('J - L', -menu => $jl);
$tb->entryconfigure('M - O', -menu => $mo);
$tb->entryconfigure('P - R', -menu => $pr);
$tb->entryconfigure('S - U', -menu => $su);
$tb->entryconfigure('V - X', -menu => $vx);
$tb->entryconfigure('Y - Z', -menu => $yz);
foreach $name ($telObject->telNames()) {
if ($name =~ /^[A-Ca-c]/) {
$ac->command(-label =>$name, -command=> sub{
$tb->configure( -text => $name );
$telEnt = $name;
});
} elsif ($name =~ /^[D-Fd-f]/) {
$df->command(-label =>$name, -command=> sub{
$tb->configure( -text => $name );
$telEnt = $name;
});
} elsif ($name =~ /^[G-Ig-i]/) {
$gi->command(-label =>$name, -command=> sub{
$tb->configure( -text => $name );
$telEnt = $name;
});
}elsif ($name =~ /^[J-Lj-l]/) {
$jl->command(-label =>$name, -command=> sub{
$tb->configure( -text => $name );
$telEnt = $name;
});
}elsif ($name =~ /^[M-Om-o]/) {
$mo->command(-label =>$name, -command=> sub{
$tb->configure( -text => $name );
$telEnt = $name;
});
}elsif ($name =~ /^[P-Rp-r]/) {
$pr->command(-label =>$name, -command=> sub{
$tb->configure( -text => $name );
$telEnt = $name;
});
} elsif ($name =~ /^[S-Us-u]/) {
$su->command(-label =>$name, -command=> sub{
$tb->configure( -text => $name );
$telEnt = $name;
});
} elsif ($name =~ /^[V-Xv-x]/) {
$vx->command(-label =>$name, -command=> sub{
$tb->configure( -text => $name );
$telEnt = $name;
});
} elsif ($name =~ /^[Y-Zy-z]/) {
$yz->command(-label =>$name, -command=> sub{
$tb->configure( -text => $name );
$telEnt = $name;
});
}
}
# place the middle time field
$topFrame->Label (
-text => "Center Time (in HST):"
)->grid(-column=>0, -row=>2, -padx =>5, -pady=>5);
my $timeEnt = $TIME;
$topFrame->Entry(-relief=>'sunken',
-textvariable=> \$timeEnt,
-width=>10)->grid(-column=>1, -row=>2, -padx =>10, -pady=>5);
# place the y-axis menuButton
$topFrame->Label (
-text => "y-axis:"
)->grid(-column=>0, -row=>4, -padx =>5, -pady=>5);
my $yEnt = $Y_AXIS;
my $yb = $topFrame->Menubutton( -text=>$yEnt, -relief => 'raised', -width => 15);
foreach $name (@axes) {
$yb->command(-label =>$name, -command=> sub{
$yb->configure( -text => $name );
$yEnt = $name;
});
}
$yb->grid(-column=>1, -row=>4, -padx =>10, -pady=>5, -sticky=>'w');
# inserting cascading choices
# place the x-axis menuButton
$topFrame->Label (
-text => "x-axis:"
)->grid(-column=>0, -row=>3, -padx =>5, -pady=>5);
my $xEnt = $X_AXIS;
my $xb = $topFrame->Menubutton( -text=>$xEnt, -relief => 'raised', -width => 15);
foreach $name (@axes) {
$xb->command(-label =>$name, -command=> sub{
$xb->configure( -text => $name );
$xEnt = $name;
});
}
$xb->grid(-column=>1, -row=>3, -padx =>10, -pady=>5, -sticky=>'w');
# place the Apply button
my $F = $Top->Frame->pack();
my $buttonF = $F->Frame->pack(-side => 'left', -padx => 5, -pady => 10);
$buttonF->Button(
-text => 'Apply',
-command => sub{
$X_AXIS = $xEnt;
$Y_AXIS = $yEnt;
if ($TEL ne $telEnt || $TIME ne $timeEnt) {
foreach (@SOURCE_LIST) {$_->erasePoints()};
} else {
foreach (@SOURCE_LIST) {$_->eraseTimeDot()};
}
$TIME = $timeEnt;
$TIME =~ s/^\s+//;
$TIME =~ s/\s+$//;
$TIME =~ s/\s+/:/g;
$TEL = $telEnt;
$telObject->name($TEL);
plot();
}
)->pack(-side=>'left');
# place the ok button
$buttonF = $F->Frame->pack(-side => 'right', -padx => 5, -pady => 10);
$buttonF->Button(
-text => 'Accept',
-command => sub{
$X_AXIS = $xEnt;
$Y_AXIS = $yEnt;
if ($TEL ne $telEnt || $TIME ne $timeEnt) {
foreach (@SOURCE_LIST) {$_->erasePoints()};
} else {
foreach (@SOURCE_LIST) {$_->eraseTimeDot()};
}
$TIME = $timeEnt;
$TIME =~ s/^\s+//;
$TIME =~ s/\s+$//;
$TIME =~ s/\s+/:/g;
$TEL = $telEnt;
$telObject->name($TEL);
plot();
destroy $Top;
}
)->pack(-side=>'right');
# place the cancel button
$buttonF->Button(
-text => 'Cancel',
-command => sub{
destroy $Top;
}
)->pack(-side=>'right');
# Closing the window should reset $optBut.
$Top->bind('<Destroy>', sub {
my $widget = shift;
return unless $widget == $Top;
$optBut->configure(-state => 'normal');
});
$buttonF->Button(
-text => 'Save Setting',
-command => sub{
$X_AXIS = $xEnt;
$Y_AXIS = $yEnt;
$TEL = $telEnt;
$TIME = $timeEnt;
$defaults->newval('Options', 'TEL', $TEL);
$defaults->newval('Options', 'XAXIS', $X_AXIS);
$defaults->newval('Options', 'YAXIS', $Y_AXIS);
$defaults->newval('Options', 'TIME', $TIME);
$defaults->RewriteConfig();
}
)->pack(-side=>'right');
$MW->update;
}
=item B<addPlanetSource>
Adds a planet source into the plotting list. The planet number
is currently unused.
addPlanetSource($name, $number);
=cut
sub addPlanetSource {
my $name = shift;
my $number = shift;
my $source = new App::SourcePlot::Source($name);
if (!isWithin ($source, @SOURCE_LIST)) {
$source->color(getColor());
push (@SOURCE_LIST, $source);
}
}
=item B<getSource>
Prompts the user to enter source coords and name.
Can specify previous source object to edit.
Returns a Source object.
=cut
sub getSource {
my $source = shift;
my @Epocs = qw/RJ RB GA AZ/;
my $name;
my $Top = $MW->Toplevel;
$Top->title('Source Plot');
$Top->resizable(0,0);
my $topFrame = $Top->Frame(-relief=>'groove', -borderwidth =>2, -width =>50)->pack(-padx=>10, -fill => 'x', -ipady=>10, -pady => 10);
$topFrame->Label (
-text => "Name:"
)->grid(-column=>0, -row=>0);
my $nameEnt = $topFrame->Entry(-relief=>'sunken',
-width=>15)->grid(-column=>1, -row=>0, -padx =>10, -pady=>3);
$nameEnt->insert(0,$source->name());
$topFrame->Label (
-text => "Ra:"
)->grid(-column=>0, -row=>1);
my $raEnt = $topFrame->Entry(-relief=>'sunken',
-width=>15)->grid(-column=>1, -row=>1, -padx =>10, -pady=>3);
$raEnt->insert(0,$source->ra()) unless $source->is_blank();
$topFrame->Label (
-text => "Dec:"
)->grid(-column=>0, -row=>2);
my $decEnt = $topFrame->Entry(-relief=>'sunken',
-width=>15)->grid(-column=>1, -row=>2, -padx =>10, -pady=>3);
$decEnt->insert(0,$source->dec()) unless $source->is_blank();
$topFrame->Label (
-text => "Epoc:"
)->grid(-column=>0, -row=>3, -padx =>5, -pady=>5);
my $epocEnt = $source->is_blank() ? $Epocs[0] : $source->epoc();
my $epocB = $topFrame->Menubutton( -text=>$epocEnt, -relief => 'raised', -width => 15);
foreach $name (@Epocs) {
$epocB->command(-label =>$name, -command=> sub{
$epocB->configure( -text => $name );
$epocEnt = $name;
});
}
$epocB->grid(-column=>1, -row=>3, -padx =>10, -pady=>5, -sticky=>'w');
my $buttonF = $Top->Frame->pack(-padx=>10, -pady=>10);
$buttonF->Button(
-text => 'Ok',
-command => sub{
$source->configure($nameEnt->get(),
$raEnt->get(),
$decEnt->get(),
$epocEnt);
destroy $Top;
$RESPONSE = 1;
}
)->pack(-side=>'right');
$buttonF->Button(
-text => 'Cancel',
-command => sub{
destroy $Top;
$RESPONSE = -1;
}
)->pack(-side=>'right');
$Top->update;
$Top->grab;
}
=item B<editSource>
Edits the existing source list.
=cut
sub editSource {
my (@selected) = ();
my $udeb = 0; #debugger
$eBut->configure(-state=>'disabled');
my $Top = $MW->Toplevel ();
$Top->geometry('+600+90');
$Top->title('Source Plot: Edit Window');
$Top->resizable(0,0);
my $topFrame = $Top->Frame(-relief=>'groove', -borderwidth =>2, -width =>50)->pack(-padx=>10, -fill => 'x', -ipady=>3, -pady => 10);
# create the header
my $headFrame = $topFrame->Frame(-relief=>'flat', -borderwidth =>2)->grid(-row=>0, -sticky=>'nsew', -ipadx => 3);
my $head = $topFrame->Text(
-wrap => 'none',
-relief => 'flat',
-foreground => 'midnightblue',
-width => 50,
-height => 1,
-font => '-*-Courier-Medium-R-Normal--*-120-*-*-*-*-*-*',
-takefocus => 0
)->grid (-sticky=>'ew', -row =>0);
my $title = sprintf "%5s %-16s %-12s %-13s %-4s", 'Index', 'Name', 'Ra', 'Dec', 'Epoc';
$head->insert ('end', $title);
$head->configure(-state=>'disabled');
# create the scollable text
my $T = $topFrame->Scrolled('Text',
-scrollbars => 'e',
-background => '#333333',
-wrap => 'none',
-width => 60,
-height => 15,
-font => '-*-Courier-Medium-R-Normal--*-120-*-*-*-*-*-*',
-setgrid => 1,
)->grid(qw/-sticky nsew/, -row=>1);
$T->bindtags(qw/widget_demo/); # remove all bindings but dummy "widget_demo"
# create the done button
my $buttonF = $Top->Frame->pack(-padx=>10, -pady=>10, -fill =>'x');
my $doneBut = $buttonF->Button(
-text => 'Done',
-width => 4,
-command => sub{
destroy $Top;
}
)->pack(-side=>'right');
$balloon->attach($doneBut,
-balloonmsg => "Press to close the Edit Window");
# Closing the window should reset $eBut.
$Top->bind('<Destroy>', sub {
my $widget = shift;
return unless $widget == $Top;
$eBut->configure(-state=>'normal');
$EDIT_OPEN = 0;
});
# create the On / Off button
my $OOBut = $buttonF->Button( -text=>'ON/Off',
-width => 5,
-command=> sub{
my $source;
foreach $source (@selected) {
if ($source->active()) {
tagOnOff ($source->name(), 0);
$source->active(0);
$T->tag('configure', 'd'.$source->index(), -foreground => 'black');
} else {
tagOnOff ($source->name(), 1, $source->color(), 'yellow');
$source->active(1);
$T->tag('configure', 'd'.$source->index(), -foreground => $source->color());
}
}
@selected = ();
})->pack(-side => 'left');
$balloon->attach($OOBut,
-balloonmsg => "Press to turn the selected\nsources On/Off");
# create the delete all button
my $manipF = $buttonF->Frame->pack(-side=>'left', -padx=>30);
my $DABut = $manipF->Button(
-text => 'Delete All',
-width => 8,
-command => sub{
$LAST_COMMAND = 'Delete';
foreach my $source (@SOURCE_LIST) {
tagOnOff ($source->name(), 1, $source->color(), 'yellow');
$source->active(1);
$EditWin->{Text}->tagDelete('d'.$source->index());
}
@UNDO_LIST = @SOURCE_LIST;
$undoBut->configure(-state => 'normal') if $EDIT_OPEN;
@SOURCE_LIST = ();
@selected = ();
fillWithSourceList ($EditWin, 'full');
plot();
}
)->pack(-side=>'left');
$balloon->attach($DABut,
-balloonmsg => "Press to delete all sources\nfrom the Edit Window");
# create the delete button
my $delBut = $manipF->Button(
-text => 'Delete',
-width => 6,
-command => sub{
my $len = @selected;
$LAST_COMMAND = 'Delete';
@UNDO_LIST = ();
if ($len > 0) {
$undoBut->configure(-state => 'normal');
} else {
$undoBut->configure(-state => 'disabled');
}
foreach my $source (@selected) {
tagOnOff ($source->name(), 1, $source->color(), 'yellow');
$source->active(1);
push (@UNDO_LIST, $source);
$plotter->delete('l'.$source->name());
$plotter->delete('fo'.$source->name());
$plotter->delete('t'.$source->name());
$plotter->delete('time'.$source->name());
$plotter->delete('ttimeDot'.$source->name());
$plotter->delete('ltimeDot'.$source->name());
$plotter->delete('timeDot'.$source->name());
$EditWin->{Text}->tagDelete('d'.$source->index());
remove ($source, \@SOURCE_LIST);
}
@selected = ();
fillWithSourceList ($EditWin, 'full');
}
)->pack(-side=>'left');
$balloon->attach($delBut,
-balloonmsg => "Press to delete the\nselected sources");
# create the delete button
$undoBut = $manipF->Button(
-text => 'Undo',
-command => sub{
if ($LAST_COMMAND eq 'Add') {
print "Last command was add\n" if $udeb;
foreach my $source (@UNDO_LIST) {
remove ($source, \@SOURCE_LIST);
print "removing ".$source->name."\n" if $udeb;
}
} elsif ($LAST_COMMAND eq 'Delete') {
print "Last command was removed\n" if $udeb;
foreach my $source (@UNDO_LIST) {
push (@SOURCE_LIST, $source);
print "Adding ".$source->name."\n" if $udeb;
}
} elsif ($LAST_COMMAND eq 'Change') {
print "Last command was changed\n" if $udeb;
my $source = pop (@UNDO_LIST);
my $orig = pop (@UNDO_LIST);
print "Reconfiguring ".$source->name."\n" if $udeb;
$source->configure($orig);
print "Clearing ".$source->name."\n" if $udeb;
$source->erasePoints();
}
$LAST_COMMAND = '';
@UNDO_LIST = ();
@selected = ();
$undoBut->configure(-state => 'disabled') if $EDIT_OPEN;
fillWithSourceList ($EditWin, 'full');
plot();
}
)->pack(-side=>'left');
$undoBut->configure(-state => 'disabled') if $EDIT_OPEN;
$balloon->attach($undoBut,
-balloonmsg => "Press to undo the last\nmajor command");
# create the add button
my $newBut = $manipF->Button( -text=>'New',
-width => 4,
-command=> sub{
my $source = new App::SourcePlot::Source();
&getSource ($source);
my $res = &waitForResponse();
if ($res && !isWithin ($source, @SOURCE_LIST)) {
$source->color(getColor());
push (@SOURCE_LIST, $source);
fillWithSourceList ($EditWin, 'full');
$LAST_COMMAND = 'Add';
@UNDO_LIST = ();
$undoBut->configure(-state => 'normal') if $EDIT_OPEN;
push (@UNDO_LIST, $source);
plot();
}
})->pack(-side => 'left');
$balloon->attach($newBut,
-balloonmsg => "Press to manually\nadd a source");
# create the planet button
my $planBut = $manipF->Menubutton(
-text =>'Planets',
-width => 7,
-relief => 'raised');
{
my $k = 0;
foreach my $plan (@planets) {
my $c = $k;
$planBut->command(-label=>$plan, -command=> sub{
&addPlanetSource ($plan, $c);
print "$TEL is tel \n" if $locateBug;
fillWithSourceList ($EditWin, 'full');
plot();
});
$k++;
}
}
$planBut->command(-label=>'All', -command=> sub{
my $j = 0;
foreach (@planets) {
&addPlanetSource ($_, $j);
$j++;
}
fillWithSourceList ($EditWin, 'full');
plot();
});
$planBut->pack(-side => 'left');
$balloon->attach($planBut,
-balloonmsg => "Press to add a planet\nto the plot");
$EditWin->{'Window'} = $Top;
$EditWin->{'Text'} = $T;
$EditWin->{'Selected'} = \@selected;
$EditWin->{'Sources'} = \@SOURCE_LIST;
fillWithSourceList ($EditWin, 'full');
$MW->update;
}
=item B<fillWithSourceList>
Fills a Text box with the list of current sources.
=cut
sub fillWithSourceList {
my(@bold, @normal);
my $Window = shift;
my $T = $Window->{'Text'};
my $selected = $Window->{'Selected'};
my $arrayOfSources = $Window->{'Sources'};
my $task = shift;
my $sort = shift;
my $index = shift;
my (@sources);
my ($source);
my @entered = ();
my($line,$itag);
@sources = @$arrayOfSources;
#@sources = sort by_name @$arrayOfSources;
if(defined $sort)
{
for($sort)
{
/unsorted/ and do { @sources = @$arrayOfSources;last;};
/name/ and do { @sources = sort by_name @$arrayOfSources;last;};
/ra/ and do { @sources = sort by_ra @$arrayOfSources;last;};
/dec/ and do { @sources = sort by_dec @$arrayOfSources;last;};
die "Unknown value for form variable \$sort \n";
}
}
# Enable infobox for access
$T->configure(-state=>'normal');
# Clear the existing widgets
if (defined $task && $task eq 'full') {
$T->delete('1.0','end');
foreach my $source (@sources) {
$T->tagDelete('d'.$source->index());
}
}
# Set up display styles
if ($T->depth > 1) {
@bold = (-background => "#eeeeee", qw/-relief raised -borderwidth 1/);
@normal = (-background => undef, qw/-relief flat/);
} else {
@bold = (qw/-foreground white -background black/);
@normal = (-foreground => undef, -background => undef);
}
$T->tag(qw/configure normal -foreground blue/);
$T->tag(qw/configure inactive -foreground black/);
$T->tag(qw/configure selected -foreground red/);
foreach ( @COLOR_LIST ){
$T->tag('configure',$_, -foreground => $_);
}
# Insert the current values
if (defined $task && $task eq 'full') {
my $len = @sources;
for ($index=0; $index < $len; $index++) {
$source = $sources[$index];
$source->index($index);
$line = $source->dispLine();
if (&isWithin ($source, @$selected)) {
inswt($T,"$line\n","d$index",'selected');
} else {
if ($source->active()) {
if ($source->color() ne '') {
inswt($T,"$line\n","d$index",$source->color());
} else {
inswt($T,"$line\n","d$index",'normal');
}
} else {
inswt($T,"$line\n","d$index",'inactive');
}
}
}
$len = @sources;
for ($itag=0; $itag < $len; $itag++) {
my $dtag = "d$itag";
$T->tag('bind', $dtag, '<Any-Enter>' =>
sub {
shift->tag('configure', $dtag, @bold);
if ($Window->{'Changeable'} && $sources[substr($dtag,1,99)]->active()) {
$plotter->configureTag ('l'.$sources[substr($dtag,1,99)]->name(), -width => 3);
}
}
);
$T->tag('bind', $dtag, '<Any-Leave>' =>
sub {
shift->tag('configure', $dtag, @normal);
if ($Window->{'Changeable'} && $sources[substr($dtag,1,99)]->active()) {
$plotter->configureTag ('l'.$sources[substr($dtag,1,99)]->name(), -width => 1);
}
}
);
$T->tag('bind', $dtag, '<ButtonRelease-1>' =>
sub {
if (!$BUSY){
if (! &isWithin ($sources[substr($dtag,1,99)], @$selected) ) {
shift->tag('configure', $dtag, -foreground => 'red');
push (@$selected, $sources[substr($dtag,1,99)]);
} else {
if ($sources[substr($dtag,1,99)]->color() ne '') {
shift->tag('configure', $dtag, -foreground => $sources[substr($dtag,1,99)]->color());
} else {
shift->tag('configure', $dtag, -foreground => 'blue');
}
&remove ($sources[substr($dtag,1,99)], $selected);
}
}
}
);
if ($Window->{'Changeable'}) {
$T->tag('bind', $dtag, '<Double-1>' => sub {
my $source = $sources[substr($dtag,1,99)];
$LAST_COMMAND = 'Change';
@UNDO_LIST = ($source->coords());
&getSource ($source);
my $res = &waitForResponse();
if( $res ) {
$undoBut->configure(-state => 'normal') if $EDIT_OPEN;
} else {
$undoBut->configure(-state => 'disabled') if $EDIT_OPEN;
}
push (@UNDO_LIST, $source);
$source->erasePoints();
fillWithSourceList ($EditWin, 'full');
plot();
});
}else {
$T->tag('bind', $dtag, '<Double-1>' => sub {
$BUSY = 1;
my $source = $sources[substr($dtag,1,99)];
push (@$selected, $source);
$LAST_COMMAND = 'Add';
@UNDO_LIST = ();
$undoBut->configure(-state => 'normal') if $EDIT_OPEN;
$MW->update();
my $T = shift;
foreach $source (@$selected) {
if (!isWithin ($source, @SOURCE_LIST)) {
my $s = $source->copy();
$s->color(getColor());
push (@SOURCE_LIST, $s);
push (@UNDO_LIST, $s);
}
$T->tag('configure', 'd'.$source->index(), -foreground => 'blue');
}
@$selected = ();
if ($EDIT_OPEN) {
fillWithSourceList ($EditWin, 'full');
}
plot();
$BUSY = 0;
});
}
}
}
$T->mark(qw/set insert 1.0/);
# Disable access to infobox
$T->configure(-state=>'disabled');
}
=item B<plot>
Plots the graphs, including axis.
=cut
sub plot {
print "Entered plot\n" if $locateBug;
my ($xworldRatio, $yworldRatio);
my $debug = 0;
$TIME =~ s/^\s+//;
$TIME =~ s/\s+$//;
$TIME =~ s/\s+/:/g;
my ($timeH, $min, $sec) = split (/:/, $TIME,3);
my $XSpaceForTime = 0;
my $YSpaceForTime = 0;
my $special = 0;
my $lstDiffX;
my $lstDiffY;
$XSpaceForTime = $xborder if $Y_AXIS =~ /time/i;
$YSpaceForTime = $yborder if $X_AXIS =~ /time/i;
$timeH += $min/60 + $sec/3600;
#clear the plot
$plotter->clean();
# calc the world coords
print " Calculating the world coordinate system\n" if $locateBug;
if ($X_AXIS =~ /time/i) {
my $strp = new DateTime::Format::Strptime(
pattern => '%Y/%m/%d %H:%M:%S',
on_error => 'croak');
my $dt = $strp->parse_datetime($DATE . ' ' . $TIME);
$dt->subtract(hours => 2);
my ($lst, $mjd) = ut2lst($dt->year(), $dt->month(), $dt->day(),
$dt->hour(), $dt->minute(), $dt->second(),
$telObject->long_by_rad());
$dt->add(hours => 1);
my ($lst2, $mjd2) = ut2lst($dt->year(), $dt->month(), $dt->day(),
$dt->hour(), $dt->minute(), $dt->second(),
$telObject->long_by_rad());
if ($lst2 < $lst) {
$lst2 += 2*pi;
}
$lstDiffX = $lst2 - $lst;
$xworldRatio = (24*$lstDiffX)/($xplot - $xborder*2 - $XSpaceForTime);
$minX = $lst;
$maxX = $lst + 24*$lstDiffX;
print "x-axis is time\n" if $debug;
} elsif ($X_AXIS =~ /air mass/i) {
$xworldRatio = 90/($xplot - $xborder*2 - $XSpaceForTime);
$minX = 0;
$maxX = 90;
print "x-axis is ele\n" if $debug;
} elsif ($X_AXIS =~ /elevation/i) {
$xworldRatio = 90/($xplot - $xborder*2 - $XSpaceForTime);
$minX = 0;
$maxX = 90;
print "x-axis is ele\n" if $debug;
} elsif ($X_AXIS =~ /azimuth/i) {
$xworldRatio = 360/($xplot - $xborder*2 - $XSpaceForTime);
$minX = 0;
$maxX = 360;
print "x-axis is pa\n" if $debug;
} elsif ($X_AXIS =~ /parallactic angle/i) {
$xworldRatio = 360/($xplot - $xborder*2 - $XSpaceForTime);
$minX = -180;
$maxX = 180;
print "x-axis is az\n" if $debug;
} else {
print "ERROR: X axis undefined!!\n\n";
$minX = 0;
$maxX = 0;
}
if ($Y_AXIS =~ /time/i) {
my $strp = new DateTime::Format::Strptime(
pattern => '%Y/%m/%d %H:%M:%S',
on_error => 'croak');
my $dt = $strp->parse_datetime($DATE . ' ' . $TIME);
$dt->subtract(hours => 2);
my ($lst, $mjd) = ut2lst($dt->year(), $dt->month(), $dt->day(),
$dt->hour(), $dt->minute(), $dt->second(),
$telObject->long_by_rad());
$dt->add(hours => 1);
my ($lst2, $mjd2) = ut2lst($dt->year(), $dt->month(), $dt->day(),
$dt->hour(), $dt->minute(), $dt->second(),
$telObject->long_by_rad());
if ($lst2 < $lst) {
$lst2 += 2*pi;
}
$lstDiffY = $lst2 - $lst;
$yworldRatio = (-24*$lstDiffY)/($yplot - $yborder*2 - $YSpaceForTime);
$maxY = $lst + 24 * $lstDiffY;
$minY = $lst;
print "y-axis is time\n" if $debug;
} elsif ($Y_AXIS =~ /air mass/i) {
$yworldRatio = -90/($yplot - $yborder*2- $YSpaceForTime);
$maxY = 90;
$minY = 0;
print "y-axis is ele\n" if $debug;
} elsif ($Y_AXIS =~ /elevation/i) {
$yworldRatio = -90/($yplot - $yborder*2- $YSpaceForTime);
$maxY = 90;
$minY = 0;
print "y-axis is ele\n" if $debug;
} elsif ($Y_AXIS =~ /azimuth/i) {
$yworldRatio = -360/($yplot - $yborder*2- $YSpaceForTime);
$minY = 0;
$maxY = 360;
print "y-axis is az\n" if $debug;
} elsif ($Y_AXIS =~ /parallactic angle/i) {
$yworldRatio = -360/($yplot - $yborder*2- $YSpaceForTime);
$minY = -180;
$maxY = 180;
print "y-axis is pa\n" if $debug;
} else {
print "ERROR: Y axis undefined!!\n\n";
$minY = 0;
$maxY = 0;
}
if ((($Y_AXIS =~ /azimuth/i) && ($X_AXIS =~ /elevation/i)) || (($X_AXIS =~ /azimuth/i) && ($Y_AXIS =~ /elevation/i))|| (($Y_AXIS =~ /azimuth/i) && ($X_AXIS =~ /air mass/i)) || (($X_AXIS =~ /azimuth/i) && ($Y_AXIS =~ /air mass/i))) {
$yworldRatio = -180/($yplot - $yborder*2);
$xworldRatio = 180/($xplot - $xborder*2);
$minY = 0;
$maxY = 180;
$minX = 0;
$maxX = 180;
$plotter->usingWorld(0);
$plotter->worldCenter($xplot/2, $yplot/2);
$plotter->usingWorld(1);
if(($X_AXIS =~ /elevation/i) || ($Y_AXIS =~ /elevation/i))
{
$special = 1;
}
else # air mass case
{
$special = 2;
}
print "axes is special\n" if $debug;
}
print " Setting the world coordinate system\n" if $locateBug;
$plotter->worldToPixRatio($xworldRatio, $yworldRatio);
$plotter->worldAtZero($minX, $minY);
# draw the grid lines - x axis
$plotter->fontColor('White');
$plotter->drawColor('dark green');
if (($special == 1) ||($special == 2)) {
# draw the angle grid lines
for (my $z = 0; $z < 360; $z += 15) {
my $tz = (90 - $z) * pi / 180;
my $x = 92 * cos $tz;
my $y = 92 * sin $tz;
$plotter->drawLine (0,0,$x, $y, 'grid');
}
#label the angle grid lines
for (my $z = 0; $z < 360; $z += 15) {
my $tz = (90 - $z) * pi / 180;
my $x = 95 * cos $tz;
my $y = 95 * sin $tz;
if ($z == 0) {
$plotter->drawText ($x, $y, 'N', 'twhite');
} elsif ($z == 90) {
$plotter->drawText ($x, $y, 'E', 'twhite');
} elsif ($z == 180) {
$plotter->drawText ($x, $y, 'S', 'twhite');
} elsif ($z == 270) {
$plotter->drawText ($x, $y, 'W', 'twhite');
} else {
$plotter->drawText ($x, $y, $z, 'twhite');
}
}
# draw the outside circle
$plotter->drawColor('white');
$plotter->penWidth(2);
$plotter->drawOval(-90,-90,90,90, 'owhite');
$plotter->penWidth(1);
if($special == 1) # elevation -azimuth case
{
#draw the elevation labels
$plotter->drawColor ('dark green');
my $texty = $plotter->toWy(5) - $plotter->toWy(0);
for (my $i = 0; $i<90; $i+=10) {
$plotter->drawOval(-$i, -$i, $i, $i, 'grid');
$plotter->drawText(0,$i-$texty,90 - $i, 'twhite');
}
$plotter->usingWorld(0);
$plotter->worldCenter($xborder, $yplot-$yborder);
$plotter->usingWorld(1);
}
if ($special == 2) #azimuth - air mass case
{
#draw the elevation labels
$plotter->drawColor ('dark green');
my $ytext = $plotter->toWx(5) - $plotter->toWx(0);
for (qw/1.0 1.25 1.5 1.75 2.0 3.0 4.0 5.0/)
{
my $new_el = 90 - am_to_deg($_);
$plotter->drawOval(-$new_el, -$new_el, $new_el, $new_el, 'grid');
}
$plotter->drawText (0, 0,1.0,'twhite');
$plotter->drawText (0, 39,1.25,'twhite');
$plotter->drawText (0, 50,1.5,'twhite');
$plotter->drawText (0, 57,1.75,'twhite');
$plotter->drawText (0, 61,2.0,'twhite');
$plotter->drawText (0, 71,3.0,'twhite');
$plotter->drawText (0, 76,4.0,'twhite');
$plotter->drawText (0, 80,5.0,'twhite');
$plotter->usingWorld(0);
$plotter->worldCenter($xborder, $yplot-$yborder);
$plotter->usingWorld(1);
}
}elsif ($X_AXIS =~ /time/i) {
my $ho = $timeH - 8;
for ( my $hour = $minX+$lstDiffX*4; $hour < $maxX; $hour += $lstDiffX*4 ) {
my $tel = $telObject;
$plotter->drawLine($hour, $minY, $hour, $maxY, 'grid');
# calculate the HST time
my $ti = $ho;
$ti -= 24 if $ti >= 24;
$ti += 24 if $ti < 0;
my ($h, $m) = split (/\./, $ti);
$m = ('.'.$m) * 60;
$m = '0'.$m if $m < 10;
my $y = $plotter->toWy ($yplot - $yborder*3/2);
$plotter->drawText($hour, $y, sprintf("%2s:%2s", $h, $m), 'twhite');
# calculate the ut time
my $uh = $h + 10;
$uh -= 24 if $uh >= 24;
my $uy = $plotter->toWy ($yborder/2);
$plotter->drawText($hour, $uy, sprintf("%2s:%2s", $uh, $m), 'twhite');
#calculate the lst time
my $lst;
my $si = '+';
my @parts;
$si = '-' if ($hour < 0);
$lst = $hour;
while ($lst > 2 * pi) {
$lst -= 2*pi;
}
my ($si, @hms) = palDr2tf(2, $lst);
my ($lh, $lm, $ls, undef) = @hms;
$lm += 1 if $ls >= 30;
if ($lm >= 60) {
$lm -= 60;
$lh += 1;
$lh -= 24 if $lh > 24;
}
$lm = '0'.$lm if $lm < 10;
my $ly = $plotter->toWy ($yborder*3/2);
$plotter->drawText($hour, $ly, sprintf("%2s:%2s", $lh, $lm), 'twhite');
$ho += 4;
}
} elsif ($X_AXIS =~ /elevation/i) {
my $y = $plotter->toWy ($yplot - $yborder - 10);
for ( my $deg = $minX+20; $deg < $maxX; $deg += 10 ) {
$plotter->drawLine($deg, $minY, $deg, $maxY, 'grid');
$plotter->drawText($deg, $y, $deg, 'twhite');
}
} elsif ($X_AXIS =~ /air mass/i) {
my $dAt12 = 90 - (180 / pi * acos (1/1.25));
my $dAt15 = 90 - (180 / pi * acos (1/1.5));
my $dAt17 = 90 - (180 / pi * acos (1/1.75));
my $dAt3 = 90 - (180 / pi * acos (1/3));
my $dAt4 = 90 - (180 / pi * acos (1/4));
my $dAt5 = 90 - (180 / pi * acos (1/5));
$plotter->drawLine(30, $minY,30, $maxY, 'grid');
$plotter->drawLine($dAt12, $minY, $dAt12, $maxY, 'grid');
$plotter->drawLine($dAt15, $minY, $dAt15, $maxY, 'grid');
$plotter->drawLine($dAt17, $minY, $dAt17, $maxY, 'grid');
$plotter->drawLine($dAt3, $minY, $dAt3, $maxY, 'grid');
$plotter->drawLine($dAt4, $minY, $dAt4, $maxY, 'grid');
$plotter->drawLine($dAt5, $minY, $dAt5, $maxY, 'grid');
my $y = $plotter->toWy ($yplot - $yborder - 10);
$plotter->drawText($dAt12, $y, '1.25', 'twhite');
$plotter->drawText($dAt15, $y, '1.5', 'twhite');
$plotter->drawText($dAt17, $y, '1.75', 'twhite');
$plotter->drawText(30, $y, '2', 'twhite');
$plotter->drawText($dAt3, $y, '3', 'twhite');
$plotter->drawText($dAt4, $y, '4', 'twhite');
$plotter->drawText($dAt5, $y, '5', 'twhite');
} elsif ($X_AXIS =~ /azimuth/i) {
for ( my $deg = $minX+40; $deg < $maxX; $deg += 40 ) {
$plotter->drawLine($deg, $minY, $deg, $maxY, 'grid');
my $y = $plotter->toWy ($yplot - $yborder - 10);
$plotter->drawText($deg, $y, $deg, 'twhite');
}
} elsif ($X_AXIS =~ /parallactic angle/i) {
for ( my $deg = $minX+40; $deg < $maxX; $deg += 40 ) {
$plotter->drawLine($deg, $minY, $deg, $maxY, 'grid');
my $y = $plotter->toWy ($yplot - $yborder - 10);
$plotter->drawText($deg, $y, $deg, 'twhite');
}
}
# draw the grid lines - y axis
if (($special == 1) || ($special == 2)) {}
elsif ($Y_AXIS =~ /time/i) {
my $ho = $timeH - 8;
for ( my $hour = $minY+$lstDiffY*4; $hour < $maxY; $hour += $lstDiffY*4 ) {
# draw the time lines
my $tel = $telObject;
$plotter->drawLine($minX, $hour, $maxX, $hour, 'grid');
# calculate the HST time
my $ti = $ho;
$ti -= 24 if $ti >= 24;
$ti += 24 if $ti < 0;
my ($h, $m) = split (/\./, $ti);
$m = ('.'.$m) * 60;
$m = '0'.$m if $m < 10;
$m = '0'.$m if $m < 10;
my $x = $plotter->toWx ($xborder*3/2);
$plotter->drawText($x, $hour, sprintf("%2s:%2s", $h, $m), 'twhite');
# calculate the ut time
my $uh = $h + 10;
$uh -= 24 if $uh >= 24;
my $ux = $plotter->toWx ($xplot - $xborder/2);
$plotter->drawText($ux, $hour, sprintf("%2s:%2s", $uh, $m), 'twhite');
#calculate the lst time
my $lst;
my $si = '+';
my @parts;
$si = '-' if ($hour < 0);
$lst = $hour;
while ($lst > 2 * pi) {
$lst -= 2*pi;
}
my ($si, @hms) = palDr2tf(2, $lst);
my ($lh, $lm, $ls, undef) = @hms;
$lm += 1 if $ls >= 30;
if ($lm >= 60) {
$lm -= 60;
$lh += 1;
$lh -= 24 if $lh > 24;
}
$lm = '0'.$lm if $lm < 10;
my $lx = $plotter->toWx ($xplot - $xborder*3/2);
$plotter->drawText($lx, $hour, sprintf("%2s:%2s", $lh, $lm), 'twhite');
$ho += 4;
}
} elsif ($Y_AXIS =~ /elevation/i) {
for ( my $deg = $minY+20; $deg < $maxY; $deg += 10 ) {
$plotter->drawLine($minX, $deg, $maxX, $deg, 'grid');
my $x = $plotter->toWx ($xborder + 10);
$plotter->drawText($x, $deg, $deg, 'twhite');
}
} elsif ($Y_AXIS =~ /air mass/i) {
my $dAt12 = 90 - (180 / pi * acos (1/1.25));
my $dAt15 = 90 - (180 / pi * acos (1/1.5));
my $dAt17 = 90 - (180 / pi * acos (1/1.75));
my $x = $plotter->toWx ($xborder +15);
for my $am (qw/1 2 3 4 5/)
{
my $dAt = 90 - (180 / pi * acos (1/$am));
$plotter->drawLine($minX, $dAt, $maxX, $dAt, 'grid');
$plotter->drawText ($x, $dAt, $am, 'twhite');
}
$plotter->drawLine($minX, $dAt12, $maxX, $dAt12, 'grid');
$plotter->drawLine($minX, $dAt15, $maxX, $dAt15, 'grid');
$plotter->drawLine($minX, $dAt17, $maxX, $dAt17, 'grid');
$plotter->drawText ($x, $dAt12, '1.25', 'twhite');
$plotter->drawText ($x, $dAt15, '1.5', 'twhite');
$plotter->drawText ($x, $dAt17, '1.75', 'twhite');
} elsif ($Y_AXIS =~ /azimuth/i) {
for ( my $deg = $minY+40; $deg < $maxY; $deg += 40 ) {
$plotter->drawLine($minX, $deg, $maxX, $deg, 'grid');
my $x = $plotter->toWx ($xborder + 10);
$plotter->drawText($x, $deg, $deg, 'twhite');
}
} elsif ($Y_AXIS =~ /parallactic angle/i) {
for ( my $deg = $minY+40; $deg < $maxY; $deg += 40 ) {
$plotter->drawLine($minX, $deg, $maxX, $deg, 'grid');
my $x = $plotter->toWx ($xborder + 10);
$plotter->drawText($x, $deg, $deg, 'twhite');
}
}
if (!$special) {
# erase any mess under the y axis or to the left of the x axis
my $bottom = $plotter->toWy($yplot+10);
my $top = $plotter->toWy(-10);
my $west = $plotter->toWx(-10);
my $east = $plotter->toWx($xplot + 10);
$plotter->drawColor('black');
$plotter->drawBox($west, $maxY, $east, $top, 'foblack');
$plotter->drawBox($west, $minY, $east, $bottom, 'foblack');
$plotter->drawBox($west, $bottom, $minX, $top, 'foblack');
$plotter->drawBox($maxX, $bottom, $east, $top, 'foblack');
# plot the axes and graphs - last to make look neat
print " Plotting the grid border and labels\n" if $locateBug;
$plotter->drawColor('White');
$plotter->fontColor('White');
$plotter->drawLine($minX, $minY, $minX, $maxY, 'lwhite');
$plotter->drawLine($minX, $minY, $maxX, $minY, 'lwhite');
$plotter->drawLine($maxX, $minY, $maxX, $maxY, 'lwhite');
$plotter->drawLine($minX, $maxY, $maxX, $maxY, 'lwhite');
$plotter->raiseAbove('lwhite', 'foblack');
$plotter->raiseAbove('twhite', 'foblack');
}
# label the axes
$plotter->usingWorld(0);
if ($special) {}
elsif ($X_AXIS =~ /time/i) {
$plotter->drawText($xplot/2,$yplot - $yborder/2, 'HST', 'twhite');
$plotter->drawText($xplot/2,$yborder, 'UT', 'twhite');
$plotter->drawText($xplot/2,$yborder*2, 'LST', 'twhite');
$plotter->drawLine($xborder, $yborder, $xplot - $xborder, $yborder, 'lwhite');
$plotter->drawLine($xborder, $yborder, $xborder, $yborder + $YSpaceForTime, 'lwhite');
$plotter->drawLine($xplot - $xborder, $yborder, $xplot - $xborder, $yborder + $YSpaceForTime, 'lwhite');
} elsif ($X_AXIS =~ /elevation/i) {
$plotter->usingWorld(1);
my $yoff = $maxY + $plotter->toWy (5) - $plotter->toWy(0);
my $toff = $plotter->toWy (10) - $plotter->toWy(0);
$plotter->drawLine (30, $maxY, 30, $yoff, 'lwhite');
my $dAt12 = 90 - (180 / pi * acos (1/1.25));
my $dAt15 = 90 - (180 / pi * acos (1/1.5));
my $dAt17 = 90 - (180 / pi * acos (1/1.75));
my $dAt3 = 90 - (180 / pi * acos (1/3));
my $dAt4 = 90 - (180 / pi * acos (1/4));
my $dAt5 = 90 - (180 / pi * acos (1/5));
$plotter->drawLine ($dAt12, $maxY, $dAt12, $yoff, 'lwhite');
$plotter->drawLine ($dAt15, $maxY, $dAt15, $yoff, 'lwhite');
$plotter->drawLine ($dAt17, $maxY, $dAt17, $yoff, 'lwhite');
$plotter->drawLine ($dAt3, $maxY, $dAt3, $yoff, 'lwhite');
$plotter->drawLine ($dAt4, $maxY, $dAt4, $yoff, 'lwhite');
$plotter->drawLine ($dAt5, $maxY, $dAt5, $yoff, 'lwhite');
$plotter->drawText (90, $maxY + $toff, '1', 'twhite');
$plotter->drawText (30, $maxY + $toff, '2', 'twhite');
$plotter->drawText ($dAt12, $maxY + $toff, '1.25', 'twhite');
$plotter->drawText ($dAt15, $maxY + $toff, '1.5', 'twhite');
$plotter->drawText ($dAt17, $maxY + $toff, '1.75', 'twhite');
$plotter->drawText ($dAt3, $maxY + $toff, '3', 'twhite');
$plotter->drawText ($dAt4, $maxY + $toff, '4', 'twhite');
$plotter->drawText ($dAt5, $maxY + $toff, '5', 'twhite');
$plotter->usingWorld(0);
$plotter->drawText($xplot/2,$yborder/2, 'Air Mass', 'twhite');
$plotter->drawText($xplot/2,$yplot - $yborder/2, 'Elevation', 'twhite');
} elsif ($X_AXIS =~ /air mass/i) {
$plotter->usingWorld(1);
my $yoff = $maxY + $plotter->toWy (5) - $plotter->toWy(0);
my $toff = $plotter->toWy (10) - $plotter->toWy(0);
for ( my $deg = $minX+20; $deg < $maxX; $deg += 20 ) {
$plotter->drawLine ($deg, $maxY, $deg, $yoff, 'lwhite');
$plotter->drawText ($deg, $maxY + $toff, $deg, 'twhite');
}
$plotter->usingWorld(0);
$plotter->drawText($xplot/2,$yborder/2, 'Elevation', 'twhite');
$plotter->drawText($xplot/2,$yplot - $yborder/2, 'Air Mass', 'twhite');
} elsif ($X_AXIS =~ /azimuth/i) {
$plotter->drawText($xplot/2,$yplot - $yborder/2, 'Azimuth', 'twhite');
} elsif ($X_AXIS =~ /parallactic angle/i) {
$plotter->drawText($xplot/2,$yplot - $yborder/2, 'Parallactic Angle', 'twhite');
}
if ($special) {}
elsif ($Y_AXIS =~ /time/i) {
$plotter->usingWorld(1);
my $y = $plotter->toPy(4) - $plotter->toPy(0);
$plotter->usingWorld(0);
$plotter->drawText($xborder/2, $yplot/2 + $y/2, 'HST', 'twhite');
$plotter->drawText($xplot - $xborder, $yplot/2 + $y/2,'UT', 'twhite');
$plotter->drawText($xplot - $xborder*2, $yplot/2 + $y/2, 'LST', 'twhite');
$plotter->drawLine($xplot - $xborder, $yborder, $xplot - $xborder, $yplot- $yborder, 'lwhite');
$plotter->drawLine($xplot - $xborder, $yborder, $xplot - $xborder - $XSpaceForTime, $yborder, 'lwhite');
$plotter->drawLine($xplot - $xborder, $yplot - $yborder, $xplot - $xborder - $XSpaceForTime, $yplot - $yborder, 'lwhite');
} elsif ($Y_AXIS =~ /elevation/i) {
$plotter->usingWorld(1);
my $xoff = $maxX - ($plotter->toWx (5) - $plotter->toWx(0));
my $toff = $plotter->toWx (15) - $plotter->toWx(0);
my $dAt12 = 90 - (180 / pi * acos (1/1.25));
my $dAt15 = 90 - (180 / pi * acos (1/1.5));
my $dAt17 = 90 - (180 / pi * acos (1/1.75));
my $dAt3 = 90 - (180 / pi * acos (1/3));
my $dAt4 = 90 - (180 / pi * acos (1/4));
my $dAt5 = 90 - (180 / pi * acos (1/5));
my $dAt1 = 90 - (180 / pi * acos (1));
my $dAt2 = 90 - (180 / pi * acos (1/2));
$plotter->drawLine ($maxX, $dAt1, $xoff, $dAt1, 'lwhite');
$plotter->drawLine ($maxX, $dAt2, $xoff, $dAt2, 'lwhite');
$plotter->drawLine ($maxX, $dAt12, $xoff, $dAt12, 'lwhite');
$plotter->drawLine ($maxX, $dAt15, $xoff, $dAt15, 'lwhite');
$plotter->drawLine ($maxX, $dAt17, $xoff, $dAt17, 'lwhite');
$plotter->drawLine ($maxX, $dAt3, $xoff, $dAt3, 'lwhite');
$plotter->drawLine ($maxX, $dAt4, $xoff, $dAt4, 'lwhite');
$plotter->drawLine ($maxX, $dAt5, $xoff, $dAt5, 'lwhite');
$plotter->drawText ($maxX-$toff, 90, '1', 'twhite');
$plotter->drawText ($maxX-$toff, $dAt2, '2', 'twhite');
$plotter->drawText ($maxX-$toff, $dAt12, '1.25', 'twhite');
$plotter->drawText ($maxX-$toff, $dAt15, '1.5', 'twhite');
$plotter->drawText ($maxX-$toff, $dAt17, '1.75', 'twhite');
$plotter->drawText ($maxX-$toff, $dAt3, '3', 'twhite');
$plotter->drawText ($maxX-$toff, $dAt4, '4', 'twhite');
$plotter->drawText ($maxX-$toff, $dAt5, '5', 'twhite');
$plotter->usingWorld(0);
$plotter->drawTextVert($xplot - $xborder/2,$yplot/2, 'Air Mass', 'twhite');
$plotter->drawTextVert($xborder/2,$yplot/2, 'Elevation', 'twhite');
} elsif ($Y_AXIS =~ /air mass/i) {
$plotter->usingWorld(1);
my $xoff = $maxX - ($plotter->toWx (5) - $plotter->toWx(0));
my $x = $plotter->toWx ($xplot - $xborder -10);
for ( my $deg = $minY+20; $deg < $maxY; $deg += 20 ) {
$plotter->drawText($x, $deg, $deg, 'twhite');
$plotter->drawLine ($maxX, $deg, $xoff, $deg, 'twhite');
}
$plotter->usingWorld(0);
$plotter->drawTextVert($xplot - $xborder/2,$yplot/2, 'Elevation', 'twhite');
$plotter->drawTextVert($xborder/2,$yplot/2, 'Air Mass', 'twhite');
} elsif ($Y_AXIS =~ /azimuth/i) {
$plotter->drawTextVert($xborder/2,$yplot/2, 'Azimuth', 'twhite');
} elsif ($Y_AXIS =~ /parallactic angle/i) {
$plotter->drawTextVert($xborder/2,$yplot/2, 'Parallactic Angle', 'twhite');
}
$plotter->usingWorld(1);
# now interact them
print " Plotting the grid lines\n" if $locateBug;
if (($X_AXIS eq $Y_AXIS) || (($X_AXIS =~ /elevation/i) &&
($Y_AXIS =~ /air mass/i)) || (($X_AXIS =~ /air mass/i) &&
($Y_AXIS =~ /elevation/i))) {
$plotter->drawColor('purple');
$plotter->drawLine ($minX, $minY, $maxX, $maxY);
} else {
my ($y, $mo, $d) = split (/\//, $DATE, 3);
my @points = ();
my $Top;
foreach my $source (@SOURCE_LIST) {
if ($source->active()) {
$source->calcPoints($DATE, $TIME, $NUM_POINTS, $MW, $telObject);
if ($X_AXIS =~ /time/i) {
if (($Y_AXIS =~ /elevation/i) || ($Y_AXIS =~ /air mass/i)){
@points = $source->time_ele_points();
} elsif ($Y_AXIS =~ /azimuth/i) {
@points = $source->time_az_points();
} elsif ($Y_AXIS =~ /parallactic angle/i) {
@points = $source->time_pa_points();
}
} elsif (($X_AXIS =~ /elevation/i) || ($X_AXIS =~ /air mass/i)) {
if ($Y_AXIS =~ /time/i) {
@points = $source->ele_time_points();
} elsif ($Y_AXIS =~ /azimuth/i) {
@points = $source->ele_az_points();
} elsif ($Y_AXIS =~ /parallactic angle/i) {
@points = $source->ele_pa_points();
}
} elsif ($X_AXIS =~ /azimuth/i) {
if ($Y_AXIS =~ /time/i) {
@points = $source->az_time_points();
} elsif (($Y_AXIS =~ /elevation/i) || ($Y_AXIS =~ /air mass/i)){
@points = $source->ele_az_points();
} elsif ($Y_AXIS =~ /parallactic angle/i) {
@points = $source->az_pa_points();
}
} elsif ($X_AXIS =~ /parallactic angle/i) {
if ($Y_AXIS =~ /time/i) {
@points = $source->pa_time_points();
} elsif (($Y_AXIS =~ /elevation/i) || ($Y_AXIS =~ /air mass/i)){
@points = $source->pa_ele_points();
} elsif ($Y_AXIS =~ /azimuth/i) {
@points = $source->pa_az_points();
}
}
$plotter->penWidth($source->lineWidth());
$plotter->penWidth($H_WIDTH) if ($H_LIGHT == $source);
$plotter->drawColor($source->color());
print "color plotted is ".$source->color()."\n" if $debug;
my @times = $source->time_ele_points();
if (($special == 1) || ($special == 2)) {
my $i;
my $len = @points;
my @newpoints = ();
my $prevr;
my $prevt;
my $dotColor = "#55ffff";
for ($i = 0; $i < $len; $i += 2) {
my $r = $points[$i];
$r = 90 - $r;
my $theta = $points[$i+1];
$theta = (90 - $theta)*pi/180;
my $x = 90 + $r * cos $theta;
my $y = 90 + $r * sin $theta;
my $len2 = @newpoints;
if ($r < 90 ) {
if ($len2 == 0 && defined $prevr && $prevr > 90) {
my $dr = 90 - $r;
my $pdr = $prevr - 90;
my $nt = $dr / ($dr + $pdr) * ($prevt - $theta) + $theta;
my $x2 = 90+90 * cos $nt;
my $y2 = 90+90 * sin $nt;
push (@newpoints, $x2);
push (@newpoints, $y2);
}
if ($r > 60 && defined $prevr && $prevr < 60){
my $time = $timeH - 12 + int($i/2) * ( 24 / ($NUM_POINTS - 1));
$time += 24 if $time < 0;
$time -= 24 if $time > 24;
my ($h, $m) = split (/\./, $time);
$m = ('.'.$m) * 60;
$m = '0'.$m if $m < 10;
$h = '0'.$h if $h < 10;
$time = sprintf ("%2s:%2s", $h, $m);
$plotter->fontColor('#ffff00');
$plotter->drawText($x, $y+3, "$time", 'time'.$source->name());
$plotter->drawFillOval ($x-1, $y-1, $x+1, $y+1, 'fo'.$source->name());
} elsif ($r < 60 && defined $prevr && $prevr > 60) {
my $time = $timeH - 12 + int($i/2) * ( 24 / ($NUM_POINTS - 1));
my $di = 24 / ($NUM_POINTS - 1);
$time -= $di;
$time += 24 if $time < 0;
$time -= 24 if $time > 24;
my ($h, $m) = split (/\./, $time);
$m = ('.'.$m) * 60;
$m = '0'.$m if $m < 10;
$h = '0'.$h if $h < 10;
$time = sprintf ("%2s:%2s", $h, $m);
my $x2 = 90 + $prevr * cos $prevt;
my $y2 = 90 + $prevr * sin $prevt;
$plotter->fontColor('#ffff00');
$plotter->drawText($x2, $y2+3, "$time", 'time'.$source->name());
$plotter->drawFillOval ($x2-1, $y2-1, $x2+1, $y2+1, 'fo'.$source->name());
}
$plotter->drawFillOval ($x-.5, $y-.5, $x+.5, $y+.5, 'fo'.$source->name());
push (@newpoints, $x);
push (@newpoints, $y);
} elsif ($len2 > 0 && defined $prevr && $prevr < 90) {
my $dr = $r - 90;
my $pdr = 90 - $prevr;
my $nt = $dr / ($dr + $pdr) * ($prevt - $theta) + $theta;
$x = 90+90 * cos $nt;
$y = 90+90 * sin $nt;
$plotter->fontColor($source->color());
$plotter->drawTextFromLeft($x+5, $y, $source->name(), 't'.$source->name());
push (@newpoints, $x);
push (@newpoints, $y);
$plotter->drawSmoothLine (@newpoints, 'l'.$source->name());
@newpoints = ();
} elsif (defined $len2 and $len2 >= 4) {
$plotter->drawSmoothLine (@newpoints, 'l'.$source->name());
@newpoints = ();
}
$prevr = $r;
$prevt = $theta;
}
@points = @newpoints;
} else {
#give the points some labels
my $max = $points[1];
my $xm = $points[0];
my $len = @points;
my $xw = $plotter->toWx(1) - $plotter->toWx(0);
my $yw = $plotter->toWy(1) - $plotter->toWy(0);
for (my $i = 1; $i < $len; $i+=2) {
my $x = $points[$i-1];
my $y = $points[$i];
$plotter->drawFillOval ($x-$xw, $y-$yw, $x+$xw, $y+$yw, 'fo'.$source->name());
if ($points[$i] > $max) {
$max = $points[$i];
$xm = $points[$i-1];
}
}
my $ya = $plotter->toWy(5) - $plotter->toWy(0);
$plotter->fontColor($source->color());
$plotter->drawTextFromLeft($xm, $max-$ya, $source->name(), 't'.$source->name());
my @second = ();
my $prevx = $points[0];
my $prevy = $points[1];
my $plotNow = 0;
for (my $i = 3; $i < $len; $i+=2) {
my $x = $points[$i-1];
my $y = $points[$i];
my $plotx = $x;
my $ploty = $y;
my $plotxn = $x;
my $plotyn = $y;
push (@second, $prevx);
push (@second, $prevy);
if ($Y_AXIS =~ /parallactic angle/i) {
if ($prevy < -80 && $y > 80) {
my $len2 = @second;
if ($len2 > 0) {
$ploty = $y - 360;
$plotyn = $prevy + 360;
$plotNow = 1;
}
} elsif ($y < -80 && $prevy > 80) {
my $len2 = @second;
if ($len2 > 0) {
$ploty = $y + 360;
$plotyn = $prevy - 360;
$plotNow = 1;
}
}
} elsif ($X_AXIS =~ /parallactic angle/i) {
if ($prevx < -80 && $x > 80) {
my $len2 = @second;
if ($len2 > 0) {
$plotx = $x - 360;
$plotxn = $prevx + 360;
$plotNow = 1;
}
} elsif ($x < -80 && $prevx > 80) {
my $len2 = @second;
if ($len2 > 0) {
$plotx = $x + 360;
$plotxn = $prevx - 360;
$plotNow = 1;
}
}
}
if ($Y_AXIS =~ /azimuth/i) {
if ($prevy < 100 && $y > 260) {
my $len2 = @second;
if ($len2 > 0) {
$ploty = $y - 360;
$plotyn = $prevy + 360;
$plotNow = 1;
}
} elsif ($y < 100 && $prevy > 260) {
my $len2 = @second;
if ($len2 > 0) {
$ploty = $y + 360;
$plotyn = $prevy - 360;
$plotNow = 1;
}
}
} elsif ($X_AXIS =~ /azimuth/i) {
if ($prevx < 100 && $x > 260) {
my $len2 = @second;
if ($len2 > 0) {
$plotx = $x - 360;
$plotxn = $prevx + 360;
$plotNow = 1;
}
} elsif ($x < 100 && $prevx > 260) {
my $len2 = @second;
if ($len2 > 0) {
$plotx = $x + 360;
$plotxn = $prevx - 360;
$plotNow = 1;
}
}
}
if ($plotNow) {
$plotNow = 0;
push (@second, $plotx);
push (@second, $ploty);
$plotter->drawSmoothLine (@second, 'l'.$source->name());
@second = ();
push (@second, $plotxn);
push (@second, $plotyn);
} elsif ($i == $len - 1) {
push (@second, $x);
push (@second, $y);
}
$prevx = $x;
$prevy = $y;
}
@points = @second;
}
my $len2 = @points;
$plotter->drawSmoothLine (@points, 'l'.$source->name()) if ($len2 > 0);
my @bold = (-background => "#bbbbbb", -foreground => 'black', qw/-relief raised -borderwidth 1/);
my @normal = (-background => undef, -foreground => $source->color(), qw/-relief flat/);
$plotter->bindTag( 'l'.$source->name(), '<Any-Enter>' =>
sub {
my $s = $source;
if ($s->active()) {
$EditWin->{'Text'}->tag('configure', 'd'.$s->index(), @bold) if $EDIT_OPEN;
$plotter->configureTag ('l'.$s->name(), -width => 3);
}
});
$plotter->bindTag( 'l'.$source->name(), '<Any-Leave>' =>
sub {
my $s = $source;
if ($s->active()) {
$EditWin->{'Text'}->tag('configure', 'd'.$s->index(), @normal) if $EDIT_OPEN;
$plotter->configureTag ('l'.$s->name(), -width => 1);
}
});
}
$plotter->penWidth(1);
#calculate where the time dots go.
calcTime($source);
if ($plotter->existTag('foblack')) {
if ($plotter->existTag('timeDot'.$source->name())){
$plotter->raiseAbove('timeDot'.$source->name(),'l'.$source->name());
$plotter->raiseAbove('foblack', 'timeDot'.$source->name());
} elsif ($plotter->existTag('l'.$source->name())){
$plotter->raiseAbove('foblack', 'l'.$source->name());
}
$plotter->raiseAbove('t'.$source->name(), 'foblack') if $plotter->existTag('t'.$source->name());
$plotter->raiseAbove('lwhite', 'foblack');
$plotter->raiseAbove('twhite', 'foblack');
} else {
$plotter->raiseAbove('lwhite', 'l'.$source->name()) if $plotter->existTag('l'.$source->name());
$plotter->raiseAbove('twhite', 'l'.$source->name()) if $plotter->existTag('l'.$source->name());
}
}
if ($Top != undef) {
destroy $Top;
}
$plotter->penWidth(1);
}
print "Exit plot\n" if $locateBug;
}
=item B<am_to_deg>
Converts air mass to degrees.
=cut
sub am_to_deg {
(180/pi)*asin(1/$_[0]);
}
=item B<getColor>
Returns a color.
=cut
sub getColor {
my $color = $COLOR_LIST[$COLOR_INDEX];
my $len = @COLOR_LIST;
$COLOR_INDEX++;
$COLOR_INDEX = $COLOR_INDEX % $len;
return $color;
}
=item B<calcTime>
Draws a dot at the current time on each source.
=cut
sub calcTime {
$TIMER->cancel if defined $TIMER;
my $so = shift;
my $timeBug = 0;
my ($ss, $mm, $hh, $md, $mo, $yr, $wd, $yd, $isdst) = gmtime(time);
$mo++; # this catches the month up to the current date
$mo = '0'.$mo if length($mo) < 2;
$md = '0'.$md if length($md) < 2;
$mm = '0'.$mm if length($mm) < 2;
$ss = '0'.$ss if length($ss) < 2;
$yr += 1900;
my ($sety, $setm, $setd) = split (/\//, $DATE);
if ((!$yr =~ /$sety/) || $setm != $mo || $setd != $md) {
return;
}
print "The gm time is $hh:$mm:$ss and date is $yr\/$mo\/$md\n" if $timeBug;
#calculate the local time
my $strp = new DateTime::Format::Strptime(
pattern => '%Y/%m/%d %H:%M:%S',
on_error => 'croak');
my $dt = $strp->parse_datetime("$yr\/$mo\/$md $hh:$mm:$ss");
$dt->subtract(hours => 10);
my $t = $dt->strftime('%H:%M:%S');
my $d = $dt->strftime('%Y/%m/%d');
#my $time = $hh2+$mm2/60+$ss2/3600;
$plotter->drawColor ('white');
foreach my $source (@SOURCE_LIST) {
next if (!$source->active());
next if (defined $so && $source != $so);
if (($so == undef) || ($so->timeDotX() eq '')) {
$plotter->delete ('timeDot'.$source->name());
print "Real time = $TIME and real date = $DATE\n" if $timeBug;
print "Date = $d and time = $t before calcpoint\n" if $timeBug;
my ($lst, $ele, $az, $pa, $elex, $eley, $azx, $azy) = $source->calcPoint($d, $t, $telObject);
if ($lst < $minX) {
$lst += 2*pi;
} elsif ($lst > $maxX) {
$lst -= 2*pi;
}
my ($x, $y);
if ($X_AXIS =~ /time/i) {
$x = $lst;
} elsif (($X_AXIS =~ /elevation/i) || ($X_AXIS =~ /air mass/i)) {
$x = $ele;
} elsif ($X_AXIS =~ /azimuth/i) {
$x = $az;
} elsif ($X_AXIS =~ /parallactic angle/i) {
$x = $pa;
}
if ($Y_AXIS =~ /time/i) {
$y = $lst;
} elsif (($Y_AXIS =~ /elevation/i) || ($Y_AXIS =~ /air mass/i)) {
$y = $ele;
} elsif ($Y_AXIS =~ /azimuth/i) {
$y = $az;
} elsif ($Y_AXIS =~ /parallactic angle/i) {
$y = $pa;
}
if ((($Y_AXIS =~ /azimuth/i) && ($X_AXIS =~ /elevation/i)) || (($X_AXIS =~ /azimuth/i) && ($Y_AXIS =~ /elevation/i)) || (($Y_AXIS =~ /azimuth/i) &&
($X_AXIS =~ /air mass/i)) || (($Y_AXIS =~ /air mass/i) &&
($X_AXIS =~ /azimuth/i))) {
$x = $ele;
$y = $az;
}
$source->timeDotX($x);
$source->timeDotY($y);
$elex = $plotter->toWx($elex) - $plotter->toWx(0);
$eley = $plotter->toWy($eley) - $plotter->toWy(0);
$azx = $plotter->toWx($azx) - $plotter->toWx(0);
$azy = $plotter->toWy($azy) - $plotter->toWy(0);
# draw the time dot
if ($ele > 0) {
$plotter->drawColor($source->color());
if ((($Y_AXIS =~ /azimuth/i) && ($X_AXIS =~ /elevation/i)) || (($X_AXIS =~ /azimuth/i) && ($Y_AXIS =~ /elevation/i)) || (($Y_AXIS =~ /azimuth/i) &&
($X_AXIS =~ /air mass/i)) || (($Y_AXIS =~ /air mass/i) &&
($X_AXIS =~ /azimuth/i))) {
my $r = 90 - $source->timeDotX();
if ($r < 90) {
my $theta = (90 - $source->timeDotY()) * pi / 180;
my $x2 = 90 + $r * cos $theta;
my $y2 = 90 + $r * sin $theta;
$plotter->drawFillOval($x2 - 2, $y2 - 2, $x2 + 2, $y2+2, 'timeDot'.$source->name());
$plotter->bindTag( 'timeDot'.$source->name(), '<Any-Enter>' => sub {});
$plotter->bindTag('ltimeDot'.$source->name(), '<Any-Leave>' => sub {});
}
}
else {
my ($sX, $sY) = $plotter->toW($dotSizeX, $dotSizeY);
$sX -= $plotter->toWx(0);
$sY -= $plotter->toWy(0);
my $x2 = $source->timeDotX();
my $y2 = $source->timeDotY();
if ($x2 ne '' && $y2 ne '') {
$plotter->drawFillOval($x2 - $sX, $y2 - $sY, $x2 + $sX, $y2+$sY, 'timeDot'.$source->name());
$plotter->bindTag( 'timeDot'.$source->name(), '<Any-Enter>' =>
sub {
$plotter->drawColor($source->color());
$plotter->drawFillOval($x2 - $sX, $y2 - $sY, $x2 + $sX, $y2+$sY, 'ltimeDot'.$source->name());
# draw AZ-EL pointers
$plotter->drawColor('red');
$plotter->fontColor('red');
$plotter->penWidth(2);
$plotter->drawLine ($x2 - $elex, $y2 - $eley, $x2, $y2, 'ltimeDot'.$source->name());
$plotter->drawLine ($x2 - $azx, $y2 - $azy, $x2, $y2, 'ltimeDot'.$source->name());
$plotter->penWidth(1);
my ($xoff, $yoff) = $plotter->toW(7,7);
$xoff -= $plotter->toWx(0);
$yoff -= $plotter->toWy(0);
$plotter->drawText ($x2 - $elex-$xoff, $y2 - $eley-$yoff, 'El', 'ltimeDot'. $source->name());
# and RA-Dec box
my $bx = 1.5*sqrt($azx*$azx+$elex*$elex);
my $by = 1.5*sqrt($azy*$azy+$eley*$eley);
$plotter->drawLine ($x2 - $bx, $y2 - $by, $x2 - $bx, $y2 + $by, 'ltimeDot'.$source->name());
$plotter->drawLine ($x2 - $bx, $y2 + $by, $x2 + $bx, $y2 + $by, 'ltimeDot'.$source->name());
$plotter->drawLine ($x2 + $bx, $y2 + $by, $x2 + $bx, $y2 - $by, 'ltimeDot'.$source->name());
$plotter->drawLine ($x2 + $bx, $y2 - $by, $x2 - $bx, $y2 - $by, 'ltimeDot'.$source->name());
$plotter->drawText ($x2, $y2 - 1.2*$by, 'R.A.', 'ltimeDot'. $source->name());
$plotter->drawText ($x2 - 1.12*$bx, $y2+0.2*$by, 'D', 'ltimeDot'. $source->name());
$plotter->drawText ($x2 - 1.1*$bx, $y2, 'e', 'ltimeDot'. $source->name());
$plotter->drawText ($x2 - 1.1*$bx, $y2-0.2*$by, 'c', 'ltimeDot'. $source->name());
$plotter->fontColor($source->color());
});
$plotter->bindTag('ltimeDot'.$source->name(), '<Any-Leave>' =>
sub {
$plotter->delete ('ltimeDot'.$source->name());
});
$source->AzElOffsets($elex, $eley, $azx, $azy);
}
}
}
}
else {
# draw the time dot
if ($source->elevation() > 0) {
$plotter->drawColor($source->color());
if ((($Y_AXIS =~ /azimuth/i) && ($X_AXIS =~ /elevation/i)) || (($X_AXIS =~ /azimuth/i) && ($Y_AXIS =~ /elevation/i)) || (($Y_AXIS =~ /azimuth/i) &&
($X_AXIS =~ /air mass/i)) || (($Y_AXIS =~ /air mass/i) &&
($X_AXIS =~ /azimuth/i))) {
my $r = 90 - $source->timeDotX();
if ($r < 90) {
my $theta = (90 - $source->timeDotY()) * pi / 180;
my $x2 = 90 + $r * cos $theta;
my $y2 = 90 + $r * sin $theta;
$plotter->drawFillOval($x2 - 2, $y2 - 2, $x2 + 2, $y2+2, 'timeDot'.$source->name());
$plotter->bindTag( 'timeDot'.$source->name(), '<Any-Enter>' => sub {});
$plotter->bindTag('ltimeDot'.$source->name(), '<Any-Leave>' => sub {});
}
}
else {
my ($sX, $sY) = $plotter->toW($dotSizeX, $dotSizeY);
$sX -= $plotter->toWx(0);
$sY -= $plotter->toWy(0);
my $x2 = $source->timeDotX();
my $y2 = $source->timeDotY();
if ($x2 ne '' && $y2 ne '') {
my ($elex, $eley, $azx, $azy) = $source->AzElOffsets();
$plotter->drawColor($source->color());
$plotter->drawFillOval($x2 - $sX, $y2 - $sY, $x2 + $sX, $y2+$sY, 'timeDot'.$source->name());
$plotter->bindTag( 'timeDot'.$source->name(), '<Any-Enter>' =>
sub {
$plotter->drawColor($source->color());
$plotter->drawFillOval($x2 - $sX, $y2 - $sY, $x2 + $sX, $y2+$sY, 'ltimeDot'.$source->name());
$plotter->drawColor('red');
$plotter->fontColor('red');
$plotter->penWidth(2);
# draw AZ-EL pointers
$plotter->drawLine ($x2 - $elex, $y2 - $eley, $x2, $y2, 'ltimeDot'.$source->name());
$plotter->drawLine ($x2 - $azx, $y2 - $azy, $x2, $y2, 'ltimeDot'.$source->name());
$plotter->penWidth(1);
my ($xoff, $yoff) = $plotter->toW(7,7);
$xoff -= $plotter->toWx(0);
$yoff -= $plotter->toWy(0);
$plotter->drawText ($x2 - $elex-$xoff, $y2 - $eley-$yoff, 'El', 'ltimeDot'. $source->name());
# and RA-Dec box
my $bx = 1.5*sqrt($azx*$azx+$elex*$elex);
my $by = 1.5*sqrt($azy*$azy+$eley*$eley);
$plotter->drawLine ($x2 - $bx, $y2 - $by, $x2 - $bx, $y2 + $by, 'ltimeDot'.$source->name());
$plotter->drawLine ($x2 - $bx, $y2 + $by, $x2 + $bx, $y2 + $by, 'ltimeDot'.$source->name());
$plotter->drawLine ($x2 + $bx, $y2 + $by, $x2 + $bx, $y2 - $by, 'ltimeDot'.$source->name());
$plotter->drawLine ($x2 + $bx, $y2 - $by, $x2 - $bx, $y2 - $by, 'ltimeDot'.$source->name());
$plotter->drawText ($x2, $y2 - 1.2*$by, 'R.A.', 'ltimeDot'. $source->name());
$plotter->drawText ($x2 - 1.12*$bx, $y2+0.2*$by, 'D', 'ltimeDot'. $source->name());
$plotter->drawText ($x2 - 1.1*$bx, $y2, 'e', 'ltimeDot'. $source->name());
$plotter->drawText ($x2 - 1.1*$bx, $y2-0.2*$by, 'c', 'ltimeDot'. $source->name());
$plotter->fontColor($source->color());
});
$plotter->bindTag('ltimeDot'.$source->name(), '<Any-Leave>' =>
sub {
$plotter->delete ('ltimeDot'.$source->name());
});
}
}
}
}
}
$TIMER = $MW->after($TimeLap, \&calcTime );
}
=item B<tagOnOff>
Configures all object with the tag name off or on.
=cut
sub tagOnOff {
my $tag = shift;
my $turnOn = shift;
my $color = shift;
my $color2 = shift;
if (!$turnOn) {
$plotter->configureTag('l'.$tag, -fill=>'black');
$plotter->configureTag('fo'.$tag, -fill=>'black', outline=>'black');
$plotter->configureTag('o'.$tag, outline=>'black');
$plotter->configureTag('t'.$tag, -fill=>'black');
$plotter->configureTag('time'.$tag, -fill=>'black');
$plotter->configureTag('timeDot'.$tag, -fill=>'black', -outline => 'black');
if (!($tag =~/(white)/i)) {
$plotter->lowerBelow('l'.$tag, 'grid');
$plotter->lowerBelow('o'.$tag, 'grid');
$plotter->lowerBelow('fo'.$tag, 'grid');
$plotter->lowerBelow('t'.$tag, 'grid');
$plotter->lowerBelow('time'.$tag, 'grid');
$plotter->lowerBelow('timeDot'.$tag, 'grid');
}
} else {
$plotter->configureTag('l'.$tag, -fill=>$color);
$plotter->configureTag('fo'.$tag, -fill=>$color, outline=>$color);
$plotter->configureTag('o'.$tag, outline=>$color);
$plotter->configureTag('timeDot'.$tag, -fill=>$color, outline=>$color);
$plotter->configureTag('t'.$tag, -fill=>$color);
if (defined $color2) {
$plotter->configureTag('time'.$tag, -fill=>$color2);
} else {
$plotter->configureTag('time'.$tag, -fill=>$color);
}
$plotter->raiseAbove('l'.$tag, 'grid');
$plotter->raiseAbove('o'.$tag, 'grid');
$plotter->raiseAbove('fo'.$tag, 'grid');
$plotter->raiseAbove('t'.$tag, 'grid');
$plotter->raiseAbove('time'.$tag, 'grid');
$plotter->raiseAbove('timeDot'.$tag, 'grid');
}
}
=item B<error>
Displays an error message in Tk.
=cut
sub error {
my $MW = shift;
my $errWin = $MW->Toplevel(-borderwidth=>10);
$errWin->title('Observation Log Error!');
$errWin->resizable(0,0);
$errWin->Button(
-text => 'Ok',
-command => sub{
destroy $errWin;
})->pack(-side=>'bottom');
my $message = shift;
$errWin->Label (
-text => "\nError!\n\n ".$message." \n",
-relief=>'sunken'
)->pack(-side=>'bottom', -pady => 10);
$errWin->title(shift) if @_;
$MW->update;
$errWin->grab;
}
=item B<inswt>
The "Insert With Tags" procedure inserts text into a given text widget
and applies one or more tags to that text.
Parameters:
$w - Window in which to insert
$text - Text to insert (it's inserted at the "insert" mark)
$args - One or more tags to apply to text. If this is empty
then all tags are removed from the text.
Returns: Nothing
=cut
sub inswt {
my($w, $text, @args) = @_;
my $start = $w->index('insert');
$w->insert('insert', $text);
foreach my $tag ($w->tag('names', $start)) {
$w->tag('remove', $tag, $start, 'insert');
}
foreach my $i (@args) {
$w->tag('add', $i, $start, 'insert');
}
}
1;
__END__
=back
=head1 SEE ALSO
L<Astro::Coords>
L<Astro::Catalog>
L<Tk::AstroCatalog>
=head1 AUTHORS
Casey Best (University of Victoria),
Pam Shimek (University of Victoria),
Tim Jenness (Joint Astronomy Centre),
Remo Tilanus (Joint Astronomy Centre),
Graham Bell (Joint Astronomy Centre).
=head1 COPYRIGHT
Copyright (C) 2012, 2013 Science and Technology Facilities Council.
Copyright (C) 1998, 1999 Particle Physics and Astronomy Research
Council. All Rights Reserved.
This program is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation; either version 3 of the License, or
(at your option) any later version.
This program is distributed in the hope that it will be useful,but
WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
General Public License for more details.
You should have received a copy of the GNU General Public License
along with this program. If not, see <http://www.gnu.org/licenses/>.
=cut