CN102198009A - 基于超声背散射信号参量的松质骨诊断系统 - Google Patents

基于超声背散射信号参量的松质骨诊断系统 Download PDF

Info

Publication number
CN102198009A
CN102198009A CN2011101589700A CN201110158970A CN102198009A CN 102198009 A CN102198009 A CN 102198009A CN 2011101589700 A CN2011101589700 A CN 2011101589700A CN 201110158970 A CN201110158970 A CN 201110158970A CN 102198009 A CN102198009 A CN 102198009A
Authority
CN
China
Prior art keywords
ultrasonic
cancellous bone
backscattering
signals
signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN2011101589700A
Other languages
English (en)
Inventor
他得安
王文斌
韩海杰
杨白丁
王威琪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fudan University
Original Assignee
Fudan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fudan University filed Critical Fudan University
Priority to CN2011101589700A priority Critical patent/CN102198009A/zh
Publication of CN102198009A publication Critical patent/CN102198009A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明属于医学超声技术领域,具体为一种基于超声背散射信号参量的松质骨诊断系统。该系统包括:超声背散射信号的获取模块、预处理模块、背散射参量计算模块和松质骨状况评价模块。本发明利用超声收发一体探头从人体跟骨处获取超声背散射信号,并提取出有效的松质骨超声背散射信号,再利用计算模块计算出背散射系数、表观积分背散射系数、频谱质心偏移量和平均骨小梁间距四个参量,最后用这四个参量通过与存储在系统内部标准数据库中的标准值比较来分析松质骨的健康状况。相比于传统的使用宽带超声衰减和超声传播速度的超声透射法诊断系统,本发明诊断系统能得到松质骨结构的完整信息,从而更好地检测人体松质骨的健康状况。

Description

基于超声背散射信号参量的松质骨诊断系统
技术领域
本发明属于医学超声技术领域,具体涉及一种松质骨超声诊断系统。
背景技术
随着人口老龄化的加剧,骨质疏松已经成为影响公众健康的重要问题之一。患有骨质疏松时,会引起患者骨矿密度的下降和骨微结构的退化。目前常用的诊断骨质疏松的方法是双能X射线测量骨密度,但是由于它具有电离辐射且仪器成本比较昂贵,因此不能作为日常的检测手段。超声诊断具有无电离辐射、便捷、速度快、价格低廉等优势而日益受到人们的广泛关注。
目前常用的方法是超声透射法,它使用双探头获取人体跟骨处的超声透射信号,以宽带超声衰减(BUA)和超声传播速度(SOS)反映骨矿密度信息。超声背散射法相比于这种传统的超声透射法具有以下优点:(1)超声背散射法使用收发一体探头,因此相对于超声透射法需要2个探头,它的测量位置更加灵活,易于探测跟骨以外的骨组织(例如脊骨、腕骨等)。(2)研究表明,骨质疏松的引起不仅是骨矿密度的下降,也和骨微结构的变化密切相关,相比于超声透射法只能检测骨矿密度的变化,超声背散射信号中含有骨微结构的全部信息,因此理论上能更加准确地诊断骨质疏松。
目前理论上常用的用超声背散射信号定征松质骨结构状况的参数为背散射系数(BC)、表观积分背散射系数(AIB)、频谱质心偏移参数(SCS)和平均骨小梁间距(MTBS),它们和松质骨的骨矿密度和微结构状况都有较高的相关性。
发明内容
本发明的目的在于提出一种客观性好、诊断正确性高的人体松质骨诊断系统。
本发明提供的人体松质骨诊断系统,是基于超声背散射信号参量的松质骨诊断软件系统。该系统包括:超声背散射信号的获取模块、预处理模块、背散射参量的计算模块、松质骨状况评价模块4个部分。
下面对本发明内容做进一步介绍:
1.超声背散射信号的获取模块:通过使用松质骨超声诊断系统控制超声收发一体探头,在探头表面均匀涂敷医用超声耦合剂后将其置于待测者脚跟内侧较平整的位置,然后进行测量。测量过程中,每隔1ms超声探头发射一组超声脉冲到松质骨并接收超声背散射信号,每次测量发射128组超声脉冲并将所有得到的信号平均,以减小噪声;
2.预处理模块:为得到有效的背散射信号,首先通过阈值判断法来确定背散射信号的开始位置,然后根据待测者输入的体重/身高比来估计他的软组织和皮质骨的厚度,根据采样率和超声传播速度计算出背散射信号中混合的反射信号的长度,然后去除这部分反射信号。最后通过选取感兴趣的松质骨厚度,用同样的方法计算信号长度,最后加时域窗截取出这部分有效的超声背散射信号;
3.背散射参量的计算模块:根据得到的有效背散射信号,再结合松质骨诊断系统中预存的参考信号,先利用FFT分别得到他们的幅度谱                                               
Figure 2011101589700100002DEST_PATH_IMAGE002
Figure 2011101589700100002DEST_PATH_IMAGE004
(其中
Figure 622568DEST_PATH_IMAGE002
是有效的背散射信号的幅度谱,
Figure 114729DEST_PATH_IMAGE004
是参考信号的幅度谱)。再利用以下公式计算背散射参量:
Figure 2011101589700100002DEST_PATH_IMAGE006
                                   (1)
