python - 在函数中实现脚本。一些建议？

Question

首先，我对 Python(一个编程领域)还很陌生，但我想学习和转换 jwpat7 开发的一个函数。 .给定一组从凸包导出的点

hull= [(560023.44957588764,6362057.3904932579), 
      (560023.44957588764,6362060.3904932579), 
      (560024.44957588764,6362063.3904932579), 
      (560026.94957588764,6362068.3904932579), 
      (560028.44957588764,6362069.8904932579), 
      (560034.94957588764,6362071.8904932579), 
      (560036.44957588764,6362071.8904932579), 
      (560037.44957588764,6362070.3904932579), 
      (560037.44957588764,6362064.8904932579), 
      (560036.44957588764,6362063.3904932579), 
      (560034.94957588764,6362061.3904932579), 
      (560026.94957588764,6362057.8904932579), 
      (560025.44957588764,6362057.3904932579), 
      (560023.44957588764,6362057.3904932579)]

此脚本返回此 post problem 之后所有可能区域的打印.代码开发者 jwpat7 是:

import math

def mostfar(j, n, s, c, mx, my): # advance j to extreme point
    xn, yn = hull[j][0], hull[j][1]
    rx, ry = xn*c - yn*s, xn*s + yn*c
    best = mx*rx + my*ry
    while True:
        x, y = rx, ry
        xn, yn = hull[(j+1)%n][0], hull[(j+1)%n][1]
        rx, ry = xn*c - yn*s, xn*s + yn*c
        if mx*rx + my*ry >= best:
            j = (j+1)%n
            best = mx*rx + my*ry
        else:
            return (x, y, j)

n = len(hull)
iL = iR = iP = 1                # indexes left, right, opposite
pi = 4*math.atan(1)
for i in range(n-1):
    dx = hull[i+1][0] - hull[i][0]
    dy = hull[i+1][1] - hull[i][1]
    theta = pi-math.atan2(dy, dx)
    s, c = math.sin(theta), math.cos(theta)
    yC = hull[i][0]*s + hull[i][1]*c    
    xP, yP, iP = mostfar(iP, n, s, c, 0, 1)
    if i==0: iR = iP
    xR, yR, iR = mostfar(iR, n, s, c,  1, 0)
    xL, yL, iL = mostfar(iL, n, s, c, -1, 0)
    area = (yP-yC)*(xR-xL) 
    print '    {:2d} {:2d} {:2d} {:2d} {:9.3f}'.format(i, iL, iP, iR, area)

结果是:

i iL iP iR    Area
 0  6  8  0   203.000
 1  6  8  0   211.875
 2  6  8  0   205.800
 3  6 10  0   206.250
 4  7 12  0   190.362
 5  8  0  1   203.000
 6 10  0  4   201.385
 7  0  1  6   203.000
 8  0  3  6   205.827
 9  0  3  6   205.640
10  0  4  7   187.451
11  0  4  7   189.750
12  1  6  8   203.000

我希望创建一个返回最小矩形的长度、宽度和面积的函数。例如:

Length, Width, Area = get_minimum_area_rectangle(hull)
print Length, Width, Area
18.036, 10.392, 187.451

我的问题是:

1. 我需要创建一个函数还是两个函数？例如:def mostfar 和 get_minimum_area_rectangle
2. hull 是一个值列表。这是最好的格式吗？
   1. 遵循单一功能方法，我有一个问题要集成到最里面

提前致谢

1)解决方案:一个函数 按照 Scott Hunter 建议的第一个解决方案，我在将 mostfar() 集成到 get_minimum_area_rectangle() 中时遇到了问题。非常感谢任何建议或帮助，因为我可以学习。

#!/usr/bin/python
import math

def get_minimum_area_rectangle(hull):
    # get pi greek
    pi = 4*math.atan(1)
    # number of points
    n = len(hull)
     # indexes left, right, opposite
    iL = iR = iP = 1
    # work clockwise direction
    for i in range(n-1):
        # distance on x axis
        dx = hull[i+1][0] - hull[i][0]
        # distance on y axis
        dy = hull[i+1][1] - hull[i][1]
        # get orientation angle of the edge
        theta = pi-math.atan2(dy, dx)
        s, c = math.sin(theta), math.cos(theta)
        yC = hull[i][0]*s + hull[i][1]*c

从上面的 jwpat7 示例开始，我需要使用 mostfar()。我很难理解在这一点上如何整合(抱歉用词不对)

Answer 1

1. 您可以使用单个函数或两个函数，但使用两个函数可能更简洁、更容易。您可以按原样保留 mostfar 函数。然后，通过添加函数定义行将代码的后半部分转换为函数:

   def get_minimum_area_rectangle(hull):
   

   …然后缩进其余代码(以 n = len(hull) 开头)以形成函数的主体。您还需要更改函数以返回您想要获取的值(长度、宽度和面积)。这将使您的代码保持模块化和干净，并且需要很少的更改。

2. 使用 hull 的值列表似乎可以达到此目的。另一种方法是使用数组(如 NumPy 数组)，但在这种情况下，您将迭代地处理数据，一次一个项目，而不是同时对多个数据点进行任何计算。所以一个列表应该没问题。访问列表中的项目很快，与您必须进行的数学运算相比，它不应该成为瓶颈。