<?php if (!defined('PmWiki')) exit();
/* NZTopo 
Use http://www.topomap.co.nz/ API to display an excerpt from a map in PmWiki
See http://www.pmwiki.org/wiki/Cookbook/NZTopo and http://kiwiwiki.co.nz/pmwiki/pmwiki.php/Cookbook/NZTopo
+
  Copyright 2013-2015 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.
+
# 2015-07-01 Attachments corrected (gpx, kml requires extra processing)
# 2014-08-09 Markup_e PHP 5.5 compatible, add mapref parameter
# 2013-07-15 Initial version
*/
# Version date
$RecipeInfo['NZTopo']['Version'] = '2015-07-01';
$FmtPV['$NZTopoVersion'] = "'NZTopo version {$RecipeInfo['NZTopo']['Version']}'"; // return version as a custom page variable

SDV($NZTopoZoom, 13); # set default zoom factor
SDV($NZTopoNewwin, 1); # set default new window setting
SDV($NZTopoDebug, false); # set default debug setting

# declare $NZTopo for (:if enabled NZTopo:) recipe installation check
global $NZTopo; $NZTopo = 1;
# set debug flag
  $debugon = $NZTopoDebug; # if on writes input and output to web page

## Add a PmWiki 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:)
## the following builds the regex to parse the NZTopo PmWiki directive input parameters
# 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?:\/\/'; # 
  $vurl     = "(?:(?:'" . $vscheme . "[^']+')" . # provide for single and double quoted URL strings, or no quotes 
                '|(?:"' . $vscheme . '[^"]+")' . 
                '|(?:'  . $vscheme . '[^.?$#\s].[^\s\r\n]*))'; #from https://mathiasbynens.be/demo/url-regex
#   kml= -- URL to kml file
  $pkml     = 'kml=' . $vurl;
  if( $debugon ) { # add extra test parameters
    $pkml2  = 'kml2=' . $vurl; # added for debug and testing
    $pkml3  = 'kml3=' . $vurl; # added for debug and testing
    $pkml4  = 'kml4=' . $vurl; # added for debug and testing
  } # debugon
#   gpx= -- URL to gpx 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-+!?()]+))'; # provide for single and double quoted strings, and unquoted string (no spaces)
#   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 . ')';
  if( $debugon ) { # add extra test parameters
	$qlocale .= '|(?:' . $pkml2 . ')|(?:' . $pkml3 . ')|(?:' . $pkml4 . ')';
  } # debugon
# display modifications
  $qmods    = "(" . $pheight . ")\s*|(" . $pwidth . ")\s*|(" . $ppin . ")\s*|(" . $plabel . ")\s*|(" . $pzoom . ")\s*|(" . $pfloat . ")\s*|(" . $pclear .  ")\s*";