Figure 2011101589700100002DEST_PATH_IMAGE008
                          (2)
                                     (3)
Figure 2011101589700100002DEST_PATH_IMAGE012
                              (4)
其中
Figure 2011101589700100002DEST_PATH_IMAGE014
分别是探头-6dB有效带宽的最小值和最大值;为超声探头的中心频率,为128次测量的背散射信号的频谱均值;
Figure 2011101589700100002DEST_PATH_IMAGE022
为背散射信号的基频。利用公式(1)计算背散射系数(BC),利用公式(2)计算表观积分背散射系数(AIB),利用公式(3)计算频谱质心偏移量(SCS),利用公式(4)计算平均骨小梁间距(MTBS)。
4.松质骨状况评价模块:系统中的数据库保存了不同年龄、不同性别和不同体重/身高比下的各个背散射参量的标准值。在每次测量前,输入待测者的编号、姓名、性别、年龄、身高和体重等个人信息,系统利用计算出来的背散射参量和对应的标准值进行比较,最后综合评价待测用户的松质骨健康状况。
本发明中,背散射信号的基频
Figure 670869DEST_PATH_IMAGE022
的具体计算方法:首先利用简易反向滤波器最小化超声换能器脉冲响应和散射响应的影响,然后利用自相关技术加强背散射信号的周期性,最后用阈值判断和平均的方法求出
Figure 8309DEST_PATH_IMAGE022
。见图2所示。
本发明系统的具体操作步骤为:首先利用松质骨超声诊断硬件系统,控制超声收发一体探头得到人体的松质骨超声背散射信号,然后预处理,得到有效的背散射信号,在根据预存储在诊断系统中的参考信号,分别计算BC、AIB、SCS及MTBS等四个背散射参量,最后将这四个参量与诊断系统内数据库中的标准值进行对照,来判断松质骨的健康状况。
传统的基于超声透射法的诊断系统都是使用宽带超声衰减(BUA)和超声传播速度(SOS)这两个参数来评价人体的骨质状况,与这种方法不同,本发明提出了将背散射系数(BC)、表观积分背散射系数(AIB)、频谱质心偏移量(SCS)和平均骨小梁间距(MTBS)这四个背散射参量用于松质骨诊断系统,从而来更简单有效地评价松质骨的健康状况。
附图说明
图1 为松质骨超声诊断系统结构框图。
图2 为基频
Figure 714097DEST_PATH_IMAGE022
的计算框图。
图3为松质骨超声背散射诊断软件系统使用背散射参量诊断松质骨健康状况系统框图。
图4为待测者的人体松质骨超声背散射时域信号,其中矩形框中的是提取出来的有效的松质骨超声背散射信号。
图5为参考信号。
图6 为不同年龄情况下,背散射系数(BC)的分布图。
图7为不同年龄情况下,表观积分背散射系数(AIB)的分布图。
图8为不同年龄情况下,频谱质心偏移量(SCS)的分布图。
图9为不同年龄情况下,平均骨小梁间距(MTBS)的分布图。
具体实施方式
实施例:
下面以一名待测者为例,介绍整个测量过程,并在最后对本系统和传统的超声透射法进行比较。
具体步骤为:
(1)将超声探头固定于人体的脚跟位置,在松质骨超声诊断系统的控制下,向人体松质骨垂直发射超声脉冲波,超声波的传播受到松质骨微结构的作用,背向散射回来的超声信号由同一超声探头接收。
(2)将步骤(1)中接收到的信号通过并行数据总线传输到预处理模块中,进行信号预处理,再根据软组织和皮质骨的厚度去除其中的反射信号,最后通过加时域窗的方法来提取有效的背散射信号。
(3)对步骤(2)中得到的有效背散射信号,结合系统中预存的参考信号,采用软件模块分别计算背散射系数(BC),表观积分背散射系数(AIB)、频谱质心偏移量(SCS)和平均骨小梁间距(MTBS)四个参量。
(4)将步骤(3)中得到的四个参量与已经建立并储存在系统数据库中的标准值进行比较,来诊断并显示松质骨的健康状况。
步骤(1)中,使用松质骨超声诊断系统控制超声探头,在探头表面均匀涂敷医用超声耦合剂后,将探头固定在待测者的脚跟内侧平整处进行测量;测量过程中,探头向脚跟松质骨发射超声脉冲,背散射信号由同一探头接收,每一次测量经过128次平均以减小噪声干扰。因此要求在测量过程中,探头的位置要固定不动。
步骤(2)中,由超声收发探头模块接收到松质骨的超声背散射信号之后,通过并行数据总线传输到预处理模块,进行数据预处理,预处理的方法是首先去除人体软组织和皮质骨的反射信号,具体方法为首先用阈值判断法得到背散射信号的起始位置,然后根据待测者的身高和体重估计其软组织和皮质骨的厚度,再利用超声的传播速度和ADC的采样率计算背散射信号中软组织和皮质骨反射信号的部分并去除它们。最后设定感兴趣的松质骨的厚度,用同样的方法计算出信号长度,通过简单的矩形窗截取得到有效的背散射信号。
松质骨超声诊断系统使用超声背散射参量诊断松质骨健康状况的过程由图2给出:首先使用超声收发一体探头在系统的控制下获取待测者跟骨处的超声背散射信号(如图3示),然后通过预处理截取有效的背散射信号(如图3矩形框中所示),最后结合参考信号(如图4矩形框所示),计算各个背散射参量(如图5—8所示),再根据系统内部的标准值判断松质骨的健康状况。
从图1可以看出,与传统的使用两个探头的超声透射法相比,使用超声背散射法评价松质骨健康状况的方法只需要使用单一的超声探头,因此不仅节省成本,也能够测量更多的部位。从测量的背散射参量可以看出,与透射法参数只能反映松质骨骨矿密度相比,背散射法能够估计骨小梁间距,即能够反映松质骨微结构的变化,因此能够更加全面的评价松质骨的健康状况。

