CN110179498A - 一种用于淋巴水肿显像的热声成像系统及其控制方法 - Google Patents

一种用于淋巴水肿显像的热声成像系统及其控制方法 Download PDF

Info

Publication number
CN110179498A
CN110179498A CN201910589582.4A CN201910589582A CN110179498A CN 110179498 A CN110179498 A CN 110179498A CN 201910589582 A CN201910589582 A CN 201910589582A CN 110179498 A CN110179498 A CN 110179498A
Authority
CN
China
Prior art keywords
thermoacoustic
signal
lymphedema
image
pulse
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
CN201910589582.4A
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.)
West China Hospital of Sichuan University
Original Assignee
West China Hospital of Sichuan 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 West China Hospital of Sichuan University filed Critical West China Hospital of Sichuan University
Priority to CN201910589582.4A priority Critical patent/CN110179498A/zh
Publication of CN110179498A publication Critical patent/CN110179498A/zh
Pending legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0093Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明公开了一种用于淋巴水肿显像的热声成像系统及其控制方法,属于医疗领域。该用于淋巴水肿显像的热声成像系统包括:脉冲微波源,用于激发淋巴水肿区域的热声效应,产生热声信号;超声换能器,用于接收热声信号;热声信号处理单元包括分别用于对热声信号进行放大和滤波的放大器和滤波器;数据采集单元,用于对经过放大和滤波的热声信号进行A/D转换,得到数字化的热声信号;计算机,用于控制脉冲微波源,以及对数字化热声信号进行图像重建得到热声图像,并对热声图像进行颜色编码实现淋巴水肿的突出显像。

Description

