《一种基于时空调制的血流速度测量装置及方法.pdf》由会员分享，可在线阅读，更多相关《一种基于时空调制的血流速度测量装置及方法.pdf（11页珍藏版）》请在专利查询网上搜索。
1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201611033407.X (22)申请日 2016.11.14 (71)申请人 佛山科学技术学院 地址 528231 广东省佛山市禅城区江湾一 路18号 (72)发明人 曾亚光 韩定安 王茗祎 张新 (74)专利代理机构 广州市红荔专利代理有限公 司 44214 代理人 吝秀梅 李彦孚 (51)Int.Cl. A61B 5/026(2006.01) A61B 5/0285(2006.01) A61B 5/00(2006.01) (54)发明名称 一种基于时空调制的血流速度测量。
2、装置及 方法 (57)摘要 本发明涉及生物组织血流成像技术领域， 更 具体地， 涉及一种基于时空调制的血流速度测量 装置及方法， 该方法包括以下步骤： 将光源通过 端口照射在置于显微镜载物台上的被测样品， 对 显微镜进行聚焦， 聚焦结束后用成像装置采集显 微镜目镜中被测样品的虚像作为图像数据传输 至电脑， 该图像数据的信噪比高， 电脑对采集到 的图像数据进行时间调制后再进行空间调制， 以 进一步消除信号中的噪声， 获取更准确的血红细 胞信号， 进而利用自相关算法及图像分割方法分 别计算血红细胞的渡越时间及血红细胞的长度， 并将血红细胞的长度与血红细胞的渡越时间的 比值作为被测样品的血流速度， 。
3、计算的血流速度 的精度高。 权利要求书2页 说明书6页 附图2页 CN 106419890 A 2017.02.22 CN 106419890 A 1.一种基于时空调制的血流速度测量装置， 其特征在于， 包括电脑(1)、 成像装置(2)、 显微镜(3)、 端口(4)及光源(5)； 光源(5)通过端口(4)照射在置于显微镜(3)载物台(7)上的 被测样品(6)， 成像装置(2)固定于显微镜(3)目镜上方， 成像装置(2)用于采集被测样品(6) 的图像数据并将所述图像数据传输至电脑(1)， 电脑(1)对接收的图像数据进行处理以获取 被测样品(6)的血流速度。 2.根据权利要求1所述的一种基于时空调。
4、制的血流速度测量装置， 其特征在于， 所述成 像装置(2)为COMS相机。 3.根据权利要求1至2任一项所述的一种基于时空调制的血流速度测量装置， 其特征在 于， 所述光源(5)发出的光为白光。 4.一种基于时空调制的血流速度测量方法， 其特征在于， 使用权利要求3所述的基于时 空调制的血流速度测量装置， 包括以下步骤： S1.采集被测样品(6)的图像数据： 将被测样品(6)放在显微镜(3)的载物台(7)上， 并使 被测样品(6)位于显微镜(3)的物镜的正下方， 将光源(5)发出的光通过端口(4)照射在被测 样品(6)上， 对显微镜(3)进行聚焦， 聚焦结束后用成像装置(2)采集被测样品(6)。
5、在显微镜 (3)目镜中的成像以作为采集的被测样品(6)的图像数据， 成像装置(2)采集的图像总帧数 为F， 成像装置(2)的帧率为f； S2.成像装置(2)将采集到的被测样品(6)的图像数据传输至电脑(1)， 电脑(1)对接收 的图像数据进行处理以得到被测样品(6)的血流速度， 包括如下步骤： S21.对采集到的F帧图像数据逐一像素点进行时间调制后再进行空间调制， 以将背景 组织信号与血红细胞信号进行分离， 得到血红细胞信号Ist(x， y， t)； S22.对采集到的血红细胞信号Ist(x， y， t)进行自相关运算， 以获取血红细胞的渡越时 间； S23.对采集到的血红细胞信号Ist(x，。
6、 y， t)进行图像分割， 以获取血红细胞的长度； S24.用血红细胞的长度与血红细胞的渡越时间的比值作为血红细胞的流动速度， 即被 测样品(6)的血流速度。 5.根据权利要求4所述的一种基于时空调制的血流速度测量方法， 其特征在于， 步骤 S21中， 对F帧图像数据上每个像素点处的原始信号I(x， y， t)进行快速傅里叶变换， 对快速 傅里叶变换后的信号FFTI(x,y,t)进行第一次带通滤波， 对第一次带通滤波后的信号IRBC (f)进行快速逆傅里叶变换， 以得到经过时间调制后的血红细胞信号It(x,y,t)； 对经过时 间调制后的血红细胞信号It(x,y,t)进行二维空间傅里叶变换， 。
7、对二维空间傅里叶变换后 的信号Ist(fx,fy,t)进行第二次带通滤波， 对第二次带通滤波后的信号进行二维空间逆傅 里叶变换， 以得到经过时空调制后的血红细胞信号Ist(x， y， t)； I(x,y,t)I0+IN(x,y,t)+IRBC(x,y,t) (1) 其中， I(x， y， t)为像素点(x， y)位置处的原始信号， I0为背景组织散射产生的散射光强 度， IN(x,y,t)为噪声信号强度， IRBC(x,y,t)为运动的血红细胞散射产生的散射光强度， 信 号为非周期的脉冲信号； 权 利 要 求 书 1/2 页 2 CN 106419890 A 2 其中， N为采集过程中通过的血。
8、红细胞个数， M为血红细胞散射光强幅值， rect(t-ti- 0/2)/ 0为单位矩形函数， 0为矩形函数的宽度， ti是第i个血红细胞通过该像素位置的起 始时间； 其中，为执行第一次带通滤波， 带通窗口范围为f1f2， f1为f/F， f2f/2， f为成 像装置(2)的帧率， F为采集的图像总帧数， FFTI(x,y,t)为进行快速傅里叶变换； It(x,y,t)FFT-1IRBC(x,y,f) (4) 其中， FFT-1为进行快速逆傅里叶变换； Ist(fx,fy,t)FFT2It(x,y,t) (5) 其中， FFT2为进行二维空间傅里叶变换； 第二次带通滤波的带通窗口范围为f3f4。
9、， f3 0， f4f/2-10Hz； Ist(x,y,t)IFFT2Ist(fx,fy,t) (6) 其中， IFFT2为进行二维空间逆傅里叶变换。 6.根据权利要求4所述的一种基于时空调制的血流速度测量方法， 其特征在于， 步骤 S22中利用公式(7)得到血红细胞的自相关函数； 其中， E为平均值运算， T为信号的周期， F为成像装置(2)采集的图像总帧数为， 为移 动时间变量； 对自相关函数按公式(8)进行归一化处理， 得到血红细胞的归一化的自相关函数； 所述血红细胞的渡越时间为自相关曲线初始时刻点与曲线纵坐标首次为0时刻点两点 间曲线的斜率的倒数。 7.根据权利要求4所述的一种基于时空。
10、调制的血流速度测量方法， 其特征在于， 步骤 S23中使用单阈值图像分割方法对采集到的血红细胞信号Ist(x， y， t)进行图像分割， 以提取 血红细胞的边界， 计算多个血红细胞沿流经方向的最大长度并取平均值作为血红细胞的长 度。 权 利 要 求 书 2/2 页 3 CN 106419890 A 3 一种基于时空调制的血流速度测量装置及方法 技术领域 0001 本发明涉及生物组织血流成像技术领域， 更具体地， 涉及一种基于时空调制的血 流速度测量装置及方法。 背景技术 0002 当今许多病情的分析都通过血液来分析， 而对血流速度的分析可以判断组织器官 等状态的好坏， 可以分析高压病、 冠心病。
11、等病情， 对治疗有很大的帮助， 所以市场上急需要 血流速度测量方法和测量装置。 0003 现有的血流速度测量方法主要有超声多普勒法、 激光多普勒法和激光散斑衬比成 像法。 超声和激光多普勒的方法都是激光散斑衬比血流成像法。 上述方法容易受多种成像 环境因素影响， 获取的图像数据信噪比较低， 因此基于图像数据计算的血流速度精度较低。 发明内容 0004 本发明的目的在于克服现有技术的不足， 提供一种基于时空调制的血流速度测量 装置， 该装置不易受成像环境的影响， 获取的图像数据信噪比高， 基于图像数据计算的血流 速度精度高。 0005 为解决上述技术问题， 本发明采用的技术方案是： 0006 提。
12、供一种基于时空调制的血流速度测量装置， 包括电脑、 成像装置、 显微镜、 端口 及光源； 光源通过端口照射在置于显微镜载物台上的被测样品， 成像装置固定于显微镜目 镜上方， 成像装置用于采集被测样品的图像数据并将所述图像数据传输至电脑， 电脑对接 收的图像数据进行处理以获取被测样品的血流速度。 0007 上述方案中， 通过设置电脑、 成像装置、 显微镜、 端口及光源， 使得当光源通过端口 照射置于显微镜载物台上的被测样品时， 显微镜能对被测样品造影成像， 成像装置采集被 测样品的图像数据后将其传输至电脑进行处理， 即可获得被测样品的血流速度。 本发明一 种基于时空调制的血流速度测量装置， 不易。
13、受成像环境的影响， 获取的图像数据信噪比高， 基于图像数据计算的血流速度精度高。 0008 优选地， 所述成像装置为COMS相机。 COMS相机不但能获取高空间分辨率的图像数 据， 进而提高计算的血流速度的精度， 而且能节省成本。 0009 优选地， 为了便于使用者判定毛细血管的位置， 将光源发出的光设置为白光。 0010 本发明的另一个目的是提供一种基于时空调制的血流速度测量方法， 该方法使用 上述基于时空调制的血流速度测量装置， 包括如下步骤： 0011 S1.采集被测样品的图像数据： 将被测样品放在显微镜的载物台上， 并使被测样品 位于显微镜的物镜的正下方， 将光源发出的光通过端口照射在。
14、被测样品上， 对显微镜进行 聚焦， 聚焦结束后用成像装置采集被测样品在显微镜目镜中的成像以作为采集的被测样品 的图像数据， 成像装置采集的图像总帧数为F， 成像装置的帧率为f； 0012 S2.成像装置将采集到的被测样品的图像数据传输至电脑， 电脑对接收的图像数 说 明 书 1/6 页 4 CN 106419890 A 4 据进行处理以得到被测样品的血流速度， 包括如下步骤： 0013 S21.对采集到的F帧图像数据逐一像素点进行时间调制后再进行空间调制， 以将 背景组织信号与血红细胞信号进行分离， 得到血红细胞信号Ist(x， y， t)； 0014 S22.对采集到的血红细胞信号Ist(x。
15、， y， t)进行自相关运算， 以获取血红细胞的渡 越时间； 0015 S23.对采集到的血红细胞信号Ist(x， y， t)进行图像分割， 以获取血红细胞的长度； 0016 S24.用血红细胞的长度与血红细胞的渡越时间的比值作为血红细胞的流动速度， 即被测样品的血流速度。 0017 本发明一种基于时空调制的血流速度测量方法， 对采集到的图像数据进行时间调 制后再进行空间调制， 能进一步消除信号中的噪声， 获取更准确的血红细胞信号， 进而利用 自相关算法及图像分割方法分别计算血红细胞的渡越时间及血红细胞的长度， 并将血红细 胞的长度与血红细胞的渡越时间的比值作为被测样品的血流速度， 计算的血流。
16、速度的精度 高。 0018 优选地， 步骤S21中， 对F帧图像数据上每个像素点处的原始信号I(x， y， t)进行快 速傅里叶变换， 对快速傅里叶变换后的信号FFTI(x,y,t)进行第一次带通滤波， 对第一次 带通滤波后的信号IRBC(f)进行快速逆傅里叶变换， 以得到经过时间调制后的血红细胞信号 It(x,y,t)； 对经过时间调制后的血红细胞信号It(x,y,t)进行二维空间傅里叶变换， 对二 维空间傅里叶变换后的信号Ist(fx,fy,t)进行第二次带通滤波， 对第二次带通滤波后的信 号进行二维空间逆傅里叶变换， 以得到经过时空调制后的血红细胞信号Ist(x， y， t)； 0019。
17、 I(x,y,t)I0+IN(x,y,t)+IRBC(x,y,t) (1) 0020 其中， I(x， y， t)为像素点(x， y)位置处的原始信号， I0为背景组织散射产生的散射 光强度， IN(x,y,t)为噪声信号强度， IRBC(x,y,t)为运动的血红细胞散射产生的散射光强 度， 信号为非周期的脉冲信号； 0021 0022 其中， N为采集过程中通过的血红细胞个数， M为血红细胞散射光强幅值， rect(t- ti- 0/2)/ 0为单位矩形函数，0为矩形函数的宽度， ti是第i个血红细胞通过该像素位置的 起始时间， IRBC(x,y,t)远远小于I0， 以达到将背景组织信号和血。
18、红细胞信号进行分离的目 的， 从而得到血红细胞信号； 0023 0024其中，为执行第一次带通滤波， 带通窗口范围为f1f2， f1为f/F， f2f/2， f 为成像装置(2)的帧率， F为采集的图像总帧数， FFTI(x,y,t)为进行快速傅里叶变换； 0025 It(x,y,t)FFT-1IRBC(x,y,f) (4) 0026 其中， FFT-1为进行快速逆傅里叶变换； 0027 Ist(fx,fy,t)FFT2It(x,y,t) (5) 0028 其中， FFT2为进行二维空间傅里叶变换； 第二次带通滤波的带通窗口范围为f3 f4， f30， f4f/2-10Hz； 0029 Ist。
19、(x,y,t)IFFT2Ist(fx,fy,t) (6) 说 明 书 2/6 页 5 CN 106419890 A 5 0030 其中， IFFT2为进行二维空间逆傅里叶变换。 经过时间调制后的图像数据会出现类 似于 “散斑” 的噪声点， 对时间调制后的信号再进行空间调制， 能获取更准确的血红细胞信 号， 有利于提高计算的血流速度的精度。 0031 优选地， 步骤S22中利用公式(7)得到血红细胞的自相关函数； 0032 0033 其中， E为平均值运算， T为信号的周期， F为成像装置(2)采集的图像总帧数为， 为移动时间变量； 0034 对自相关函数按公式(8)进行归一化处理， 得到血红细。
20、胞的归一化的自相关函数； 0035 0036 所述血红细胞的渡越时间为自相关曲线初始时刻点与曲线纵坐标首次为0时刻点 两点间曲线的斜率的倒数。 0037 优选地， 步骤S23中使用单阈值图像分割方法对采集到的血红细胞信号Ist(x， y， t) 进行图像分割， 以提取血红细胞的边界， 计算多个血红细胞沿流经方向的最大长度并取平 均值作为血红细胞的长度。 取平均值能有效去除细胞形状不规则的影响， 提高获取的血红 细胞的长度。 单阈值图像分割方法具体做法为： 对一幅大小为CD像素的灰度图像， 灰度级 数为L,像素灰度级为Ist(x,y,t)， 阈值为l， 对所有像素的灰度级做映射： 0038 00。
21、39 公式(9)将灰度图像变为二值图像， 此时血红细胞所在区域赋值为0， 背景组织区 域赋值为1， 此时可以通过0和1的分解来获取血红细胞的边界。 0040 与现有技术相比， 本发明的有益效果是： 0041 本发明一种基于时空调制的血流速度测量装置， 通过设置电脑、 成像装置、 显微 镜、 端口及光源， 使得当光源通过端口照射置于显微镜载物台上的被测样品时， 显微镜能对 被测样品造影成像， 成像装置采集被测样品的图像数据后将其传输至电脑进行处理， 即可 获得被测样品的血流速度， 该装置不易受成像环境的影响， 获取的图像数据信噪比高， 基于 图像数据计算的血流速度精度高； 本发明一种基于时空调制。
22、的血流速度测量方法， 对采集 到的图像数据进行时间调制后再进行空间调制， 能进一步消除信号中的噪声， 获取更准确 的血红细胞信号， 进而利用自相关算法及图像分割方法分别计算血红细胞的渡越时间及血 红细胞的长度， 并将血红细胞的长度与血红细胞的渡越时间的比值作为被测样品的血流速 度， 计算的血流速度的精度高。 附图说明 0042 图1为本实施例一种基于时空调制的血流速度测量装置的示意图。 0043 图2为本实施例一种基于时空调制的血流速度测量方法的流程图， 其中矩形框内 为处理内容或结果， 圆角矩形框内为处理方法。 0044 图3为本实施例一种基于时空调制的血流速度测量方法中确定血红细胞的渡越时。
23、 间的示意图。 说 明 书 3/6 页 6 CN 106419890 A 6 具体实施方式 0045 下面结合具体实施方式对本发明作进一步的说明。 其中， 附图仅用于示例性说明， 表示的仅是示意图， 而非实物图， 不能理解为对本专利的限制； 为了更好地说明本发明的实 施例， 附图某些部件会有省略、 放大或缩小， 并不代表实际产品的尺寸； 对本领域技术人员 来说， 附图中某些公知结构及其说明可能省略是可以理解的。 0046 本发明实施例的附图中相同或相似的标号对应相同或相似的部件； 在本发明的描 述中， 需要理解的是， 若有术语 “上” 、“下” 、“左” 、“右” 等指示的方位或位置关系为基于。
24、附图 所示的方位或位置关系， 仅是为了便于描述本发明和简化描述， 而不是指示或暗示所指的 装置或元件必须具有特定的方位、 以特定的方位构造和操作， 因此附图中描述位置关系的 用语仅用于示例性说明， 不能理解为对本专利的限制， 对于本领域的普通技术人员而言， 可 以根据具体情况理解上述术语的具体含义。 0047 实施例 0048 本实施例一种基于时空调制的血流速度测量装置的示意图如图1所示， 包括电脑 1、 成像装置2、 显微镜3、 端口4及光源5； 光源5通过端口4照射在置于显微镜3的载物台7上的 被测样品6， 成像装置2固定于显微镜3目镜上方， 成像装置2用于采集被测样品6的图像数据 并将所。
25、述图像数据传输至电脑1， 电脑1对接收的图像数据进行处理以获取被测样品6的血 流速度。 0049 使用该血流速度测量装置测量被测样品6的毛细血管内的血流速度时， 将被测样 品6置于显微镜3的载物台7上， 并使被测样品6位于显微镜3的物镜下， 将光源5发出的光通 过端口4照射在被测样品6上， 对显微镜3进行聚焦， 聚焦结束后用成像装置2采集被测样品6 在显微镜3目镜中的成像， 成像装置2将采集的图像数据传输至电脑处理， 即可获得被测样 品6的血流速度。 本发明一种基于时空调制的血流速度测量装置， 不易受成像环境的影响， 获取的图像数据信噪比高， 基于图像数据计算的血流速度精度高。 0050 其中。
26、， 所述成像装置2为COMS相机。 COMS相机不但能获取高空间分辨率的图像数 据， 进而提高计算的血流速度的精度， 而且能节省成本。 0051 另外， 所述光源5发出的光为白光。 白光的设置便于使用者判定被测样品6的毛细 血管的位置， 有利于成像装置2获取高质量的的图像数据。 0052 本发明还提供了一种基于时空调制的血流速度测量方法， 该测量方法使用上述基 于时空调制的血流速度测量装置， 其流程图如图2所示， 包括以下步骤： 0053 S1.采集被测样品6的图像数据： 将被测样品6放在显微镜3的载物台7上， 并使被测 样品6位于显微镜3的物镜的正下方， 将光源5发出的光通过端口4照射在被测。
27、样品6上， 对显 微镜3进行聚焦， 聚焦结束后用成像装置2采集被测样品6在显微镜3目镜中的成像以作为采 集的被测样品6的图像数据， 成像装置2采集的图像总帧数为F， 成像装置2的帧率为f； 0054 S2.成像装置2将采集到的被测样品6的图像数据传输至电脑1， 电脑1对接收的图 像数据进行处理以得到被测样品6的血流速度， 包括如下步骤： 0055 S21.对采集到的F帧图像数据逐一像素点进行时间调制后再进行空间调制， 以将 背景组织信号与血红细胞信号进行分离， 得到血红细胞信号Ist(x， y， t)； 0056 S22.对采集到的血红细胞信号Ist(x， y， t)进行自相关运算， 以获取血。
28、红细胞的渡 说 明 书 4/6 页 7 CN 106419890 A 7 越时间； 0057 S23.对采集到的血红细胞信号Ist(x， y， t)进行图像分割， 以获取血红细胞的长度； 0058 S24.用血红细胞的长度与血红细胞的渡越时间的比值作为血红细胞的流动速度， 即被测样品6的血流速度。 0059 本发明一种基于时空调制的血流速度测量方法， 对采集到的图像数据进行时间调 制后再进行空间调制， 能进一步消除信号中的噪声， 获取更准确的血红细胞信号， 进而利用 自相关算法及图像分割方法分别计算血红细胞的渡越时间及血红细胞的长度， 并将血红细 胞的长度与血红细胞的渡越时间的比值作为被测样品。
29、的血流速度， 计算的血流速度的精度 高。 0060 其中， 步骤S21中， 对F帧图像数据上每个像素点处的原始信号I(x， y， t)进行快速 傅里叶变换， 对快速傅里叶变换后的信号FFTI(x,y,t)进行第一次带通滤波， 对第一次带 通滤波后的信号IRBC(f)进行快速逆傅里叶变换， 以得到经过时间调制后的血红细胞信号 It(x,y,t)； 对经过时间调制后的血红细胞信号It(x,y,t)进行二维空间傅里叶变换， 对二 维空间傅里叶变换后的信号Ist(fx,fy,t)进行第二次带通滤波， 对第二次带通滤波后的信 号进行二维空间逆傅里叶变换， 以得到经过时空调制后的血红细胞信号Ist(x， 。
30、y， t)； 0061 I(x,y,t)I0+IN(x,y,t)+IRBC(x,y,t) (1) 0062 其中， I(x， y， t)为像素点(x， y)位置处的原始信号， I0为背景组织散射产生的散射 光强度， IN(x,y,t)为噪声信号强度， IRBC(x,y,t)为运动的血红细胞散射产生的散射光强 度， 信号为非周期的脉冲信号； 0063 0064 其中， N为采集过程中通过的血红细胞个数， M为血红细胞散射光强幅值， 取决于血 红细胞的吸收系数， 可以近似认为其为常数， rect(t-ti- 0/2)/ 0为单位矩形函数， 0为 矩形函数的宽度， ti是第i个血红细胞通过该像素位置。
31、的起始时间， IRBC(x,y,t)远远小于I0， 以达到将背景组织信号和血红细胞信号进行分离的目的， 从而得到血红细胞信号； 0065 0066其中，为执行第一次带通滤波， 带通窗口范围为f1f2， f1为f/F， f2f/2， f 为成像装置2的帧率， F为采集的图像总帧数， FFTI(x,y,t)为进行快速傅里叶变换； 0067 It(x,y,t)FFT-1IRBC(x,y,f) (4) 0068 其中， FFT-1为进行快速逆傅里叶变换； 0069 Ist(fx,fy,t)FFT2It(x,y,t) (5) 0070 其中， FFT2为进行二维空间傅里叶变换； 第二次带通滤波的带通窗口。
32、范围为f3 f4， f30， f4f/2-10Hz； 0071 Ist(x,y,t)IFFT2Ist(fx,fy,t) (6) 0072 其中， IFFT2为进行二维空间逆傅里叶变换。 经过时间调制后的图像数据会出现类 似于 “散斑” 的噪声点， 对时间调制后的信号再进行空间调制， 能获取更准确的血红细胞信 号， 有利于提高计算的血流速度的精度。 0073 本实施例中， 步骤S22中利用公式(7)得到血红细胞的自相关函数； 说 明 书 5/6 页 8 CN 106419890 A 8 0074 0075 自相关函数描述了信号Ist(x， y， t)在两个不同时刻t、 t+ 取值之间的相关程度，。
33、 其 中， E为平均值运算， T为信号的周期， F为成像装置(2)采集的图像总帧数为， 为移动时 间变量； 0076 对自相关函数按公式(8)进行归一化处理， 得到血红细胞的归一化的自相关函数； 0077 0078 假设血红细胞的自相关曲线如图3所示， 则血红细胞的渡越时间为自相关曲线初 始时刻点A与曲线纵坐标首次为0时刻点B两点间曲线的斜率的倒数。 0079 其中， 步骤S23中使用单阈值图像分割方法对采集到的血红细胞信号Ist(x， y， t)进 行图像分割， 以提取血红细胞的边界， 计算多个血红细胞沿流经方向的最大长度并取平均 值作为血红细胞的长度。 取平均值能有效去除细胞形状不规则的影。
34、响， 提高获取的血红细 胞的长度。 单阈值图像分割方法具体做法为： 对一幅大小为CD像素的灰度图像， 灰度级数 为L,像素灰度级为Ist(x,y,t)， 阈值为l， 对所有像素的灰度级做映射： 0080 0081 公式(9)将灰度图像变为二值图像， 此时血红细胞所在区域赋值为0， 背景组织区 域赋值为1， 此时可以通过0和1的分解来获取血红细胞的边界。 0082 显然， 本发明的上述实施例仅仅是为清楚地说明本发明所作的举例， 而并非是对 本发明的实施方式的限定。 对于所属领域的普通技术人员来说， 在上述说明的基础上还可 以做出其它不同形式的变化或变动。 这里无需也无法对所有的实施方式予以穷举。 凡在本 发明的精神和原则之内所作的任何修改、 等同替换和改进等， 均应包含在本发明权利要求 的保护范围之内。 说 明 书 6/6 页 9 CN 106419890 A 9 图1 图2 说 明 书 附 图 1/2 页 10 CN 106419890 A 10 图3 说 明 书 附 图 2/2 页 11 CN 106419890 A 11 。