传单-使用图块图层时如何保持原始坐标系有效?

时间:2019-03-19 10:01:12

标签: javascript leaflet gdal tiling cartesian-coordinates

在我的Leaflet应用程序中,我曾经显示带有ImageOverlay的背景层。但是,由于图像太大,导致Leaflet处理变慢,因此我改用平铺方法。 我使用gdal2tiles-leaflet来生成我的图块。工作正常。

但是现在投影背景层的笛卡尔坐标系(我使用Leaflet CRS Simple)不再有效。

这是我的图像规格:

  • 宽度:21002像素
  • 高度:14694像素
  • 分辨率:0.02(1米= 50像素)

这是使用 ImageOverlay 时我的图像的边界: enter image description here

这是使用 TileLayer 时我的图像的边界: enter image description here

这是我一起显示它们时的样子: enter image description here

有人知道发生了什么事吗?

1 个答案:

答案 0 :(得分:0)

如@TomazicM所建议,我研究了leaflet-rastercoords插件,该插件是 gdal2tiles-leaflet 的补充。基于此,我实现了一个API,该API允许以栅格笛卡尔方式转换坐标。

为此,我首先将 leaflet-rastercoords 代码提供给一个方法:

export let loadLeafletRastercoordsPlugin = (L: any) => {
    // rastercoords.js source code
}

然后我写了一个类来处理坐标转换:

import * as L from 'leaflet';
import { loadLeafletRastercoordsPlugin } from './leaflet-rastercoords';

export class RasterCoords {

    public rc: any;
    public map: L.Map;
    public rasterWidth: number;
    public rasterHeight: number;
    public resolution: number;

    constructor(map: L.Map, rasterWidth: number, rasterHeight: number, resolution: number) {
        loadLeafletRastercoordsPlugin(L);
        this.rc = new L['RasterCoords'](map, [rasterWidth, rasterHeight]);
        this.map = map;
        this.rasterWidth = rasterWidth;
        this.rasterHeight = rasterHeight;
        this.resolution = resolution;
    }

}

使用一种方法将光栅坐标根据其宽度,高度和分辨率并以自下而上的方式将Y轴投影到图像的原始正交平面中:

public project(coordinates: L.LatLngTuple): L.LatLngTuple {
    coordinates = this.applyResolution(coordinates);
    const projectedCoordinates = this.rc.project(
        coordinates
    );

    return this.applyCartesianProjection([projectedCoordinates.y, projectedCoordinates.x] as L.LatLngTuple);
}

private applyResolution(coordinates: L.LatLngTuple): L.LatLngTuple {
    return coordinates.map((v: number) => v * this.resolution) as L.LatLngTuple;
}

private applyCartesianProjection(coordinates: L.LatLngTuple): L.LatLngTuple {
    return [(this.rasterHeight * this.resolution) - coordinates[0], coordinates[1]];
}

并使用一种方法“ 取消投影”笛卡尔坐标(即逐点对 project 方法执行的操作进行回溯):

public unproject(coordinates: L.LatLngTuple): L.LatLngTuple {
    coordinates = this.unapplyResolution(this.unapplyCartesianProjection(coordinates));

    return this.rc.unproject([coordinates[1], coordinates[0]]);
}

private unapplyResolution(coordinates: L.LatLngTuple): L.LatLngTuple {
    return coordinates.map((v: number) => v / this.resolution) as L.LatLngTuple;
}

private unapplyCartesianProjection(coordinates: L.LatLngTuple): L.LatLngTuple {
    return [Math.abs(coordinates[0] - (this.rasterHeight * this.resolution)), coordinates[1]];
}

然后,该API可帮助我根据其笛卡尔坐标将对象有效地添加到地图中:

const imageWidth = 21002;
const imageHeight = 14694;
const imageResolution = 0.02;
const map = L.map('map');
const rc = new RasterCoords(map, imageWidth, imageHeight, imageResolution);
map.setView(rc.unproject([imageWidth, imageHeight]), 2);
L.tileLayer('./image/{z}/{x}/{y}.png', {
    noWrap: true
}).addTo(map);
new L.CircleMarker(this.rc.unproject([293, 420])).addTo(map);