全局Threshold和动态阈值分割Dyn_Threshold的应用场景

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

处理图片和结果图片：

以下是代码：

 

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

dev_display (Image)
*显示图像  

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

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

dev_display (Image)
dev_set_draw ('margin')
dev_set_line_width ()
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, , )
*平滑处理图像  

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

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

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

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

dev_set_colored ()
dev_display (SmallConnection)  

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