Площадь Площадь Плоских Квадратов

Question

Я пытаюсь вычислить площадь общей площади рассеянных квадратов, как изображение ниже.

Я знаю, что есть функция, с помощью которой мы могли бы оценить размер области, сделав всеобъемлющую выпуклую, но я не уверен, могу ли я использовать какой-либо подобный метод для решения этой проблемы.

Есть ли у кого-то есть хорошая идея?

Изображение:

Edit [01/01/2016]: Код я попытался, как это. Он работает, но каждый раз он должен сохранять и загружать изображение .png. Я хотел спросить, можно ли вычислить размер без процесса сохранения изображений.

clear 

b=zeros(127,2); %matrix 

for i=1:127 
    rnd=randn(1,2); 
    b(i,:)=30.*rnd; 
end 

BW=scatter(b(:,1),b(:,2),15,[],'r','s','filled'); 

view(2) %view from Z+ 
axis off 
saveas(gcf,'scatter.png') 

close all 

BWbase = imread('scatter.png'); %import 
BW = im2bw(BWbase,0.5); %convert to binary data 
imshow(BW); 
bwarea(BW) 

Answer 1

Вам не нужно, чтобы сохранить и загрузить изображение если вы знаете, соответствующие функции. У вашего кода есть некоторые проблемы, которые я исправил, и я объясню здесь:

1- Вам не нужно инициализировать b нулями, а затем заполнить его, используя цикл for по случайным значениям. Вместо этого вы можете делать все, что в одной строке кода:

b = 30 * randn(127,2); 

2- Тогда вы построить точки с помощью scatter. Я удалил [], чтобы исправить эту команду.

BW = scatter(b(:,1),b(:,2),15,'r','s','filled'); 

Обратите внимание, что это не изображение, а фигура, которую нужно преобразовать в изображение. Функция, которую вы ищете, - getframe. Эта функция записывает кадр из рисунка. Вы можете позвонить ему, предоставив ему текущую цифру gcf.

F = getframe(gcf); 

4- Теперь ваш кадр в хранится в F, но это не образ еще. Чтобы преобразовать его в изображение, вам нужно использовать функцию frame2im.

[X, Map] = frame2im(F); 

5- Тогда, как вы делали это сами, вы преобразовать изображение в двоичную и рассчитать площадь:

BW2 = im2bw(X,0.5); 
bwarea(BW2) 

Весь код можно увидеть здесь:

b = 30 * randn(127,2); 
BW = scatter(b(:,1),b(:,2),15,'r','s','filled'); 
axis off 
F = getframe(gcf); 
[X, Map] = frame2im(F); 
BW2 = im2bw(X,0.5); 
imshow(BW2); 
bwarea(BW2) 

Для получения дополнительной информации ознакомьтесь со следующими документами: getframe, frame2im.