"Избавиться" от объектов соприкасающихся с границей изображения

% Пример демонстрирует построение гистограммы площадей объектов, не касающихся границ изображения.
% Чтение исходного изображения и вывод его на экран.
BW1=imread(split_bw.tif);
imshow(BWI);
% Поиск объектов.
[L num]=bwlabel(BW1, 8);
% Вычисление признаков объектов:
% площадей и координат центров масс.
feats=imfeature(L, 'Area', 'Centroid', 8);
% Переписывание значений признаков
% из массива структур feats в отдельные массивы.
Areas=zeros(num);
CentX=zeros(num);
CentY=zeros(num);
for i=1:1:num
Areas(i)=feats(l).Area;
CentX(i)=feats(i).Centroid(1);
CentY(i)=feats(i).Centroid(2);
end;
% Формирование бинарного изображения с объектами,
% площадь которых более 30 пикселов, а центр масс объекта удален
% от границ изображения более чем на 20 пикселов.
[rows cols]=size(BW1);
idx=find(Areas>30&CentX>20&CentY>20...
        &CentX<(cols-20)&CentY<(rows-20));
BW1=ismember(Lidx);
figure, imshow(BW1);
% Снова поиск объектов.
[L num]=bwlabel(BW1, 8);
figure, imshow(L+1, jet(num+1));
% Вычисление площадей объектов.
feats=imfeature(L, 'Area', 8);
% Построение гистограммы площадей объектов.
figure, hist([feats, Area]);