var D2R = 0.01745329251994329577;
|
var R2D = 57.29577951308232088;
|
var PJD_3PARAM = 1;
|
var PJD_7PARAM = 2;
|
var datum_transform = require('./datum_transform');
|
var adjust_axis = require('./adjust_axis');
|
var proj = require('./Proj');
|
var toPoint = require('./common/toPoint');
|
module.exports = function transform(source, dest, point) {
|
var wgs84;
|
if (Array.isArray(point)) {
|
point = toPoint(point);
|
}
|
function checkNotWGS(source, dest) {
|
return ((source.datum.datum_type === PJD_3PARAM || source.datum.datum_type === PJD_7PARAM) && dest.datumCode !== "WGS84");
|
}
|
|
// Workaround for datum shifts towgs84, if either source or destination projection is not wgs84
|
if (source.datum && dest.datum && (checkNotWGS(source, dest) || checkNotWGS(dest, source))) {
|
wgs84 = new proj('WGS84');
|
transform(source, wgs84, point);
|
source = wgs84;
|
}
|
// DGR, 2010/11/12
|
if (source.axis !== "enu") {
|
adjust_axis(source, false, point);
|
}
|
// Transform source points to long/lat, if they aren't already.
|
if (source.projName === "longlat") {
|
point.x *= D2R; // convert degrees to radians
|
point.y *= D2R;
|
}
|
else {
|
if (source.to_meter) {
|
point.x *= source.to_meter;
|
point.y *= source.to_meter;
|
}
|
source.inverse(point); // Convert Cartesian to longlat
|
}
|
// Adjust for the prime meridian if necessary
|
if (source.from_greenwich) {
|
point.x += source.from_greenwich;
|
}
|
|
// Convert datums if needed, and if possible.
|
point = datum_transform(source.datum, dest.datum, point);
|
|
// Adjust for the prime meridian if necessary
|
if (dest.from_greenwich) {
|
point.x -= dest.from_greenwich;
|
}
|
|
if (dest.projName === "longlat") {
|
// convert radians to decimal degrees
|
point.x *= R2D;
|
point.y *= R2D;
|
}
|
else { // else project
|
dest.forward(point);
|
if (dest.to_meter) {
|
point.x /= dest.to_meter;
|
point.y /= dest.to_meter;
|
}
|
}
|
|
// DGR, 2010/11/12
|
if (dest.axis !== "enu") {
|
adjust_axis(dest, true, point);
|
}
|
|
return point;
|
};
|