Claims (3)

1.一种基于超声背散射信号参量的松质骨诊断系统,其特征在于包括如下模块:
(1)超声背散射信号的获取模块:该模块通过控制超声收发一体探头,测量背散射信号;测量过程中,每隔1ms超声探头发射一组超声脉冲到松质骨并接收超声背散射信号,每次测量发射128组超声脉冲并将所有得到的信号平均,以减小噪声;  
(2)预处理模块:该模块首先通过阈值判断法来确定背散射信号的开始位置,然后根据待测者输入的体重/身高比来估计他的软组织和皮质骨的厚度,根据采样率和超声传播速度计算出背散射信号中混合的反射信号的长度,然后去除这部分反射信号;最后选取感兴趣的松质骨厚度,用同样的方法计算信号长度,最后加时域窗截取出有效的超声背散射信号;
(3)背散射参量的计算模块:该模块根据得到的有效背散射信号,再结合松质骨诊断系统中预存的参考信号,利用FFT分别得到背散射信号的幅度谱                                                
Figure 437963DEST_PATH_IMAGE002
和参考信号的幅度谱
Figure 833172DEST_PATH_IMAGE004
,再由以下公式计算4个背散射参量:
 
Figure 211326DEST_PATH_IMAGE006
                                     (1)
   
Figure 96105DEST_PATH_IMAGE008
                            (2)
 
