#!/usr/bin/perl
#
# FixMyStreet:Map::OSM
# OSM maps on FixMyStreet.
#
# Copyright (c) 2010 UK Citizens Online Democracy. All rights reserved.
# Email: matthew@mysociety.org; WWW: http://www.mysociety.org/
package FixMyStreet::Map::OSM;
use strict;
use Math::Trig;
use mySociety::Web qw(ent NewURL);
use Utils;
sub header_js {
return '
';
}
sub map_type {
return 'OpenLayers.Layer.OSM.Mapnik';
}
# display_map Q PARAMS
# PARAMS include:
# latitude, longitude for the centre point of the map
# TYPE is 1 if the map is clickable, 2 if clickable and has a form upload,
# 0 if not clickable
# PINS is array of pins to show, location and colour
# PRE/POST are HTML to show above/below map
sub display_map {
my ($self, $q, %params) = @_;
$params{pre} ||= '';
# Map centre may be overridden in the query string
$params{latitude} = Utils::truncate_coordinate($q->param('lat')+0)
if defined $q->param('lat');
$params{longitude} = Utils::truncate_coordinate($q->param('lon')+0)
if defined $q->param('lon');
my $zoom = defined $q->param('zoom') ? $q->param('zoom') : 2;
my $zoom_act = 14 + $zoom;
my ($x_tile, $y_tile) = latlon_to_tile_with_adjust($params{latitude}, $params{longitude}, $zoom_act);
my $tl = ($x_tile-1) . "/" . ($y_tile-1);
my $tr = "$x_tile/" . ($y_tile-1);
my $bl = ($x_tile-1) . "/$y_tile";
my $br = "$x_tile/$y_tile";
my $tl_src = "http://a.tile.openstreetmap.org/$zoom_act/$tl.png";
my $tr_src = "http://b.tile.openstreetmap.org/$zoom_act/$tr.png";
my $bl_src = "http://c.tile.openstreetmap.org/$zoom_act/$bl.png";
my $br_src = "http://tile.openstreetmap.org/$zoom_act/$br.png";
map { s{/}{.} } ($tl, $tr, $bl, $br);
my @pins;
my $pins = '';
foreach my $pin (@{$params{pins}}) {
$pin->[3] ||= '';
push @pins, "[ $pin->[0], $pin->[1], '$pin->[2]', '$pin->[3]' ]";
$pins .= display_pin($q, $pin, $x_tile, $y_tile, $zoom_act);
}
my $pins_js = join(",\n", @pins);
my $img_type;
if ($params{type}) {
$img_type = 'OpenStreetMap and contributors, CC-BY-SA');
my $compass = compass($q, $x_tile, $y_tile, $zoom);
my $map_type = $self->map_type();
$out .= <
$params{pre}
$copyright
EOF
return $out;
}
sub display_pin {
my ($q, $pin, $x_tile, $y_tile, $zoom) = @_;
my ($px, $py) = latlon_to_px($pin->[0], $pin->[1], $x_tile, $y_tile, $zoom);
my $num = '';
my $host = Page::base_url_with_lang($q, undef);
my %cols = (red=>'R', green=>'G', blue=>'B', purple=>'P');
my $out = '';
return $out unless $pin->[3];
my $cobrand = Page::get_cobrand($q);
my $url = Cobrand::url($cobrand, NewURL($q, -url => '/report/' . $pin->[3]), $q);
# XXX Would like to include title here in title=""
$out = '' . $out . '';
return $out;
}
# Given a lat/lon, convert it to OSM tile co-ordinates (precise).
sub latlon_to_tile($$$) {
my ($lat, $lon, $zoom) = @_;
my $x_tile = ($lon + 180) / 360 * 2**$zoom;
my $y_tile = (1 - log(tan(deg2rad($lat)) + sec(deg2rad($lat))) / pi) / 2 * 2**$zoom;
return ( $x_tile, $y_tile );
}
# Given a lat/lon, convert it to OSM tile co-ordinates (nearest actual tile,
# adjusted so the point will be near the centre of a 2x2 tiled map).
sub latlon_to_tile_with_adjust($$$) {
my ($lat, $lon, $zoom) = @_;
my ($x_tile, $y_tile) = latlon_to_tile($lat, $lon, $zoom);
# Try and have point near centre of map
if ($x_tile - int($x_tile) > 0.5) {
$x_tile += 1;
}
if ($y_tile - int($y_tile) > 0.5) {
$y_tile += 1;
}
return ( int($x_tile), int($y_tile) );
}
sub tile_to_latlon {
my ($x, $y, $zoom) = @_;
my $n = 2 ** $zoom;
my $lon = $x / $n * 360 - 180;
my $lat = rad2deg(atan(sinh(pi * (1 - 2 * $y / $n))));
return ( $lat, $lon );
}
# Given a lat/lon, convert it to pixel co-ordinates from the top left of the map
sub latlon_to_px($$$$$) {
my ($lat, $lon, $x_tile, $y_tile, $zoom) = @_;
my ($pin_x_tile, $pin_y_tile) = latlon_to_tile($lat, $lon, $zoom);
my $pin_x = tile_to_px($pin_x_tile, $x_tile);
my $pin_y = tile_to_px($pin_y_tile, $y_tile);
return ($pin_x, $pin_y);
}
# Convert tile co-ordinates to pixel co-ordinates from top left of map
# C is centre tile reference of displayed map
sub tile_to_px {
my ($p, $c) = @_;
$p = 256 * ($p - $c + 1);
$p = int($p + .5 * ($p <=> 0));
return $p;
}
sub click_to_tile {
my ($pin_tile, $pin) = @_;
$pin -= 256 while $pin > 256;
$pin += 256 while $pin < 0;
return $pin_tile + $pin / 256;
}
# Given some click co-ords (the tile they were on, and where in the
# tile they were), convert to WGS84 and return.
sub click_to_wgs84 {
my ($self, $q, $pin_tile_x, $pin_x, $pin_tile_y, $pin_y) = @_;
my $tile_x = click_to_tile($pin_tile_x, $pin_x);
my $tile_y = click_to_tile($pin_tile_y, $pin_y);
my $zoom = 14 + (defined $q->param('zoom') ? $q->param('zoom') : 2);
my ($lat, $lon) = tile_to_latlon($tile_x, $tile_y, $zoom);
return ( $lat, $lon );
}
sub compass ($$$$) {
my ( $q, $x, $y, $zoom ) = @_;
my ($lat, $lon) = map { Utils::truncate_coordinate($_) } tile_to_latlon($x, $y-1, $zoom+14);
my $north = NewURL( $q, lat => $lat, lon => $lon, zoom => $zoom );
($lat, $lon) = map { Utils::truncate_coordinate($_) } tile_to_latlon($x, $y+1, $zoom+14);
my $south = NewURL( $q, lat => $lat, lon => $lon, zoom => $zoom );
($lat, $lon) = map { Utils::truncate_coordinate($_) } tile_to_latlon($x-1, $y, $zoom+14);
my $west = NewURL( $q, lat => $lat, lon => $lon, zoom => $zoom );
($lat, $lon) = map { Utils::truncate_coordinate($_) } tile_to_latlon($x+1, $y, $zoom+14);
my $east = NewURL( $q, lat => $lat, lon => $lon, zoom => $zoom );
($lat, $lon) = map { Utils::truncate_coordinate($_) } tile_to_latlon($x, $y, $zoom+14);
my $zoom_in = $zoom < 3 ? NewURL( $q, lat => $lat, lon => $lon, zoom => $zoom+1 ) : '#';
my $zoom_out = $zoom > 0 ? NewURL( $q, lat => $lat, lon => $lon, zoom => $zoom-1 ) : '#';
my $world = NewURL( $q, lat => $lat, lon => $lon, zoom => 0 );
#my $host = Page::base_url_with_lang( $q, undef );
my $dir = "/jslib/OpenLayers-2.10/img";
return <