<p>我也在为一个类似的问题而挣扎。
我还想画一个菲涅耳区,但我想把它画在视线范围内,这条线连接着a点和B点</p>
<p>使用ewcz提供的代码,我添加了这条线并绘制了所有内容。在</p>
<p>结果表明,旋转直线与椭圆轴线不对应,因而与视线不一致。在</p>
<pre><code>#!/usr/bin/env python
import math
from shapely.geometry import Point, LineString
from shapely.affinity import scale, rotate
from matplotlib import pyplot as plt
#input parameters
A = Point(0, 0)
B = Point(400, 10)
R = 5
d = A.distance(B)
#first, rotate B to B' around A so that |AB'| = |AB| and B'.y = A.y
#and then take S as midpoint of AB'
S = Point(A.x + d/2, A.y)
#Make a straight line
LOS = LineString([(A.x, A.y), (B.x, A.y)])
#alpha represents the angle of this rotation
alpha = math.atan2(B.y - A.y, B.x - A.x)
#create a circle with center at S passing through A and B'
C = S.buffer(d/2)
#rescale this circle in y-direction so that the corresponding
#axis is R units long
C = scale(C, 1, R/(d/2))
#rotate the ellipse obtained in previous step around A into the
#original position (positive angles represent counter-clockwise rotation)
C = rotate(C, alpha, origin=A, use_radians=True)
f_x, f_y = C.exterior.xy
#plot the ellipse
plt.plot(f_x, f_y)
#rotate the line in the same way as the ellipse
LOS_R = rotate(LOS, alpha, origin=A, use_radians=True)
f_x, f_y = LOS_R.xy
#plot the line
plt.plot(f_x, f_y)
plt.show()
</code></pre>
<p>结果图为:
<a href="https://i.stack.imgur.com/lAXaZ.png" rel="nofollow noreferrer">Plotted image with matplot</a></p>