Generated: Sat Oct 12 17:22:43 2013 from astro.pl 2013/03/05 8.4 KB. text copy
#!/usr/bin/perl -w # NAME: astro.pl # AIM: SPECIFIC - Load xplane earth_astro.dat and output to csv file use strict; use warnings; use File::Basename; # split path ($name,$dir,$ext) = fileparse($file [, qr/\.[^.]*/] ) use Cwd; my $os = $^O; my $perl_dir = '/home/geoff/bin'; my $PATH_SEP = '/'; my $temp_dir = '/tmp'; if ($os =~ /win/i) { $perl_dir = 'C:\GTools\perl'; $temp_dir = $perl_dir; $PATH_SEP = "\\"; } unshift(@INC, $perl_dir); require 'lib_utils.pl' or die "Unable to load 'lib_utils.pl' Check paths in \@INC...\n"; # log file stuff our ($LF); my $pgmname = $0; if ($pgmname =~ /(\\|\/)/) { my @tmpsp = split(/(\\|\/)/,$pgmname); $pgmname = $tmpsp[-1]; } my $outfile = $temp_dir.$PATH_SEP."temp.$pgmname.txt"; open_log($outfile); # user variables my $VERS = "0.0.1 2013-03-02"; my $load_log = 0; my $in_file = ''; my $verbosity = 0; my $out_file = $temp_dir.$PATH_SEP.'temp-astro.csv'; my $out_file2 = $temp_dir.$PATH_SEP.'temp-visible.csv'; # ### DEBUG ### my $debug_on = 1; my $def_file = 'D:\FG\xplane\1000\earth_astro.dat'; ### program variables my @warnings = (); my $cwd = cwd(); sub VERB1() { return $verbosity >= 1; } sub VERB2() { return $verbosity >= 2; } sub VERB5() { return $verbosity >= 5; } sub VERB9() { return $verbosity >= 9; } sub show_warnings($) { my ($val) = @_; if (@warnings) { prt( "\nGot ".scalar @warnings." WARNINGS...\n" ); foreach my $itm (@warnings) { prt("$itm\n"); } prt("\n"); } else { prt( "\nNo warnings issued.\n\n" ) if (VERB9()); } } sub pgm_exit($$) { my ($val,$msg) = @_; if (length($msg)) { $msg .= "\n" if (!($msg =~ /\n$/)); prt($msg); } show_warnings($val); close_log($outfile,$load_log); exit($val); } sub prtw($) { my ($tx) = shift; $tx =~ s/\n$//; prt("$tx\n"); push(@warnings,$tx); } sub process_in_file($) { my ($inf) = @_; if (! open INF, "<$inf") { pgm_exit(1,"ERROR: Unable to open file [$inf]\n"); } my @lines = <INF>; close INF; my $lncnt = scalar @lines; prt("Processing $lncnt lines, from [$inf]...\n"); my ($line,$i,$lnn,$indata,$len,@arr,$acnt,$type); my ($rasc,$dec,$lum,$name,$total); $lnn = 0; $indata = 0; my @astro = (); my @all = (); $total = 0; for ($i = 0; $i < $lncnt; $i++) { $line = $lines[$i]; chomp $line; $line = trim_all($line); $lnn++; $len = length($line); next if ($len == 0); if ($indata) { @arr = split(/\s+/,$line); $acnt = scalar @arr; $type = $arr[0]; if ($type == 99) { prt("$lnn: End of file 99\n"); last; } if ($acnt < 3) { prt("$lnn: Unknown! [$line]\n"); next; } $rasc = $arr[0]; $dec = $arr[1]; $lum = $arr[2]; $name = ''; $name = join(' ',splice(@arr,3)) if ($acnt > 3); #$name = ($acnt > 3) ? $arr[3] : ''; #prt("$rasc,$dec,$lum,$name\n"); if ($lum <= 6.5 ) { push(@astro,[$rasc,$dec,$lum,$name]); } push(@all,[$rasc,$dec,$lum,$name]); $total++; } else { # wait for version line # 740 Version - build 1795 if ($line =~ /\d+\s+Version\s+-\s+build/) { prt(substr($line,0,24)."\n"); $indata = 1; } } } $acnt = scalar @astro; prt("Loaded $acnt of $total BRIGHT stars from [$inf]\n"); my ($ofile,$ra); $line = '"Right Ascension",Declination,"Visible Magnitude",name'."\n"; $ofile = $out_file2; for ($i = 0; $i < $acnt; $i++) { $ra = $astro[$i]; $rasc = ${$ra}[0]; $dec = ${$ra}[1]; $lum = ${$ra}[2]; $name = ${$ra}[3]; $line .= "$rasc,$dec,$lum,$name\n"; } write2file($line,$ofile); prt("$acnt written to $ofile\n"); $acnt = scalar @all; $line = '"Right Ascension",Declination,"Visible Magnitude",name'."\n"; $ofile = $out_file; for ($i = 0; $i < $acnt; $i++) { $ra = $all[$i]; $rasc = ${$ra}[0]; $dec = ${$ra}[1]; $lum = ${$ra}[2]; $name = ${$ra}[3]; $line .= "$rasc,$dec,$lum,$name\n"; } write2file($line,$ofile); prt("$acnt written to $ofile\n"); } ######################################### ### MAIN ### parse_args(@ARGV); process_in_file($in_file); pgm_exit(0,""); ######################################## sub need_arg { my ($arg,@av) = @_; pgm_exit(1,"ERROR: [$arg] must have a following argument!\n") if (!@av); } sub parse_args { my (@av) = @_; my ($arg,$sarg); while (@av) { $arg = $av[0]; if ($arg =~ /^-/) { $sarg = substr($arg,1); $sarg = substr($sarg,1) while ($sarg =~ /^-/); if (($sarg =~ /^h/i)||($sarg eq '?')) { give_help(); pgm_exit(0,"Help exit(0)"); } elsif ($sarg =~ /^v/) { if ($sarg =~ /^v.*(\d+)$/) { $verbosity = $1; } else { while ($sarg =~ /^v/) { $verbosity++; $sarg = substr($sarg,1); } } prt("Verbosity = $verbosity\n") if (VERB1()); } elsif ($sarg =~ /^l/) { if ($sarg =~ /^ll/) { $load_log = 2; } else { $load_log = 1; } prt("Set to load log at end. ($load_log)\n") if (VERB1()); } elsif ($sarg =~ /^o/) { need_arg(@av); shift @av; $sarg = $av[0]; $out_file = $sarg; prt("Set out file to [$out_file].\n") if (VERB1()); } else { pgm_exit(1,"ERROR: Invalid argument [$arg]! Try -?\n"); } } else { $in_file = $arg; prt("Set input to [$in_file]\n") if (VERB1()); } shift @av; } if ((length($in_file) == 0) && $debug_on) { $in_file = $def_file; prt("Set DEFAULT input to [$in_file]\n"); } if (length($in_file) == 0) { pgm_exit(1,"ERROR: No input files found in command!\n"); } if (! -f $in_file) { pgm_exit(1,"ERROR: Unable to find in file [$in_file]! Check name, location...\n"); } } sub give_help { prt("$pgmname: version $VERS\n"); prt("Usage: $pgmname [options] in-file\n"); prt("Options:\n"); prt(" --help (-h or -?) = This help, and exit 0.\n"); prt(" --verb[n] (-v) = Bump [or set] verbosity. def=$verbosity\n"); prt(" --load (-l) = Load LOG at end. ($outfile)\n"); prt(" --out <file> (-o) = Write output to this file.\n"); } # X-Plane apt.dat codes my $x_code = <<EOF FROM : http://data.x-plane.com/designers.html#Formats Astronomical: 740 Version (latest) =================================== from : http://data.x-plane.com/file_specs/Astro740.htm sample 6.752569 -16.713143 -1.43 Sirius 19.846301 8.867385 0.77 Altair 2.529743 89.264138 1.97 Polaris 6.752569 Right Ascension in decimal hours. Always a positive number. -16.713143 Declination in decimal degrees. Positive declinations are north of the celestial equator (eg. the pole star, Polaris, is at a declination of 89.264138 degrees). -1.43 Visible magnitude of the star. This is a weird logarithmic scale (low numbers are brightest), and stars to a magnitude of +6.5 are considered visible to the naked eye (though this will vary hugely with your local seeing conditions, light pollution, altitude, etc.). Sirius (the brightest star in the night sky) has a negative magnitude (-1.43) because it is very, very bright. Sirius Star name (optional - not used by X-Plane). Any units of angular measure can be used for right ascension, but it is customarily measured in hours ( h ), minutes ( m ), and seconds ( s ), with 24h being equivalent to a full circle. EOF # ================================ # eof - astro.pl