BW = zeros(9,10);
BW(4:6,4:7) = 1
SE = strel('square',3)
BW2 = imdilate(BW,SE)