排序坐标 [已关闭]

问：

你能帮我正确地对点（x 和 y 对谱）列表（数组）进行排序吗？ 我使用了以下功能，但它并不总是给我想要的。 示例是所附的图片。

我希望按以下顺序连接点：1 -> 2 -> 3 -> 4 -> 5。

但是，它们按拧紧顺序连接：1 -> 4 -> 2 -> 3 -> 5。

如何编辑代码以制作我想要的东西？

感谢您的帮助！

我想要这个特定的顺序，以便 voronoi 区域正确连接，如下图所示（蓝线）。我拥有的点是定义特定区域的 voronoi 图中的点。 简而言之，我将尝试在示例中解释我正在尝试做什么：

1. 假设我有 15 组点。在每组点中，可以有不同数量的点。
2. 15 个点中的第一组点是我感兴趣的集合，因为考虑到所有其他 14 组点，我想在第一组点周围获得沃罗诺伊区域。假设在第一组点中有 10 个点。
3. 因此，为了做到这一点，我合并了所有 15 组点的点，创建了一个 voronoi 图。
4. 从那个 voronoi 图中，我知道前 10 个点来自第一组点周围的区域，所以我可以得到它们并合并在一起。
5. “combined_coords”是我的第一组点周围的沃罗诺伊地区的点。
6. 现在我想对它们进行排序，以便以正确的顺序从它们定义多边形，就像在 voronoi 图上一样。
7. 我无法访问 shapely、geopandas 等库。只有 numpy 和 scipy。

combined_coords = np.array(
[[ 5766.28476155, 19516.99300136],
 [ 5765.3326747 , 19516.2217852 ],
 [ 5764.33255982, 19515.68693966],
 [ 5762.69492707, 19515.17292775],
 [ 5762.05045578, 19515.29517553],
 [ 5762.02505098, 19516.45055273],
 [ 5762.28121558, 19517.78513446],
 [ 5762.62958104, 19518.74867786],
 [ 5763.31238375, 19519.86920443],
 [ 5763.72620097, 19520.19250031],
 [ 5764.24999132, 19520.88744217],
 [ 5764.65547718, 19521.18117486],
 [ 5765.21676807, 19522.12631202],
 [ 5765.67501467, 19522.40253271],
 [ 5766.79878351, 19528.51536924],
 [ 5766.81357218, 19528.52972979],
 [ 5766.54317078, 19523.59611998],
 [ 5766.92157601, 19525.04371133],
 [ 5767.67457849, 19529.12581027],
 [ 5771.5263795 , 19524.43240175],
 [ 5770.59556598, 19522.45791971],
 [ 5770.45102759, 19521.34643935],
 [ 5769.91919814, 19521.21760643],
 [ 5768.52296931, 19520.39572883],
 [ 5767.63542881, 19519.1100208 ],
 [ 5767.57293753, 19519.06926975],
 [ 5766.99422734, 19518.06749945],
 [ 5766.76994398, 19517.86007797]]    
)
new_points = sorting_function(combined_coords)
def sorting_function(xy_e):
    #step 1
    xl=xy_e[:,0]
    yl=xy_e[:,1]
    xc=sum(xl)/len(xl)
    yc=sum(yl)/len(yl)

    xa=np.array(xl)-xc
    ya=np.array(yl)-yc

    #step 3
    q=[{},{},{},{}]
    for i in range(len(xa)):
        if (xa[i]<0 and ya[i]<0):
            q[0][i]=(xa[i],ya[i])
        elif (xa[i]<0 and ya[i]>0):
            q[1][i]=(xa[i],ya[i])
        elif (xa[i]>0 and ya[i]>0):
            q[2][i]=(xa[i],ya[i])
        else:
            q[3][i]=(xa[i],ya[i])

    #step 4        
    alpha={}
    for i in range(len(q)):
        for j in q[i].keys():
            beta=np.degrees(np.arctan(abs(q[i][j][1])/abs(q[i][j][0])))
            if i==0:
                s=90-beta
            elif i==1:
                s=90+beta
            elif i==2:
                s=270-beta
            else:
                s=270+beta
            alpha[j]=s

    #step 5  

    xy = xy_e.tolist()

    re_xy=[]
    sorted_alpha=sorted(alpha.values())
    for i in sorted_alpha:
        for k,l in alpha.items():
            if i==l:
                check_p=xy[k] in re_xy
                if not check_p:
                    re_xy.append(xy[k])
    return re_xy