#!/usr/bin/perl 
 
my $filename = 'R001.txt'; 
open (my $fh, '>', $filename) or die "Could not open file '$filename' $!"; 
 
while ( <> ) 
{ 
	if ( /R   1/ ) 
	{ 
		$easting = substr($_,6,8); 
		$northing = substr($_,14,9); 
		 
		print $fh "$easting,$northing\n"; 
	} 
} 
close $fh

R   1 
R 564 
R 565 
R1128 
R1129

#!/usr/bin/perl 
use strict; 
use warnings; 
# use Devel::Size 'total_size'; 
 
my $format = "x a4 x a8 a9 x4 a4 x a8 a9 x4 a4 x a8 a9"; 
            #   ID   E  N     ID   E  N     ID   E  N 
 
my @data; 
 
# reads from 'file2' 
while (<>) { 
    last unless /^R/; 
    my @temp = unpack $format; 
     
    for my $i (0 .. $#temp) { 
        if ($i % 3 == 0) { 
            $data[ $temp[$i] ] = [ @temp[ $i+1, $i+2 ] ]; 
        } 
    } 
} 
 
printf "%4s %10s %s\n", qw/ID East North/; 
 
for my $id (1, 564, 565, 1128, 1129) { 
    printf "%4s %10s %s\n", $id, $data[$id][0], $data[$id][1]; 
} 
 
#print total_size \@data; # prints '1_897_456'

 
C:\Old_Data\perlp>perl test2.pl file2 
  ID       East North 
   1   453448.0 7998331.8 
 564   447738.7 8002432.1 
 565   453478.2 7998372.2 
1128   447767.3 8002471.8 
1129   453507.3 7998412.9 
 
C:\Old_Data\perlp>

$data[$id][0], $data[$id][1];

453448.0,7998331.8 
453454.4,7998327.5 
453461.6,7998322.9

grep ^"R   1" filename > output

C:\Users\Bill\forums\guru>type marky9074.pl 
use strict; 
use warnings; 
use List::MoreUtils qw( any ); 
use Data::Dumper; 
my @required_ids = ('R   1', 'R 564', 'R 565', 'R1128', 'R1129'); 
open my $POS, '<', 'sample3.txt'; 
my %data; 
while (<$POS>) { 
    my ($prefix, $id0, $e0, $n0, $id1, $e1, $n1, $id2, $e2, $n2) 
        = unpack 'A (A4xA8A9x4)2 A4xA8A9', $_; 
    $id0 = $prefix.$id0; 
    if (any {$id0 eq $_} @required_ids) { 
        push @{$data{$id0}}, [$e0, $n0]; 
    } 
    $id1 = $prefix.$id1; 
    if (any {$id1 eq $_} @required_ids) { 
        push @{$data{$id1}}, [$e1, $n1]; 
    } 
    $id2 = $prefix.$id2; 
    if (any {$id2 eq $_} @required_ids) { 
        push @{$data{$id2}}, [$e2, $n2]; 
    } 
} 
print Dumper(\%data) 
 
 
C:\Users\Bill\forums\guru>perl marky9074.pl 
$VAR1 = { 
          'R1129' => [ 
                       [ 
                         '453507.3', 
                         '7998412.9' 
                       ], 
                       [ 
                         '453513.7', 
                         '7998408.6' 
                       ], 
                       [ 
                         '453520.7', 
                         '7998404.1' 
                       ] 
                     ], 
          'R1128' => [ 
                       [ 
                         '447767.3', 
                         '8002471.8' 
                       ], 
                       [ 
                         '447773.2', 
                         '8002466.7' 
                       ], 
                       [ 
                         '447779.8', 
                         '8002461.3' 
                       ] 
                     ], 
          'R 565' => [ 
                       [ 
                         '453478.2', 
                         '7998372.2' 
                       ], 
                       [ 
                         '453484.6', 
                         '7998367.9' 
                       ], 
                       [ 
                         '453491.7', 
                         '7998363.3' 
                       ] 
                     ], 
          'R   1' => [ 
                       [ 
                         '453448.0', 
                         '7998331.8' 
                       ], 
                       [ 
                         '453454.4', 
                         '7998327.5' 
                       ], 
                       [ 
                         '453461.6', 
                         '7998322.9' 
                       ] 
                     ], 
          'R 564' => [ 
                       [ 
                         '447738.7', 
                         '8002432.1' 
                       ], 
                       [ 
                         '447744.7', 
                         '8002427.0' 
                       ], 
                       [ 
                         '447751.3', 
                         '8002421.5' 
                       ] 
                     ] 
        }; 
 
