Page MenuHomePhabricator
Files F2706

emaps-bg.php
Public
Actions

Authored By
• bzimport
Nov 21 2014, 9:10 PM
Size
8 KB
Referenced Files
None
Subscribers
None

emaps-bg.php
View Options

<?php
/**
This file has two parts:
* a heavily stripped-down version of the "gPoint.php class"
written by Brenor Brophy http://www.brenorbrophy.com
(original found at http://www.phpclasses.org)
* a program to connect to www.emaps.bg, converting geographic coordinates
to UTM coordinates, made for the purposes of bg.wikipedia.org
written by Petko Yotov http://5ko.free.fr/
This file may be copied and modified under the terms of the
*** GNU Ceneral Public License version 2 or newer ***
*/
/*------------------------------------------------------------------------------
** PHP class to convert Latitude & Longitude coordinates into
** UTM & Lambert Conic Conformal Northing/Easting coordinates.
** Version: 1.0
** Author: Brenor Brophy
** Email & URL: brenor@sbcglobal.net www.brenorbrophy.com
**
** Code for datum and UTM conversion was converted from C++ code written by
** Chuck Gantz (chuck.gantz@globalstar.com) from http://www.gpsy.com/gpsinfo/geotoutm/
** This URL has many other references to useful information concerning conversion
** of coordinates.
*/
define ("meter2nm", (1/1852));
define ("nm2meter", 1852);
/*------------------------------------------------------------------------------
** class gPoint ... for Geographic Point
**
** This class encapsulates the methods for representing a geographic point on the
** earth in three different coordinate systema. Lat/Long, UTM and Lambert Conic
** Conformal.
*/
class gPoint
{
/* Reference ellipsoids derived from Peter H. Dana's website-
** http://www.utexas.edu/depts/grg/gcraft/notes/datum/elist.html
** Department of Geography, University of Texas at Austin
** Internet: pdana@mail.utexas.edu 3/22/95
**
** Source:
** Defense Mapping Agency. 1987b. DMA Technical Report: Supplement to Department
** of Defense World Geodetic System 1984 Technical Report. Part I and II.
** Washington, DC: Defense Mapping Agency
*/
var $ellipsoid = array(//Ellipsoid name, Equatorial Radius, square of eccentricity
"Airy" =>array (6377563, 0.00667054),
"Australian National" =>array (6378160, 0.006694542),
"Bessel 1841" =>array (6377397, 0.006674372),
"Bessel 1841 Nambia" =>array (6377484, 0.006674372),
"Clarke 1866" =>array (6378206, 0.006768658),
"Clarke 1880" =>array (6378249, 0.006803511),
"Everest" =>array (6377276, 0.006637847),
"Fischer 1960 Mercury" =>array (6378166, 0.006693422),
"Fischer 1968" =>array (6378150, 0.006693422),
"GRS 1967" =>array (6378160, 0.006694605),
"GRS 1980" =>array (6378137, 0.00669438),
"Helmert 1906" =>array (6378200, 0.006693422),
"Hough" =>array (6378270, 0.00672267),
"International" =>array (6378388, 0.00672267),
"Krassovsky" =>array (6378245, 0.006693422),
"Modified Airy" =>array (6377340, 0.00667054),
"Modified Everest" =>array (6377304, 0.006637847),
"Modified Fischer 1960" =>array (6378155, 0.006693422),
"South American 1969" =>array (6378160, 0.006694542),
"WGS 60" =>array (6378165, 0.006693422),
"WGS 66" =>array (6378145, 0.006694542),
"WGS 72" =>array (6378135, 0.006694318),
"WGS 84" =>array (6378137, 0.00669438));
// Properties
var $a; // Equatorial Radius
var $e2; // Square of eccentricity
var $datum; // Selected datum
var $Xp, $Yp; // X,Y pixel location
var $lat, $long; // Latitude & Longitude of the point
var $utmNorthing, $utmEasting, $utmZone; // UTM Coordinates of the point
var $lccNorthing, $lccEasting; // Lambert coordinates of the point
var $falseNorthing, $falseEasting; // Origin coordinates for Lambert Projection
var $latOfOrigin; // For Lambert Projection
var $longOfOrigin; // For Lambert Projection
var $firstStdParallel; // For lambert Projection
var $secondStdParallel; // For lambert Projection
// constructor
function gPoint($datum='WGS 84') // Default datum is WGS 84
{
$this->a = $this->ellipsoid[$datum][0]; // Set datum Equatorial Radius
$this->e2 = $this->ellipsoid[$datum][1]; // Set datum Square of eccentricity
$this->datum = $datum; // Save the datum
}
//
// Set/Get/Output Longitude & Latitude of the point
//
function setLongLat($long, $lat)
{
$this->long = $long; $this->lat = $lat;
}
function convertLLtoTM($LongOrigin='')
{
$k0 = 0.9996;
$falseEasting = 500000.0;
//Make sure the longitude is between -180.00 .. 179.9
$LongTemp = ($this->long+180)-(integer)(($this->long+180)/360)*360-180; // -180.00 .. 179.9;
$LatRad = deg2rad($this->lat);
$LongRad = deg2rad($LongTemp);
$ZoneNumber=35; // Added by Petko Yotov
$LongOrigin = ($ZoneNumber - 1)*6 - 180 + 3; //+3 puts origin in middle of zone
$LongOriginRad = deg2rad($LongOrigin);
$eccPrimeSquared = ($this->e2)/(1-$this->e2);
$N = $this->a/sqrt(1-$this->e2*sin($LatRad)*sin($LatRad));
$T = tan($LatRad)*tan($LatRad);
$C = $eccPrimeSquared*cos($LatRad)*cos($LatRad);
$A = cos($LatRad)*($LongRad-$LongOriginRad);
$M = $this->a*((1 - $this->e2/4 - 3*$this->e2*$this->e2/64 - 5*$this->e2*$this->e2*$this->e2/256)*$LatRad
- (3*$this->e2/8 + 3*$this->e2*$this->e2/32 + 45*$this->e2*$this->e2*$this->e2/1024)*sin(2*$LatRad)
+ (15*$this->e2*$this->e2/256 + 45*$this->e2*$this->e2*$this->e2/1024)*sin(4*$LatRad)
- (35*$this->e2*$this->e2*$this->e2/3072)*sin(6*$LatRad));
$this->utmEasting = ($k0*$N*($A+(1-$T+$C)*$A*$A*$A/6
+ (5-18*$T+$T*$T+72*$C-58*$eccPrimeSquared)*$A*$A*$A*$A*$A/120)
+ $falseEasting);
$this->utmNorthing = ($k0*($M+$N*tan($LatRad)*($A*$A/2+(5-$T+9*$C+4*$C*$C)*$A*$A*$A*$A/24
+ (61-58*$T+$T*$T+600*$C-330*$eccPrimeSquared)*$A*$A*$A*$A*$A*$A/720)));
if($this->lat < 0)
$this->utmNorthing += 10000000.0; //10000000 meter offset for southern hemisphere
}
} // end of class gPoint
// Next code by Petko Yotov
function create_link($lat, $lon, $ell)
{
$gpoint =& new gPoint($ell);
$gpoint->setLongLat($lon, $lat);
$gpoint->convertLLtoTM('');
$mapX = round($gpoint->utmEasting);
$mapY = round($gpoint->utmNorthing);
$minx = $mapX-14425;
$maxx = $mapX+14425;
$miny = $mapY-10000;
$maxy = $mapY+10000;
$base_url = "http://www.emaps.bg/emaps/content.asp?mapSize=1&cboMaps=0";
$url = "$base_url&minx=$minx&maxx=$maxx&miny=$miny&maxy=$maxy";
return $url;
}
function message($en, $bg)
{
header("Content-type: text/html; charset=UTF-8");
echo <<<EOF
<html><title>Error... Грешка...</title><body>
<h1>Error... Грешка...</h1>
<strong><big>$en</big></strong><br/>
You can:
<ul>
<li><a href='javascript:history.go(-1);'>Go back</a>.</li>
<li>Go to <a href='http://www.emaps.bg/'>www.emaps.bg</a>.</li>
<li>Go to <a href='http://bg.wikipedia.org/'>bg.wikipedia.org</a>.</li>
<li>View the <a href='{$_SERVER["PHP_SELF"]}?q=source'>source code of this script</a>.</li>
</ul>
<br/>
<br/>
<strong><big>$bg</big></strong><br/>
Можете:
<ul>
<li><a href='javascript:history.go(-1);'>Да се върнете</a>.</li>
<li>Да посетите <a href='http://www.emaps.bg/'>www.emaps.bg</a>.</li>
<li>Да посетите <a href='http://bg.wikipedia.org/'>bg.wikipedia.org</a>.</li>
<li>Да прегледате <a href='{$_SERVER["PHP_SELF"]}?q=source'>изходния код на този скрипт</a>.</li>
</ul>
<br/>
<br/>
<hr/>
<address>
<small>Script by <a href="http://5ko.free.fr/">Petko Yotov</a>.</small>
</address>
</body></html>
EOF;
exit;
}
if(strtolower($_GET['q']) == "source")
{
header("Content-type: text/html; charset=UTF-8");
highlight_file(__FILE__);
exit;
}
if(!isset($_GET['q']) || preg_match("/[^0-9,\.]/", trim($_GET['q']) ) )
{
message("No coordinates requested or illegal character", "Не са посочени координати или е въведен неразрешен символ");
}
list($lat, $lon) = explode(",", $_GET['q'], 2);
$lat = floatval($lat); // 41..45
$lon = floatval($lon); // 22..29
if($lat<41 || $lat>45 || $lon<22 || $lon>29 )
{
message("Coordinates are out of range. Lattitude '$lat' (must be 41..45), Longitude '$lon' (must be 22..29)", "Координатите са извън картата. Ширина: '$lat' (трябва да е 41..45), дължина: '$lon' (трябва да е 22..29)");
}
$url = create_link($lat, $lon, "WGS 84");
header("Location: $url");
exit;
?>

File Metadata

Mime Type
text/x-php
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
2281
Default Alt Text
emaps-bg.php (8 KB)

Event Timeline

Content licensed under Creative Commons Attribution-ShareAlike (CC BY-SA) 4.0 unless otherwise noted; code licensed under GNU General Public License (GPL) 2.0 or later and other open source licenses. By using this site, you agree to the Terms of Use, Privacy Policy, and Code of Conduct. · Wikimedia Foundation · Privacy Policy · Code of Conduct · Terms of Use · Disclaimer · CC-BY-SA · GPL