查找Lat\Lon数组中的最大斜率

def find_max_slope(points):
    max_slope = float('-inf')
    point1_index = -1
    point2_index = -1
    
    for i in range(len(points)):
        for j in range(i+1, len(points)):
            lat1, lon1 = points[i]
            lat2, lon2 = points[j]
            
            if lat2 - lat1 != 0:
                slope = (lon2 - lon1) / (lat2 - lat1)
                
                if slope > max_slope:
                    max_slope = slope
                    point1_index = i
                    point2_index = j
    
    return point1_index, point2_index

# 示例数据
points = [(39.9042, 116.4074), (31.2304, 121.4737), (23.1291, 113.2644)]
point1_index, point2_index = find_max_slope(points)

print("最大斜率对应的两个点的索引：", point1_index, point2_index)
print("最大斜率：", (points[point2_index][1] - points[point1_index][1]) / (points[point2_index][0] - points[point1_index][0]))