<?php if (!defined('PmWiki')) exit();
/*
Use http://www.topomap.co.nz/ to display an excerpt from a map in PmWiki
+
  Copyright 2013-2014 Simon Davis
  This text is written for PmWiki; you can redistribute it and/or modify it 
  under the terms of the GNU General Public License as published by the 
  Free Software Foundation; either version 3 of the License, or (at your option)
  any later version. See pmwiki.php for full details and lack of warranty.
+
# 2014-08-09 Markup_e PHP 5.5 compatible, add mapref parameter
# 2013-07-15 Initial version
*/
# Version date
$RecipeInfo['NZTopo']['Version'] = '2014-08-10';
$FmtPV['$NZTopoVersion'] = "'NZTopo version {$RecipeInfo['NZTopo']['Version']}'"; // return version as a custom page variable

SDV($NZTopoZoom, 12); # set default zoom factor
SDV($NZTopoNewwin, 1); # set default new window setting
SDV($NZTopoDebug, true); # set default debug setting

# declare $NZTopo for (:if enabled NZTopo:) recipe installation check
global $NZTopo; $NZTopo = 1;

## Add a custom markup 
# (:nztopo ll=-41.293722,174.871482 :)
# (:nztopo topo50=BP33912708,BP33876687 height=300 width=400 pin=1 label='destination' float=right clear=both zoom=12:)
# (:nztopo mapref=BN33991890:)
# directive arguments are
#  ll= -- decimal latitude,longitude
#  llbs = -- -- decimal latitude,longitude;latitude,longitude
  $vlat     = '[-+]?\d{1,2}[.]\d+'; # -90 .. 90
  $vlong    = '[-+]?[1]?\d{1,2}[.]\d+'; # -180 .. 180
  $vlatlong = $vlat . '[,]' . $vlong; # latitude,longitude
  $pll      = 'll=' . $vlatlong;
  $pllbs    = 'llbs=' . $vlatlong . '(?:[;]' . $vlatlong . ')+'; # two or more pairs of lat long co-ordinates
# nztm
  $vnztm    = '\d{7}(?:[.]\d{0,3})?'; # single NZTM co-ordinate
  $vnztm2   = $vnztm . '[,]' . $vnztm; # pair of NZTM co-ordinates, easting, northing
# mapref
  $mapref   ='[ABC][A-Z][0-4]\d{7}'; # topo50 map reference e.g. BN33991890
## parameters
  $pnzne    = 'nzne=' . $vnztm2;
  $pnzbs    = 'nzbs=' . $vnztm2 . '(?:[;]' . $vnztm2 . ')+'; # two or more pairs of NZTM co-ordinates
# topo50= --  topo50 grid coordinate, or two grid coordinates (not currently implemented in nztopomap)
  $ptopo50  = 'topo50=' . $mapref . '(?:[,]' . $mapref .')?';
# -- only one of topo50, ll, llbs, nzne, nzbs,, kml or gpx can be supplied
  $vscheme  = '(?:https?|ftp):\/\/'; # 
  $vurl     = '(?:(?:\'' . $vscheme . '[^\']+\')|(?:"' . $vscheme . '[^"]+"))'; # provide for single and double quoted strings
# kml= -- URL to kml file
  $pkml     = 'kml=' . $vurl;
# kml= -- URL to kml file
  $pgpx     = 'gpx=' . $vurl;
# map reference
  $pmapref  = 'mapref=' . $mapref;
#  height= -- image height in pixels
  $pheight  = 'height=\d{1,5}';
#  width= -- image width in pixels
  $pwidth   = 'width=\d{1,5}';
#  pin= -- show pin
  $ppin     = 'pin=[01]';
#  label -- tool tip label for pin
  $plabel   = 'label=(?:(?:\'[^\']+\')|(?:"[^"]+")|(?:[\w\s]+))'; # provide for single and double quoted strings, and unquoted string
# zoom= -- scale factor for map
  $pzoom    = 'zoom=\d{1,2}';
# float= -- left or right
  $pfloat   = 'float=(?:left|right)';
# clear= -- left, right, both
  $pclear   = 'clear=(?:left|right|both)';
# only one location can be supplied
  $qlocale  = '(?:' . $pll . ')|(?:' . $pllbs . ')|(?:' . $pkml . ')|(?:' . $pgpx . ')|(?:' . $pnzne . ')|(?:' . $pnzbs . ')|(?:' . $pmapref . ')';
# display modifications
  $qmods    = "(" . $pheight . ")\s*|(" . $pwidth . ")\s*|(" . $ppin . ")\s*|(" . $plabel . ")\s*|(" . $pzoom . ")\s*|(" . $pfloat . ")\s*|(" . $pclear .  ")\s*";
