如何在给定矩阵中找到最大区域并返回其位置/边界

def find_largest_region(matrix):
    max_region_size = 0
    max_region_coords = [(0, 0), (0, 0)]

    def dfs(i, j, curr_region_size, coords):
        if i < 0 or i >= len(matrix) or j < 0 or j >= len(matrix[0]) or matrix[i][j] != 1:
            return
        matrix[i][j] = -1  # 标记已访问过的元素
        curr_region_size += 1
        coords[1] = (i, j)  # 更新最大区域的右下角坐标

        dfs(i - 1, j, curr_region_size, coords)  # 上
        dfs(i + 1, j, curr_region_size, coords)  # 下
        dfs(i, j - 1, curr_region_size, coords)  # 左
        dfs(i, j + 1, curr_region_size, coords)  # 右

    for i in range(len(matrix)):
        for j in range(len(matrix[0])):
            if matrix[i][j] == 1:
                coords = [(i, j), (i, j)]
                dfs(i, j, 0, coords)
                if coords[1][0] - coords[0][0] + 1 > max_region_size:  # 比较当前区域和最大区域的大小
                    max_region_size = coords[1][0] - coords[0][0] + 1
                    max_region_coords = coords

    return max_region_coords

# 示例用法
matrix = [
    [1, 1, 0, 0, 0],
    [1, 1, 1, 1, 0],
    [0, 0, 1, 1, 0],
    [0, 0, 0, 0, 0]
]
result = find_largest_region(matrix)
print("最大区域的位置/边界:", result)