以锯齿形顺序处理完整圆的所有角度的函数

如何解决以锯齿形顺序处理完整圆的所有角度的函数

我需要以度为单位计算很多角度的值。为了快速建立粗略的形状，以及稍后的精细位，我想按这个顺序计算形状（0°、180°、90°、270°、45°、135°...）

下面的代码做我想要的。我想知道：有没有办法以更直接的方式做到这一点？它需要适用于任何（整数）数字（例如 72 或 7465）

感谢您的帮助。

import numpy as np

def evenly_spaced_star_order(number):
    
    Total=np.linspace(0,360,number,endpoint=False)

    Res=[]
    
    for devider in [2**_ for _ in range(1000)]:
        for counter in range(devider):
            Number=(counter*len(Total))//devider
            if np.isfinite(Total[Number]):
                Res.append(Total[Number])
                Total[Number]=np.nan
        if np.all(np.isnan(Total)):
            break
    return(Res)
print(evenly_spaced_star_order(16))

解决方法

我的解决方案递归地分离偶数和奇数索引。然后将奇数行放在最终列表的末尾（按顺序），偶数行再次递归拆分。

我的顺序与您的原始函数一致，但速度要快得多（一个数量级或更多），并且确实适用于任何整数。

# recursive evenly_spaced_star_order()
def esso(number):
    def interleave(arr):
        return arr if len(arr) <= 1 else np.append(interleave(arr[0::2]),arr[1::2])
    
    return interleave(np.linspace(0,360,number,endpoint=False))

print(esso(16))

我的时间：

%timeit evenly_spaced_star_order(16)
885 µs ± 8.68 µs per loop (mean ± std. dev. of 7 runs,1000 loops each)

%timeit esso(16)
60.1 µs ± 998 ns per loop (mean ± std. dev. of 7 runs,10000 loops each)

%timeit evenly_spaced_star_order(1000)
5.88 ms ± 192 µs per loop (mean ± std. dev. of 7 runs,100 loops each)

%timeit esso(1000)
111 µs ± 10.9 µs per loop (mean ± std. dev. of 7 runs,10000 loops each)

随着点数的增加，我的性能会越来越好（与原始代码相比）。

第二种解决方案

它几乎没有那么漂亮，但顺序更接近，速度仍然更快。

def esso2(number):
    def interleave(arr):
        if arr.shape[0] <= 1:
            return arr
        mid = arr.shape[0] // 2
        it1 = iter(interleave(arr[0:mid]))
        it2 = iter(interleave(arr[mid:]))
        return sum(zip(it1,it2),()) + tuple(it2)
    return np.array(interleave(np.linspace(0,endpoint=False)))
print(esso2(72))