This solves canvas antialiasing issues and also simplifies the code handling wrapX in the source projection.
267 lines
7.7 KiB
JavaScript
267 lines
7.7 KiB
JavaScript
goog.provide('ol.reproj.Triangulation');
|
|
|
|
goog.require('goog.array');
|
|
goog.require('goog.math');
|
|
goog.require('ol.extent');
|
|
goog.require('ol.proj');
|
|
|
|
|
|
/**
|
|
* Single triangle; consists of 3 source points and 3 target points.
|
|
*
|
|
* @typedef {{source: Array.<ol.Coordinate>,
|
|
* target: Array.<ol.Coordinate>}}
|
|
*/
|
|
ol.reproj.Triangle;
|
|
|
|
|
|
|
|
/**
|
|
* @param {ol.proj.Projection} sourceProj
|
|
* @param {ol.proj.Projection} targetProj
|
|
* @param {ol.Extent} targetExtent
|
|
* @param {ol.Extent=} opt_maxSourceExtent
|
|
* @param {number=} opt_maxSubdiv Maximal subdivision.
|
|
* @param {number=} opt_errorThreshold Acceptable error (in source units).
|
|
* @constructor
|
|
*/
|
|
ol.reproj.Triangulation = function(sourceProj, targetProj, targetExtent,
|
|
opt_maxSourceExtent, opt_maxSubdiv, opt_errorThreshold) {
|
|
|
|
/**
|
|
* @type {ol.proj.Projection}
|
|
* @private
|
|
*/
|
|
this.sourceProj_ = sourceProj;
|
|
|
|
/**
|
|
* @type {ol.proj.Projection}
|
|
* @private
|
|
*/
|
|
this.targetProj_ = targetProj;
|
|
|
|
/**
|
|
* @type {ol.TransformFunction}
|
|
* @private
|
|
*/
|
|
this.transformInv_ = ol.proj.getTransform(this.targetProj_, this.sourceProj_);
|
|
|
|
/**
|
|
* @type {ol.Extent}
|
|
* @private
|
|
*/
|
|
this.maxSourceExtent_ = goog.isDef(opt_maxSourceExtent) ?
|
|
opt_maxSourceExtent : null;
|
|
|
|
var errorThreshold = goog.isDef(opt_errorThreshold) ?
|
|
opt_errorThreshold : 0; //TODO: define
|
|
/**
|
|
* @type {number}
|
|
* @private
|
|
*/
|
|
this.errorThresholdSquared_ = errorThreshold * errorThreshold;
|
|
|
|
/**
|
|
* @type {Array.<ol.reproj.Triangle>}
|
|
* @private
|
|
*/
|
|
this.triangles_ = [];
|
|
|
|
/**
|
|
* @type {ol.Extent}
|
|
* @private
|
|
*/
|
|
this.trianglesSourceExtent_ = null;
|
|
|
|
/**
|
|
* Indicates that source coordinates has to be shifted during reprojection.
|
|
* This is needed when the triangulation crosses
|
|
* edge of the source projection (dateline).
|
|
* @type {boolean}
|
|
* @private
|
|
*/
|
|
this.wrapsXInSource_ = false;
|
|
|
|
/**
|
|
* @type {number}
|
|
* @private
|
|
*/
|
|
this.sourceWorldWidth_ = ol.extent.getWidth(this.sourceProj_.getExtent());
|
|
|
|
var tlDst = ol.extent.getTopLeft(targetExtent);
|
|
var trDst = ol.extent.getTopRight(targetExtent);
|
|
var brDst = ol.extent.getBottomRight(targetExtent);
|
|
var blDst = ol.extent.getBottomLeft(targetExtent);
|
|
var tlDstSrc = this.transformInv_(tlDst);
|
|
var trDstSrc = this.transformInv_(trDst);
|
|
var brDstSrc = this.transformInv_(brDst);
|
|
var blDstSrc = this.transformInv_(blDst);
|
|
|
|
this.addQuadIfValid_(tlDst, trDst, brDst, blDst,
|
|
tlDstSrc, trDstSrc, brDstSrc, blDstSrc,
|
|
opt_maxSubdiv || 0);
|
|
};
|
|
|
|
|
|
/**
|
|
* Adds triangle to the triangulation.
|
|
* @param {ol.Coordinate} a
|
|
* @param {ol.Coordinate} b
|
|
* @param {ol.Coordinate} c
|
|
* @param {ol.Coordinate} aSrc
|
|
* @param {ol.Coordinate} bSrc
|
|
* @param {ol.Coordinate} cSrc
|
|
* @private
|
|
*/
|
|
ol.reproj.Triangulation.prototype.addTriangle_ = function(a, b, c,
|
|
aSrc, bSrc, cSrc) {
|
|
this.triangles_.push({
|
|
source: [aSrc, bSrc, cSrc],
|
|
target: [a, b, c]
|
|
});
|
|
};
|
|
|
|
|
|
/**
|
|
* Adds quad (points in clock-wise order) to the triangulation
|
|
* (and reprojects the vertices) if valid.
|
|
* @param {ol.Coordinate} a
|
|
* @param {ol.Coordinate} b
|
|
* @param {ol.Coordinate} c
|
|
* @param {ol.Coordinate} d
|
|
* @param {ol.Coordinate} aSrc
|
|
* @param {ol.Coordinate} bSrc
|
|
* @param {ol.Coordinate} cSrc
|
|
* @param {ol.Coordinate} dSrc
|
|
* @param {number} maxSubdiv Maximal allowed subdivision of the quad.
|
|
* @private
|
|
*/
|
|
ol.reproj.Triangulation.prototype.addQuadIfValid_ = function(a, b, c, d,
|
|
aSrc, bSrc, cSrc, dSrc, maxSubdiv) {
|
|
|
|
var srcQuadExtent = ol.extent.boundingExtent([aSrc, bSrc, cSrc, dSrc]);
|
|
if (!goog.isNull(this.maxSourceExtent_)) {
|
|
if (!ol.extent.intersects(srcQuadExtent, this.maxSourceExtent_)) {
|
|
// whole quad outside source projection extent -> ignore
|
|
return;
|
|
}
|
|
}
|
|
var srcCoverageX = ol.extent.getWidth(srcQuadExtent) / this.sourceWorldWidth_;
|
|
|
|
// when the quad is wrapped in the source projection
|
|
// it covers most of the projection extent, but not fully
|
|
var wrapsX = this.sourceProj_.canWrapX() &&
|
|
srcCoverageX > 0.5 && srcCoverageX < 1;
|
|
|
|
if (maxSubdiv > 0) {
|
|
var needsSubdivision = !wrapsX && this.sourceProj_.isGlobal() &&
|
|
srcCoverageX > 0.25; //TODO: define
|
|
|
|
var center = [(a[0] + c[0]) / 2, (a[1] + c[1]) / 2];
|
|
var centerSrc = this.transformInv_(center);
|
|
|
|
if (!needsSubdivision) {
|
|
var dx;
|
|
if (wrapsX) {
|
|
var centerSrcEstimX =
|
|
(goog.math.modulo(aSrc[0], this.sourceWorldWidth_) +
|
|
goog.math.modulo(cSrc[0], this.sourceWorldWidth_)) / 2;
|
|
dx = centerSrcEstimX -
|
|
goog.math.modulo(centerSrc[0], this.sourceWorldWidth_);
|
|
} else {
|
|
dx = (aSrc[0] + cSrc[0]) / 2 - centerSrc[0];
|
|
}
|
|
var dy = (aSrc[1] + cSrc[1]) / 2 - centerSrc[1];
|
|
var centerSrcErrorSquared = dx * dx + dy * dy;
|
|
needsSubdivision = centerSrcErrorSquared > this.errorThresholdSquared_;
|
|
}
|
|
if (needsSubdivision) {
|
|
var ab = [(a[0] + b[0]) / 2, (a[1] + b[1]) / 2];
|
|
var abSrc = this.transformInv_(ab);
|
|
var bc = [(b[0] + c[0]) / 2, (b[1] + c[1]) / 2];
|
|
var bcSrc = this.transformInv_(bc);
|
|
var cd = [(c[0] + d[0]) / 2, (c[1] + d[1]) / 2];
|
|
var cdSrc = this.transformInv_(cd);
|
|
var da = [(d[0] + a[0]) / 2, (d[1] + a[1]) / 2];
|
|
var daSrc = this.transformInv_(da);
|
|
|
|
this.addQuadIfValid_(a, ab, center, da,
|
|
aSrc, abSrc, centerSrc, daSrc, maxSubdiv - 1);
|
|
this.addQuadIfValid_(ab, b, bc, center,
|
|
abSrc, bSrc, bcSrc, centerSrc, maxSubdiv - 1);
|
|
this.addQuadIfValid_(center, bc, c, cd,
|
|
centerSrc, bcSrc, cSrc, cdSrc, maxSubdiv - 1);
|
|
this.addQuadIfValid_(da, center, cd, d,
|
|
daSrc, centerSrc, cdSrc, dSrc, maxSubdiv - 1);
|
|
|
|
return;
|
|
}
|
|
}
|
|
|
|
if (wrapsX) {
|
|
this.wrapsXInSource_ = true;
|
|
}
|
|
|
|
this.addTriangle_(a, c, d, aSrc, cSrc, dSrc);
|
|
this.addTriangle_(a, b, c, aSrc, bSrc, cSrc);
|
|
};
|
|
|
|
|
|
/**
|
|
* @return {ol.Extent}
|
|
*/
|
|
ol.reproj.Triangulation.prototype.calculateSourceExtent = function() {
|
|
if (!goog.isNull(this.trianglesSourceExtent_)) {
|
|
return this.trianglesSourceExtent_;
|
|
}
|
|
|
|
var extent = ol.extent.createEmpty();
|
|
|
|
if (this.wrapsXInSource_) {
|
|
// although only some of the triangles are crossing the dateline,
|
|
// all coordiantes need to be "shifted" to be positive
|
|
// to properly calculate the extent (and then possibly shifted back)
|
|
|
|
goog.array.forEach(this.triangles_, function(triangle, i, arr) {
|
|
var src = triangle.source;
|
|
ol.extent.extendCoordinate(extent,
|
|
[goog.math.modulo(src[0][0], this.sourceWorldWidth_), src[0][1]]);
|
|
ol.extent.extendCoordinate(extent,
|
|
[goog.math.modulo(src[1][0], this.sourceWorldWidth_), src[1][1]]);
|
|
ol.extent.extendCoordinate(extent,
|
|
[goog.math.modulo(src[2][0], this.sourceWorldWidth_), src[2][1]]);
|
|
}, this);
|
|
|
|
var sourceProjExtent = this.sourceProj_.getExtent();
|
|
var right = sourceProjExtent[2];
|
|
if (extent[0] > right) extent[0] -= this.sourceWorldWidth_;
|
|
if (extent[2] > right) extent[2] -= this.sourceWorldWidth_;
|
|
} else {
|
|
goog.array.forEach(this.triangles_, function(triangle, i, arr) {
|
|
var src = triangle.source;
|
|
ol.extent.extendCoordinate(extent, src[0]);
|
|
ol.extent.extendCoordinate(extent, src[1]);
|
|
ol.extent.extendCoordinate(extent, src[2]);
|
|
});
|
|
}
|
|
|
|
this.trianglesSourceExtent_ = extent;
|
|
return extent;
|
|
};
|
|
|
|
|
|
/**
|
|
* @return {boolean}
|
|
*/
|
|
ol.reproj.Triangulation.prototype.getWrapsXInSource = function() {
|
|
return this.wrapsXInSource_;
|
|
};
|
|
|
|
|
|
/**
|
|
* @return {Array.<ol.reproj.Triangle>}
|
|
*/
|
|
ol.reproj.Triangulation.prototype.getTriangles = function() {
|
|
return this.triangles_;
|
|
};
|