一种用于淋巴水肿显像的热声成像系统及其控制方法
技术领域
本发明涉及医疗领域,具体涉及一种用于淋巴水肿显像的热声成像系统及其控制方法。
背景技术
淋巴水肿是指因淋巴液回流障碍所引起的组织肿胀,主要是淋巴管阻塞和发育异常。常发生于小腿、上臂、生殖器和面部等处,淋巴水肿的临床表现:早期用手指按压时可呈现明显的凹陷性压窝,晚期皮肤和皮下组织增生,皮皱加深,皮肤增厚变硬粗糙,外观似大象皮肤。早期抬高患处或卧床休息后肿胀可以消失或减轻。晚期患肢肿大明显,表面角化粗糙,呈橡皮样肿。少数可有皮肤裂开、溃疡或出现疣状赘生物。
关于淋巴水肿的测量,一般采用直接测量的方式,这种测量方式具有一定的主观性,并且测量精度较低。
发明内容
针对现有技术中的上述不足,本发明旨在提供一种对淋巴水肿的测量精度更高的用于淋巴水肿显像的热声成像系统及其控制方法。
为了达到上述发明创造的目的,本发明采用的技术方案为:
提供一种用于淋巴水肿显像的热声成像系统,其包括:
脉冲微波源,用于激发淋巴水肿区域的热声效应,产生热声信号;
超声换能器,用于接收热声信号;
热声信号处理单元,包括分别用于对热声信号进行放大和滤波的放大器和滤波器;
数据采集单元,用于对经过放大和滤波的热声信号进行A/D转换,得到数字化的热声信号;
计算机,用于控制脉冲微波源,以及对数字化热声信号进行图像重建得到热声图像,并对热声图像进行颜色编码实现淋巴水肿的突出显像。
进一步地,图像重建的方法包括:有限元定量重建算法、滤波反投影或延迟叠加重建算法。
进一步地,图像重建的方法包括:
S1、对数字化的热声信号进行离散化处理,然后将数字化的热声信号传播延迟时间和空间距离按照声波在介质中的传播速度转变为网格点数;
S2、在超声换能器的探头扫描轨迹圆所在平面内,以传播延迟时间对应的网格点数代表的距离为半径,晶元中心所在位置为圆心画圆;
S3、将位于以晶元中心为圆心的圆的圆弧上且在超声换能器的孔径角以内的网格点赋予对应传播延迟时间点的信号幅值;
S4、在探头扫描轨迹圆截面内,将距探头扫描轨迹圆圆心不同距离处被赋予了信号幅值的网格点的信号强度进行延迟叠加;
S5、将延迟叠加后的结果用图像表示得到热声图像。
进一步地,信号幅值为对应晶元上接收到的信号强度,其计算公式为:
其中,A(i,j)为坐标为的晶元上接收到的信号强度;B(k,l)为第l点吸收体到达第k个晶元的信号强度;M为超声换能器的晶元个数。
进一步地,超声换能器为多元阵列超声探头,换能通道数为2~2048,频率范围为5.0~10.0MHz,形状包括线阵和相控阵。
进一步地,脉冲微波源包括依次连接的脉冲微波发生器、传输线缆和辐射器。
进一步地,脉冲微波发生器产生的脉冲微波信号的中心频率为0.5~10.0GHz,脉冲宽度为10~500ns,脉冲峰值功率为30~100kW,脉冲重复频率为1~200Hz。
进一步地,辐射器为高功率增益天线。
另一方面,提供一种基于本方案设计的热声成像系统的控制方法,其包括:
对脉冲微波源进行参数设置和初始化,并开启脉冲微波源进行预热;
触发脉冲微波源发射脉冲微波,脉冲微波源将脉冲微波辐射到淋巴水肿区域,使淋巴水肿区域与未水肿区域吸收脉冲微波能量后产生热声信号;
通过超能换声器接受热声信号,并通过热声信号处理单元对热声信号进行放大和滤波,通过数据采集单元对经过放大和滤波的热声信号进行A/D转换,得到数字化的热声信号;
通过计算机并对数字化热声信号进行图像重建得到热声图像,然后对热声图像进行颜色编码实现淋巴水肿的突出显像。
本发明的有益效果为:
利用淋巴水肿具有较高的微波吸收的特性,以及热声成像具有的微波吸收高对比度和超声成像的高分辨率特点,结合图像重建得到淋巴水肿区域的热声图像,并对热声图像进行颜色编码实现淋巴水肿的突出显像,进而使得显像结果较为准确,精度较高,从而有利于后期的准确分析。
整个热声成像过程具有无创和快速的特点。并且该方法不仅可以单独使用,也可配合超声成像、光声成像形成多模态淋巴水肿成像技术使用,用于对患者异常部位的显像。
附图说明
图1为具体实施例中热声成像系统的结构图;
图2为一种超声换能器的局部结构示意图;
图3为具体实施例中图像重建算法的原理图;
图4为图1所示的热声成像系统对人体手臂进行热声成像的结果。
其中,1、晶元;2、探头扫描轨迹圆;3、吸收体;4、目标物体。
具体实施方式
下面结合附图,对本发明的具体实施方式做详细说明,以便于本技术领域的技术人员理解本发明。但应该清楚,下文所描述的实施例仅仅是本发明的一部分实施例,而不是全部实施例。在不脱离所附的权利要求限定和确定的本发明的精神和范围内,本领域普通技术人员在没有做出任何创造性劳动所获得的所有其他实施例,都属于本发明的保护范围。
如图1所示,该用于淋巴水肿显像的热声成像系统包括:
脉冲微波源,用于激发淋巴水肿区域的热声效应,产生热声信号;
超声换能器,用于接收热声信号;
热声信号处理单元,包括分别用于对热声信号进行放大和滤波的放大器和滤波器;
数据采集单元,用于对经过放大和滤波的热声信号进行A/D转换,得到数字化的热声信号;
计算机,用于控制脉冲微波源,以及对数字化热声信号进行图像重建得到热声图像,并对热声图像进行颜色编码实现淋巴水肿的突出显像。
实施时,如图2所示,由于超声换能器自身的特点,使得其只能接收一定孔径角以内的信号。因此,传统的图像重建算法会出现图像伪影,影响淋巴水肿显像的准确性。对此,为提高淋巴水肿显像的准确性,如图3所示,其中箭头所示方向为探头扫描轨迹方向,目标物体4位于探头扫描轨迹圆2截面内,本方案优选图像重建的方法包括:
S1、对数字化的热声信号进行离散化处理,然后将数字化的热声信号传播延迟时间和空间距离按照声波在介质中的传播速度转变为网格点数;
对于第m个探头采集到的数字化的热声信号可以表示为Pm(1*L)的一维矩阵,其中L表示总的采样点数。对于网格大小为Δd的情况,Pm(1*L)对应的传播空间距离对应的网格位置点(也即网格点数)Pgrid为:
上式中Vs为组织中超声速度,SR为系统采样率。
Pm(1*L)对应的传播延迟时间对应的网格点数为
S2、在超声换能器的探头扫描轨迹圆2所在平面内,以传播延迟时间对应的网格点数代表的距离为半径,晶元1中心所在位置为圆心画圆;
S3、将位于以晶元1中心为圆心的圆的圆弧上且在超声换能器的孔径角以内的网格点赋予对应传播延迟时间点的信号幅值;
其中,孔径角可以通过以下计算公式得到:Sin(α/2)=0.514λ/D,其中,α为孔径角,λ为超声换能器中心频率对应的波长,D为晶元1直径(对于矩形超声换能器,D为成像平面所在边的长度)。
信号幅值为对应晶元1上接收到的信号强度,其计算公式为:
其中,A(i,j)为坐标为i,j的晶元1上接收到的信号强度;B(k,l)为第l点吸收体3到达第k个晶元1的信号强度;M为超声换能器的晶元个数,
S4、在探头扫描轨迹圆2截面内,将距探头扫描轨迹圆2圆心不同距离处被赋予了信号幅值的网格点的信号强度进行延迟叠加;
S5、将延迟叠加后的结果用图像表示得到热声图像。
图4为采用上述图像重建算法对人体手臂进行热声成像的结果,其中箭头所示为手臂皮肤,虚线框所示为血管。
其中,超声换能器为多元阵列超声探头,换能通道数为2~2048,频率范围为5.0~10.0MHz,形状包括线阵和相控阵。
其中,脉冲微波源包括依次连接的脉冲微波发生器、传输线缆和辐射器。
其中,脉冲微波发生器产生的脉冲微波信号的中心频率为0.5~10.0GHz,脉冲宽度为10~500ns,脉冲峰值功率为30~100kW,脉冲重复频率为1~200Hz。
其中,辐射器为高功率增益天线,其包括开口波导天线、喇叭天线、偶极子和贴片天线等。
在另一实施例中,图像重建的方法包括:有限元定量重建算法、滤波反投影或延迟叠加重建算法。
另一方面,本方案还提供一种基于本方案设计的热声成像系统的控制方法,其包括:
对脉冲微波源进行参数设置和初始化,并开启脉冲微波源进行预热;
触发脉冲微波源发射脉冲微波,脉冲微波源将脉冲微波辐射到淋巴水肿区域,使淋巴水肿区域与未水肿区域吸收脉冲微波能量后产生热声信号;
通过超能换声器接受热声信号,并通过热声信号处理单元对热声信号进行放大和滤波,通过数据采集单元对经过放大和滤波的热声信号进行A/D转换,得到数字化的热声信号;
通过计算机并对数字化热声信号进行图像重建得到热声图像,然后对热声图像进行颜色编码实现淋巴水肿的突出显像。

