项目开发过程中有绘制历史轨迹的功能，上面说直线绘制效果不是太好，想看看曲线绘制的效果。
由于leaflet没有提供绘制曲线的方法，只能自己算出曲线的座标点以直线相连，绘制成曲线
算法部分https://blog.csdn.net/TarkuNi/article/details/99054036
<!DOCTYPE html>
<html>
<head>
<meta charset=utf-8 />
<meta name='viewport' content='width=device-width, initial-scale=1.0, maximum-scale=1.0, user-scalable=no, minimal-ui' />
<title>Leaflet Side-by-side</title>
<script src="https://www.jq22.com/jquery/jquery-3.3.1.js"></script>
<script src='http://cdn.leafletjs.com/leaflet/v1.3.1/leaflet.js'></script>
<link href='http://cdn.leafletjs.com/leaflet/v1.3.1/leaflet.css' rel='stylesheet' />
<style>
body {
margin: 0;
padding: 0;
}
#map {
position: absolute;
top: 0px;
bottom: 0;
width: 100%;
}
</style>
</head>
<body>
<!-- <input type="checkbox" name="async" checked id="" onclick="checkClick()">同步 -->
<div id='map'></div>
<script>
var map = L.map('map').setView([51.505, -0.09], 13);
L.tileLayer('https://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png', {
attribution: '© <a href="https://www.openstreetmap.org/copyright">OpenStreetMap</a> contributors'
}).addTo(map);
L.marker([51.5, -0.09]).addTo(map)
.bindPopup('A pretty CSS3 popup.<br> Easily customizable.')
.openPopup();
calcCoorArr=(point_start,point_end,num)=>{// 入参两个座标（字符串格式："114.30911,30.600052"）,num是返回座标点个数
//第一步把座标字符串转为对象,为了方便计算转为了数字格式
var p_start={x:parseFloat(point_start.split(",")[0]),y:parseFloat(point_start.split(",")[1])};
var p_end={x:parseFloat(point_end.split(",")[0]),y:parseFloat(point_end.split(",")[1])};
//此处敲黑板，是任务的第二大难点，求出相对的第三个点，用于固定曲线的弯曲度，下面公式是已知三角形两点座标和两个座标点的夹角求第三点座标，两个夹角我们是自定义任意值，不过不要加起来超过180度
// 已知两点p1(x1,y1)、p2(x2,y2)和两点所对应的角度A和B,x3、y3是对应第三点的座标，cot25°=2.14
//x3 = (x1*cotB+x2*cotA+y2-y1)/(cotA+cotB)
//y3 = (y1*cotB+y2*cotA+x1-x2)/(cotA+cotB)
let x3=(p_start.x*2.14+p_end.x*2.14-p_start.y+p_end.y)/(2*2.14)
let y3=(p_start.y*2.14+p_end.y*2.14-p_end.x+p_start.x)/(2*2.14)
var p_crt1={x:x3,y:y3};
var p_crt2={x:x3,y:y3};
//下面计算贝叶斯曲线，不是几个字能说清，直接拿去用没毛病
/**
* 计算公式：
* | 1 0 0 0| |P0|
* [1 t t*t t*t*t] |-3 3 0 0| |P1|
* |3 -6 3 0| |P2|
* |-1 3 -3 1| |p3|
*
* **/
let paths=[];
for(let i=0;i<num+1;i++){
let t=i/num;
var _matrix1=[1,t,t*t,t*t*t];
var _matrix2=[
[1,0,0,0]
,[-3,3,0,0]
,[3,-6,3,0]
,[-1,3,-3,1]
];
var _matrix3=[
[p_start.x,p_start.y]
,[p_crt1.x,p_crt1.y]
,[p_crt2.x,p_crt2.y]
,[p_end.x,p_end.y]
];
var _matrix_tmp=[
_matrix1[0]*_matrix2[0][0]+_matrix1[1]*_matrix2[1][0]+_matrix1[2]*_matrix2[2][0]+_matrix1[3]*_matrix2[3][0]
,_matrix1[0]*_matrix2[0][1]+_matrix1[1]*_matrix2[1][1]+_matrix1[2]*_matrix2[2][1]+_matrix1[3]*_matrix2[3][1]
,_matrix1[0]*_matrix2[0][2]+_matrix1[1]*_matrix2[1][2]+_matrix1[2]*_matrix2[2][2]+_matrix1[3]*_matrix2[3][2]
,_matrix1[0]*_matrix2[0][3]+_matrix1[1]*_matrix2[1][3]+_matrix1[2]*_matrix2[2][3]+_matrix1[3]*_matrix2[3][3]
];
var _matrix_final=[
_matrix_tmp[0]*_matrix3[0][0]+_matrix_tmp[1]*_matrix3[1][0]+_matrix_tmp[2]*_matrix3[2][0]+_matrix_tmp[3]*_matrix3[3][0]
,_matrix_tmp[0]*_matrix3[0][1]+_matrix_tmp[1]*_matrix3[1][1]+_matrix_tmp[2]*_matrix3[2][1]+_matrix_tmp[3]*_matrix3[3][1]
];
//下面注释掉的原因是入参是经纬度，但leaflet渲染需要的是纬度在前经度在后的数组，然后你懂的
var _res_point=[
_matrix_final[0]
,_matrix_final[1]
];
// var _res_point=[_matrix_final[1],_matrix_final[0]];
paths.push(_res_point);
}
return paths;
}
var path = calcCoorArr('51.5,-0.09','51.48,-0.041',20)
console.log(path)
L.polyline(path,{color:'blue'}).addTo(map).bindPopup("green to red");
L.marker([51.48, -0.041]).addTo(map)
var path2 = calcCoorArr('51.48,-0.041','51.48,-0.087',20)
L.polyline(path2,{color:'blue'}).addTo(map).bindPopup("green to red");
L.marker([51.48, -0.087]).addTo(map)
var path3 = calcCoorArr('51.48,-0.087','51.49,-0.014',20)
L.polyline(path3,{color:'blue'}).addTo(map).bindPopup("green to red");
L.marker([51.49, -0.014]).addTo(map)
map.on('click',function(e){
// 用于调试 获取点击的座标和地图级别
console.log(e.latlng)
console.log(map.getZoom())
})
</script>
</body>
</html>