Figure 363139DEST_PATH_IMAGE010
                                    (3)
     
Figure 2011101589700100001DEST_PATH_IMAGE011
                               (4)
其中,
Figure 2011101589700100001DEST_PATH_IMAGE013
分别是超声探头-6dB有效带宽的最小值和最大值;
Figure 919333DEST_PATH_IMAGE014
为超声探头的中心频率,
Figure 975014DEST_PATH_IMAGE016
为128次测量的背散射信号的频谱均值;
Figure 2011101589700100001DEST_PATH_IMAGE017
为背散射信号的基频;BC为背散射系数,AIB为表观积分背散射系数,SCS频谱质心偏移量,MTBS为平均骨小梁间距;
(4)松质骨状况评价模块:系统中的数据库保存有不同年龄、不同性别和不同体重/身高比下的各个背散射参量的标准值;在每次测量前,输入待测者的编号、姓名、性别、年龄、身高和体重等个人信息,该模块根据上述模块计算得到的背散射参量和对应的标准值进行比较,综合评价待测者的松质骨健康状况。
2.根据权利要求1所述的基于超声背散射信号参量的松质骨诊断系统,其特征在于所述背散射信号的基频的具体计算方法为:首先利用简易反向滤波器最小化超声换能器脉冲响应和散射响应的影响,然后利用自相关技术加强背散射信号的周期性,最后用阈值判断和平均的方法求出
Figure 528672DEST_PATH_IMAGE017
3.一种基于权利要求1所述松质骨诊断系统的松质骨诊断方法,其特征在于具体步骤为:
(1)将超声探头固定于人体的脚跟位置,在松质骨超声诊断系统的控制下,向人体松质骨垂直发射超声脉冲波,超声波的传播受到松质骨微结构的作用,背向散射回来的超声信号由同一超声探头接收;
(2)将步骤(1)中接收到的信号通过并行数据总线传输到预处理模块中,进行信号预处理,再根据软组织和皮质骨的厚度去除其中的反射信号,最后通过加时域窗的方法来提取有效的背散射信号;
    (3)对步骤(2)中得到的有效背散射信号,结合系统中预存的参考信号,采用软件模块分别计算背散射系数(BC),表观积分背散射系数(AIB)、频谱质心偏移量(SCS)和平均骨小梁间距(MTBS)四个参量;
    (4)将步骤(3)中得到的四个参量与已经建立并储存在系统数据库中的标准值进行比较,来评价并显示松质骨的健康状况。
CN2011101589700A 2011-06-14 2011-06-14 基于超声背散射信号参量的松质骨诊断系统 Pending CN102198009A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011101589700A CN102198009A (zh) 2011-06-14 2011-06-14 基于超声背散射信号参量的松质骨诊断系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011101589700A CN102198009A (zh) 2011-06-14 2011-06-14 基于超声背散射信号参量的松质骨诊断系统

Publications (1)

Publication Number Publication Date
CN102198009A true CN102198009A (zh) 2011-09-28

Family

ID=44659174

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011101589700A Pending CN102198009A (zh) 2011-06-14 2011-06-14 基于超声背散射信号参量的松质骨诊断系统

Country Status (1)