##
Markup_e( 'NZTopo', # an internal PmWiki function that defines the custom markup for the wiki (see http://www.pmwiki.org/wiki/PmWiki/CustomMarkup)
  '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
#
/** Main NZTopo parser
 *   /param   arguments 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; # configuration variables from configuration file, e.g. config.php
global $debugon; # for debug
// Initialise variables
  $retval    = ''; # return value
  $opt       = ParseArgs($p1 . ' ' . $p2 . ' ' . $p3 . ' ' . $p4 . ' ' . $p5 . ' ' . $p6 . ' ' . $p7 . ' ' . $p8 . ' ' . $p9); # see pmWiki documentation
  $src_url   = 'http://www.topomap.co.nz/NZTopoMap'; # URL for browser service
  $src_urlemb = 'http://www.topomap.co.nz/NZTopoMapEmbedded'; # URL for embedded service
  $src_proxy = 'http://www.topomap.co.nz/proxy.ashx?'; # URL for gpx proxy
  $version   = 'v=2'; # fixed version of nztopomap interface
  $newwin    = $NZTopoNewwin; # open map in new window
  //
  ## need dpi to calculate size of map to ask for, have not been told what topomap generates
  $ppcm      = 80; # 80 pixels per centimetre to provide base scale for maps (guess)
  //
  $width     = $opt['width']; # iframe parameter
  $height    = $opt['height']; # iframe parameter
  $float     = (bool) $opt['float']  ? 'float:' . $opt['float'] . ';' : ''; # div parameter
  $clear     = (bool) $opt['clear']  ? 'clear:' . $opt['clear'] . ';' : ''; # div parameter
  $pdiv      = (bool) $clear . $float ? ' style="' . $clear . $float . '"' : ''; # div positioning

  $debugval  = $debugon ? '<div style="font-size:smaller;' . $float . $clear . '">' : ''; # initialise
  $dbgval    = ''; # initialise
  
## 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
  }
## build query string (i.e. topomap API)
  $query  = ''; # initialise
  $query .= (bool) $opt['ll']     ? '&ll='     . $opt['ll']     : '';
  $query .= (bool) $opt['llbs']   ? '&llbs='   . $opt['llbs']   : '';
  $query .= (bool) $opt['topo50'] ? '&topo50=' . $opt['topo50'] : '';
  $query .= (bool) $opt['nzbs']   ? '&nzbs='   . $opt['nzbs']   : '';
  $query .= (bool) $opt['nzne']   ? '&nzne='   . $opt['nzne']   : '';
  // 
  # trim URLs for spaces and quotes; rawurldecode URLs before rawurlencode them for nztopomap
  # rawurldecode resets any pasted url to a 'known' state
//$query .= (bool) $opt['kml']    ? '&kml='    . rawurlencode(htmlentities(rawurldecode(trim($opt['kml'], ' \'\"')))) : ''; # does not comply with topomap API
  $query .= (bool) $opt['kml']    ? '&kml='    . encode_path($opt['kml'], $dbgval) : ''; # call function to comply with topomap API
  $debugval .= (bool) $opt['kml'] ?  $dbgval : ''; # debug output
  if( $debugon ) { # apply different combinations of encoding for match with topomap API
  # kml filename double rawurlencode (for debug and testing only)
  $query .= (bool) $opt['kml2']   ? '&kml='    . rawurlencode(rawurlencode(htmlentities(rawurldecode(trim($opt['kml2'], ' \'\"'))))) : '';
  # kml filename single rawurlencode with no rawurldecode (for debug and testing only)
  $query .= (bool) $opt['kml3']   ? '&kml='    . rawurlencode(htmlentities(trim($opt['kml3'], ' \'\"'))) : '';
  # kml filename double rawurlencode with no rawurldecode (for debug and testing only)
  $query .= (bool) $opt['kml4']   ? '&kml='    . rawurlencode(rawurlencode(htmlentities(trim($opt['kml4'], ' \'\"')))) : '';
  } # debugon
  // note gpx works fine
  $query .= (bool) $opt['gpx']    ? '&gpx='    . rawurlencode($src_proxy . htmlentities(rawurldecode(trim($opt['gpx'], ' \'\"')))) : '';
  //
  $query .= (bool) $opt['mapref'] ? '&mapref=' . $opt['mapref'] : '';
  $query .= (bool) $opt['pin']    ? '&pin='    . $opt['pin']    : '';
  $query .= (bool) $opt['label']  ? '&lbl='    . htmlentities(trim($opt['label'])) : ''; # trin spaces, encode special chars for them to work with the topo map
  $query .= (bool) $opt['zoom']   ? '&z='      . $opt['zoom']   : '&z=' . $NZTopoZoom;
  $query .= '&new=' . $newwin; # set open in new window option
  $awidth  = (bool) $width  ? ' width="'  . $width  . '"' : ''; # iframe parameter
  $aheight = (bool) $height ? ' height="' . $height . '"' : ''; # iframe parameter
  $src    = $src_urlemb . '?' . $version . $query;

  $retval .= '<div class="nztopo"' . $pdiv . '>' . "\n";
  $retval .= '<iframe class="nztopoframe" seamless frameborder="0" scrolling="no" marginheight="0" marginwidth="0"' . $awidth . $aheight . ' src="' . $src . '">';
  $retval .= "\n" . '</iframe>';
  $retval .= '</div>' . "\n";
/* example query string for nztopomap http://www.topomap.co.nz/NZTopoMap or http://www.topomap.co.nz/NZTopoMapEmbedded
?v=2&ll=-41.293722,174.871482&z=15&pin=1
?v=2&kml=http%3A%2F%2Fkiwiwiki.co.nz%2Fpmwiki%2Fuploads%2FTest%2FNZTopo-Attach%2F29%2520May%25202015%252020_01_35.kml
?v=2&gpx=http%3A%2F%2Fwww.topomap.co.nz%2Fproxy.ashx%3Fhttp%3A%2F%2Fkiwiwiki.co.nz%2Fpmwiki%2Fuploads%2FTest%2FNZTopo-Attach%2F20150531.gpx
*/
  if( $debugon ) { # display inputs and outputs to wiki page
    $debugval .= '<i>input</i>: p1=' . $p1 . ' p2=' . $p2 . ' p3=' . $p3 . ' p4=' . $p4 . ' p5=' . $p5 . ' p6=' . $p6 . ' p7=' . $p7 . ' p8=' . $p8 . ' p9=' . $p9 . '</br>' . "\n";
    $debugval .= '<i>output</i>: src="' . '<a href="' . $src_url . '?' . $version . $query . '" target="_blank">' . $version . $query . '</a>';
    $debugval .= '"</div>' . "\n";
  } # debugon
  return $retval . $debugval;
} # NZTopo_Parse

function encode_path($kml_url, &$dbgval) {
  global $debugon; # for debug
# for the kml parameter the topomap API is inconsistent. The kml path is urlencoded separately to the full url string
  // specifically rawurlencode($src_proxy . rawurldecode(trim(filename.gpx, ' \'\"'))) works,
  // but          rawurlencode(            rawurldecode(trim(filename.kml, ' \'\"'))) does not
	$temp_url = # decode for parse_url and to start with known state (user may for example provide spaces in the URL string)
	    htmlentities(         # protect against injection attacks
		rawurldecode(trim($kml_url, ' \'\"'))); # trim single and double quotes and spaces
	$parsed_url = parse_url($temp_url); # extract url components
	$dbgval .= $debugon ? ' path="' . $parsed_url['path'] . '" </br>' . "\n" : ''; # debug output
	$parsed_url['path'] = rawurlencode(substr($parsed_url['path'], 1)); # encode path separately again,  dropping leading '/'
	$dbgval .= $debugon ? ' enc_path="' . $parsed_url['path'] . '" </br>' . "\n" : ''; # debug output
	$enc_url = ''; # reconstruct URL
	$enc_url .= (bool) $parsed_url['scheme'] ? $parsed_url['scheme'] . '://' : ''; # aka http or https
	$enc_url .= (bool) $parsed_url['host']   ? $parsed_url['host']           : ''; # domain name
	$enc_url .= (bool) $parsed_url['port']   ? ':' . $parsed_url['port']     : ''; # port
	$enc_url .= (bool) $parsed_url['path']   ? '/' . $parsed_url['path']     : ''; # add '/' dropped before as topomap considers it part of the host
	# query and fragment omitted
	return rawurlencode ($enc_url); # finally encode url and return it
} # encode_path

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";
  } # debugon
  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";
  } # debugon
  return array ($widthkm, $heightkm, $dbgval);
}

# http://snipplr.com/view/2531/
function distance($lat1, $lng1, $lat2, $lng2)
# convert distance in km between two lat longs
{
	$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