求点集的最小面积外接矩形

刀么克瑟拉莫

于 2019-12-17 09:59:08 发布

阅读量1.6k

点赞数

CC 4.0 BY-SA版权

分类专栏：小功能

本文链接：https://blog.youkuaiyun.com/random_repick/article/details/103574844

小功能专栏收录该内容

22 篇文章

订阅专栏

本文介绍了一种求解点集最小面积外接矩形的算法，通过旋转坐标系并遍历所有可能的角度来寻找最紧凑的包围矩形。算法首先定义了一个旋转函数，用于将点绕原点旋转特定角度，然后通过遍历0到90度的所有角度，计算每个角度下包围所有点的矩形面积，并保留最小面积及其对应的矩形尺寸。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

求点集的最小面积外接矩形

void Rotate(float& x, float& y, float angle) {
  float a = x * cos(angle) - y * sin(angle);
  float b = x * sin(angle) + y * cos(angle);
  x = a;
  y = b;
}

void Rect(const vector<float>& x,
          const vector<float>& y,
          float& length,
          float& width) {
  float area = 1024;
  int size = x.size();
  // rotate 0~90 degree
  for (int i = 0; i < 91; ++i) {
    float tmpx = x[0], tmpy = y[0];
    Rotate(tmpx, tmpy, i);
    float xl = tmpx, xr = tmpx, yt = tmpy, yb = tmpy;
    // traverse all points
    for (int j = 1; j < size; ++j) {
      tmpx = x[j];
      tmpy = y[j];
      Rotate(tmpx, tmpy, i);
      if (tmpx < xl)
        xl = tmpx;
      if (tmpx > xr)
        xr = tmpx;
      if (tmpy < yb)
        yb = tmpy;
      if (tmpy > yt)
        yt = tmpy;
    }
    float xx = xr - xl, yy = yt - yb;
    if (area > xx * yy) {
      area = xx * yy;
      length = xx;
      width = yy;
    }
  }
  //
  if (length < width) {
    float tmp = length;
    length = width;
    width = tmp;
  }
}