Думаю так:
for i:=1 to 16 do
for j:=1 to 16 do
Begin
xStart := 1 + (i-1)*16;
yStart := 1 + (j-1)*16;
kb:=0;
for ImX:=xStart to xStart + 15 do
for ImY:=yStart to yStart +15 do
if a[ImX,ImY]=1 then kb:=kb+1;//kb -количество черных пикселей
if kb>porog then b[i,j]:=1 else b[i,j]:=0; {porog макс знач черных точек в масиве}
End;