在特殊情況下，坡度為無(wú)窮大/零，但基本方法是使用插值計(jì)算某個(gè)點(diǎn)的坡度，然后找到垂直坡度，然后計(jì)算該距離處的點(diǎn)。在

我修改了來(lái)自here的示例，以添加第二個(gè)圖。它與data file you provided一起工作，但您可能需要更改另一個(gè)信封的符號(hào)計(jì)算。在

編輯根據(jù)您關(guān)于希望信封連續(xù)的評(píng)論，我在結(jié)尾處添加了一個(gè)低俗的半圓，使非常接近。基本上，在創(chuàng)建封套時(shí)，越圓越凸，效果越好。另外，你需要把開(kāi)頭和結(jié)尾重疊，否則你會(huì)有一個(gè)間隙。在

而且，幾乎可以肯定的是，它可以變得更有效率我不是一個(gè)numpy專(zhuān)家，所以這只是純Python。在def offset(coordinates, distance):

coordinates = iter(coordinates)

x1, y1 = coordinates.next()

z = distance

points = []

for x2, y2 in coordinates:

# tangential slope approximation

try:

slope = (y2 - y1) / (x2 - x1)

# perpendicular slope

pslope = -1/slope # (might be 1/slope depending on direction of travel)

except ZeroDivisionError:

continue

mid_x = (x1 + x2) / 2

mid_y = (y1 + y2) / 2

sign = ((pslope > 0) == (x1 > x2)) * 2 - 1

# if z is the distance to your parallel curve,

# then your delta-x and delta-y calculations are:

# z**2 = x**2 + y**2

# y = pslope * x

# z**2 = x**2 + (pslope * x)**2

# z**2 = x**2 + pslope**2 * x**2

# z**2 = (1 + pslope**2) * x**2

# z**2 / (1 + pslope**2) = x**2

# z / (1 + pslope**2)**0.5 = x

delta_x = sign * z / ((1 + pslope**2)**0.5)

delta_y = pslope * delta_x

points.append((mid_x + delta_x, mid_y + delta_y))

x1, y1 = x2, y2

return points

def add_semicircle(x_origin, y_origin, radius, num_x = 50):

points = []

for index in range(num_x):

x = radius * index / num_x

y = (radius ** 2 - x ** 2) ** 0.5

points.append((x, -y))

points += [(x, -y) for x, y in reversed(points)]

return [(x + x_origin, y + y_origin) for x, y in points]

def round_data(data):

# Add infinitesimal rounding of the envelope

assert data[-1] == data[0]

x0, y0 = data[0]

x1, y1 = data[1]

xe, ye = data[-2]

x = x0 - (x0 - x1) * .01

y = y0 - (y0 - y1) * .01

yn = (x - xe) / (x0 - xe) * (y0 - ye) + ye

data[0] = x, y

data[-1] = x, yn

data.extend(add_semicircle(x, (y + yn) / 2, abs((y - yn) / 2)))

del data[-18:]

from pylab import *

with open('ah79100c.dat', 'rb') as f:

f.next()

data = [[float(x) for x in line.split()] for line in f if line.strip()]

t = [x[0] for x in data]

s = [x[1] for x in data]

round_data(data)

parallel = offset(data, 0.1)

t2 = [x[0] for x in parallel]

s2 = [x[1] for x in parallel]

plot(t, s, 'g', t2, s2, 'b', lw=1)

title('Wing with envelope')

grid(True)

axes().set_aspect('equal', 'datalim')

savefig("test.png")

show()

總結(jié)

以上是生活随笔為你收集整理的python求点到曲线距离_Python。如何从点和偏移距离的x，y列表中获取偏移样条曲线的x，y坐标...的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。