#!/usr/bin/perl
$FILE= "ShoreLines-14959.dat"; # http://shoreline.noaa.gov/data/datasheets/wvs.html -> WCL 1:5 000 000

sub transfo {
	my ($lon, $lat) = @_;

	my $lo = ($lon + 180);
	my $la = (90 - $lat);
	return("$lo $la");
	}
print "<!DOCTYPE svg>\n";
print "<svg xmlns=\"http://www.w3.org/2000/svg\" version=\"1.1\" width=\"1080\" height=\"540\" viewBox=\"0 0 360 180\">\n";

open FILE or die "Can't open $FILE\n";
$eop = ""; # end of polyline
while ($line = <FILE>)
	{
	if ($line =~ /^>/)
		{ # end previous polyline and start a new one
		print $eop;
		$eol = "";
		print "\n<polyline points=\"";
		$eop = "\" style=\"stroke:black; stroke-width:.1; fill:none;\"/>\n";
		}
	else	{ # put new point
		chomp $line;
		my ($lon, $lat) = split /\t/,$line;
		print $eol . transfo($lon, $lat);
		$eol = ",\n";
		}
	}
print "$eop";
for ($i = 0; $i < 180; $i += 10)
	{ # draw parallels
	print "<line x1=\"0\" y1=\"$i\" x2=\"360\" y2=\"$i\" style=\"stroke:blue; stroke-width:.1; fill:none;\"/>\n";
	}
# draw Equator with double width
print "<line x1=\"0\" y1=\"90\" x2=\"360\" y2=\"90\" style=\"stroke:blue; stroke-width:.2; fill:none;\"/>\n";
for ($i = 0; $i < 360; $i += 15)
	{ # Draw meridians
	print "<line x1=\"$i\" y1=\"0\" x2=\"$i\" y2=\"180\" style=\"stroke:green; stroke-width:.1; fill:none;\"/>\n";
	}
# Draw Prime meridian with double width
print "<line x1=\"180\" y1=\"0\" x2=\"180\" y2=\"180\" style=\"stroke:green; stroke-width:.2; fill:none;\"/>\n";
print "</svg>\n";