我正在寻找一种简短的智能方法来查找线段上的所有整数点。 2个点也是整数,线可以是0,45,90,135等角度。
这是我的长代码(到目前为止90度案例):
def getPoints(p1,p2)
if p1[0] == p2[0]:
if p1[1] < p2[1]:
return [(p1[0],x) for x in range(p1[1],p2[1])]
else:
return [(p1[0],x) for x in range(p1[1],p2[1],-1)]
if p2[1] == p2[1]:
if p1[0] < p2[0]:
return [(x,p1[1]) for x in range(p1[0],p2[0])]
else:
return [(x,p1[1]) for x in range(p1[0],p2[0],-1)]
编辑:我没有提到它,但斜率将始终为整数-1,0或1,有8种情况需要检查。
答案 0 :(得分:1)
将斜率减小到最低项(p / q),然后从线段的一个端点到另一个端点以p垂直和q水平增量步进。如果缩减到最低项的代码将5/0减少到1/0,则相同的代码可以用于垂直线段。
答案 1 :(得分:0)
对每对点做一点数学计算m&amp; c表示mx + c,并将其与您正在考虑的行的公式进行比较。 (N.B.你 将 通过零来区分以应对。)
答案 2 :(得分:0)
我可以编写有效的代码,但重复代码的数量是 把我扔了,这就是我转向你们的原因
这可能会有很大的改进,但也许它会让你走上正轨。 (抱歉,没有时间让它变得更好!)
def points(p1,p2):
slope = (p2[1]-p1[1])/float(p2[0]-p1[0])
[(x,x*slope) for x in range (p1[0], p2[0]) if int(x*slope) == x*slope)]
答案 3 :(得分:0)
扩展@Jon Kiparsky的答案。
def points_on_line(p1, p2):
fx, fy = p1
sx, sy = p2
if fx == sx and fy == sy:
return []
elif fx == sx:
return [(fx, y) for y in range(fy+1, sy)]
elif fy == sy:
return [(x, fy) for x in range(fx+1, sx)]
elif fx > sx and fy > sy:
p1, p2 = p2, p1
slope = (p2[1] - p1[1]) / float(p2[0] - p1[0])
return [(x, int(x*slope)) for x in range(p1[0], p2[0]) if int(x*slope) == x*slope and (x, int(x*slope)) != p1]
答案 4 :(得分:-1)
def getpoints(p1, p2):
# Sort both points first.
(x1, y1), (x2, y2) = sorted([p1, p2])
a = b = 0.0
# Not interesting case.
if x1 == x2:
yield p1
# First point is in (0, y).
if x1 == 0.0:
b = y1
a = (y2 - y1) / x2
elif x2 == 0.0:
# Second point is in (0, y).
b = y2
a = (y1 - y2) / x1
else:
# Both points are valid.
b = (y2 - (y1 * x2) / x1) / (1 - (x2 / x1))
a = (y1 - b) / x1
for x in xrange(int(x1), int(x2) + 1):
y = a * float(x) + b
# Delta could be increased for lower precision.
if abs(y - round(y)) == 0:
yield (x, y)