aze*_*q6d 8 python geometry center radius
我写了这个函数:
def get_x_y_co(circles):
xc = circles[0] #x-co of circle (center)
yc = circles[1] #y-co of circle (center)
r = circles[2] #radius of circle
arr=[]
for i in range(360):
y = yc + r*math.cos(i)
x = xc+ r*math.cos(i)
x=int(x)
y=int(y)
#Create array with all the x-co and y-co of the circle
arr.append([x,y])
return arr
Run Code Online (Sandbox Code Playgroud)
'circles' 是一个数组 [X-center, Y-center, Radius]
我想提取圆圈中存在的整数分辨率的所有点。现在,我意识到我正在创建一组位于圆边界上的点,但我无法访问圆内的点。
我想只减小半径,然后对半径的所有值进行迭代,直到半径为 0
但我觉得有一种更有效的方法。欢迎任何帮助
Maa*_*bré 10
from itertools import product
def points_in_circle(radius):
for x, y in product(range(int(radius) + 1), repeat=2):
if x**2 + y**2 <= radius**2:
yield from set(((x, y), (x, -y), (-x, y), (-x, -y),))
list(points_in_circle(2))
Run Code Online (Sandbox Code Playgroud)
Run Code Online (Sandbox Code Playgroud)[(0, 0), (0, 1), (0, -1), (0, -2), (0, 2), (1, 0), (-1, 0), (-1, 1), (1, -1), (1, 1), (-1, -1), (2, 0), (-2, 0)]
def points_in_circle_np(radius):
a = np.arange(radius + 1)
for x, y in zip(*np.where(a[:,np.newaxis]**2 + a**2 <= radius**2)):
yield from set(((x, y), (x, -y), (-x, y), (-x, -y),))
Run Code Online (Sandbox Code Playgroud)
def points_in_circle_np(radius, x0=0, y0=0, ):
x_ = np.arange(x0 - radius - 1, x0 + radius + 1, dtype=int)
y_ = np.arange(y0 - radius - 1, y0 + radius + 1, dtype=int)
x, y = np.where((x_[:,np.newaxis] - x0)**2 + (y_ - y0)**2 <= radius**2)
# x, y = np.where((np.hypot((x_-x0)[:,np.newaxis], y_-y0)<= radius)) # alternative implementation
for x, y in zip(x_[x], y_[y]):
yield x, y
Run Code Online (Sandbox Code Playgroud)
| 归档时间: |
|
| 查看次数: |
12117 次 |
| 最近记录: |