Method: CED
C:      1.000000
sigma:  0.500000
rho:    5.000000
alpha:  0.001000
ht:     0.250000
kmax:   10000
dir_v:  20.000000
dir_h:  0.000000

bwthreshold: 0.41
im2bw(result_ced, bwthreshold)

bwmorph(result_im2bw, 'clean', Inf)

Original file
alpha: 410.2162 beta: 128.7239 gamma: 3.1868

Modified file
alpha: 259 beta: 3.2086 gamma: 80.7217