Claims (9)

1.一种用于淋巴水肿显像的热声成像系统,其特征在于,包括:
脉冲微波源,用于激发淋巴水肿区域的热声效应,产生热声信号;
超声换能器,用于接收热声信号;
热声信号处理单元,包括分别用于对热声信号进行放大和滤波的放大器和滤波器;
数据采集单元,用于对经过放大和滤波的热声信号进行A/D转换,得到数字化的热声信号;
计算机,用于控制脉冲微波源,以及对数字化热声信号进行图像重建得到热声图像,并对热声图像进行颜色编码实现淋巴水肿的突出显像。
2.根据权利要求1所述的热声成像系统,其特征在于,所述图像重建的方法包括:有限元定量重建算法、滤波反投影或延迟叠加重建算法。
3.根据权利要求1所述的热声成像系统,其特征在于,所述图像重建的方法包括:
S1、对数字化的热声信号进行离散化处理,然后将数字化的热声信号传播延迟时间和空间距离按照声波在介质中的传播速度转变为网格点数;
S2、在超声换能器的探头扫描轨迹圆(2)所在平面内,以传播延迟时间对应的网格点数代表的距离为半径,晶元(1)中心所在位置为圆心画圆;
S3、将位于以晶元(1)中心为圆心的圆的圆弧上且在超声换能器的孔径角以内的网格点赋予对应传播延迟时间点的信号幅值;
S4、在所述探头扫描轨迹圆(2)截面内,将距探头扫描轨迹圆(2)圆心不同距离处被赋予了信号幅值的网格点的信号强度进行延迟叠加;
S5、将延迟叠加后的结果用图像表示得到所述热声图像。
4.根据权利要求3所述的热声成像方法,其特征在于,所述信号幅值为对应晶元(1)上接收到的信号强度,其计算公式为:
其中,A(i,j)为坐标为(i,j)的晶元(1)上接收到的信号强度;B(k,l)为第l点吸收体(3)到达第k个晶元(1)的信号强度;M为超声换能器的晶元个数。
5.根据权利要求1所述的热声成像系统,其特征在于,所述超声换能器为多元阵列超声探头,换能通道数为2~2048,频率范围为5.0~10.0MHz,形状包括线阵和相控阵。
6.根据权利要求1-5任一所述的热声成像系统,其特征在于,所述脉冲微波源包括依次连接的脉冲微波发生器、传输线缆和辐射器。
7.根据权利要求6所述的热声成像系统,其特征在于,所述脉冲微波发生器产生的脉冲微波信号的中心频率为0.5~10.0GHz,脉冲宽度为10~500ns,脉冲峰值功率为30~100kW,脉冲重复频率为1~200Hz。
8.根据权利要求6所述的热声成像系统,其特征在于,所述辐射器为高功率增益天线。
9.一种基于权利要求1-8任一热声成像系统的控制方法,其特征在于,包括:
对脉冲微波源进行参数设置和初始化,并开启脉冲微波源进行预热;
触发脉冲微波源发射脉冲微波,脉冲微波源将脉冲微波辐射到淋巴水肿区域,使淋巴水肿区域与未水肿区域吸收脉冲微波能量后产生热声信号;
通过超能换声器接受热声信号,并通过热声信号处理单元对热声信号进行放大和滤波,通过数据采集单元对经过放大和滤波的热声信号进行A/D转换,得到数字化的热声信号;
通过计算机并对数字化热声信号进行图像重建得到热声图像,然后对热声图像进行颜色编码实现淋巴水肿的突出显像。
CN201910589582.4A 2019-07-02 2019-07-02 一种用于淋巴水肿显像的热声成像系统及其控制方法 Pending CN110179498A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910589582.4A CN110179498A (zh) 2019-07-02 2019-07-02 一种用于淋巴水肿显像的热声成像系统及其控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910589582.4A CN110179498A (zh) 2019-07-02 2019-07-02 一种用于淋巴水肿显像的热声成像系统及其控制方法

