1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
|
#!/usr/bin/perl -w
use strict;
use lib qw(./lib ../lib);
use Test;
plan tests => 46;
# grid->lat/long test for GB module
# $Revision: 1.2 $
use Geography::NationalGrid;
ok(1);
my $point;
# 1. worked example from online refs
$point = new Geography::NationalGrid( 'GB',
Easting => 651409.903,
Northing => 313177.270,
Userdata => { Name => 'Worked Example from east/north' },
);
ok( qeq($point->latitude, [52, 39, 27.25312]) );
ok( qeq($point->longitude, [1, 43, 4.51772]) );
ok( $point->gridReference() eq 'TG 514131' );
ok( $point->gridReference(1) eq 'TG 5140913177' );
ok( $point->gridReference(10) eq 'TG 51401317' );
ok( $point->gridReference(100) eq 'TG 514131' );
ok( $point->gridReference(1000) eq 'TG 5113' );
# 2. broadcasting house
comp( 'TQ 289006', '50d 47m 24.19s', '-0d 10m 15.69s');
# 3. origin
comp( 'sv000000', '49d 45m 58.27s', '-7d 33m 23.21s' );
# 4 Llanfair
$point = new Geography::NationalGrid( 'GB',
Easting => 300100,
Northing => 171600,
);
ok( qeq($point->latitude, '51d 26m 2.73s') );
ok( qeq($point->longitude, '-3d 26m 14.34s') );
comp( 'TQ 389773', '51d 28m 38.40s', '0d 0m 0s' ); # Greenwich
comp( 'NT 261745', '55d 57m 28.00s', '-3d 11m 2.00s' ); # Edinburgh
comp( 'SP 510071', '51d 45m 35.40s', '-1d 15m 39.00s' ); #Oxford
comp( 'TL 431595', '52d 12m 51.60s', '0d 5m 42.00s' ); #Cambridge
comp( 'NZ 269416', '54d 46m 6.20s', '-1d 34m 56.00s' ); #Durham
comp( 'SJ 287900', '53d 24m 5.40s', '-3d 4m 20.50s' ); #Liverpool
# arbitrary points follow
comp('TV 500998','50d 46m 40.03s', '0d 7m 40.20s');
comp('SM 600300', '51d 54m 58.61s', '-5d 29m 25.12s' );
comp('SC 200850', '54d 13m 43.90s', '-4d 45m 42.43s' );
comp('NR 700150', '55d 22m 29.85s', '-5d 37m 50.36s');
comp('NB 050050', '57d 56m 11.91s', '-6d 59m 6.71s' );
comp('HP 600010', '60d 41m 17.87s', '-0d 54m 4.53s' );
comp('NK 002543', '57d 34m 44.37s', '-1d 59m 47.96s' );
comp('NZ 003005', '54d 23m 59.37s', '-1d 59m 43.36s' );
comp('SP 500500', '52d 8m 44.21s', '-1d 16m 9.21s' );
comp('TG 450023', '52d 33m 46.45s ', '1d 36m 55.01s' );
sub qeq {
my $rv = 0;
my $i = Geography::NationalGrid->deg2rad($_[0]);
my $j = Geography::NationalGrid->deg2rad($_[1]);
if (abs($i - $j) < 0.000007) {
$rv = 1;
} else {
warn "# got $i, expecting $j, difference " . ($i - $j);
}
return $rv;
}
sub comp {
my ($grid, $lat, $long) = @_;
my $point = new Geography::NationalGrid( 'GB',
GridReference => $grid
);
ok( qeq($point->latitude, $lat) );
ok( qeq($point->longitude, $long) );
}
|
