以前ご質問頂いた内容のつづきと想定して、以下のような処理でしょうか?
I = imread('image.jpeg'); % 画像読み込み
BW = imbinarize(rgb2gray(I)); % 2値化
BW = ~BW; % 白黒反転
H = filter2(ones(3),~BW);
BWext = H == 7;
se = strel('square',2);
BWext2 = imdilate(BWext, se); % 値が7のピクセルを膨張
BWout = BW | BWext2; % 元の2値画像とマージ
cMap = jet(10);
% 結果を表示
tiledlayout('flow')
nexttile
imshow(BW)
title('[1] 元のバイナリ画像','FontSize',12)
ax = nexttile;
imagesc(H,[-0.5 9.5])
colormap(ax,cMap)
colorbar
axis tight
axis equal
title('[2] 近傍の1の個数','FontSize',12)
nexttile
imshow(BWext)
title('[3] 値が7のピクセル','FontSize',12)
nexttile
imshow(BWext2)
title('[4] [3]を膨張','FontSize',12)
nexttile
imshow(BWout)
title('[5] 元画像とマージ([1]+[4])','FontSize',12)