archive-gh-me/WebWorldWind
1
0
Fork 0
mirror of https://github.com/NASAWorldWind/WebWorldWind.git synced 2026-02-01 16:38:15 +00:00
Code Issues Projects Releases Wiki Activity

Added handling of EPSG:3857 in WmtsLayer. Filled out OwsServiceProvider and OwsServiceIdentification.

Browse Source
This commit is contained in:
Tom Gaskins 2015-09-14 10:24:07 -07:00
parent 47bf732cce
commit db50c16de5
4 changed files with 92 additions and 28 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

67
src/layer/WmtsLayer.js
View File

@ -34,7 +34,6 @@ define([
// TODO: Test Mercator layers.
// TODO: Support tile matrix limits.
// TODO: Extensibility for other projections.
// TODO: Determine where to get the layer's bounding box from if WGS84 bounding box not specified.
// TODO: Finish parsing capabilities document (ServiceIdentification and ServiceProvider).
// TODO: Time dimensions.
@ -173,19 +172,26 @@ define([
"No supported Tile Matrix Set could be found."));
}
// Determine the layer's sector.
// TODO: The WGS84 bounding box is optional. Determine how to determine the sector if it is not specified.
// Determine the layer's sector if possible. Mandatory for EPSG:4326 tile matrix sets. (Others compute
// it from tile Matrix Set metadata.)
if (layerCaps.wgs84BoundingBox) {
this.sector = new Sector(
layerCaps.wgs84BoundingBox.lowerCorner[1],
layerCaps.wgs84BoundingBox.upperCorner[1],
layerCaps.wgs84BoundingBox.lowerCorner[0],
layerCaps.wgs84BoundingBox.upperCorner[0]);
} else {
// For now, throw an exception if there is no WGS84 bounding box.
} else if (this.tileMatrixSet.boundingBox &&
WmtsLayerCapabilities.isEpsg4326Crs(this.tileMatrixSet.boundingBox.crs)) {
this.sector = new Sector(
this.tileMatrixSet.boundingBox.lowerCorner[1],
this.tileMatrixSet.boundingBox.upperCorner[1],
this.tileMatrixSet.boundingBox.lowerCorner[0],
this.tileMatrixSet.boundingBox.upperCorner[0]);
} else if (WmtsLayerCapabilities.isEpsg4326Crs(this.tileMatrixSet.supportedCRS)) {
// Throw an exception if there is no 4326 bounding box.
throw new ArgumentError(
Logger.logMessage(Logger.LEVEL_SEVERE, "WmtsLayer", "constructor",
"No WGS84 bounding box was specified in the layer capabilities."));
"No EPSG:4326 bounding box was specified in the layer or tile matrix set capabilities."));
}
// Form a unique string to identify cache entries.
@ -470,29 +476,48 @@ define([
};
WmtsLayer.prototype.createTile3857 = function (tileMatrix, row, column) {
var mapWidth = tileMatrix.tileWidth * tileMatrix.matrixWidth,
mapHeight = tileMatrix.tileHeight * tileMatrix.matrixHeight,
swX = WWMath.clamp(column * tileMatrix.tileWidth - 0.5, 0, mapWidth),
neY = WWMath.clamp(row * tileMatrix.tileHeight - 0.5, 0, mapHeight),
neX = WWMath.clamp(swX + (tileMatrix.tileWidth) + 0.5, 0, mapWidth),
swY = WWMath.clamp(neY + (tileMatrix.tileHeight) + 0.5, 0, mapHeight),
if (!tileMatrix.mapWidth) {
this.computeTileMatrixValues3857(tileMatrix);
}
var swX = WWMath.clamp(column * tileMatrix.tileWidth - 0.5, 0, tileMatrix.mapWidth),
neY = WWMath.clamp(row * tileMatrix.tileHeight - 0.5, 0, tileMatrix.mapHeight),
neX = WWMath.clamp(swX + (tileMatrix.tileWidth) + 0.5, 0, tileMatrix.mapWidth),
swY = WWMath.clamp(neY + (tileMatrix.tileHeight) + 0.5, 0, tileMatrix.mapHeight),
x, y, swLat, swLon, neLat, neLon;
x = (swX / mapWidth) - 0.5;
y = 0.5 - (swY / mapHeight);
swLat = 90 - 360 * Math.atan(Math.exp(-y * 2 * Math.PI)) / Math.PI; // TODO: Use top-left corner metadata
swLon = 360 * x;
x = swX / tileMatrix.mapWidth;
y = swY / tileMatrix.mapHeight;
swLon = tileMatrix.topLeftCorner[0] + x * tileMatrix.tileMatrixDeltaX;
swLat = tileMatrix.topLeftCorner[1] - y * tileMatrix.tileMatrixDeltaY;
var swDegrees = WWMath.epsg3857ToEpsg4326(swLon, swLat);
x = (neX / mapWidth) - 0.5;
y = 0.5 - (neY / mapHeight);
neLat = 90 - 360 * Math.atan(Math.exp(-y * 2 * Math.PI)) / Math.PI; // TODO: Use top-left corner metadata
neLon = 360 * x;
x = neX / tileMatrix.mapWidth;
y = neY / tileMatrix.mapHeight;
neLon = tileMatrix.topLeftCorner[0] + x * tileMatrix.tileMatrixDeltaX;
neLat = tileMatrix.topLeftCorner[1] - y * tileMatrix.tileMatrixDeltaY;
var neDegrees = WWMath.epsg3857ToEpsg4326(neLon, neLat);
var sector = new Sector(swLat, neLat, swLon, neLon);
var sector = new Sector(swDegrees[0], neDegrees[0], swDegrees[1], neDegrees[1]);
return this.makeTile(sector, tileMatrix, row, column);
};
WmtsLayer.prototype.computeTileMatrixValues3857 = function (tileMatrix) {
var pixelSpan = tileMatrix.scaleDenominator * 0.28e-3,
tileSpanX = tileMatrix.tileWidth * pixelSpan,
tileSpanY = tileMatrix.tileHeight * pixelSpan,
tileMatrixMaxX = tileMatrix.topLeftCorner[0] + tileSpanX * tileMatrix.matrixWidth,
tileMatrixMinY = tileMatrix.topLeftCorner[1] - tileSpanY * tileMatrix.matrixHeight,
bottomRightCorner = [tileMatrixMaxX, tileMatrixMinY],
topLeftCorner = tileMatrix.topLeftCorner;
tileMatrix.tileMatrixDeltaX = bottomRightCorner[0] - topLeftCorner[0];
tileMatrix.tileMatrixDeltaY = topLeftCorner[1] - bottomRightCorner[1];
tileMatrix.mapWidth = tileMatrix.tileWidth * tileMatrix.matrixWidth;
tileMatrix.mapHeight = tileMatrix.tileHeight * tileMatrix.matrixHeight;
};
WmtsLayer.prototype.makeTile = function (sector, tileMatrix, row, column) {
var path = this.cachePath + "-layer/" + tileMatrix.identifier + "/" + row + "/" + column + "."
+ WWUtil.suffixForMimeType(this.imageFormat);

29
src/ogc/OwsServiceIdentification.js
View File

@ -20,7 +20,7 @@ define([
* @classdesc Represents an OWS Service Identification section of an OGC capabilities document.
* This object holds as properties all the fields specified in the OWS Service Identification.
* Fields can be accessed as properties named according to their document names converted to camel case.
* For example, "serviceType" and "titles".
* For example, "serviceType" and "title".
* Note that fields with multiple possible values are returned as arrays, such as "titles" and "abstracts".
* @param {Element} element An XML DOM element representing the OWS Service Identification section.
* @throws {ArgumentError} If the specified XML DOM element is null or undefined.
@ -31,10 +31,31 @@ define([
Logger.logMessage(Logger.LEVEL_SEVERE, "OwsServiceIdentification", "constructor", "missingDomElement"));
}
this.assembleDocument(element);
};
var children = element.children;
for (var c = 0; c < children.length; c++) {
var child = children[c];
OwsServiceIdentification.prototype.assembleDocument = function (element) {
if (child.localName === "ServiceType") {
this.serviceType = child.textContent;
} else if (child.localName === "ServiceTypeVersion") {
this.serviceTypeVersion = child.textContent;
} else if (child.localName === "Profile") {
this.profile = this.profiles || [];
this.profile.push(child.textContent);
} else if (child.localName === "Title") {
this.title = this.title|| [];
this.title.push(child.textContent);
} else if (child.localName === "Abstract") {
this.abstract = this.title|| [];
this.abstract.push(child.textContent);
} else if (child.localName === "Fees") {
this.fees = child.textContent;
} else if (child.localName === "AccessConstraints") {
this.accessConstraints = this.accessConstraints || [];
this.accessConstraints.push(child.textContent);
}
// TODO: Keywords
}
};
return OwsServiceIdentification;

13
src/ogc/OwsServiceProvider.js
View File

@ -30,10 +30,17 @@ define([
Logger.logMessage(Logger.LEVEL_SEVERE, "OwsServiceProvider", "constructor", "missingDomElement"));
}
this.assembleDocument(element);
};
var children = element.children;
for (var c = 0; c < children.length; c++) {
var child = children[c];
OwsServiceProvider.prototype.assembleDocument = function (element) {
if (child.localName === "ProviderName") {
this.providerName = child.textContent;
} else if (child.localName === "ProviderSite") {
this.providerSite = child.getAttribute("xlink:href");
}
// TODO: Service Contact
}
};
return OwsServiceProvider;

11
src/util/WWMath.js
View File

@ -706,6 +706,17 @@ define([
*/
gudermannianInverse: function (latitude) {
return Math.log(Math.tan(Math.PI / 4 + (latitude * Angle.DEGREES_TO_RADIANS) / 2)) / Math.PI;
},
epsg3857ToEpsg4326: function (easting, northing) {
var r = 6.3781e6,
latRadians = (Math.PI / 2) - 2 * Math.atan(Math.exp(-northing / r)),
lonRadians = easting / r;
return [
WWMath.clamp(latRadians * Angle.RADIANS_TO_DEGREES, -90, 90),
WWMath.clamp(lonRadians * Angle.RADIANS_TO_DEGREES, -180, 180)
];
}
};