Publications (1)

Publication Number Publication Date
CN110179498A true CN110179498A (zh) 2019-08-30

Family

ID=67724605

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910589582.4A Pending CN110179498A (zh) 2019-07-02 2019-07-02 一种用于淋巴水肿显像的热声成像系统及其控制方法

Country Status (1)

Country Link
CN (1) CN110179498A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110559016A (zh) * 2019-10-21 2019-12-13 四川瑞利浦医疗科技有限公司 基于微波热声和超声成像的甲状腺检测装置及方法
CN110742588A (zh) * 2019-10-28 2020-02-04 四川大学华西医院 微波热声、光声和超声三模态肠道组织成像方法及系统
CN112782279A (zh) * 2019-11-07 2021-05-11 华南师范大学 非接触式热声成像方法与装置
CN112773351A (zh) * 2019-11-07 2021-05-11 华南师范大学 一种基于无油耦合与阻抗匹配紧凑型探头的热声成像装置与方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015162215A1 (en) * 2014-04-25 2015-10-29 Helmholtz Zentrum München Deutsches Forschungszentrum Für Gesundheit Und Umwelt (Gmbh) Device and method for frequency-domain thermoacoustic sensing
CN105259426A (zh) * 2014-07-18 2016-01-20 中国科学院沈阳自动化研究所 一种热声效应辐射场空间分布测量装置及方法
US20190038220A1 (en) * 2017-08-01 2019-02-07 Endra Life Sciences Inc. Method and system for determining fractional fat content of tissue
CN109480777A (zh) * 2018-11-05 2019-03-19 四川华微康医疗科技有限公司 一种用于热消融边界显像的热声成像系统及方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015162215A1 (en) * 2014-04-25 2015-10-29 Helmholtz Zentrum München Deutsches Forschungszentrum Für Gesundheit Und Umwelt (Gmbh) Device and method for frequency-domain thermoacoustic sensing
CN105259426A (zh) * 2014-07-18 2016-01-20 中国科学院沈阳自动化研究所 一种热声效应辐射场空间分布测量装置及方法
US20190038220A1 (en) * 2017-08-01 2019-02-07 Endra Life Sciences Inc. Method and system for determining fractional fat content of tissue
CN109480777A (zh) * 2018-11-05 2019-03-19 四川华微康医疗科技有限公司 一种用于热消融边界显像的热声成像系统及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
黄林,刘立新,鲁奎等: "3GHz微波热声成像系统", 《西南交通大学学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110559016A (zh) * 2019-10-21 2019-12-13 四川瑞利浦医疗科技有限公司 基于微波热声和超声成像的甲状腺检测装置及方法
CN110742588A (zh) * 2019-10-28 2020-02-04 四川大学华西医院 微波热声、光声和超声三模态肠道组织成像方法及系统
CN112782279A (zh) * 2019-11-07 2021-05-11 华南师范大学 非接触式热声成像方法与装置
CN112773351A (zh) * 2019-11-07 2021-05-11 华南师范大学 一种基于无油耦合与阻抗匹配紧凑型探头的热声成像装置与方法
CN112782279B (zh) * 2019-11-07 2023-04-07 华南师范大学 非接触式热声成像方法与装置
CN112773351B (zh) * 2019-11-07 2024-03-22 华南师范大学 一种基于无油耦合与阻抗匹配紧凑型探头的热声成像装置与方法

