光栅重投影 OpenLayers 3

Question

我正在尝试将 EPSG:3395 中的 WMS 源重新投影到 OpenLayers 3 中的 EPSG:32661。我有 proj4js 并加载了投影定义，但没有显示任何内容。

奇怪的是，甚至http://openlayers.org/en/v3.12.1/examples/reprojection-by-code.html 的示例似乎都不起作用。

Answer 1

我正在调试为什么会中断并发现将以下坐标从 32661 转换为 3857 时会引发第一个错误。假设 EPSG:32661 proj4js def 已加载，您可以通过在控制台中运行以下 JS 来单独重新创建此错误：

proj4("EPSG:32661", "EPSG:3857", [2000000.0000000005, 2000000])

-> 休息

然而，对于 WGS84，

proj4("EPSG:32661", "EPSG:4326", [2000000.0000000005, 2000000])

-> [90, 90]

我通过以下查询验证了 proj.4 本身是相同的：

https://mygeodata.cloud/cs2cs/

我不知道这是正常还是错误。有趣的是，它基本上发生在北极，所以在这个极端点可能存在一些舍入错误或数学错误？

Answer 2

当预期包含 x 和 y 属性的对象时，由 OpenLayers 调用的 proj4 代码处理空值时会发生错误。根据使用的 proj4.js 的版本，它可能发生在任何转换中，其中一个坐标在另一个投影中无效，但在转换到 EPSG 或从 EPSG 转换时似乎不太可能：4326 OpenLayers 示例试图通过设置范围来防止错误在投影上，但是如果您想在其国家投影中显示一个国家，同时也显示周边国家（尽管有一些失真），或者在 UTM 投影中包括相邻区域，那么这是不可取的。我使用ol.proj.addCoordinateTransforms 来确保所有转换都通过 EPSG:4326 消除大部分错误并添加 try/catch 以覆盖任何剩余的错误（尽管如果使用 proj4.js 的最新版本 2.5.0，try/catch 似乎冗余，即使该版本在其他投影之间直接转换也会产生错误）。然而，ol.proj.getPointResolution() 所需的从 EPSG:2163 到 EPSG:4326 的一些转换即使使用最新的 proj4.js 也会失败，这使得在某些位置无法居中或缩小地图。幸运的是 EPSG:4269 是在 proj4 中预定义的，当用作中间而不是 EPSG:4326 时，它允许对涉及 EPSG:4326 的错误进行 try/catch 处理（除非应该使用 EPSG:4326 时重新投影到 EPSG:4269 或从 EPSG:4269 重新投影！ ) 这是对 OpenLayers 示例的重写，所有转换均通过 EPSG:4269 进行，并且对投影范围没有限制

