Matlab描绘堆积角图像边界并拟合函数求休止角

休止角（堆积角）

• 在工程应用中，堆积角（休止角）是一个经常听到的名词，目前测量堆积角普遍使用到堆积角测量仪，通过读取堆积物料高度和堆积物料底基直径对堆积角进行估算，缺点在于最后得到的数据并不是十分的准确，特别是针对于黏性较大的物料。
• 因此我们可以采用对堆积物料进行拍照，读取边界的形式进行测量，使测量结果更加精确。目前此方法虽已经被广泛使用，但因为其涉及到Matlab，故更多人怕麻烦不愿采用。因此我将一系列工作内容进行一个细致说明，希望帮助到大家。

图像

拍照

• 如何进行拍照？当拍摄时镜头角度有一个较小偏移时，也会导致最后结果有很大的误差，这时更多的人想到的是固定调平镜头。我想说这样岂不是太麻烦了，在这里我给大家推荐一个方法，采用软件对图像进行垂直水平调整。
• 我们现在拍摄背景后面放置一张绘有标准矩形的白色打印纸，当然也可以提前将矩形打印在纸上，然后再拍摄，拍摄完成通过软件调试，将图像上的矩形刚好铺满整个画面。该类软件比较多，我就不做过多的赘述。

图像处理

• 这里我们用到了PhotoShop，怎么处理的我先放几张图片，网上教程太多了，我也不细说，今天我们的任务在于代码。
  第一张图片是原始图片，我这里用仿真图片代理代替，第二张是处理完之后的图片，图像大小是804x804像素（可以在PS里面调节，若不是此像素，MATLAB会报错）。

Matlab处理代码

话不多说，上代码！！！
这一段单纯是处理PS里面不想去干的，对图像进行简单的二值化处理并裁剪一半（原因不用多说，如果不裁剪，最后拟合函数就是个二次函数）。

%读取一张图片，并显示（C:\Users\ASUS\Desktop\Matlab-LIJunru\15.jpg这部分内容改为处理完的照片位置）
original_picture=imread('D:\Matlab-LIJunru\15.jpg');%横竖像素点为804
% %把图像转换成灰度图像
% GrayPic=rgb2gray(original_picture);%把RGB图像转化成灰度图像
% figure(2)
% imshow(GrayPic);
% title('RGB图像转化为灰度图像')
% 
%对图像进行二值化处理
figure(1);
subplot(2,4,1);
imshow(original_picture);
title('原始图像')

%thresh=graythresh(original_picture);
thresh = 0.95;
Pic=im2bw(original_picture,thresh);
figure(1);
subplot(2,4,2);
imshow(Pic);
title('RGB图像转化为二值化图像')
imwrite(Pic,'Pic.jpg');
%%
BW1 = Pic;
BW2 = bwperim(BW1,8);
BW3 = imcomplement(BW2);
BW3(:,1) = 1;
BW3(:,804) = 1;
BW3(1,:) = 1;

figure(1);
subplot(2,4,3);
imshow(BW2);
title('边界')
subplot(2,4,4);
imshow(BW3);
title('反相')
imwrite(BW3,'边界.jpg');

BW3(:,[403:804]) = 1;
figure(1);
subplot(2,4,5);
imshow(BW3);
title('取一半'<