动态阈值分割 dyn_threshold

手册里面的particle例子，例子的任务是分析颗粒在液体中。在这个应用程序的主要困难：存在两种类型的对象:大明亮物体和较低的小物体的对比。此外噪音使分割的存在困难；无法使用全局灰度阈值threshold进行分割；所以采用先将大块不需要检测的部分去除掉，再通过灰度动态阈值dyn_threshold分割图像得到想要的内容。

处理图片和结果图片：

以下是代码：

 read_image (Image, 'particle')
 *获取图像

 dev_display (Image)
 *显示图像

 threshold (Image, Large, 110, 255)
 *灰度阈值分割图像

 dilation_circle (Large, LargeDilation, 7.5)
 *圆角膨胀

 dev_display (Image)
 dev_set_draw ('margin')
 dev_set_line_width (3)
 dev_set_color ('red')
 dev_display (LargeDilation)
 dev_set_draw ('fill')
 *显示图像

 complement (LargeDilation, NotLarge)
 *返回补充图像，即获得去除大斑点后的图像NotLarge

 reduce_domain (Image, NotLarge, ParticlesRed)
 *减去除了NotLarge图像，即去除大斑点后的图像，减少运算

 mean_image (ParticlesRed, Mean, 31, 31)
 *平滑处理图像

 dyn_threshold (ParticlesRed, Mean, SmallRaw, 3, 'light')
 *选择灰度阈值；网友详解：当前背景之间差异明显时，可以设定全局阈值进行threshold，但很多情况下由于背景不均一，
 *目标体经常表现为比背景局部亮一些或暗一些，无法确定全局阈值操作，需要通过其邻域找到一个合适的阈值进
 *行分割。ThresholdImage是参考图像，通过与OrigImage对比找到领域确定阈值，一般采用平滑滤波算子
 *（如mean_image）获取参考图像。offset设定邻域比较的区间范围，灰度值变化在offset范围内均是可以接受的。

 opening_circle (SmallRaw, Small, 2.5)
 *消除小区域(小于圆形结构元素)和光滑的边界地区

 connection (Small, SmallConnection)
 *显示联通区域

 dev_display (Image)
 *这句不加窗口显示效果会有雪花

 dev_set_colored (12)
 dev_display (SmallConnection)
 *显示结果图像

下图是只用threshold时候的实验效果，无法分割出小斑点：