function reprojectionErrorHandler (projections, opt_intermediate) {

var intermediate = opt_intermediate || 'EPSG:4269';

function transform (projA, projB) {

    return function (input, opt_output, opt_dimension) {
        var length = input.length;
        var dimension = opt_dimension !== undefined ? opt_dimension : 2;
        var output = opt_output !== undefined ? opt_output : new Array(length);
        var ll, point, i, j;
        try {
            for (i = 0; i < length; i += dimension) {
                ll = ol.proj.transform([input[i], input[i + 1]], projA, intermediate);
                point = ol.proj.transform([ll[i], ll[i + 1]], intermediate, projB);
                output[i] = point[0];
                output[i + 1] = point[1];
                for (j = dimension - 1; j >= 2; --j) {
                    output[i + j] = input[i + j];
                }
            }
        } catch (e) {}
        return output;
    };

}

if (Array.isArray(projections)) {
    for (i = 0; i < projections.length-1; i++) {
        for (j = i+1; j < projections.length; j++) {
            if (ol.proj.get(projections[i]).getCode() != ol.proj.get(projections[j]).getCode() &&
                ol.proj.get(projections[i]).getCode() != ol.proj.get(intermediate).getCode() &&
                ol.proj.get(projections[j]).getCode() != ol.proj.get(intermediate).getCode() ) {

                ol.proj.addCoordinateTransforms(
                    projections[i],
                    projections[j],
                    transform(projections[i], projections[j]),
                    transform(projections[j], projections[i])
                );

                ol.proj.addCoordinateTransforms(
                    projections[j],
                    projections[i],
                    transform(projections[j], projections[i]),
                    transform(projections[i], projections[j])
                );

            }
        }
    }
}

}

  var map = new ol.Map({
    layers: [
      new ol.layer.Tile({
        source: new ol.source.OSM()
      })
    ],
    target: 'map',
    view: new ol.View({
      projection: 'EPSG:3857',
      center: [0, 0],
      zoom: 1
    })
  });


  var queryInput = document.getElementById('epsg-query');
  var searchButton = document.getElementById('epsg-search');
  var resultSpan = document.getElementById('epsg-result');
  var renderEdgesCheckbox = document.getElementById('render-edges');

  function setProjection(code, name, proj4def, bbox) {
    if (code === null || name === null || proj4def === null || bbox === null) {
      resultSpan.innerHTML = 'Nothing usable found, using EPSG:3857...';
      map.setView(new ol.View({
        projection: 'EPSG:3857',
        center: [0, 0],
        zoom: 1
      }));
      return;
    }

    resultSpan.innerHTML = '(' + code + ') ' + name;

    var newProjCode = 'EPSG:' + code;
    if (newProjCode != 'EPSG:4269') {
      // 'EPSG:4269' is predefined in proj4 but not OL so proj4 must be registered if using OL5
      proj4.defs(newProjCode, proj4def);
    }
    if (ol.proj.proj4 && ol.proj.proj4.register) { ol.proj.proj4.register(proj4); };
    var newProj = ol.proj.get(newProjCode);
   // var fromLonLat = ol.proj.getTransform('EPSG:4326', newProj);

    if (newProjCode != 'EPSG:4326' && newProjCode != 'EPSG:3857') {
      reprojectionErrorHandler(['EPSG:4326', newProj]);
      reprojectionErrorHandler(['EPSG:3857', newProj]);
    }

    // very approximate calculation of projection extent
   // var extent = ol.extent.applyTransform(
   //     [bbox[1], bbox[2], bbox[3], bbox[0]], fromLonLat);
   // newProj.setExtent(extent);
    var newView = new ol.View({
      projection: newProj,
      center: [0, 0],
      zoom: 4
    });
    map.setView(newView);
   // newView.fit(extent);
  }


  function search(query) {
    resultSpan.innerHTML = 'Searching ...';
    fetch('https://epsg.io/?format=json&q=' + query).then(function(response) {
      return response.json();
    }).then(function(json) {
      var results = json['results'];
      if (results && results.length > 0) {
        for (var i = 0, ii = results.length; i < ii; i++) {
          var result = results[i];
          if (result) {
            var code = result['code'], name = result['name'],
                proj4def = result['proj4'], bbox = result['bbox'];
           // if (code && code.length > 0 && proj4def && proj4def.length > 0 &&
           //     bbox && bbox.length == 4) {
              setProjection(code, name, proj4def, bbox);
              return;
           // }
          }
        }
      }
      setProjection(null, null, null, null);
    });
  }


  /**
   * Handle click event.
   * @param {Event} event The event.
   */
  searchButton.onclick = function(event) {
    search(queryInput.value);
    event.preventDefault();
  };


  /**
   * Handle change event.
   */
  renderEdgesCheckbox.onchange = function() {
    map.getLayers().forEach(function(layer) {
      if (layer instanceof ol.layer.Tile) {
        var source = layer.getSource();
        if (source instanceof ol.source.TileImage) {
          source.setRenderReprojectionEdges(renderEdgesCheckbox.checked);
        }
      }
    });
  };

html,
body {
  height: 100%;
  width: 100%;
  padding: 0px;
  margin: 0px;
}

.map {
  width: 100%;
}

<link rel="stylesheet" href="https://openlayers.org/en/v4.6.5/css/ol.css" type="text/css">
<!-- The line below is only needed for old environments like Internet Explorer and Android 4.x -->
<script src="https://cdn.polyfill.io/v2/polyfill.min.js?features=requestAnimationFrame,Element.prototype.classList,URL"></script>
<script src="https://openlayers.org/en/v4.6.5/build/ol.js"></script>
<script src="https://cdnjs.cloudflare.com/ajax/libs/proj4js/2.5.0/proj4.js"></script>
<div id="map" class="map"></div>
<form class="form-inline">
  <label for="epsg-query">Search projection:</label>
  <input type="text" id="epsg-query" placeholder="4326, 27700, US National Atlas, Swiss, France, ..." class="form-control" size="50" />
  <button id="epsg-search" class="btn">Search</button>
  <span id="epsg-result"></span>
  <div>
    <label for="render-edges">
      Render reprojection edges
      <input type="checkbox" id="render-edges">
    </label>
  </div>
</form>