[英]Finding equally spaced points given a set of points with line segments in between
我目前正在编程 map,并试图将街道段分成相等的部分。 如果街道段是直的,只需将街道段的长度除以您想要的任何因素即可。 然而,弯曲的街道更难分成相等的部分,因为它们由多个部分组成。 我想要做的是找出一种方法将街道段分成相等的点,不管它有多弯曲,或者每段有多长。 我尝试参数化曲线以使其工作,但它仍然不起作用。 要了解我所说的参数化是什么意思,请查看这个。
目前,这就是我目前正在实施的方式。
//street_seg_length is length of the ith street segment
for (double j = 0.0; j < street_seg_length[i]; j+=inc) {
double rem = j - int(j);
//A street segment can be split up into multiple curve points.
//The more curve points, the curvier a road is
LatLon from = curve_points_pos[i][int(j)];
LatLon to = curve_points_pos[i][int(j)+1];
//Returns cartesian coordinates from latitude and longitude of nth and n+1th curve point
double x1 = x_from_lon(from.longitude());
double y1 = y_from_lat(from.latitude());
double x2 = x_from_lon(to.longitude());
double y2 = y_from_lat(to.latitude());
//arrowX, arrowY are supposed to be the coordinates of a point between 2 curve points
double arrowX = (1 - rem)*x1 + rem*x2;
double arrowY = (1 - rem) * y1 + rem * y2;
}
}
rem 是上篇文章中讨论的余数,j 是与上篇文章中的 t 相同的值,p 是第 n 个点,q 是上篇文章中讨论的第 n+1 个点。
有人可以解释我能做些什么来实现这一目标或我做错了什么吗? 我想在街段上找到等距的点,不管它有多弯曲。 我是否正确遵循链接帖子中的算法? 我相信我是,但很明显算法是错误的,或者我没有正确实现它,后者更有可能。
您生成的路径+
将具有等距点,但您的原始路径o
具有不同长度的段。 你有两个不同的事情要处理:
curr
和next
以及t
。 这是curr / next
值的示例:
o + 0 / 1 t = 0.0
| |
| + 0 / 1 t = 0.667
o |
| + 1 / 2 t = 0.2
| |
| + 1 / 2 t = 0.6
| |
o + 1 / 2 t = 1.0
所以你的算法是这样的:
curr = 0
和next = 1
开始。z
curr
和next
,使curr ≤ z ≤ next
。t
并进行插值这是一个将路径seg
拆分为n
个相同长度的段的实现。 (它不在 C++ 中,而是在 Javascript 中,它直接处理(x,y)坐标而不是 lon/lat,但我想你可以理解它。)
function split(seg, n) {
let acc = []; // n + 1 accumulated lengths
let len = 0; // overall length
let p = seg[0];
let res = []; // array of n + 1 result points
// find segemnt and overall lengths
for (let i = 0; i < seg.length; i++) {
let q = seg[i];
len += Math.hypot(q.x - p.x, q.y - p.y);
acc.push(len);
p = q;
}
acc.push(2 * len); // sentinel
let curr = 0;
let next = 1;
// create equidistant result points
for (let i = 0; i < n; i++) {
let z = len * i / n; // running length of point i
// advance to current segment
while (z > acc[next]) {
curr++;
next++;
}
// interpolate in segment
let p = seg[curr];
let q = seg[next];
let t = (z - acc[curr]) / (acc[next] - acc[curr]);
res.push(new Point(p.x * (1 - t) + q.x * t,
p.y * (1 - t) + q.y * t));
}
// push end point (leave out when joining consecutive segments.)
res.push(seg[seg.length - 1]);
return res;
}
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.