hofig6s6
Detecting cells of a FISH image using open top-hat

from morph import *

f=mmreadgray('fish_m.tif')

th = mmareaopen(mmopenth(f,mmsedisk(4)),2)

arrow=mmseunion(mmseunion(mmseline(6,45),mmseline(6,-45)),mmseline(10))

fbin = mmbinary(th,50)

fcc  = mmblob(mmlabel(fbin,mmsebox()),'centroid')

mmshow(f)

mmshow(th)

mmshow(fbin)

aux = mmunion(f, mmgray(mmdil(fcc,mmsetrans(arrow,int32([0,5])))))

mmshow(aux)

zoom = aux[9:140,29:200]

mmshow(zoom)