def iathick(f, Iab=None, n=-1, theta=45, DIRECTION="CLOCKWISE"): from iaisbinary import iaisbinary from iaintersec import iaintersec from iasupgen import iasupgen from iainterot import iainterot from iaunion import iaunion from iaisequal import iaisequal from iahomothick import iahomothick if Iab is None: Iab = iahomothick() DIRECTION = upper(DIRECTION) assert iaisbinary(f), 'f must be binary image' if n == -1: n = product(f.shape) y = f zero = iaintersec(f, 0) for i in range(n): aux = zero for t in range(0, 360, theta): sup = iasupgen(y, iainterot(Iab, t, DIRECTION)) aux = iaunion(aux, sup) y = iaunion(y, sup) if iaisequal(aux, zero): break return y
def iathick(f, Iab=None, n=-1, theta=45, DIRECTION="CLOCKWISE"): from iaisbinary import iaisbinary from iaintersec import iaintersec from iasupgen import iasupgen from iainterot import iainterot from iaunion import iaunion from iaisequal import iaisequal from iahomothick import iahomothick if Iab is None: Iab = iahomothick() DIRECTION = upper(DIRECTION) assert iaisbinary(f),'f must be binary image' if n == -1: n = product(f.shape) y = f zero = iaintersec(f,0) for i in range(n): aux = zero for t in range(0,360,theta): sup = iasupgen( y, iainterot(Iab, t, DIRECTION)) aux = iaunion( aux, sup) y = iaunion( y, sup) if iaisequal(aux,zero): break return y
def iasupcanon(f, Iab, theta=45, DIRECTION="CLOCKWISE"): from iaintersec import iaintersec from iainterot import iainterot from iaunion import iaunion from iasupgen import iasupgen DIRECTION = upper(DIRECTION) y = iaintersec(f, 0) for t in range(0, 360, theta): Irot = iainterot(Iab, t, DIRECTION) y = iaunion(y, iasupgen(f, Irot)) return y
def iasupcanon(f, Iab, theta=45, DIRECTION="CLOCKWISE"): from iaintersec import iaintersec from iainterot import iainterot from iaunion import iaunion from iasupgen import iasupgen DIRECTION = upper(DIRECTION) y = iaintersec(f,0) for t in range(0,360,theta): Irot = iainterot( Iab, t, DIRECTION ) y = iaunion( y, iasupgen(f, Irot)) return y
def iacthin(f, g, Iab=None, n=-1, theta=45, DIRECTION="CLOCKWISE"): from iaisbinary import iaisbinary from iasupgen import iasupgen from iainterot import iainterot from iaunion import iaunion from iasubm import iasubm from iaisequal import iaisequal from iahomothin import iahomothin if Iab is None: Iab = iahomothin() DIRECTION = upper(DIRECTION) assert iaisbinary(f),'f must be binary image' if n == -1: n = product(f.shape) y = f old = y for i in range(n): for t in range(0,360,theta): sup = iasupgen( y, iainterot(Iab, t, DIRECTION)) y = iaunion( iasubm( y, sup),g) if iaisequal(old,y): break old = y return y