地图中坐标系是非常多的,最常用的是EPSG:3857(等于谷歌的900913,等于esri的102100)的web墨卡托投影和GPS的EPSG:4326的WGS84坐标系。业务场景中,很多时候并不是这些常用坐标系,比如要使用北京54,西安80等坐标系怎么办咧?本文仅仅抛砖引玉,介绍方法,阐述如何应用自定义坐标系。
一 定义坐标系
每个坐标系都是有一个EPSG定义,本文准备以EPSG:3395举例。
首先在http://epsg.io/网站,查找坐标系定义。
这里我们是定义ol3的,于是选择proj4js的,复制粘贴下来。
var projection_3395 = new ol.proj.Projection({
code: 'EPSG:3395',
extent:[-20026376.39,-15496570.74,20026376.39,18764656.23],
units: 'm',
axisOrientation: 'neu'
});
//定义3395坐标系,且与其他4326,3857的互相转换
proj4.defs("EPSG:3395","+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs");
二 坐标系转换
说白了,其实是通过proj4对坐标系转换,然后作为匿名回调重写ol.proj的坐标转换功能,以下代码定义了3395转4326,3857方法。
//结合proj4在ol3中自定义坐标系
ol.proj.addProjection(projection_3395);
ol.proj.addCoordinateTransforms("EPSG:4326", "EPSG:3395",
function(coordinate) {
return proj4("EPSG:4326","EPSG:3395",coordinate);
},
function(coordinate) {
return proj4("EPSG:3395","EPSG:4326",coordinate);;
}
);
ol.proj.addCoordinateTransforms("EPSG:3857", "EPSG:3395",
function(coordinate) {
return proj4("EPSG:3857","EPSG:3395",coordinate);
},
function(coordinate) {
return proj4("EPSG:3395","EPSG:3857",coordinate);;
}
);
三 测试坐标系转换
ol.proj.transform([118,32],'EPSG:4326','EPSG:3395');
//这种如果测试成功,代表定义成功了。
四 加载自定义坐标系图层
//定义自定义坐标系
proj4.defs("EPSG:2384","+proj=tmerc +lat_0=0 +lon_0=117 +k=1 +x_0=500000 +y_0=0 +a=6378140 +b=6356755.288157528 +units=m +no_defs");
var projection = new ol.proj.Projection({
code: 'EPSG:2384',
extent : [345754.23,2501018.29,607809.05,5528581.54],
units: 'm',
axisOrientation: 'neu',
global: false
});
//结合proj4在ol3中自定义坐标系
ol.proj.addProjection(projection);
ol.proj.addCoordinateTransforms("EPSG:4326", "EPSG:2384",
function (coordinate) {
return proj4("EPSG:4326", "EPSG:2384", coordinate);
},
function (coordinate) {
return proj4("EPSG:2384", "EPSG:4326", coordinate);
}
);
//测试坐标系转换
var coor1 = ol.proj.transform([118.8260199966696, 32.03620907291769], 'EPSG:4326', 'EPSG:2384');
console.log(coor1);
var format = 'image/png';
var bounds = [628556.2744487218, 3452829.882138792,
716571.1992278139, 3618292.411380542];
//定义一个坐标系图层
var tiled = new ol.layer.Tile({
source: new ol.source.TileWMS({
url: 'http://localhost:8090/geoserver/cite/wms',
params: {
'FORMAT': format,
'VERSION': '1.1.1',
tiled: true,
STYLES: '',
LAYERS: 'cite:roads',
//tilesOrigin: 628556.2744487218 + "," + 3452829.882138792
}
})
});
var map = new ol.Map({
target : 'map',
layers : [
tiled
],
view : new ol.View({
projection: projection
})
});
map.getView().fit(bounds, map.getSize());