#
Markup_e( 'NZTopo', 
  'directives',
  "/\\(:nztopo (" . $qlocale . ")\s*(?:" . $qmods . "){0,7}\s*:\\)/i", # 
  "Keep(NZTopo_Parse(\$m[1], \$m[2], \$m[3], \$m[4], \$m[5], \$m[6], \$m[7], \$m[8], \$m[9]))" );
# s = dot matches all chars including newline
# i = case insensitive
# m = multiline
# uses lazy evaluation, preserves leading and trailing white space
# Keep prevents PmWiki markup being applied
#
  $debugon = $NZTopoDebug; # debug
/** Main NZTopo parser
 *   /param   arguements as documented above
 *   /return  The HTML-formatted NZTopo markup wrapped in a <div> of class "nztopo", see http://www.topomap.co.nz/.
 */
function NZTopo_Parse($p1, $p2, $p3, $p4, $p5, $p6, $p7, $p8, $p9) {
#
global $NZTopoZoom, $NZTopoNewwin;
global $debugon; # for debug
// Initialise variables
  $retval  = ''; # return value
  $opt     = ParseArgs($p1 . ' ' . $p2 . ' ' . $p3 . ' ' . $p4 . ' ' . $p5 . ' ' . $p6 . ' ' . $p7 . ' ' . $p8 . ' ' . $p9);
  $version = 2; # fixed version of nztopomap interface
  $newwin  = $NZTopoNewwin; # open map in new window
  $ppcm    = 80; # 80 pixels per centrimetre to provide base scale for maps
  $width   = $opt['width']; # iframe parameter
  $height  = $opt['height']; # iframe parameter
  $float   = $opt['float']  ? 'float:' . $opt['float'] . ';' : ""; # div parameter
  $clear   = $opt['clear']  ? 'clear:' . $opt['clear'] . ';' : ""; # div parameter
  $pdiv    = $clear . $float ? ' style="' . $clear . $float . '"' : ""; # div positioning

  $debugval = $debugon ? '<div style="font-size:smaller;' . $float . $clear . '">' : '';
  
// query string parameters
  $bllbs = (boolean) $opt['llbs'];
  if( $bllbs ) { # get size of map in km
    list ($widthkm, $heightkm, $dbgval) = llbskm($opt['llbs']);
    $debugval .= $dbgval;
    $width  = max ($widthkm * $ppcm, $width); # override height
    $height = max ($heightkm * $ppcm, $width); # override height
  }
  $bnzbs = (boolean) $opt['nzbs'];
  if( $bnzbs ) { # get size of map in km
    list ($widthkm, $heightkm, $dbgval) = nzbskm($opt['nzbs']);
    $debugval .= $dbgval;
    $width  = max ($widthkm * $ppcm, $width); # override height
    $height = max ($heightkm * $ppcm, $width); # override height
  }
  $query  = ''; 
  $query .= '&new=' . $newwin;
  $query .= $opt['ll']     ? '&ll='     . $opt['ll']     : "";
  $query .= $opt['llbs']   ? '&llbs='   . $opt['llbs']   : "";
  $query .= $opt['topo50'] ? '&topo50=' . $opt['topo50'] : "";
  $query .= $opt['nzbs']   ? '&nzbs='   . $opt['nzbs']   : "";
  $query .= $opt['nzne']   ? '&nzne='   . $opt['nzne']   : "";
  $query .= $opt['kml']    ? '&kml='    . rawurlencode($opt['kml']) : "";
  $query .= $opt['gpx']    ? '&gpx='    . rawurlencode($opt['gpx']) : "";
  $query .= $opt['mapref'] ? '&mapref=' . $opt['mapref'] : "";
  $query .= $opt['pin']    ? '&pin='    . $opt['pin']    : "";
  $query .= $opt['label']  ? '&lbl='    . rawurlencode(htmlentities($opt['label'])) : ""; # double url encode unicode chars for them to work with the topo map
  $query .= $opt['zoom']   ? '&z='      . $opt['zoom']   : '&z=' . $NZTopoZoom;
  $awidth  = $width  ? ' width="'  . $width  . '"' : ""; # iframe parameter
  $aheight = $height ? ' height="' . $height . '"' : ""; # iframe parameter
  $src    = 'http://www.topomap.co.nz/NZTopoMapEmbedded?&v=' . $version . $query;

  $retval .= '<div class="nztopo"' . $pdiv . '>' . "\n";
  $retval .= '<iframe seamless frameborder="0" scrolling="no" marginheight="0" marginwidth="0"' . $awidth . $aheight . ' src="' . $src . '">';
  $retval .= "\n" . '</iframe>';
  $retval .= '</div>' . "\n";
/* example HTML for nztopomap
<iframe width="300" height="300" frameborder="0" scrolling="no" marginheight="0" marginwidth="0" src="http://www.topomap.co.nz/NZTopoMapEmbedded?v=2&ll=-41.293722,174.871482&z=15&pin=1"></iframe>
*/
  if( $debugon ) {
    $debugval .= 'p1=' . $p1 . ' p2=' . $p2 . ' p3=' . $p3 . ' p4=' . $p4 . ' p5=' . $p5 . ' p6=' . $p6 . ' p7=' . $p7 . ' p8=' . $p8 . ' p9=' . $p9 . '</br>' . "\n";
    $debugval .= 'src="' . $version . $query . '"</div>' . "\n";
  }
  return $retval . $debugval;
} # NZTopo_Parse