Country Link
CN (1) CN102198009A (zh)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102512203A (zh) * 2011-12-27 2012-06-27 王培勇 基于超声的人体骨密度监测仪
CN102670255A (zh) * 2012-05-10 2012-09-19 中国科学院合肥物质科学研究院 自助式超声骨强度检测系统及其使用方法
CN103565468A (zh) * 2012-07-30 2014-02-12 Ge医疗系统环球技术有限公司 用于获得骨质状况信息的超声诊断设备及其方法
CN104138277A (zh) * 2014-07-29 2014-11-12 哈尔滨工业大学 定量超声系统中基于经验模态分解的组织微观结构检测方法
CN104146729A (zh) * 2014-07-29 2014-11-19 哈尔滨工业大学 一种骨骼超声系统中编码增强的聚焦超声骨组织微观结构检测方法
CN104837411A (zh) * 2012-12-06 2015-08-12 古野电气株式会社 测量装置及测量方法
CN104997534A (zh) * 2015-07-24 2015-10-28 复旦大学 基于带约束最小二乘法的超声衰减系数成像方法
CN106022369A (zh) * 2016-05-12 2016-10-12 南京大学 一种基于超声背向散射信号成像的骨质评价方法
CN106037794A (zh) * 2016-05-12 2016-10-26 南京大学 基于频谱成像与最大似然估计的骨密度测量方法
CN106037799A (zh) * 2016-06-22 2016-10-26 华南理工大学 基于超声rf背散射信号时频分析的弹性参数成像方法
CN106264616A (zh) * 2016-09-07 2017-01-04 复旦大学 一种基于安卓平台的背散射超声骨质诊断系统
CN106777950A (zh) * 2016-12-08 2017-05-31 宇龙计算机通信科技(深圳)有限公司 一种健康监控方法及装置
CN107854141A (zh) * 2017-11-24 2018-03-30 南京大学 一种基于超声背向散射信号的骨小梁间距测量系统
CN108366782A (zh) * 2015-12-08 2018-08-03 奥林巴斯株式会社 超声波诊断装置、超声波诊断装置的工作方法以及超声波诊断装置的工作程序
US10641739B2 (en) 2014-12-23 2020-05-05 Trustees Of Boston University Method and system for oblique backscattering ultrasound transmissive contrast imaging

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101401732A (zh) * 2008-10-20 2009-04-08 中国科学院合肥物质科学研究院 超声骨密度测量分析系统

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101401732A (zh) * 2008-10-20 2009-04-08 中国科学院合肥物质科学研究院 超声骨密度测量分析系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
他得安等: "评价骨质状况及骨质疏松症的超声方法", 《中国医疗器械信息》 *
黄凯: "超声背散射法评价松质骨状况的参数估计及其成像研究", 《中国博士学位论文全文数据库医药卫生科技辑(月刊)》 *

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102512203A (zh) * 2011-12-27 2012-06-27 王培勇 基于超声的人体骨密度监测仪
CN102670255A (zh) * 2012-05-10 2012-09-19 中国科学院合肥物质科学研究院 自助式超声骨强度检测系统及其使用方法
CN102670255B (zh) * 2012-05-10 2013-12-04 中国科学院合肥物质科学研究院 自助式超声骨强度检测系统及其使用方法
CN103565468A (zh) * 2012-07-30 2014-02-12 Ge医疗系统环球技术有限公司 用于获得骨质状况信息的超声诊断设备及其方法
CN103565468B (zh) * 2012-07-30 2018-01-19 Ge医疗系统环球技术有限公司 用于获得骨质状况信息的超声诊断设备及其方法
CN104837411A (zh) * 2012-12-06 2015-08-12 古野电气株式会社 测量装置及测量方法
CN104146729A (zh) * 2014-07-29 2014-11-19 哈尔滨工业大学 一种骨骼超声系统中编码增强的聚焦超声骨组织微观结构检测方法
CN104138277B (zh) * 2014-07-29 2016-03-16 哈尔滨工业大学 定量超声系统中基于经验模态分解的组织微观结构检测方法
CN104138277A (zh) * 2014-07-29 2014-11-12 哈尔滨工业大学 定量超声系统中基于经验模态分解的组织微观结构检测方法
US10641739B2 (en) 2014-12-23 2020-05-05 Trustees Of Boston University Method and system for oblique backscattering ultrasound transmissive contrast imaging
CN104997534A (zh) * 2015-07-24 2015-10-28 复旦大学 基于带约束最小二乘法的超声衰减系数成像方法
CN104997534B (zh) * 2015-07-24 2017-10-10 复旦大学 基于带约束最小二乘法的超声衰减系数成像方法
CN108366782A (zh) * 2015-12-08 2018-08-03 奥林巴斯株式会社 超声波诊断装置、超声波诊断装置的工作方法以及超声波诊断装置的工作程序
CN106037794A (zh) * 2016-05-12 2016-10-26 南京大学 基于频谱成像与最大似然估计的骨密度测量方法
CN106022369A (zh) * 2016-05-12 2016-10-12 南京大学 一种基于超声背向散射信号成像的骨质评价方法
CN106037794B (zh) * 2016-05-12 2019-04-19 南京大学 基于频谱成像与最大似然估计的骨密度测量方法
CN106037799A (zh) * 2016-06-22 2016-10-26 华南理工大学 基于超声rf背散射信号时频分析的弹性参数成像方法
CN106037799B (zh) * 2016-06-22 2019-04-09 华南理工大学 基于超声rf背散射信号时频分析的弹性参数成像方法
CN106264616A (zh) * 2016-09-07 2017-01-04 复旦大学 一种基于安卓平台的背散射超声骨质诊断系统
CN106264616B (zh) * 2016-09-07 2023-09-01 复旦大学 一种基于安卓平台的背散射超声骨质诊断系统
CN106777950A (zh) * 2016-12-08 2017-05-31 宇龙计算机通信科技(深圳)有限公司 一种健康监控方法及装置
CN107854141A (zh) * 2017-11-24 2018-03-30 南京大学 一种基于超声背向散射信号的骨小梁间距测量系统