Similar Documents

Publication Publication Date Title
CN110179498A (zh) 一种用于淋巴水肿显像的热声成像系统及其控制方法
KR101651830B1 (ko) 고강도 집속된 초음파를 위한 의료용 초음파 영상화에서의 피드백
CN104254773B (zh) 使用二维阵列数据的移动超声诊断探测设备以及使用该设备的移动超声诊断系统
Oraevsky et al. Two-dimensional optoacoustic tomography: transducer array and image reconstruction algorithm
US8905931B2 (en) Subject information processing apparatus
CN107569256B (zh) 基于热膨胀和门控算法测量生物组织温度变化的超声方法
CN104739448B (zh) 一种超声成像方法及装置
Ding et al. Microwave-excited ultrasound and thermoacoustic dual imaging
CN103202713B (zh) 一种医用超声基波和谐波融合的图像优化方法
CN102058416B (zh) 基于压缩感知的微波热声成像装置及方法
CN109876310A (zh) 质子治疗的监控方法、装置与系统
CN106846458A (zh) 基于3d打印的超声波立体模型构建方法及装置
CN112451865A (zh) 超声治疗监测装置及方法
KR20150112057A (ko) 공간 일관성 기초 초음파 신호 처리 모듈 및 그에 의한 초음파 신호 처리 방법
CN106264722A (zh) 窗口调制结合高斯多项式拟合监测射频消融设备及方法
CN106940883B (zh) 基于超声系统点扩散函数仿真和压缩感知的超声成像方法
Xu et al. Signal processing in scanning thermoacoustic tomography in biological tissues
Chen et al. The prototype of microwave-induced thermo-acoustic tomography imaging by time reversal mirror
CN207444961U (zh) 一种超声波检测治疗仪
Park et al. Photoacoustic imaging using array transducer
Kruger et al. Thermoacoustic CT scanner for breast imaging: Design considerations
CN105342567A (zh) 一种提高重建光声图像信噪比的装置和方法
Jiang et al. Plane wave imaging combined with eigenspace-based minimum variance beamforming using a ring array in ultrasound computed tomography
CN106997045B (zh) 基于超声系统点扩散函数测量和压缩感知的超声成像方法
Mo et al. Compact ultrasound scanner with simultaneous parallel channel data acquisition capabilities

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20190830

RJ01 Rejection of invention patent application after publication