function llbskm($llbs) {
  global $debugon; # for debug
  $dbgval = '';
  $latlngs = preg_split ('/[,;]/', $llbs);
  $lats = array(); # eastings
  $lngs = array(); # northings
  foreach ($latlngs as $k => $v) {
    $k % 2 == 0 ? $lats [] = $v : $lngs [] = $v;
  }
  $maxlat = max ($lats);
  $maxlng = max ($lngs);
  $minlat = min ($lats);
  $minlng = min ($lngs);
/* rectangle coordinates are (lat/long)
       left:      right:
  top: max/min    max/max
  bot: min/min    min/max
*/
  $width1  = distance ($maxlat, $minlng, $maxlat, $maxlng);
  $width2  = distance ($minlat, $minlng, $minlat, $maxlng);
  $height1 = distance ($maxlat, $minlng, $minlat, $minlng);
  $height2 = distance ($maxlat, $maxlng, $minlat, $maxlng);
  $widthkm  = ceil (max ($width1, $width2) * 10) / 10; # take largest width and round up to km
  $heightkm = ceil (max ($height1, $height2) * 10) / 10; # take largest height and round up to km
  if( $debugon ) {
    $dbgval .= 'lats="' . var_export ($lats, TRUE) . '"</br>' . "\n";
    $dbgval .= 'lngs="' . implode ('", "', $lngs) . '"</br>' . "\n";
    $dbgval .= 'minlat=' . $minlat . ', maxlat=' . $maxlat . ', minlng=' . $minlng . ', maxlng=' . $maxlng . '</br>'  . "\n";
    $dbgval .= 'w1=' . $width1 . ', w2=' . $width2 . ', h1=' . $height1 . ', h2=' . $height2 . ', w=' . $widthkm . 'km, h=' . $heightkm . 'km</br>'  . "\n";
  }
  return array ($widthkm, $heightkm, $dbgval);
}

function nzbskm($nzbs) {
  global $debugon; # for debug
  $dbgval = '';
  $nthsests = preg_split ('/[,;]/', $nzbs); # http://www.linz.govt.nz/topography/topo-maps/topo50/prototype/map-reading-references
  $nths = array(); # northings
  $ests = array(); # eastings
  foreach ($nthsests as $k => $v) {
    $k % 2 == 0 ? $nths [] = $v : $ests [] = $v;
  }
  $maxnth = max ($nths);
  $maxest = max ($ests);
  $minnth = min ($nths);
  $minest = min ($ests);
/* rectangle coordinates are (nth/est)
       left:      right:
  top: max/min    max/max
  bot: min/min    min/max
*/
  $width1  = ($maxest - $minest) / 1000; # metres to km
  $height1 = ($maxnth - $minnth) / 1000; # metres to km
  $widthkm  = ceil ($width1 * 10) / 10; # round up to 0.1 km
  $heightkm = ceil ($height1 * 10) / 10; # round up to 0.1 km
  if( $debugon ) {
    $dbgval .= 'nthsests="' . implode ('", "', $nthsests) . '"</br>' . "\n";
    $dbgval .= 'ests="' . var_export ($ests, TRUE) . '"</br>' . "\n";
    $dbgval .= 'nths="' . implode ('", "', $nths) . '"</br>' . "\n";
    $dbgval .= 'minest=' . $minest . ', maxest=' . $maxest . ', minnth=' . $minnth . ', maxnth=' . $maxnth . '</br>'  . "\n";
    $dbgval .= 'w1=' . $width1 . ', h1=' . $height1 . ', w=' . $widthkm . 'km, h=' . $heightkm . 'km</br>'  . "\n";
  }
  return array ($widthkm, $heightkm, $dbgval);
}

# http://snipplr.com/view/2531/
function distance($lat1, $lng1, $lat2, $lng2)
{
	$pi80 = M_PI / 180;
	$lat1 *= $pi80;
	$lng1 *= $pi80;
	$lat2 *= $pi80;
	$lng2 *= $pi80;
 
	$r = 6372.797; // mean radius of Earth in km
	$dlat = $lat2 - $lat1;
	$dlng = $lng2 - $lng1;
	$a = sin($dlat / 2) * sin($dlat / 2) + cos($lat1) * cos($lat2) * sin($dlng / 2) * sin($dlng / 2);
	$c = 2 * atan2(sqrt($a), sqrt(1 - $a));
	$km = $r * $c;
 
	return $km;
} # distance