形态学颗粒分析在细胞尺寸分布中应用.docVIP

形态学颗粒分析在细胞尺寸分布中应用 　　摘 要：以数学形态学中开闭运算及其相关性质为基础，研究形态学颗粒分析方法的相关性质。在实验中，先将获取的彩色细胞图像转化为二值图像,利用形态学开运算滤除小噪声斑点后,再采用形态学开颗粒分析,在有效保证细胞形状和大小的前提下,统计出细胞的大小分布情况。实验表明，该方法能有效地给出细胞体大小分布。? 　　关键词：开闭运算；A-开运算；颗粒分析；细胞分布 ? 　　中图分类号：TP399 文献标识码：A 文章编号：1672-7800（2012）003-0155-03?? 　　? 　　作者简介：卢峥（1983-），男，湖北武汉人，中南民族大学数学与统计学学院硕士研究生，研究方向为数学应用方法与图像处理。 　　 　　 　　1 理论基础? 　　在数学形态学中，开闭运算是所有形态学运算的基础：如果一个集合?X被集合A先腐蚀再膨胀，那么可以得到一个比原来更小的集合，将其称为A对X作开运算的结果，记为X?A=(XΘA)?A，常用α?A(X)表示X被A所作的开运算，即：α?A(X)=X?A=(XΘA)?A?。? 　　如果一个集合?X?被集合A先膨胀再腐蚀，那么可以得到一个比原来更大的集合，将其称为A对?X作闭运算的结果，记为X?A=(X?A)ΘA，常用β?A(X)表示X?被A所作的闭运算，即：?β?A(X)=X?A=(X?A)ΘA?。? 　　由定义可知开闭运算具有：①递增性若?X?Y，被同一集合A作开运算有：α?A(X)?α?A(Y);被同一集合A作闭运算有：β?A(X)?β?A(Y)?；②幂等性?(X?A)?A=X?A，(X?A)?A=X?A?；③开运算具有非扩张性?α?A(X)?X；闭运算具有扩张性β?A(X)?X?。? 　　文献［4］介绍了A-开运算。A-开运算作为建立形态学颗粒分析方法的基础，下面具体讨论A-开运算及其两个相关命题的证明。? 　　定义：对于给定的集合?X?，若有集合A使得?X?A=X，则称X?是A-开的。? 　　下面论证关于A-开运算的两个命题：? 　　引理1：集合?X?是A-开的充要条件是存在集合?Y使得X=Y?A?。? 　　证明必要性：若?X?是A-开的，那么?X=X×A=(XΘA)?A，令Y=XΘA，即得X=Y?A?。? 　　充分性由闭运算的扩展性知：?(Y?A)ΘA=Y?A?Y，如果X=Y?A?， 　　因为膨胀是递增运算?((Y?A)ΘA)?A=(Y?A)?A?Y?A，于是((Y?A)ΘA)?A?Y?A=X?。? 　　另一方面,由开运算的非扩展性知：?((Y?A)ΘA)?A=(Y?A)?A?Y?A=X，于是X?A=((Y?A)ΘA)?A=(Y?A)?A=Y?A=X?，即证。? 　　引理2：设集合?A,B?满足A是B-开的，那么对于每个集合?X?R?2?有： 　　??X?A?X?B（1） ? 　　(X?A)?B=(X?B)?A=X?A（2）?? 　　 证明根据膨胀满足结合律和交换律:?A?(B?C)=(A?B)?C，A?B=B?A?。? 　　因为A是B-开的，由引理1存在集合?C使得A=B?C=C?B?。又因为：? 　　?(X?A)?B=((XΘA)?A)?B=((XΘA)?(B?C))?B? 　　=(((XΘA)?C)?B)?B=［(((XΘA)?C)?B)ΘB］?B?? 　　将?(XΘA)?C看作一个整体，由引理1证明过的结论((Y?A)ΘA)?A=Y?A可得： ? 　　=［(((XΘA)?C)?B)ΘB］?B=(((XΘA)?C)?B)?? 　　又膨胀满足结合律，可得：? 　　 ?=(((XΘA)?C)?B)=(XΘA)?(C?B)=(XΘA)?A=X?A?? 　　综上得到?(X?A)?B=X?A? 　　根据开运算的非扩展性X?A?X，再由开运算的递增性：? 　　X?A=(X?A)?B?X?B?? 　　即证明了(1)成立。? 　　为了证明?(X?B)?A=X?A，因为(X?B)?A?X?A是显然成立的。由结论（1）有：? 　　X?A=(X?A)?A?(X?B)?A?? 　　即证明?(X?B)?A=X?A?，因此（2）成立。? 　　2 颗粒分析? 　　按照A-开运算的原理，选择一个基元A，令?A?1=A，A?2=A?A=A?1?A，A?3=A?A?A=A?2?A，……，结合前面关于A-开给出的两个引理，可知对所有的r，A??r+1?是A?r?-开的。结合引理2中的结论（1）得到递增的图像序列： 　　??S?A?1?S?A?2?S?A?3?…?? 　　 可以看到在?R?2中，令{rA}（r＞0）是由结构元组成的集合，因此对于（s＞r）有sA是rA-开的。假设X是形状和尺寸不一的颗粒组成的集合，用