OpenLayers 3重新投影矢量图层(v4.3.2)

时间:2017-09-05 08:59:18

标签: jquery openlayers openlayers-3 geojson proj

向GIS爱好者和WEB制造商致意。我有一个关于预测的问题,更确切地说,从EPSG重新投影:3765到EPSG:Openlayers 4中有3857个。

在下面的代码中,我试图从EPSG重新投影:3765到EPSG:3857:

    proj4.defs(
        'EPSG:3765',
        '+proj=tmerc +lat_0=0 +lon_0=16.5 +k=0.9999 +x_0=500000 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs'
    );

var map = new ol.Map({
    view: new ol.View({
        zoom: 6,
        center: ol.proj.transform([461152.65, 5108327.47], 'EPSG:3765', 'EPSG:3857')
    }),
    target: 'js-map',
    layers: [
        new ol.layer.Tile({
            source: new ol.source.OSM()
        }),
        new ol.layer.Vector({
            source: new ol.source.Vector({
                url: 'test_1.geojson',
                format: new ol.format.GeoJSON({
                    dataProjection: 'EPSG:3765',
                    featureProjection: 'EPSG:3857'
                })
            })

        })
    ]
});

geojson文件中的坐标以米为单位,如下所示:

"coordinates": [ [ [ 461117.98, 5108379.85 ], [ 461124.53, 5108368.39 ], [ 461132.37, 5108354.26 ], [ 461141.13, 5108341.08 ...

我成功渲染了geojson,但位置错误,如下所示: enter image description here

我希望它在这里,就像这样: enter image description here

我做错了什么?是geojson坐标格式的问题还是我需要以某种方式重新投影geojson就像在第二张图片中一样?

更新

从这里你可以下载geojson并测试你是否喜欢 https://www.dropbox.com/s/ih1bh8bj4zzgutc/test_1.geojson?dl=0

使用pavlos的例子我编辑了几行,或者更精确;我正在尝试使用此代码块加载本地geojson

var geojsonObject = $.ajax({
  dataType: "json",
  url: "test_1.geojson"
});

但是我一直在犯这种错误:

Uncaught TypeError: (0 , Tm[a.type]) is not a function

代码如下:

proj4.defs(
    'EPSG:3765',
    '+proj=tmerc +lat_0=0 +lon_0=16.5 +k=0.9999 +x_0=500000 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs'
);

var geojsonObject = $.ajax({
    dataType: "json",
    url: "test_1.geojson"
  });
console.log(geojsonObject);

var styles = {
    'Polygon': new ol.style.Style({
        stroke: new ol.style.Stroke({
            color: 'red',
            width: 5
        }),
        fill: new ol.style.Fill({
            color: 'rgba(255, 0, 0, 1)'
        })
    })
};

var styleFunction = function(feature) {
    return styles[feature.getGeometry().getType()];
};



var vectorSource = new ol.source.Vector({
    features: (new ol.format.GeoJSON()).readFeatures(geojsonObject, {
        defaultDataProjection: ol.proj.get('EPSG:3765'),
        featureProjection: 'EPSG:3857'
    })
});
var vectorLayer = new ol.layer.Vector({
    source: vectorSource,
    style: styleFunction
});
var map = new ol.Map({
    layers: [
        new ol.layer.Tile({
            source: new ol.source.OSM()
        }),
        vectorLayer
    ],
    target: 'map',
    controls: ol.control.defaults({
        attributionOptions: /** @type {olx.control.AttributionOptions} */ ({
            collapsible: false
        })
    }),
    view: new ol.View({
        center: ol.proj.transform([461117.98, 5108379.85], 'EPSG:3765', 'EPSG:3857'),
        zoom: 12
    })
});

更新3

我相信我更近了一步,但仍然没有运气

我用jQuery和AJAX请求编写了一些代码,并在loader属性中设置了函数。现在代码看起来像这样:

var vectorLayer = new ol.layer.Vector({
    source: new ol.source.Vector({
        loader: function() {
            $.ajax({
                type: 'GET',
                url:'test_1.geojson',
                context: this
            }).done(function(data){
                var format = new ol.format.GeoJSON();
                this.addFeatures(format.readFeatures(data, {
                    defaultDataProjection: ol.proj.get('EPSG:3765'),
                    featureProjection: 'EPSG:3857'
                }));

            });
        }
    })
});
map.addLayer(vectorLayer);

但现在我无法看到我的GeoJSON功能,它似乎已加载,但它无法在地图上呈现。

我必须指出,如果我从readFeatures方法

中评论这两行
defaultDataProjection: ol.proj.get('EPSG:3765'),
featureProjection: 'EPSG:3857'

geojson确实渲染但在错误的位置,就像在第一张图片中一样。请帮忙......我真的不再流泪......

2 个答案:

答案 0 :(得分:1)

我认为你的问题是坐标或源EPSG。

看看https://epsg.io/3765,该投影的边界位于意大利东部,而您希望您的几何形状位于意大利西部。

使用该投影,x坐标应为负。

我的猜测是你的geojson文件没有使用EPSG:3765

答案 1 :(得分:1)

我不确定您的问题是什么,但也许您的geojson文件不包含crs定义。也许proj4js版本不正确。

检查这个小提琴 - > https://jsfiddle.net/p_tsagkis/s9vtoaak/。 我正在使用你的投影,似乎在正确的位置。

更新: 您还需要稍微更改geojson文件。 投影的定义应如下:

{"type":"FeatureCollection",
"crs": {
          "type": "name",
          "properties": {
            "name": "EPSG:3765"
          }
        },
        "features":[..........

如果您无法编辑geojson文件,还有一个选项是将geojson中存在的名称注册到ol3 js部分。 像这样:

proj4.defs('urn:ogc:def:crs:EPSG::3765', proj4.defs('EPSG:3765'));

或将投影定义js行替换为:

proj4.defs("urn:ogc:def:crs:EPSG::3765","+proj=tmerc +lat_0=0 +lon_0=16.5 +k=0.9999 +x_0=500000 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs");

我没有测试它但它应该可以工作。

ol3不知道urn:ogc:def:crs:EPSG::3765 你只需要proj4.defs("EPSG:3765","+proj=tmerc +lat_0=0 +lon_0=16.5 +k=0.9999 +x_0=500000 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs"); 所以你需要让它知道文件中存在的投影名称。