C:\Users\Bill\forums\guru>

#!/usr/bin/perl 
use strict; 
use warnings; 
use DBI; 
 
# see guru_fetch.pl 
 
my $dbh = DBI->connect("dbi:SQLite:dbname=guru.lite","","", 
    {RaiseError => 1, PrintError => 0, AutoCommit => 0}) or die "Can't connect"; 
 
2 == @ARGV or die "Usage: perl $0 tablename filename\n"; 
 
my $table = shift; 
 
#$dbh->do("DROP TABLE IF EXISTS $table"); 
 
$dbh->do(qq{CREATE TABLE $table 
            (id INTEGER, 
            east TEXT, 
            north TEXT) 
            }); 
 
 
my $sth = $dbh->prepare("INSERT INTO $table VALUES(?,?,?)"); 
 
my $format = "x a4 x a8 a9 x4 a4 x a8 a9 x4 a4 x a8 a9"; 
            #   ID   E  N     ID   E  N     ID   E  N 
 
while (<>) { 
    next unless /^R\s*\d+\s/; 
    my @temp = unpack $format; 
     
    for my $i (0 .. $#temp) { 
        if ($i % 3 == 0) { 
            $sth->execute(@temp[$i .. $i+2]); 
        } 
    } 
} 
$dbh->commit; 
$dbh->disconnect;

#!/usr/bin/perl 
use strict; 
use warnings; 
use DBI; 
use Math::Trig qw/asin pi/; # import arcsin & pi 
use Getopt::Std; 
 
# see guru_create.pl 
 
my %opts;  
getopts ('t:i:',\%opts); 
2 == keys %opts or die "Usage: perl $0 -t tablename -i 1,564,565 ...\n"; 
 
my $table = $opts{t}; 
my @ids = split /,/, $opts{i}; 
 
my $dbh = DBI->connect("dbi:SQLite:dbname=guru.lite","","", 
    {RaiseError => 1, PrintError => 0}) or die "Can't connect"; 
 
my $sth = $dbh->prepare("SELECT north, east FROM $table where id = ?"); 
 
for my $id (@ids) { 
    print "\t$id\n"; 
     
    $sth->execute($id); 
    my $aref = $sth->fetchall_arrayref; 
     
    for my $i (0 .. $#$aref-1) { 
        my $n = $aref->[$i+1][0] - $aref->[$i][0]; 
        my $e = $aref->[$i+1][1] - $aref->[$i][1]; 
        print "EAST: $e\n"; 
        print "NORTH: $n\n"; 
         
        my $dAlong =  sqrt($n**2 + $e**2); 
        print "dAlong: $dAlong\n"; 
         
        my $area = .5 * abs($n * $e); 
         
        # Area = 1/2 * b * h 
        # h = 2 * Area / b 
        my $dAcross = 2 * $area / $dAlong; 
        print "dAcross: $dAcross\n"; 
         
        # arcsin in degrees 
        my $arcsin = asin($n/$dAlong) * 180 / pi; 
         
        my $azimuth = 90 - $arcsin; 
        print "azimuth: $azimuth\n"; 
        print $/; 
    } 
} 
$dbh->disconnect;

 
C:\Old_Data\perlp>perl guru_fetch.pl -t Sample3 -i 3,7800 
        3 
EAST: 6.5 
NORTH: -4.40000000037253 
dAlong: 7.84920378148499 
dAcross: 3.64368167760968 
azimuth: 124.094977957017 
 
EAST: 7.20000000001164 
NORTH: -4.5 
dAlong: 8.49058301886082 
dAcross: 3.81599236802463 
azimuth: 122.005383208042 
 
        7800 
EAST: 5.79999999998836 
NORTH: -4.39999999944121 
dAlong: 7.28010988893352 
dAcross: 3.50544159168541 
azimuth: 127.184706449785 
 
EAST: 6.40000000002328 
NORTH: -4.40000000037253 
dAlong: 7.76659513580927 
dAcross: 3.6257844666899 
azimuth: 124.508522989836