Geozarr Stretch

<!--
  说明：波段选择与对比度拉伸的中文示例，数据在页面内生成，可独立运行。
  核心逻辑：多波段数组 + 合成/灰度/色带三种映射 + 统一拉伸区间
  依赖：OpenLayers 10.3.1 · jsDelivr CDN
-->
<!DOCTYPE html>
<html lang="zh-CN">
<head>
<meta name="robots" content="noindex,follow">
<meta charset="utf-8">
<meta name="viewport" content="width=device-width, initial-scale=1">
<title>波段选择与对比度拉伸</title>
<link rel="stylesheet" href="https://cdn.jsdelivr.net/npm/ol@10.3.1/ol.css">
<style>
  html,body{height:100%;margin:0}
  body{display:flex;flex-direction:column;overflow:hidden;
       font:13px/1.6 -apple-system,"Segoe UI","PingFang SC",sans-serif}
  .bar{flex:0 0 auto;padding:8px 12px;background:#f6f7f9;border-bottom:1px solid #e3e6ea;
       display:flex;gap:10px;align-items:center;flex-wrap:wrap;font-size:12px;color:#58606e}
  .bar label{display:flex;gap:5px;align-items:center;cursor:pointer}
  .bar button{font:inherit;padding:3px 9px;border:1px solid #c9ced6;background:#fff;
       border-radius:5px;cursor:pointer;color:#333}
  .bar button:hover{background:#eef2f8}
  .bar input[type=range]{width:110px;vertical-align:middle}
  .bar input[type=text]{font:inherit;padding:2px 6px;border:1px solid #c9ced6;border-radius:4px}
  #map{flex:1 1 auto;min-height:0;background:#eef0f3;position:relative}
  .panel{flex:0 0 auto;max-height:34%;overflow:auto;padding:8px 12px;font-size:12px;
       background:#fff;border-top:1px solid #e3e6ea;color:#3d4451}
  .panel pre{margin:0;font-family:Consolas,"Courier New",monospace;font-size:11px;
       white-space:pre-wrap;word-break:break-all;line-height:1.5}
  .ol-attribution{font-size:11px}

</style>
</head>
<body>
<div class="bar">
<label>显示方式 <select id="mode">
  <option value="rgb">真彩色</option>
  <option value="single">单波段灰度</option>
  <option value="diff">归一化差异</option>
</select></label>
<label>拉伸 <input type="range" id="stretch" min="0" max="100" value="40"></label>
<span id="info">—</span>
</div>
<div id="map"></div>

<script src="https://cdn.jsdelivr.net/npm/ol@10.3.1/dist/ol.js"></script>
<script>
(function () {
  'use strict';

  if (typeof ol === 'undefined') {
    var errBox = document.getElementById('map');
    if (errBox) {
      errBox.innerHTML = '<div style="padding:20px;color:#b42318">OpenLayers 加载失败，请检查能否访问 cdn.jsdelivr.net</div>';
    }
    return;
  }

  var fromLonLat = ol.proj.fromLonLat;
  var N = 256;

  function makeBand(fx, fy, phase) {
    var arr = new Float32Array(N * N);
    for (var y = 0; y < N; y++) {
      for (var x = 0; x < N; x++) {
        var v = 0.5 + 0.3 * Math.sin(x / fx + phase) * Math.cos(y / fy) +
                0.15 * Math.sin((x + y) / 17 + phase);
        arr[y * N + x] = Math.min(1, Math.max(0, v));
      }
    }
    return arr;
  }

  var b1 = makeBand(41, 57, 0.4);
  var b2 = makeBand(29, 33, 1.1);
  var b3 = makeBand(61, 23, 2.3);

  var canvas = document.createElement('canvas');
  canvas.width = canvas.height = N;
  var ctx = canvas.getContext('2d');
  var image = ctx.createImageData(N, N);

  var extent = ol.proj.transformExtent([115.8, 39.4, 117.0, 40.4], 'EPSG:4326', 'EPSG:3857');
  var layer = new ol.layer.Image({
    source: new ol.source.ImageStatic({
      url: canvas.toDataURL('image/png'),
      imageExtent: extent,
      projection: 'EPSG:3857'
    })
  });

  var map = new ol.Map({
    target: 'map',
    layers: [new ol.layer.Tile({source: new ol.source.XYZ({url: 'https://server.arcgisonline.com/ArcGIS/rest/services/World_Street_Map/MapServer/tile/{z}/{y}/{x}', attributions: 'Tiles © Esri'})}), layer],
    view: new ol.View({center: fromLonLat([116.4, 39.9]), zoom: 11}),
    controls: ol.control.defaults.defaults({attribution: false, zoom: true, rotate: false})
  });

  var mode = document.getElementById('mode');
  var stretch = document.getElementById('stretch');
  var info = document.getElementById('info');

  function ramp(t) {
    if (t < 0.5) {
      var k = t / 0.5;
      return [0, Math.round(140 + 90 * k), Math.round(210 * (1 - k))];
    }
    var k2 = (t - 0.5) / 0.5;
    return [Math.round(255 * k2), Math.round(215 - 25 * k2), 20];
  }

  function render() {
    var half = (+stretch.value) / 200;
    var lo = 0.5 - half;
    var hi = 0.5 + half;
    var kind = mode.value;

    for (var i = 0; i < N * N; i++) {
      var r, g, b;
      if (kind === 'rgb') {
        r = Math.round(clamp((b1[i] - lo) / (hi - lo)) * 255);
        g = Math.round(clamp((b2[i] - lo) / (hi - lo)) * 255);
        b = Math.round(clamp((b3[i] - lo) / (hi - lo)) * 255);
      } else if (kind === 'single') {
        var v = Math.round(clamp((b1[i] - lo) / (hi - lo)) * 255);
        r = g = b = v;
      } else {
        var d = (b2[i] - b1[i]) / (b2[i] + b1[i] + 0.02);
        var c = ramp(clamp(d * 1.8 + 0.5));
        r = c[0]; g = c[1]; b = c[2];
      }
      image.data[i * 4] = r;
      image.data[i * 4 + 1] = g;
      image.data[i * 4 + 2] = b;
      image.data[i * 4 + 3] = 255;
    }

    ctx.putImageData(image, 0, 0);
    layer.setSource(new ol.source.ImageStatic({
      url: canvas.toDataURL('image/png'),
      imageExtent: extent,
      projection: 'EPSG:3857'
    }));

    var label = kind === 'rgb' ? '真彩色（R=b1, G=b2, B=b3）'
      : (kind === 'single' ? '单波段灰度（b1）' : '归一化差异 ((b2-b1)/(b2+b1))');
    info.textContent = label + ' · 拉伸 [' + lo.toFixed(2) + ', ' + hi.toFixed(2) + ']';
  }

  function clamp(v) {
    return v < 0 ? 0 : (v > 1 ? 1 : v);
  }

  mode.addEventListener('change', render);
  stretch.addEventListener('input', render);
  render();
})();
</script>
</body>
</html>