Beispiel #1
0
 /**
  * @expectedException \PHPUnit_Framework_Error_Warning
  */
 public function testDistanceDifferentEllipsoids()
 {
     $work = new LatLng(51.54105, -0.12319, RefEll::Airy1830());
     $charingCross = new LatLng(51.507977, -0.124588, RefEll::WGS84());
     $expected = 3678.49665;
     self::assertEquals($expected, $work->distance($charingCross));
 }
Beispiel #2
0
 public function getReferenceEllipsoid()
 {
     return RefEll::WGS84();
 }
Beispiel #3
0
 /**
  * Convert a WGS84 latitude and longitude to an UTM reference
  *
  * Reference values for transformation are taken from OS document
  * "A Guide to Coordinate Systems in Great Britain"
  * @return UTMRef
  */
 public function toUTMRef()
 {
     if ($this->refEll != RefEll::WGS84()) {
         trigger_error('Current co-ordinates are in a non-WGS84 datum', E_USER_WARNING);
     }
     $longitudeZone = (int) (($this->lng + 180) / 6) + 1;
     // Special zone for Norway
     if ($this->lat >= 56 && $this->lat < 64 && $this->lng >= 3 && $this->lng < 12) {
         $longitudeZone = 32;
     }
     // Special zones for Svalbard
     if ($this->lat >= 72 && $this->lat < 84) {
         if ($this->lng >= 0 && $this->lng < 9) {
             $longitudeZone = 31;
         } else {
             if ($this->lng >= 9 && $this->lng < 21) {
                 $longitudeZone = 33;
             } else {
                 if ($this->lng >= 21 && $this->lng < 33) {
                     $longitudeZone = 35;
                 } else {
                     if ($this->lng >= 33 && $this->lng < 42) {
                         $longitudeZone = 37;
                     }
                 }
             }
         }
     }
     $UTMZone = $this->getUTMLatitudeZoneLetter($this->lat);
     $UTM = new UTMRef(0, 0, $UTMZone, $longitudeZone);
     //dummy to get reference data
     $scale = $UTM->getScaleFactor();
     $N0 = $UTM->getOriginNorthing();
     $E0 = $UTM->getOriginEasting();
     $phi0 = $UTM->getOriginLatitude();
     $lambda0 = $UTM->getOriginLongitude();
     $coords = $this->toTransverseMercatorEastingNorthing($scale, $E0, $N0, $phi0, $lambda0);
     if ($this->lat < 0) {
         $coords['N'] += 10000000;
     }
     return new UTMRef($coords['E'], $coords['N'], $UTMZone, $longitudeZone);
 }