Andreas Hocevar
468 lines
18 KiB
HTML
468 lines
18 KiB
HTML
<!doctype html>
|
|
<html lang="en">
|
|
<head>
|
|
<meta charset="utf-8">
|
|
<meta http-equiv="X-UA-Compatible" content="chrome=1">
|
|
<meta name="viewport" content="initial-scale=1.0, user-scalable=no, width=device-width">
|
|
<link rel="stylesheet" href="https://maxcdn.bootstrapcdn.com/font-awesome/4.4.0/css/font-awesome.min.css" type="text/css">
|
|
<link rel="stylesheet" href="https://maxcdn.bootstrapcdn.com/bootstrap/3.3.5/css/bootstrap.min.css" type="text/css">
|
|
<link rel="stylesheet" href="./resources/prism/prism.css" type="text/css">
|
|
<link rel="stylesheet" href="../css/ol.css" type="text/css">
|
|
<link rel="stylesheet" href="./resources/layout.css" type="text/css">
|
|
|
|
<link rel="stylesheet" href="shaded-relief.css">
|
|
<script src="./resources/zeroclipboard/ZeroClipboard.min.js"></script>
|
|
<title>Shaded Relief</title>
|
|
</head>
|
|
<body>
|
|
|
|
<header class="navbar" role="navigation">
|
|
<div class="container" id="navbar-inner-container">
|
|
<a class="navbar-brand" href="./"><img src="./resources/logo-70x70.png"> OpenLayers 3 Examples</a>
|
|
</div>
|
|
</header>
|
|
|
|
<div class="container-fluid">
|
|
|
|
<div class="row-fluid">
|
|
<div class="span12">
|
|
<div id="map" class="map"></div>
|
|
<table class="controls">
|
|
<tr>
|
|
<td>vertical exaggeration: <span id="vertOut"></span>x</td>
|
|
<td><input id="vert" type="range" min="1" max="5" value="1"/></td>
|
|
</tr>
|
|
<tr>
|
|
<td>sun elevation: <span id="sunElOut"></span>°</td>
|
|
<td><input id="sunEl" type="range" min="0" max="90" value="45"/></td>
|
|
</tr>
|
|
<tr>
|
|
<td>sun azimuth: <span id="sunAzOut"></span>°</td>
|
|
<td><input id="sunAz" type="range" min="0" max="360" value="45"/></td>
|
|
</tr>
|
|
</table>
|
|
</div>
|
|
</div>
|
|
|
|
|
|
<div class="row-fluid">
|
|
<div class="span12">
|
|
<h4 id="title">Shaded Relief</h4>
|
|
<p id="shortdesc">Calculate shaded relief from elevation data</p>
|
|
<div id="docs"><p>
|
|
This example uses a <code>ol.source.Raster</code> to generate data
|
|
based on another source. The raster source accepts any number of
|
|
input sources (tile or image based) and runs a pipeline of
|
|
operations on the input data. The return from the final
|
|
operation is used as the data for the output source.
|
|
</p> <p>
|
|
In this case, a single tiled source of elevation data is used as input.
|
|
The shaded relief is calculated in a single "image" operation. By setting
|
|
<code>operationType: 'image'</code> on the raster source, operations are
|
|
called with an <code>ImageData</code> object for each of the input sources.
|
|
Operations are also called with a general purpose <code>data</code> object.
|
|
In this example, the sun elevation and azimuth data from the inputs above
|
|
are assigned to this <code>data</code> object and accessed in the shading
|
|
operation. The shading operation returns an array of <code>ImageData</code>
|
|
objects. When the raster source is used by an image layer, the first
|
|
<code>ImageData</code> object returned by the last operation in the pipeline
|
|
is used for rendering.
|
|
</p>
|
|
</div>
|
|
<div id="tags">raster, shaded relief</div>
|
|
<div id="api-links">Related API documentation: <ul class="inline"><li><a href="../apidoc/ol.Map.html" title="API documentation for ol.Map">ol.Map</a></li>,<li><a href="../apidoc/ol.View.html" title="API documentation for ol.View">ol.View</a></li>,<li><a href="../apidoc/ol.layer.Image.html" title="API documentation for ol.layer.Image">ol.layer.Image</a></li>,<li><a href="../apidoc/ol.layer.Tile.html" title="API documentation for ol.layer.Tile">ol.layer.Tile</a></li>,<li><a href="../apidoc/ol.source.Raster.html" title="API documentation for ol.source.Raster">ol.source.Raster</a></li>,<li><a href="../apidoc/ol.source.TileJSON.html" title="API documentation for ol.source.TileJSON">ol.source.TileJSON</a></li>,<li><a href="../apidoc/ol.source.XYZ.html" title="API documentation for ol.source.XYZ">ol.source.XYZ</a></li></ul></div>
|
|
</div>
|
|
</div>
|
|
|
|
<div class="row-fluid">
|
|
<div id="source-controls">
|
|
<a id="copy-button"><i class="fa fa-clipboard"></i> Copy</a>
|
|
<a id="jsfiddle-button"><i class="fa fa-jsfiddle"></i> Edit</a>
|
|
</div>
|
|
<form method="POST" id="jsfiddle-form" target="_blank" action="http://jsfiddle.net/api/post/jquery/1.11.0/">
|
|
<textarea class="hidden" name="js">// NOCOMPILE
|
|
/**
|
|
* Generates a shaded relief image given elevation data. Uses a 3x3
|
|
* neighborhood for determining slope and aspect.
|
|
* @param {Array.<ImageData>} inputs Array of input images.
|
|
* @param {Object} data Data added in the "beforeoperations" event.
|
|
* @return {ImageData} Output image.
|
|
*/
|
|
function shade(inputs, data) {
|
|
var elevationImage = inputs[0];
|
|
var width = elevationImage.width;
|
|
var height = elevationImage.height;
|
|
var elevationData = elevationImage.data;
|
|
var shadeData = new Uint8ClampedArray(elevationData.length);
|
|
var dp = data.resolution * 2;
|
|
var maxX = width - 1;
|
|
var maxY = height - 1;
|
|
var pixel = [0, 0, 0, 0];
|
|
var twoPi = 2 * Math.PI;
|
|
var halfPi = Math.PI / 2;
|
|
var sunEl = Math.PI * data.sunEl / 180;
|
|
var sunAz = Math.PI * data.sunAz / 180;
|
|
var cosSunEl = Math.cos(sunEl);
|
|
var sinSunEl = Math.sin(sunEl);
|
|
var pixelX, pixelY, x0, x1, y0, y1, offset,
|
|
z0, z1, dzdx, dzdy, slope, aspect, cosIncidence, scaled;
|
|
for (pixelY = 0; pixelY <= maxY; ++pixelY) {
|
|
y0 = pixelY === 0 ? 0 : pixelY - 1;
|
|
y1 = pixelY === maxY ? maxY : pixelY + 1;
|
|
for (pixelX = 0; pixelX <= maxX; ++pixelX) {
|
|
x0 = pixelX === 0 ? 0 : pixelX - 1;
|
|
x1 = pixelX === maxX ? maxX : pixelX + 1;
|
|
|
|
// determine elevation for (x0, pixelY)
|
|
offset = (pixelY * width + x0) * 4;
|
|
pixel[0] = elevationData[offset];
|
|
pixel[1] = elevationData[offset + 1];
|
|
pixel[2] = elevationData[offset + 2];
|
|
pixel[3] = elevationData[offset + 3];
|
|
z0 = data.vert * (pixel[0] + pixel[1] * 2 + pixel[2] * 3);
|
|
|
|
// determine elevation for (x1, pixelY)
|
|
offset = (pixelY * width + x1) * 4;
|
|
pixel[0] = elevationData[offset];
|
|
pixel[1] = elevationData[offset + 1];
|
|
pixel[2] = elevationData[offset + 2];
|
|
pixel[3] = elevationData[offset + 3];
|
|
z1 = data.vert * (pixel[0] + pixel[1] * 2 + pixel[2] * 3);
|
|
|
|
dzdx = (z1 - z0) / dp;
|
|
|
|
// determine elevation for (pixelX, y0)
|
|
offset = (y0 * width + pixelX) * 4;
|
|
pixel[0] = elevationData[offset];
|
|
pixel[1] = elevationData[offset + 1];
|
|
pixel[2] = elevationData[offset + 2];
|
|
pixel[3] = elevationData[offset + 3];
|
|
z0 = data.vert * (pixel[0] + pixel[1] * 2 + pixel[2] * 3);
|
|
|
|
// determine elevation for (pixelX, y1)
|
|
offset = (y1 * width + pixelX) * 4;
|
|
pixel[0] = elevationData[offset];
|
|
pixel[1] = elevationData[offset + 1];
|
|
pixel[2] = elevationData[offset + 2];
|
|
pixel[3] = elevationData[offset + 3];
|
|
z1 = data.vert * (pixel[0] + pixel[1] * 2 + pixel[2] * 3);
|
|
|
|
dzdy = (z1 - z0) / dp;
|
|
|
|
slope = Math.atan(Math.sqrt(dzdx * dzdx + dzdy * dzdy));
|
|
|
|
aspect = Math.atan2(dzdy, -dzdx);
|
|
if (aspect < 0) {
|
|
aspect = halfPi - aspect;
|
|
} else if (aspect > halfPi) {
|
|
aspect = twoPi - aspect + halfPi;
|
|
} else {
|
|
aspect = halfPi - aspect;
|
|
}
|
|
|
|
cosIncidence = sinSunEl * Math.cos(slope) +
|
|
cosSunEl * Math.sin(slope) * Math.cos(sunAz - aspect);
|
|
|
|
offset = (pixelY * width + pixelX) * 4;
|
|
scaled = 255 * cosIncidence;
|
|
shadeData[offset] = scaled;
|
|
shadeData[offset + 1] = scaled;
|
|
shadeData[offset + 2] = scaled;
|
|
shadeData[offset + 3] = elevationData[offset + 3];
|
|
}
|
|
}
|
|
|
|
return new ImageData(shadeData, width, height);
|
|
}
|
|
|
|
var elevation = new ol.source.XYZ({
|
|
url: 'https://{a-d}.tiles.mapbox.com/v3/aj.sf-dem/{z}/{x}/{y}.png',
|
|
crossOrigin: 'anonymous'
|
|
});
|
|
|
|
var raster = new ol.source.Raster({
|
|
sources: [elevation],
|
|
operationType: 'image',
|
|
operation: shade
|
|
});
|
|
|
|
var map = new ol.Map({
|
|
target: 'map',
|
|
layers: [
|
|
new ol.layer.Tile({
|
|
source: new ol.source.TileJSON({
|
|
url: 'http://api.tiles.mapbox.com/v3/tschaub.miapgppd.jsonp'
|
|
})
|
|
}),
|
|
new ol.layer.Image({
|
|
opacity: 0.3,
|
|
source: raster
|
|
})
|
|
],
|
|
view: new ol.View({
|
|
extent: [-13675026, 4439648, -13580856, 4580292],
|
|
center: [-13615645, 4497969],
|
|
minZoom: 10,
|
|
maxZoom: 16,
|
|
zoom: 13
|
|
})
|
|
});
|
|
|
|
var controlIds = ['vert', 'sunEl', 'sunAz'];
|
|
var controls = {};
|
|
controlIds.forEach(function(id) {
|
|
var control = document.getElementById(id);
|
|
var output = document.getElementById(id + 'Out');
|
|
control.addEventListener('input', function() {
|
|
output.innerText = control.value;
|
|
raster.changed();
|
|
});
|
|
output.innerText = control.value;
|
|
controls[id] = control;
|
|
});
|
|
|
|
raster.on('beforeoperations', function(event) {
|
|
// the event.data object will be passed to operations
|
|
var data = event.data;
|
|
data.resolution = event.resolution;
|
|
for (var id in controls) {
|
|
data[id] = Number(controls[id].value);
|
|
}
|
|
});
|
|
</textarea>
|
|
<textarea class="hidden" name="css">table.controls td {
|
|
text-align: center;
|
|
padding: 2px 5px;
|
|
}
|
|
</textarea>
|
|
<textarea class="hidden" name="html"><div class="row-fluid">
|
|
<div class="span12">
|
|
<div id="map" class="map"></div>
|
|
<table class="controls">
|
|
<tr>
|
|
<td>vertical exaggeration: <span id="vertOut"></span>x</td>
|
|
<td><input id="vert" type="range" min="1" max="5" value="1"/></td>
|
|
</tr>
|
|
<tr>
|
|
<td>sun elevation: <span id="sunElOut"></span>°</td>
|
|
<td><input id="sunEl" type="range" min="0" max="90" value="45"/></td>
|
|
</tr>
|
|
<tr>
|
|
<td>sun azimuth: <span id="sunAzOut"></span>°</td>
|
|
<td><input id="sunAz" type="range" min="0" max="360" value="45"/></td>
|
|
</tr>
|
|
</table>
|
|
</div>
|
|
</div>
|
|
</textarea>
|
|
<input type="hidden" name="wrap" value="l">
|
|
<input type="hidden" name="resources" value="https://maxcdn.bootstrapcdn.com/bootstrap/3.3.4/css/bootstrap.min.css,https://maxcdn.bootstrapcdn.com/bootstrap/3.3.4/js/bootstrap.min.js,http://openlayers.org/en/v3.9.0/css/ol.css,http://openlayers.org/en/v3.9.0/build/ol.js">
|
|
</form>
|
|
<pre><code id="example-source" class="language-markup"><!DOCTYPE html>
|
|
<html>
|
|
<head>
|
|
<title>Shaded Relief</title>
|
|
<script src="https://code.jquery.com/jquery-1.11.2.min.js"></script>
|
|
<link rel="stylesheet" href="https://maxcdn.bootstrapcdn.com/bootstrap/3.3.5/css/bootstrap.min.css">
|
|
<script src="https://maxcdn.bootstrapcdn.com/bootstrap/3.3.5/js/bootstrap.min.js"></script>
|
|
<link rel="stylesheet" href="http://openlayers.org/en/v3.9.0/css/ol.css" type="text/css">
|
|
<script src="http://openlayers.org/en/v3.9.0/build/ol.js"></script>
|
|
|
|
<style>
|
|
table.controls td {
|
|
text-align: center;
|
|
padding: 2px 5px;
|
|
}
|
|
|
|
</style>
|
|
</head>
|
|
<body>
|
|
<div class="container-fluid">
|
|
|
|
<div class="row-fluid">
|
|
<div class="span12">
|
|
<div id="map" class="map"></div>
|
|
<table class="controls">
|
|
<tr>
|
|
<td>vertical exaggeration: <span id="vertOut"></span>x</td>
|
|
<td><input id="vert" type="range" min="1" max="5" value="1"/></td>
|
|
</tr>
|
|
<tr>
|
|
<td>sun elevation: <span id="sunElOut"></span>°</td>
|
|
<td><input id="sunEl" type="range" min="0" max="90" value="45"/></td>
|
|
</tr>
|
|
<tr>
|
|
<td>sun azimuth: <span id="sunAzOut"></span>°</td>
|
|
<td><input id="sunAz" type="range" min="0" max="360" value="45"/></td>
|
|
</tr>
|
|
</table>
|
|
</div>
|
|
</div>
|
|
|
|
</div>
|
|
<script>
|
|
// NOCOMPILE
|
|
/**
|
|
* Generates a shaded relief image given elevation data. Uses a 3x3
|
|
* neighborhood for determining slope and aspect.
|
|
* @param {Array.<ImageData>} inputs Array of input images.
|
|
* @param {Object} data Data added in the "beforeoperations" event.
|
|
* @return {ImageData} Output image.
|
|
*/
|
|
function shade(inputs, data) {
|
|
var elevationImage = inputs[0];
|
|
var width = elevationImage.width;
|
|
var height = elevationImage.height;
|
|
var elevationData = elevationImage.data;
|
|
var shadeData = new Uint8ClampedArray(elevationData.length);
|
|
var dp = data.resolution * 2;
|
|
var maxX = width - 1;
|
|
var maxY = height - 1;
|
|
var pixel = [0, 0, 0, 0];
|
|
var twoPi = 2 * Math.PI;
|
|
var halfPi = Math.PI / 2;
|
|
var sunEl = Math.PI * data.sunEl / 180;
|
|
var sunAz = Math.PI * data.sunAz / 180;
|
|
var cosSunEl = Math.cos(sunEl);
|
|
var sinSunEl = Math.sin(sunEl);
|
|
var pixelX, pixelY, x0, x1, y0, y1, offset,
|
|
z0, z1, dzdx, dzdy, slope, aspect, cosIncidence, scaled;
|
|
for (pixelY = 0; pixelY <= maxY; ++pixelY) {
|
|
y0 = pixelY === 0 ? 0 : pixelY - 1;
|
|
y1 = pixelY === maxY ? maxY : pixelY + 1;
|
|
for (pixelX = 0; pixelX <= maxX; ++pixelX) {
|
|
x0 = pixelX === 0 ? 0 : pixelX - 1;
|
|
x1 = pixelX === maxX ? maxX : pixelX + 1;
|
|
|
|
// determine elevation for (x0, pixelY)
|
|
offset = (pixelY * width + x0) * 4;
|
|
pixel[0] = elevationData[offset];
|
|
pixel[1] = elevationData[offset + 1];
|
|
pixel[2] = elevationData[offset + 2];
|
|
pixel[3] = elevationData[offset + 3];
|
|
z0 = data.vert * (pixel[0] + pixel[1] * 2 + pixel[2] * 3);
|
|
|
|
// determine elevation for (x1, pixelY)
|
|
offset = (pixelY * width + x1) * 4;
|
|
pixel[0] = elevationData[offset];
|
|
pixel[1] = elevationData[offset + 1];
|
|
pixel[2] = elevationData[offset + 2];
|
|
pixel[3] = elevationData[offset + 3];
|
|
z1 = data.vert * (pixel[0] + pixel[1] * 2 + pixel[2] * 3);
|
|
|
|
dzdx = (z1 - z0) / dp;
|
|
|
|
// determine elevation for (pixelX, y0)
|
|
offset = (y0 * width + pixelX) * 4;
|
|
pixel[0] = elevationData[offset];
|
|
pixel[1] = elevationData[offset + 1];
|
|
pixel[2] = elevationData[offset + 2];
|
|
pixel[3] = elevationData[offset + 3];
|
|
z0 = data.vert * (pixel[0] + pixel[1] * 2 + pixel[2] * 3);
|
|
|
|
// determine elevation for (pixelX, y1)
|
|
offset = (y1 * width + pixelX) * 4;
|
|
pixel[0] = elevationData[offset];
|
|
pixel[1] = elevationData[offset + 1];
|
|
pixel[2] = elevationData[offset + 2];
|
|
pixel[3] = elevationData[offset + 3];
|
|
z1 = data.vert * (pixel[0] + pixel[1] * 2 + pixel[2] * 3);
|
|
|
|
dzdy = (z1 - z0) / dp;
|
|
|
|
slope = Math.atan(Math.sqrt(dzdx * dzdx + dzdy * dzdy));
|
|
|
|
aspect = Math.atan2(dzdy, -dzdx);
|
|
if (aspect < 0) {
|
|
aspect = halfPi - aspect;
|
|
} else if (aspect > halfPi) {
|
|
aspect = twoPi - aspect + halfPi;
|
|
} else {
|
|
aspect = halfPi - aspect;
|
|
}
|
|
|
|
cosIncidence = sinSunEl * Math.cos(slope) +
|
|
cosSunEl * Math.sin(slope) * Math.cos(sunAz - aspect);
|
|
|
|
offset = (pixelY * width + pixelX) * 4;
|
|
scaled = 255 * cosIncidence;
|
|
shadeData[offset] = scaled;
|
|
shadeData[offset + 1] = scaled;
|
|
shadeData[offset + 2] = scaled;
|
|
shadeData[offset + 3] = elevationData[offset + 3];
|
|
}
|
|
}
|
|
|
|
return new ImageData(shadeData, width, height);
|
|
}
|
|
|
|
var elevation = new ol.source.XYZ({
|
|
url: 'https://{a-d}.tiles.mapbox.com/v3/aj.sf-dem/{z}/{x}/{y}.png',
|
|
crossOrigin: 'anonymous'
|
|
});
|
|
|
|
var raster = new ol.source.Raster({
|
|
sources: [elevation],
|
|
operationType: 'image',
|
|
operation: shade
|
|
});
|
|
|
|
var map = new ol.Map({
|
|
target: 'map',
|
|
layers: [
|
|
new ol.layer.Tile({
|
|
source: new ol.source.TileJSON({
|
|
url: 'http://api.tiles.mapbox.com/v3/tschaub.miapgppd.jsonp'
|
|
})
|
|
}),
|
|
new ol.layer.Image({
|
|
opacity: 0.3,
|
|
source: raster
|
|
})
|
|
],
|
|
view: new ol.View({
|
|
extent: [-13675026, 4439648, -13580856, 4580292],
|
|
center: [-13615645, 4497969],
|
|
minZoom: 10,
|
|
maxZoom: 16,
|
|
zoom: 13
|
|
})
|
|
});
|
|
|
|
var controlIds = ['vert', 'sunEl', 'sunAz'];
|
|
var controls = {};
|
|
controlIds.forEach(function(id) {
|
|
var control = document.getElementById(id);
|
|
var output = document.getElementById(id + 'Out');
|
|
control.addEventListener('input', function() {
|
|
output.innerText = control.value;
|
|
raster.changed();
|
|
});
|
|
output.innerText = control.value;
|
|
controls[id] = control;
|
|
});
|
|
|
|
raster.on('beforeoperations', function(event) {
|
|
// the event.data object will be passed to operations
|
|
var data = event.data;
|
|
data.resolution = event.resolution;
|
|
for (var id in controls) {
|
|
data[id] = Number(controls[id].value);
|
|
}
|
|
});
|
|
|
|
</script>
|
|
</body>
|
|
</html></code></pre>
|
|
</div>
|
|
</div>
|
|
|
|
<script src="https://code.jquery.com/jquery-1.11.2.min.js"></script>
|
|
<script src="https://maxcdn.bootstrapcdn.com/bootstrap/3.3.5/js/bootstrap.min.js"></script>
|
|
<script src="./resources/common.js"></script>
|
|
<script src="./resources/prism/prism.min.js"></script>
|
|
<script src="loader.js?id=shaded-relief"></script>
|
|
|
|
</body>
|
|
</html>
|