Similar Documents

Publication Publication Date Title
CN102198009A (zh) 基于超声背散射信号参量的松质骨诊断系统
US7535797B2 (en) High-resolution ultrasound displacement measurement apparatus and method
CN100548224C (zh) 利用超声检测内部狭窄以识别由其引起的组织振动的设备
CN101883526B (zh) 使用超声技术测量材料的厚度的方法
US6468215B1 (en) Method and device for multi-parametric ultrasonic assessment of bone conditions
Bigelow et al. Estimation of total attenuation and scatterer size from backscattered ultrasound waveforms
KR20180031711A (ko) 급성 중이염의 특성화를 위한 장치 및 방법
KR20150037689A (ko) 의료용 초음파 영상의 전단파 검출
JP2519777B2 (ja) 超音波装置、並びに超音波測定及び超音波検査を行う方法
JP4297699B2 (ja) スペクトル歪み度を描出するための方法及び装置
Deeba et al. Attenuation coefficient estimation of normal placentas
Gerig et al. Statistics of ultrasonic scatterer size estimation with a reference phantom
CN104116524B (zh) 一种超声衰减系数补偿系统及肝脏脂肪检测系统
EP2454997B1 (en) Method, device and program for analyzing cartilage using ultrasonic wave
WO2019218669A1 (zh) 一种剪切波传播速度的确定方法及装置
US20110046517A1 (en) Method and system for monitoring skeletal defects
US20190231320A1 (en) Ultrasonic shear wave imaging with background motion compensation
Omari et al. Signal to noise ratio comparisons for ultrasound attenuation slope estimation algorithms
JP2021013575A (ja) 超音波診断装置、および、超音波診断装置の制御方法
Baldwin et al. Estimating myocardial attenuation from M-mode ultrasonic backscatter
CN112075955B (zh) 粘弹性介质超声波参数测量方法与装置
Srinivasan et al. Estimating the elastographic signal-to-noise ratio using correlation coefficients
KR101840349B1 (ko) 초음파 합주파수 성분을 이용한 골밀도 예측 장치 및 방법
CN112799078B (zh) 一种剪切波传播速度的检测方法、系统和超声成像设备
CN1846628B (zh) 用时间极值法测定超声波在骨骼中传播速度的方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20110928