给定一个仅包含0和1、大小为 rows * cols的二维numpy数组，使用python numpy找出只包含1的最大矩形，并返回其面积

要找出只包含1的最大矩形并返回其面积，可以使用以下方法：
1. 首先将二维numpy数组转换为一个矩阵，其中每个元素表示从该位置向上的连续1的数量。
2. 对于每一行，使用栈来找到最大矩形。将每个元素的值与栈顶元素的值进行比较，如果当前元素的值小于栈顶元素的值，则弹出栈顶元素并计算以该元素为高度的最大矩形面积。如果当前元素的值大于或等于栈顶元素的值，则将其推入栈中。
3. 对于每一行，计算以该行为底部的最大矩形面积，并更新最大面积。
4. 返回最大面积。

下面是代码实现：

def maximalRectangle(matrix):
    if not matrix:
        return 0
    m, n = len(matrix), len(matrix[0])
    heights = [0] * n
    max_area = 0
    for i in range(m):
        for j in range(n):
            if matrix[i][j] == 0:
                heights[j] = 0
            else:
                heights[j] += 1
        stack = [-1]
        for k in range(n):
            while stack[-1] != -1 and heights[k] < heights[stack[-1]]:
                h = heights[stack.pop()]
                w = k - stack[-1] - 1
                max_area = max(max_area, h * w)
            stack.append(k)
        while stack[-1] != -1:
            h = heights[stack.pop()]
            w = n - stack[-1] - 1
            max_area = max(max_area, h * w)
    return max_area