如何通过OpenCV轮廓计算长宽比，判断谷歌地球建筑为正方形/矩形？

嘿，我来帮你搞定这个问题！其实OpenCV里有几个实用的函数，刚好能解决你算长宽比的需求，进而判断建筑是正方形还是矩形，我给你详细拆解一下：

这个方法适合轮廓没有明显旋转的情况，它会生成一个平行于坐标轴的外接矩形：

• 调用函数后会返回(x, y, width, height)四个值，分别是矩形的左上角坐标和长宽。
• 直接计算长宽比：aspect_ratio = float(width) / height（一定要转成浮点型，避免整数除法导致的精度丢失）。
• 判断逻辑：设置一个误差阈值（比如5%），如果abs(aspect_ratio - 1) < 0.05，就可以判定为正方形，否则是矩形。
• 代码示例：

import cv2

# 假设你已经提取到了目标轮廓contour
x, y, w, h = cv2.boundingRect(contour)
aspect_ratio = float(w) / h

if abs(aspect_ratio - 1) < 0.05:
    print("该建筑是正方形")
else:
    print("该建筑是矩形")

• 注意：如果谷歌地球的航拍角度导致建筑轮廓倾斜，这个方法得到的矩形会比实际建筑大，长宽比可能不准，这时候用下面的方法更合适。

这个方法会生成刚好包裹轮廓的最小矩形，哪怕轮廓是旋转的，也能贴合实际形状，精度更高：

• 调用函数后返回元组((center_x, center_y), (width, height), angle)，其中width和height是矩形的两个边长（注意这里两个值的顺序可能会因为旋转颠倒）。
• 处理步骤：
  1. 先提取边长，然后统一顺序（确保长 >= 宽，避免旋转导致的计算混乱）。
  2. 计算长宽比，再用阈值判断形状。
• 代码示例：

import cv2

rect = cv2.minAreaRect(contour)
width, height = rect[1]

# 统一边长顺序，保证width >= height
if width < height:
    width, height = height, width

aspect_ratio = width / height

if abs(aspect_ratio - 1) < 0.05:
    print("该建筑是正方形")
else:
    print("该建筑是矩形")

你已经有轮廓的面积和周长了，可以用它们来辅助验证结果，避免轮廓噪声导致误判：

• 对于正方形，理论面积应该接近(周长/4)^2；对于矩形，理论面积接近width*height（从上面方法得到的长宽）。
• 计算误差率：area_error = abs(实际轮廓面积 - 理论面积) / 实际轮廓面积，如果误差超过10%，说明轮廓可能有噪声，建议先做平滑处理（比如用cv2.GaussianBlur()模糊，或者cv2.morphologyEx()做开操作去除小噪点）。

内容的提问来源于stack exchange，提问作者eljirg