如何通过两点（直径端）画一个圆？答案

【问题标题】：How to draw a circle through two points (at diameter end)?如何通过两点（直径端）画一个圆？
【发布时间】：2017-04-26 18:58:13
【问题描述】：

如图所示在立体投影上有两个点：

这些点应该在一个直径的圆的端点上。怎么画一个经过这两点的圆？

上图的代码：

import matplotlib.pylab as plt
from mpl_toolkits.basemap import Basemap
import numpy as np
from scipy.interpolate import splev, splrep

# create instance of basemap, note we want a south polar projection to 90 = E
myMap = Basemap(projection='spstere',boundinglat=0,lon_0=180,resolution='l',round=True,suppress_ticks=True)
# set the grid up
gridX, gridY = 10.0, 15.0
parallelGrid = np.arange(-90.0,90.0,gridX)
meridianGrid = np.arange(-180.0,180.0,gridY)

# draw parallel and meridian grid, not labels are off. We have to manually create these.
myMap.drawparallels(parallelGrid,labels=[False,False,False,False])
myMap.drawmeridians(meridianGrid,labels=[False,False,False,False],labelstyle='+/-',fmt='%i')


# plot azimuth labels, with a North label.
ax = plt.gca()
ax.text(0.5,1.025,'N',transform=ax.transAxes,horizontalalignment='center',verticalalignment='bottom',size=25)
for para in np.arange(gridY,360,gridY):
    x= (1.1*0.5*np.sin(np.deg2rad(para)))+0.5
    y= (1.1*0.5*np.cos(np.deg2rad(para)))+0.5
    ax.text(x,y,u'%i\N{DEGREE SIGN}'%para,transform=ax.transAxes,horizontalalignment='center',verticalalignment='center')

summerAzi = np.array([0, 360])
summerAlt = np.array([40, 4])
summerX, summerY = myMap(summerAzi, -summerAlt)


summerX_new = np.linspace(summerX.min(), summerX.max(),30)
summerY_smooth = splev(summerX_new, splrep(summerX, summerY, k=1))

myMap.plot(summerX_new, summerY_smooth, 'g')

myMap.plot(summerX, summerY, 'go')   
plt.show()

【问题讨论】：

两个点可以定义无限多个圆，你说的是投影空间中的圆（如我们所见）还是数学空间中的圆（所表示的数据）
@NickT;是的你是对的。我更新了问题。这些点是圆直径的端点。
因为有投影，计算线的中点是不行的，用circ=plt.Circle((myX,myY), radius=myRad)
@beroe;不，它不会起作用，但如果你有任何东西，那么你可以分享。

标签： python matplotlib plot matplotlib-basemap

【解决方案1】：

内置的 tissot() 函数足以在等角投影上绘制圆圈（如本例所示）。在非等角投影上，它绘制椭圆。这里 tissot 指示线的中点是 (0, -22) 度数。它的半径 = (40-4)/2 = 18 度。点数 = 36 就可以了。

相关代码为：

myMap.tissot(0, -22, 18, 36, \
             facecolor='none', \
             edgecolor='#ff0000', \
             linewidth=1, \
             alpha=1)

【讨论】：

【解决方案2】：

此极坐标表示中的圆看起来不像矩形网格上的圆（即“圆形”）。除此之外，您可以像在笛卡尔平面上一样画一个圆，从极坐标开始，转换为笛卡尔坐标，偏移中心并使用绘图功能。

summerAzi = np.array([0, 360])
summerAlt = -np.array([40, 4])
summerX, summerY = myMap(summerAzi, summerAlt)

phi = np.linspace(0,2.*np.pi)
r = np.abs(np.diff(summerAlt))/2.
x = r*np.cos(phi)
y = -r*np.sin(phi)+summerAlt.mean()
X,Y= myMap(x,y)

myMap.plot(X,Y, color="crimson")
myMap.plot(summerX, summerY, color="gold", marker="o")

【讨论】：

有什么方法可以绘制一个圆圈，然后可以将其放置在这个极坐标上？实际情况有所不同。数据集类似于summerAzi = np.array([180, 360])。
抱歉，没有明确的问题描述，我无能为力。
我们不能这样做吗：#%matplotlib inline import matplotlib.pylab as plt import numpy as np ax = plt.subplot(111, polar=True) circle = plt.Circle((0, .3), 0.6, transform=ax.transData._b, color="red", alpha=0.1) ax.add_artist(circle) plt.show() 在问题中使用summerAzi 和summerAlt？
我试图用这个code 画圆圈，但没有成功。请问是否有机会调整圆心可能会产生一些结果？