MATLAB寻找最小矩形区域

对于一个倾斜的矩形或不规则单连通域，找出其最小矩形，并将其旋转到X,Y坐标系平行位置。如：将 下图蓝色矩形区范围内的黑色区域旋转至，蓝色矩形的边平行于X，Y坐标系。

代码1

图像以固定的步幅进行旋转，每旋转一次，获取X,Y轴上的投影长度，相乘得出投影矩形面积。找出旋转过程中最小的投影面积，此时就是最小矩形。
S=（X2-X1）*(Y2-Y1)
方案1：直接对图像进行旋转，并不断的计算投影面积。

%% 代码作用：找出图像的最小矩形，并将图像旋转到最小矩形在x,y轴方向
% 输入图像仅包含一个单连通域的图形，否则可能出错；
% 思路：
%       1.将图像在90°内旋转，找出在X,Y坐标上投影面积最小时的角度即可；
%       2.先找出最小矩形，然后判断最小矩形和X,Y轴的夹角，然后旋转一下即可。
% 下面代码采用第一种方法进行。
% 步骤：
%       1.二值化图像；2.边旋转边检查投影面积；
% 作者：cw      时间：2019.8.12
clc;
clear all;
t0 = clock;
%% 图像读入，灰度化，二值化，显示二值图
img=imread('D.jpg');
img = rgb2gray(img);  
BW=~im2bw(img);
figure;imshow(BW);

%% 计算间隔1°旋转过程中，图像的X,Y投影面积
Area =