mirror of
https://github.com/tumic0/GPXSee.git
synced 2024-10-07 07:13:21 +02:00
139 lines
3.3 KiB
C++
139 lines
3.3 KiB
C++
#include "datum.h"
|
|
#include "mercator.h"
|
|
#include "webmercator.h"
|
|
#include "transversemercator.h"
|
|
#include "lambertconic.h"
|
|
#include "albersequal.h"
|
|
#include "lambertazimuthal.h"
|
|
#include "latlon.h"
|
|
#include "gcs.h"
|
|
#include "pcs.h"
|
|
#include "projection.h"
|
|
|
|
|
|
Projection::Method::Method(int id)
|
|
{
|
|
switch (id) {
|
|
case 1024:
|
|
case 9801:
|
|
case 9802:
|
|
case 9804:
|
|
case 9807:
|
|
case 9815:
|
|
case 9820:
|
|
case 9822:
|
|
_id = id;
|
|
break;
|
|
default:
|
|
_id = 0;
|
|
}
|
|
}
|
|
|
|
Projection::Projection(const PCS *pcs) : _gcs(pcs->gcs()), _units(pcs->units()),
|
|
_cs(pcs->coordinateSystem()), _geographic(false)
|
|
{
|
|
const Ellipsoid *ellipsoid = _gcs->datum().ellipsoid();
|
|
const Projection::Setup &setup = pcs->setup();
|
|
|
|
switch (pcs->method().id()) {
|
|
case 1024:
|
|
_ct = new WebMercator();
|
|
break;
|
|
case 9801:
|
|
case 9815: // Oblique mercator aproximation using LCC1
|
|
_ct = new LambertConic1(ellipsoid, setup.latitudeOrigin(),
|
|
setup.longitudeOrigin(), setup.scale(), setup.falseEasting(),
|
|
setup.falseNorthing());
|
|
break;
|
|
case 9802:
|
|
_ct = new LambertConic2(ellipsoid, setup.standardParallel1(),
|
|
setup.standardParallel2(), setup.latitudeOrigin(),
|
|
setup.longitudeOrigin(), setup.falseEasting(),
|
|
setup.falseNorthing());
|
|
break;
|
|
case 9804:
|
|
_ct = new Mercator(ellipsoid, setup.latitudeOrigin(),
|
|
setup.longitudeOrigin(), setup.falseEasting(),
|
|
setup.falseNorthing());
|
|
break;
|
|
case 9807:
|
|
_ct = new TransverseMercator(ellipsoid, setup.latitudeOrigin(),
|
|
setup.longitudeOrigin(), setup.scale(), setup.falseEasting(),
|
|
setup.falseNorthing());
|
|
break;
|
|
case 9820:
|
|
_ct = new LambertAzimuthal(ellipsoid, setup.latitudeOrigin(),
|
|
setup.longitudeOrigin(), setup.falseEasting(),
|
|
setup.falseNorthing());
|
|
break;
|
|
case 9822:
|
|
_ct = new AlbersEqual(ellipsoid, setup.standardParallel1(),
|
|
setup.standardParallel2(), setup.latitudeOrigin(),
|
|
setup.longitudeOrigin(), setup.falseEasting(),
|
|
setup.falseNorthing());
|
|
break;
|
|
default:
|
|
_ct = 0;
|
|
}
|
|
}
|
|
|
|
Projection::Projection(const GCS *gcs, const CoordinateSystem &cs)
|
|
: _gcs(gcs), _units(LinearUnits(9001)), _cs(cs), _geographic(true)
|
|
{
|
|
_ct = new LatLon(gcs->angularUnits());
|
|
}
|
|
|
|
Projection::Projection(const Projection &p)
|
|
{
|
|
_gcs = p._gcs;
|
|
_units = p._units;
|
|
_ct = p._ct ? p._ct->clone() : 0;
|
|
_geographic = p._geographic;
|
|
_cs = p._cs;
|
|
}
|
|
|
|
Projection::~Projection()
|
|
{
|
|
delete _ct;
|
|
}
|
|
|
|
Projection &Projection::operator=(const Projection &p)
|
|
{
|
|
_gcs = p._gcs;
|
|
_units = p._units;
|
|
_ct = p._ct ? p._ct->clone() : 0;
|
|
_geographic = p._geographic;
|
|
_cs = p._cs;
|
|
|
|
return *this;
|
|
}
|
|
|
|
PointD Projection::ll2xy(const Coordinates &c) const
|
|
{
|
|
return isValid()
|
|
? _units.fromMeters(_ct->ll2xy(_gcs->fromWGS84(c))) : PointD();
|
|
}
|
|
|
|
Coordinates Projection::xy2ll(const PointD &p) const
|
|
{
|
|
return isValid()
|
|
? _gcs->toWGS84(_ct->xy2ll(_units.toMeters(p))) : Coordinates();
|
|
}
|
|
|
|
#ifndef QT_NO_DEBUG
|
|
QDebug operator<<(QDebug dbg, const Projection::Setup &setup)
|
|
{
|
|
dbg.nospace() << "Setup(" << setup.latitudeOrigin() << ", "
|
|
<< setup.longitudeOrigin() << ", " << setup.scale() << ", "
|
|
<< setup.falseEasting() << ", " << setup.falseNorthing() << ", "
|
|
<< setup.standardParallel1() << ", " << setup.standardParallel2() << ")";
|
|
return dbg.space();
|
|
}
|
|
|
|
QDebug operator<<(QDebug dbg, const Projection::Method &method)
|
|
{
|
|
dbg.nospace() << "Method(" << method.id() << ")";
|
|
return dbg.space();
|
|
}
|
|
#endif // QT_NO_DEBUG
|