public function testWGS84ToRT90() { $wgsPos = null; $rtPos = null; try { $wgsPos = new WGS84Position("N 59º 58' 55.23\" E 017º 50' 06.12\"", WGS84Format::DEGREES_MINUTES_SECONDS); $rtPos = new RT90Position($wgsPos, RT90Projection::RT90_2_5_GON_V); } catch (ParseException $e) { $this->fail($e->getMessage()); } // Conversion values from Lantmateriet.se, they convert from DMS only. // Reference: http://www.lantmateriet.se/templates/LMV_Enkelkoordinattransformation.aspx?id=11500 $xPosFromLM = 6653174.343; $yPosFromLM = 1613318.742; $lat = (double) round($rtPos->getLatitude() * 1000) / 1000; $lon = (double) round($rtPos->getLongitude() * 1000) / 1000; $this->assertEquals($lat, $xPosFromLM); //fix accuracy: 0.0001d $this->assertEquals($lon, $yPosFromLM); }
/** * Format coordinates * * @param float $lat Latitude * @param float $lon Longitude * @param string $format Format * * @return array X, Y or lat, lon */ public static function format($lat, $lon, $format) { $wgs84 = new WGS84Position($lat, $lon); switch ($format) { case 'wgs84_dms': return array($wgs84->latitudeToString(WGS84Format::DEGREES_MINUTES_SECONDS), $wgs84->longitudeToString(WGS84Format::DEGREES_MINUTES_SECONDS)); case 'wgs84_dm': return array($wgs84->latitudeToString(WGS84Format::DEGREES_MINUTES), $wgs84->longitudeToString(WGS84Format::DEGREES_MINUTES)); case 'wgs84_dd': case 'wgs84_decimal': return array($wgs84->latitudeToString(WGS84Format::DEGREES), $wgs84->longitudeToString(WGS84Format::DEGREES)); case 'rt90': $rtPos = new RT90Position($wgs84, RT90Projection::RT90_2_5_GON_V); $x = (double) round($rtPos->getLatitude() * 1000) / 1000; $y = (double) round($rtPos->getLongitude() * 1000) / 1000; return array($x, $y); case 'sweref99': $rtPos = new SWEREF99Position($wgs84, SWEREFProjection::SWEREF_99_TM); $x = (double) round($rtPos->getLatitude() * 1000) / 1000; $y = (double) round($rtPos->getLongitude() * 1000) / 1000; return array($x, $y); } }