WO2017107660A1 - 粘弹性介质的粘弹性参数检测方法和设备 - Google Patents
粘弹性介质的粘弹性参数检测方法和设备 Download PDFInfo
- Publication number
- WO2017107660A1 WO2017107660A1 PCT/CN2016/103645 CN2016103645W WO2017107660A1 WO 2017107660 A1 WO2017107660 A1 WO 2017107660A1 CN 2016103645 W CN2016103645 W CN 2016103645W WO 2017107660 A1 WO2017107660 A1 WO 2017107660A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- viscoelastic
- maximum displacement
- viscoelastic medium
- parameter
- displacement data
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000006073 displacement reaction Methods 0.000 claims abstract description 99
- 206010016654 Fibrosis Diseases 0.000 claims abstract description 12
- 230000004761 fibrosis Effects 0.000 claims abstract description 12
- 230000010355 oscillation Effects 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims description 15
- 230000002159 abnormal effect Effects 0.000 claims description 10
- 239000000523 sample Substances 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000003595 spectral effect Effects 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 3
- 230000001902 propagating effect Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 9
- 210000001519 tissue Anatomy 0.000 description 20
- 210000005228 liver tissue Anatomy 0.000 description 13
- 208000019425 cirrhosis of liver Diseases 0.000 description 10
- 238000001514 detection method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 210000004185 liver Anatomy 0.000 description 3
- 230000007882 cirrhosis Effects 0.000 description 2
- 238000013399 early diagnosis Methods 0.000 description 2
- 208000019423 liver disease Diseases 0.000 description 2
- 102000010834 Extracellular Matrix Proteins Human genes 0.000 description 1
- 108010037362 Extracellular Matrix Proteins Proteins 0.000 description 1
- 206010019668 Hepatic fibrosis Diseases 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000002091 elastography Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000002099 shear wave elastography Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
- A61B8/085—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating body or organic structures, e.g. tumours, calculi, blood vessels, nodules
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0048—Detecting, measuring or recording by applying mechanical forces or stimuli
- A61B5/0051—Detecting, measuring or recording by applying mechanical forces or stimuli by applying vibrations
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/42—Detecting, measuring or recording for evaluating the gastrointestinal, the endocrine or the exocrine systems
- A61B5/4222—Evaluating particular parts, e.g. particular organs
- A61B5/4244—Evaluating particular parts, e.g. particular organs liver
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
- G01S15/8911—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a single transducer for transmission and reception
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/899—Combination of imaging systems with ancillary equipment
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52079—Constructional features
Definitions
- the invention belongs to the technical field of medical technology, and in particular relates to a method and a device for detecting viscoelastic parameters of a viscoelastic medium.
- Hepatic fibrosis is usually caused by the deposition of excessive extracellular matrix proteins and is frequently found in most types of patients with chronic liver disease. Since early liver fibrosis or cirrhosis is reversible or controllable, accurate and effective early diagnosis of liver fibrosis is critical.
- Shear wave elastography can quantitatively assess the extent of liver fibrosis and cirrhosis by measuring liver stiffness values.
- the most widely used clinical non-invasive liver fibrosis grading detection is instantaneous elastography.
- liver viscoelastic parameters can provide very valuable information for the early diagnosis of liver fibrosis.
- the measurement of the tissue is mainly to measure the elastic parameters of the tissue, while ignoring the viscosity parameter neglects the viscoelastic parameters of the tissue, which adversely affects the early detection results of tissues such as liver fibrosis.
- the present invention provides a viscoelastic parameter detecting method and apparatus for a viscoelastic medium for obtaining viscoelastic parameters of a tissue to improve the accuracy of the measurement result of the degree of fibrosis.
- the invention provides a viscoelastic parameter detecting method for a viscoelastic medium, comprising:
- the invention provides a viscoelastic parameter detecting device for a viscoelastic medium, comprising:
- Controlling a host and a probe the probe including a vibrator and an ultrasonic transducer;
- the vibrator applies a mechanical vibration of a single predetermined frequency to the viscoelastic medium under the control of the control host, and generates a shear wave in the viscoelastic medium;
- the ultrasonic transducer emits a single-vibration source ultrasonic wave to the viscoelastic medium and receives an ultrasonic echo signal under the control of the control host during the propagation of the shear wave in the viscoelastic medium;
- the control host includes:
- a first acquiring module configured to acquire, according to the ultrasonic echo signal, maximum displacement data of the shear wave at each depth, each of the maximum displacement data characterizing that the shear wave propagates to the viscoelastic medium The maximum oscillation amplitude of the shear wave at different depths;
- a first determining module configured to determine a viscoelastic parameter of the viscoelastic medium according to the maximum displacement attenuation curve.
- the viscoelastic parameter detecting device of the viscoelastic medium applies a mechanical vibration of a predetermined frequency to the tissue to generate a shear wave of a frequency only in the viscoelastic medium. After acquiring the displacement data of the shear wave, calculating, according to the displacement data, maximum displacement data representing a maximum oscillation amplitude when the shear wave propagates to different depths, and obtaining the shear wave by fitting each maximum displacement data.
- the maximum displacement attenuation curve is used to determine the viscoelastic parameter of the viscoelastic medium based on the maximum displacement attenuation curve, which is related to both elasticity and viscosity.
- FIG. 1 is a flow chart of a first embodiment of a viscoelastic parameter detecting method for a viscoelastic medium according to the present invention
- FIG. 2 is a schematic diagram of displacement data when a shear wave propagates to a certain depth
- FIG. 3 is a flow chart of a second embodiment of a viscoelastic parameter detecting method for a viscoelastic medium according to the present invention.
- FIG. 4 is a schematic view showing a first embodiment of a viscoelastic parameter detecting device for a viscoelastic medium according to the present invention
- FIG. 5 is a schematic diagram of Embodiment 2 of a viscoelastic parameter detecting device for a viscoelastic medium according to the present invention.
- FIG. 1 is a flow chart of a first embodiment of a viscoelastic parameter detecting method for a viscoelastic medium according to the present invention.
- the method provided in this embodiment is mainly used for detecting viscoelastic parameters of liver tissue, and can be performed by a detecting device.
- the device may be an existing non-invasive liver fibrosis detector, but the processing functions required to perform the methods described in this embodiment are added to the non-invasive liver fibrosis detector.
- the detection device mainly includes a control host and a probe, and the probe includes a vibrator for generating mechanical vibration and an ultrasonic transducer for transmitting and receiving ultrasonic waves.
- the viscoelastic parameter detecting method of the viscoelastic medium may include the following steps:
- Step 101 Applying a mechanical vibration of a single predetermined frequency to the viscoelastic medium to generate a shear wave in the viscoelastic medium.
- the liver tissue is the viscoelastic medium described above.
- Applying a mechanical vibration of a single predetermined frequency to the viscoelastic medium means applying the mechanical vibration to the surface of the skin corresponding to the liver tissue.
- the vibrator exerts a sinusoidal mechanical vibration perpendicular to the surface of the skin on the surface of the skin, thereby producing a corresponding shear wave in the liver tissue, which is propagated in the liver tissue.
- the frequency of the mechanical vibration may be, for example, a low frequency such as 50 Hz.
- Step 102 During the propagation of the shear wave in the viscoelastic medium, the single-vibration source ultrasonic wave is transmitted to the viscoelastic medium, and the ultrasonic echo signal is received.
- the ultrasonic transducer emits an ultrasonic signal of a low-frequency single-vibration source to the liver tissue at a position where the vibration is applied by the vibrator, and receives the ultrasonic echo signal.
- multiple frames of ultrasound signals can be transmitted to the liver tissue at a certain time interval to track the propagation of shear waves in the liver tissue.
- the shear wave is being viscoelastic
- a series of ultrasonic signals are transmitted through the transmitter of the single-vibration source ultrasonic module integrated on the vibration probe and the ultrasonic echo signals are received, and the ultrasonic echo signal data of the time period can be processed.
- the media strain and displacement data information on the ultrasonic scanning line during this period is obtained.
- Step 103 Acquire maximum displacement data of the shear wave at each depth according to the ultrasonic echo signal.
- each maximum displacement data characterizes the maximum oscillation amplitude of the shear wave when the shear wave propagates to different depths in the viscoelastic medium.
- the ultrasonic echo signal can reflect the displacement of the shear wave in the liver tissue. Therefore, the displacement data of the shear wave can be obtained from the ultrasonic echo signal.
- a certain digital signal processing can be performed on the ultrasonic echo signal. Signal processing includes at least one of the following signal processing: time domain cross-correlation, spectral cross-correlation, squared error sum, speckle tracking, scale-invariant feature point tracking, dynamic programming, zero-crossing tracking, peak search.
- FIG. 2 illustrates a displacement result curve of a shear wave generated by mechanical vibration of the external probe as a function of time at a fixed depth in the tissue.
- mechanical vibration perpendicular to the liver tissue is applied to the liver tissue, and the ultrasonic transducer captures a displacement perpendicular to the axis of the liver tissue at which the mechanical vibration is applied, that is, a longitudinal displacement.
- the DAV in Figure 2 represents the longitudinal displacement.
- the displacement data exhibits oscillation attenuation characteristics at a fixed depth. Generally, the maximum displacement appears at the first peak. Therefore, for the displacement data corresponding to each depth obtained, it can be The maximum displacement data is extracted to obtain the maximum displacement data at different depths.
- Step 104 Fit each maximum displacement data to obtain a maximum displacement attenuation curve.
- Step 105 Determine a viscoelastic parameter of the viscoelastic medium according to a maximum displacement attenuation curve.
- the maximum The displacement data performs a certain amount of data processing.
- the maximum displacement data may be filtered in the time domain and the frequency domain, and the abnormal data in each maximum displacement data is excluded, and the abnormal data includes an average displacement value of the displacement value greater than all the maximum displacement data or an average displacement value.
- the displacement data of the multiple, or the abnormal data includes displacement data whose difference between the displacement value and the average displacement value is greater than a standard deviation of a certain multiple.
- the fitting formula is as follows:
- the fitting results can obtain three parameters a, b, and c. Since the parameters b and c affect the position of the quadratic polynomial curve and are not related to the curve attenuation trend and shape, the parameter a can be extracted to describe the attenuation trend and shape of the maximum displacement attenuation curve. This coefficient is determined by viscosity and elasticity. Co-determined, the viscoelastic parameters. That is, the coefficient corresponding to the highest power variable of the maximum displacement attenuation curve is determined as the viscoelastic parameter of the viscoelastic medium.
- the low-frequency vibration of a single frequency is used, and the viscoelastic parameter of the measured tissue can be obtained by analyzing the vibration amplitude of the shear wave.
- the specific principle is that the vibration amplitude is not only related to the elastic parameter, but also related to the viscosity, that is, related to the viscoelastic parameter, which can be described from the peak and amplitude attenuation at a specific depth.
- the drop curve formed by the decrease in peak value with the depth of propagation is affected by elasticity and viscosity. The higher the viscosity, the smaller the first wave peak is in the shallower tissue, and the deeper the depth is, the lower the viscosity of the tissue is, the more the viscous tissue is larger, and the lower the peak is more severe. .
- a mechanical vibration of a predetermined frequency is applied to the tissue to produce a shear wave of a frequency only in the viscoelastic medium.
- calculating, according to the displacement data, maximum displacement data representing a maximum oscillation amplitude when the shear wave propagates to different depths and obtaining the shear wave by fitting each maximum displacement data.
- the maximum displacement attenuation curve is used to determine the viscoelastic parameter of the viscoelastic medium based on the maximum displacement attenuation curve, which is related to both elasticity and viscosity.
- FIG. 3 is a flow chart of a second embodiment of a viscoelastic parameter detecting method for a viscoelastic medium according to the present invention. As shown in FIG. 3, on the basis of the embodiment shown in FIG. 1, after step 105, the following steps may be further included:
- Step 201 Acquire an elastic parameter of the viscoelastic medium according to the ultrasonic echo signal.
- Step 202 Determine the degree of fibrosis of the viscoelastic medium according to the elastic parameter and the viscoelastic parameter.
- the elastic parameter of the viscoelastic medium can be obtained based on the analysis processing of the received ultrasonic echo signal based on the method in the prior art.
- the degree of fibrosis of the tissue is collectively determined.
- the degree of fibrosis of tissue is generally classified as severe, general, and not severe, and each degree corresponds to a different range of elastic parameters.
- the division of the degree of fibrosis and the accuracy of the determination of the degree of fibrosis are provided to provide a usable data dimension.
- FIG. 4 is a schematic diagram of Embodiment 1 of a viscoelastic parameter detecting device for a viscoelastic medium according to the present invention. As shown in FIG. 4, the detecting device includes:
- the host 1 and the probe 2 are controlled, and the probe includes a vibrator 21 and an ultrasonic transducer 22.
- the vibrator 21 applies a mechanical vibration of a single predetermined frequency to the viscoelastic medium under the control of the control host 1, and generates shear waves in the viscoelastic medium.
- the ultrasonic transducer 22 under the control of the control host 1, emits a single-vibration source ultrasonic wave to the viscoelastic medium and receives an ultrasonic echo signal during the propagation of the shear wave in the viscoelastic medium.
- the control host 1 includes a first acquisition module 11, a calculation module 12, and a first determination module 13.
- a first obtaining module 11 configured to acquire maximum displacement data of the shear wave at each depth according to the ultrasonic echo signal, and each of the maximum displacement data characterizes the shear wave propagation The maximum amplitude of oscillation of the shear wave to different depths in the viscoelastic medium.
- the calculation module 12 is configured to fit each maximum displacement data to obtain a maximum displacement attenuation curve.
- the first determining module 13 is configured to determine a viscoelastic parameter of the viscoelastic medium according to the maximum displacement attenuation curve.
- control host further includes: a processing module 14.
- the processing module 14 is configured to perform at least one of signal processing on the ultrasonic echo signals: time domain cross correlation, spectral cross correlation, square error sum, speckle tracking, scale invariant feature point tracking, dynamic programming , zero cross tracking, peak search.
- the calculation module 12 includes: a first calculation unit 121 and a second calculation unit 122.
- the first calculating unit 121 is configured to perform filtering processing in the time domain and the frequency domain on the maximum displacement data, and reject abnormal data in each maximum displacement data.
- the second calculating unit 122 is configured to perform a polynomial fitting on each maximum displacement data after the abnormal data is removed, to obtain the maximum displacement attenuation curve.
- the first determining module 13 is specifically configured to:
- a coefficient corresponding to the highest power variable of the maximum displacement attenuation curve is determined to be a viscoelastic parameter of the viscoelastic medium.
- the detection device of this embodiment may be used to implement the technical solution of the method embodiment shown in FIG. 1 , and the implementation principle and technical effects are similar, and details are not described herein again.
- FIG. 5 is a schematic diagram of a second embodiment of a viscoelastic parameter detecting apparatus for a viscoelastic medium according to the present invention.
- the control host 1 further includes: a second acquiring module 15 a second determining module 16,
- the second obtaining module 15 is configured to acquire an elastic parameter of the viscoelastic medium according to the ultrasonic echo signal.
- the second determining module 16 is configured to determine a degree of fibrosis of the viscoelastic medium according to the elastic parameter and the viscoelastic parameter.
- the detection device of this embodiment may be used to perform the technical solution of the method embodiment shown in FIG. 3, and the implementation principle and the technical effect are similar, and details are not described herein again.
- the computer can read the storage medium, and when executed, the program executes the steps including the foregoing method embodiments; and the foregoing storage medium includes: a medium that can store program codes, such as a ROM, a RAM, a magnetic disk, or an optical disk.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Gastroenterology & Hepatology (AREA)
- Signal Processing (AREA)
- Vascular Medicine (AREA)
- Endocrinology (AREA)
- Power Engineering (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
Description
Claims (10)
- 一种粘弹性介质的粘弹性参数检测方法,其特征在于,包括:对粘弹性介质施加单一预设频率的机械振动,在所述粘弹性介质中产生剪切波;在所述剪切波在粘弹性介质传播的期间,向所述粘弹性介质发射单振源超声波,并接收超声回波信号;根据所述超声回波信号获取所述剪切波在各深度下的最大位移数据,每个所述最大位移数据表征了所述剪切波传播到所述粘弹性介质中的不同深度时所述剪切波的最大振荡幅度;拟合各最大位移数据,获得最大位移衰减曲线;根据所述最大位移衰减曲线确定所述粘弹性介质的粘弹性参数。
- 根据权利要求1所述的方法,其特征在于,所述根据所述超声回波信号获取所述剪切波的各最大位移数据之前,还包括:对所述超声回波信号进行以下信号处理中的至少一种信号处理:时域互相关、谱互相关、平方误差和、斑点跟踪、尺度不变特征点跟踪、动态规划、零交叉跟踪、峰值搜索。
- 根据权利要求1所述的方法,其特征在于,所述拟合所述各最大位移数据,获得最大位移衰减曲线,包括:对所述各最大位移数据进行时域和频域的滤波处理,剔除所述各最大位移数据中的异常数据;对剔除所述异常数据后的各最大位移数据进行多项式拟合,获得所述最大位移衰减曲线。
- 根据权利要求3所述的方法,其特征在于,所述根据所述最大位移衰减曲线确定所述粘弹性介质的粘弹性参数,包括:确定所述最大位移衰减曲线的最高次幂变量对应的系数为所述粘弹性介质的粘弹性参数。
- 根据权利要求1至4中任一项所述的方法,其特征在于,所述方法还包括:根据所述超声回波信号,获取所述粘弹性介质的弹性参数;根据所述弹性参数和所述粘弹性参数确定所述粘弹性介质的纤维化 程度。
- 一种粘弹性介质的粘弹性参数检测设备,其特征在于,包括:控制主机、探头,所述探头中包含振动器、超声换能器;所述振动器在所述控制主机的控制下,对粘弹性介质施加单一预设频率的机械振动,在所述粘弹性介质中产生剪切波;所述超声换能器在所述控制主机的控制下,在所述剪切波在粘弹性介质传播的期间,向所述粘弹性介质发射单振源超声波,并接收超声回波信号;所述控制主机,包括:第一获取模块,用于根据所述超声回波信号获取所述剪切波在各深度下的最大位移数据,每个所述最大位移数据表征了所述剪切波传播到所述粘弹性介质中的不同深度时所述剪切波的最大振荡幅度;计算模块,用于拟合各最大位移数据,获得最大位移衰减曲线;第一确定模块,用于根据所述最大位移衰减曲线确定所述粘弹性介质的粘弹性参数。
- 根据权利要求6所述的设备,其特征在于,所述控制主机还包括:处理模块,用于对所述超声回波信号进行以下信号处理中的至少一种信号处理:时域互相关、谱互相关、平方误差和、斑点跟踪、尺度不变特征点跟踪、动态规划、零交叉跟踪、峰值搜索。
- 根据权利要求6所述的设备,其特征在于,所述计算模块包括:第一计算单元,用于对所述各最大位移数据进行时域和频域的滤波处理,剔除所述各最大位移数据中的异常数据;第二计算单元,用于对剔除所述异常数据后的各最大位移数据进行多项式拟合,获得所述最大位移衰减曲线。
- 根据权利要求8所述的设备,其特征在于,所述第一确定模块具体用于:确定所述最大位移衰减曲线的最高次幂变量对应的系数为所述粘弹性介质的粘弹性参数。
- 根据权利要求6至9中任一项所述的设备,其特征在于,所述控 制主机还包括:第二获取模块,用于根据所述超声回波信号,获取所述粘弹性介质的弹性参数;第二确定模块,用于根据所述弹性参数和所述粘弹性参数确定所述粘弹性介质的纤维化程度。
Priority Applications (13)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018532567A JP6737889B2 (ja) | 2015-12-24 | 2016-10-27 | 粘弾性媒体の粘弾性パラメータ検出方法及び装置 |
KR1020187021155A KR101914021B1 (ko) | 2015-12-24 | 2016-10-27 | 점탄성 매질의 점탄성 파라미터 검출 방법 및 기기 |
ES16877471T ES2840699T3 (es) | 2015-12-24 | 2016-10-27 | Método y dispositivo para medir un parámetro viscoelástico de un medio viscoelástico |
CA3009241A CA3009241C (en) | 2015-12-24 | 2016-10-27 | Method and device for detecting viscoelastic parameter of viscoelastic medium |
AU2016378228A AU2016378228B2 (en) | 2015-12-24 | 2016-10-27 | Method and apparatus for measuring viscoelastic parameter of viscoelastic medium |
EP16877471.9A EP3395256B1 (en) | 2015-12-24 | 2016-10-27 | Method and apparatus for measuring viscoelastic parameter of viscoelastic medium |
MX2018007600A MX2018007600A (es) | 2015-12-24 | 2016-10-27 | Metodo y dispositivo para detectar el parametro viscoelastico del medio viscoelastico. |
DK16877471.9T DK3395256T3 (da) | 2015-12-24 | 2016-10-27 | Fremgangsmåde og anordning til detektering af viskoelastisk parameter for viskoelastisk medium |
PL16877471T PL3395256T3 (pl) | 2015-12-24 | 2016-10-27 | Sposób i aparat do pomiaru parametru lepkosprężystości ośrodka lepkosprężystego |
RU2018126777A RU2688299C1 (ru) | 2015-12-24 | 2016-10-27 | Способ и устройство для определения параметра вязкоупругости вязкоупругой среды |
BR112018012589-0A BR112018012589B1 (pt) | 2015-12-24 | 2016-10-27 | Método e dispositivo para detectar parâmetro viscoelástico de meio viscoelástico |
US16/012,702 US10420528B2 (en) | 2015-12-24 | 2018-06-19 | Method and device for detecting viscoelastic parameter of viscoelastic medium |
ZA2018/04119A ZA201804119B (en) | 2015-12-24 | 2018-06-20 | Method and apparatus for measuring viscoelastic parameter of viscoelastic medium |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510993421.3 | 2015-12-24 | ||
CN201510993421.3A CN105455851B (zh) | 2015-12-24 | 2015-12-24 | 粘弹性介质的粘弹性参数检测方法和设备 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US16/012,702 Continuation US10420528B2 (en) | 2015-12-24 | 2018-06-19 | Method and device for detecting viscoelastic parameter of viscoelastic medium |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2017107660A1 true WO2017107660A1 (zh) | 2017-06-29 |
Family
ID=55594354
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2016/103645 WO2017107660A1 (zh) | 2015-12-24 | 2016-10-27 | 粘弹性介质的粘弹性参数检测方法和设备 |
Country Status (15)
Country | Link |
---|---|
US (1) | US10420528B2 (zh) |
EP (1) | EP3395256B1 (zh) |
JP (1) | JP6737889B2 (zh) |
KR (1) | KR101914021B1 (zh) |
CN (1) | CN105455851B (zh) |
AU (1) | AU2016378228B2 (zh) |
BR (1) | BR112018012589B1 (zh) |
CA (1) | CA3009241C (zh) |
DK (1) | DK3395256T3 (zh) |
ES (1) | ES2840699T3 (zh) |
MX (1) | MX2018007600A (zh) |
PL (1) | PL3395256T3 (zh) |
RU (1) | RU2688299C1 (zh) |
WO (1) | WO2017107660A1 (zh) |
ZA (1) | ZA201804119B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20190049566A (ko) * | 2017-10-31 | 2019-05-09 | 지멘스 메디컬 솔루션즈 유에스에이, 인크. | 초음파 의료 이미징에서 전단 속도로부터의 조직 점탄성 추정 |
RU2744929C1 (ru) * | 2017-07-21 | 2021-03-17 | Вуси Хиски Медикал Текнолоджис Ко., Лтд. | Способ и устройство для количественной оценки вязкоупругости среды |
CN113075292A (zh) * | 2020-01-03 | 2021-07-06 | 广州汽车集团股份有限公司 | 一种汽车机油质量的测量方法、装置及存储介质 |
CN114224382A (zh) * | 2021-12-17 | 2022-03-25 | 重庆医科大学 | 粘弹性测量方法及其*** |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105455851B (zh) * | 2015-12-24 | 2018-03-13 | 无锡海斯凯尔医学技术有限公司 | 粘弹性介质的粘弹性参数检测方法和设备 |
JP6212160B1 (ja) * | 2016-04-21 | 2017-10-11 | 株式会社日立製作所 | 超音波診断装置 |
CN105726063B (zh) * | 2016-04-29 | 2019-01-29 | 深圳大学 | 一种基于圆环外部振动的瞬时剪切波激励***及方法 |
CN106383996A (zh) * | 2016-09-05 | 2017-02-08 | 四川大学 | 肝脏组织弹粘性力学新模型的建立以及基于新模型的应力应变关系和剪切波速度的获得方法 |
US11284864B2 (en) * | 2017-06-09 | 2022-03-29 | Technische Universites Eindhoven | Shear wave viscoelasticity imaging using local system identification |
CN107505393A (zh) * | 2017-07-21 | 2017-12-22 | 无锡海斯凯尔医学技术有限公司 | 介质粘弹性的测量方法和装置 |
CN107505232B (zh) | 2017-07-21 | 2019-09-03 | 无锡海斯凯尔医学技术有限公司 | 运动信息获取方法及装置 |
CN107505233A (zh) * | 2017-07-21 | 2017-12-22 | 无锡海斯凯尔医学技术有限公司 | 介质粘弹性定量方法及装置 |
GB2568946A (en) * | 2017-12-01 | 2019-06-05 | Hydramotion Ltd | Fluid property measurement by reflection of vibrational waves |
CN108095764B (zh) * | 2018-01-18 | 2024-05-03 | 北京索瑞特医学技术有限公司 | 复合探头及测量*** |
FR3078484A1 (fr) * | 2018-03-02 | 2019-09-06 | Echosens | Procede de mesure d’un parametre d’attenuation ultrasonore guide par elastographie harmonique, sonde et dispositif pour la mise en œuvre du procede |
CN116831621A (zh) * | 2018-04-28 | 2023-10-03 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声弹性测量装置及方法 |
CN110547825B (zh) * | 2019-09-03 | 2022-04-26 | 上海交通大学 | 一种基于多频剪切波幅值分析的超声弹性成像技术 |
CN113476075B (zh) * | 2020-03-16 | 2024-07-05 | 深圳市理邦精密仪器股份有限公司 | 一种超声弹性成像方法、图像数据的筛选方法及装置 |
CN114343719B (zh) * | 2022-03-17 | 2022-05-31 | 深圳华声医疗技术股份有限公司 | 超声成像控制方法、超声成像端、超声成像设备及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101699280A (zh) * | 2009-10-15 | 2010-04-28 | 北京索瑞特医学技术有限公司 | 超声无损检测粘弹性介质弹性的方法及其装置 |
CN102078205A (zh) * | 2011-03-04 | 2011-06-01 | 深圳市一体医疗科技股份有限公司 | 一种测量粘弹性介质弹性的位移估计方法及应用方法 |
WO2015116893A1 (en) * | 2014-01-30 | 2015-08-06 | General Electric Company | Methods and systems for display of shear-wave elastography and strain elastography images |
CN104825195A (zh) * | 2015-05-25 | 2015-08-12 | 无锡海斯凯尔医学技术有限公司 | 剪切波粘弹性成像方法及*** |
CN105455851A (zh) * | 2015-12-24 | 2016-04-06 | 无锡海斯凯尔医学技术有限公司 | 粘弹性介质的粘弹性参数检测方法和设备 |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6508768B1 (en) * | 2000-11-22 | 2003-01-21 | University Of Kansas Medical Center | Ultrasonic elasticity imaging |
US20030220556A1 (en) * | 2002-05-20 | 2003-11-27 | Vespro Ltd. | Method, system and device for tissue characterization |
WO2007083745A1 (ja) * | 2006-01-20 | 2007-07-26 | Hitachi Medical Corporation | 弾性画像表示方法及び弾性画像表示装置 |
US8118744B2 (en) * | 2007-02-09 | 2012-02-21 | Duke University | Methods, systems and computer program products for ultrasound shear wave velocity estimation and shear modulus reconstruction |
CN101784234B (zh) * | 2007-05-16 | 2014-07-09 | 超声成像公司 | 用于测量感兴趣区域的平均粘弹值的方法和设备 |
US9044192B2 (en) * | 2008-04-07 | 2015-06-02 | Mayo Foundation For Medical Education And Research | System and method for non-invasively measuring tissue viscoelasticity using surface waves |
US9364194B2 (en) * | 2008-09-18 | 2016-06-14 | General Electric Company | Systems and methods for detecting regions of altered stiffness |
CN103269639B (zh) * | 2010-12-22 | 2015-11-25 | 皇家飞利浦电子股份有限公司 | 利用质心估计剪切波速度 |
US8469891B2 (en) | 2011-02-17 | 2013-06-25 | Siemens Medical Solutions Usa, Inc. | Viscoelasticity measurement using amplitude-phase modulated ultrasound wave |
US9244169B2 (en) * | 2012-06-25 | 2016-01-26 | Siemens Medical Solutions Usa, Inc. | Measuring acoustic absorption or attenuation of ultrasound |
US10136876B2 (en) * | 2012-09-27 | 2018-11-27 | Mayo Foundation For Medical Education And Research | System and method for local estimation of nonlinear tissue elasticity with acoustic radiation force |
CN103054552B (zh) * | 2012-12-24 | 2014-12-10 | 深圳先进技术研究院 | 生物组织粘弹性测量方法和*** |
US10743814B2 (en) * | 2013-03-15 | 2020-08-18 | Siemens Medical Solutions Usa, Inc. | Fat fraction estimation using ultrasound with shear wave propagation |
US10512450B2 (en) * | 2013-09-24 | 2019-12-24 | Siemens Medical Solutions Usa, Inc. | Shear wave estimation from analytic data |
JP5730978B2 (ja) * | 2013-11-08 | 2015-06-10 | 日立アロカメディカル株式会社 | 超音波診断装置、及び方法 |
-
2015
- 2015-12-24 CN CN201510993421.3A patent/CN105455851B/zh active Active
-
2016
- 2016-10-27 JP JP2018532567A patent/JP6737889B2/ja active Active
- 2016-10-27 DK DK16877471.9T patent/DK3395256T3/da active
- 2016-10-27 KR KR1020187021155A patent/KR101914021B1/ko active IP Right Grant
- 2016-10-27 RU RU2018126777A patent/RU2688299C1/ru active
- 2016-10-27 ES ES16877471T patent/ES2840699T3/es active Active
- 2016-10-27 PL PL16877471T patent/PL3395256T3/pl unknown
- 2016-10-27 AU AU2016378228A patent/AU2016378228B2/en active Active
- 2016-10-27 WO PCT/CN2016/103645 patent/WO2017107660A1/zh active Application Filing
- 2016-10-27 CA CA3009241A patent/CA3009241C/en active Active
- 2016-10-27 EP EP16877471.9A patent/EP3395256B1/en active Active
- 2016-10-27 BR BR112018012589-0A patent/BR112018012589B1/pt active IP Right Grant
- 2016-10-27 MX MX2018007600A patent/MX2018007600A/es active IP Right Grant
-
2018
- 2018-06-19 US US16/012,702 patent/US10420528B2/en active Active
- 2018-06-20 ZA ZA2018/04119A patent/ZA201804119B/en unknown
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101699280A (zh) * | 2009-10-15 | 2010-04-28 | 北京索瑞特医学技术有限公司 | 超声无损检测粘弹性介质弹性的方法及其装置 |
CN102078205A (zh) * | 2011-03-04 | 2011-06-01 | 深圳市一体医疗科技股份有限公司 | 一种测量粘弹性介质弹性的位移估计方法及应用方法 |
WO2015116893A1 (en) * | 2014-01-30 | 2015-08-06 | General Electric Company | Methods and systems for display of shear-wave elastography and strain elastography images |
CN104825195A (zh) * | 2015-05-25 | 2015-08-12 | 无锡海斯凯尔医学技术有限公司 | 剪切波粘弹性成像方法及*** |
CN105455851A (zh) * | 2015-12-24 | 2016-04-06 | 无锡海斯凯尔医学技术有限公司 | 粘弹性介质的粘弹性参数检测方法和设备 |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2744929C1 (ru) * | 2017-07-21 | 2021-03-17 | Вуси Хиски Медикал Текнолоджис Ко., Лтд. | Способ и устройство для количественной оценки вязкоупругости среды |
US11719613B2 (en) | 2017-07-21 | 2023-08-08 | Wuxi Hisky Medical Technologies Co., Ltd. | Method and device for quantifying viscoelasticity of a medium |
KR20190049566A (ko) * | 2017-10-31 | 2019-05-09 | 지멘스 메디컬 솔루션즈 유에스에이, 인크. | 초음파 의료 이미징에서 전단 속도로부터의 조직 점탄성 추정 |
KR102210372B1 (ko) | 2017-10-31 | 2021-02-01 | 지멘스 메디컬 솔루션즈 유에스에이, 인크. | 초음파 의료 이미징에서 전단 속도로부터의 조직 점탄성 추정 |
US11154277B2 (en) | 2017-10-31 | 2021-10-26 | Siemens Medical Solutions Usa, Inc. | Tissue viscoelastic estimation from shear velocity in ultrasound medical imaging |
CN113075292A (zh) * | 2020-01-03 | 2021-07-06 | 广州汽车集团股份有限公司 | 一种汽车机油质量的测量方法、装置及存储介质 |
CN113075292B (zh) * | 2020-01-03 | 2023-12-22 | 广州汽车集团股份有限公司 | 一种汽车机油质量的测量方法、装置及存储介质 |
CN114224382A (zh) * | 2021-12-17 | 2022-03-25 | 重庆医科大学 | 粘弹性测量方法及其*** |
CN114224382B (zh) * | 2021-12-17 | 2023-09-15 | 重庆医科大学 | 粘弹性测量方法及其*** |
Also Published As
Publication number | Publication date |
---|---|
MX2018007600A (es) | 2019-01-30 |
CN105455851B (zh) | 2018-03-13 |
CN105455851A (zh) | 2016-04-06 |
US10420528B2 (en) | 2019-09-24 |
US20180296181A1 (en) | 2018-10-18 |
JP6737889B2 (ja) | 2020-08-12 |
AU2016378228B2 (en) | 2021-05-27 |
KR20180096768A (ko) | 2018-08-29 |
BR112018012589B1 (pt) | 2023-09-26 |
KR101914021B1 (ko) | 2018-10-31 |
RU2688299C1 (ru) | 2019-05-21 |
PL3395256T3 (pl) | 2021-04-19 |
CA3009241A1 (en) | 2017-06-29 |
JP2018538094A (ja) | 2018-12-27 |
BR112018012589A2 (zh) | 2018-12-26 |
CA3009241C (en) | 2020-05-12 |
EP3395256A4 (en) | 2019-08-14 |
EP3395256A1 (en) | 2018-10-31 |
DK3395256T3 (da) | 2021-01-04 |
ES2840699T3 (es) | 2021-07-07 |
EP3395256B1 (en) | 2020-10-14 |
ZA201804119B (en) | 2019-04-24 |
AU2016378228A1 (en) | 2018-07-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2017107660A1 (zh) | 粘弹性介质的粘弹性参数检测方法和设备 | |
Deng et al. | Ultrasonic shear wave elasticity imaging sequencing and data processing using a verasonics research scanner | |
CN101784234B (zh) | 用于测量感兴趣区域的平均粘弹值的方法和设备 | |
KR101398948B1 (ko) | 진폭―위상 변조된 초음파를 이용한 점탄성 측정 | |
AU2016347143B2 (en) | Elasticity detection method and device | |
CN106999162B (zh) | 超声波诊断装置以及弹性评价方法 | |
CN105266851B (zh) | 粘弹性介质的粘度参数检测方法和设备 | |
US8734350B2 (en) | System and method for correcting errors in shear wave measurements arising from ultrasound beam geometry | |
Urban et al. | Error in estimates of tissue material properties from shear wave dispersion ultrasound vibrometry | |
JPWO2007063619A1 (ja) | 超音波診断装置 | |
WO2017121170A1 (zh) | 组织参数检测方法和*** | |
Castaneda et al. | Application of numerical methods to elasticity imaging | |
Lu et al. | A real time displacement estimation algorithm for ultrasound elastography | |
Deng et al. | Evaluation of TT-based local PWV estimation for different propagation velocities | |
US9737283B2 (en) | Ultrasound measurement apparatus and ultrasound measurement method | |
Aristizabal et al. | Viscoelastic characterization of transverse isotropic tissue mimicking phantoms and muscle | |
BR102020011526A2 (pt) | método e dispositivo para a medição de um parâmetro de ultrassom de um meio viscoelástico | |
RU2794039C2 (ru) | Способ гибридной эластографии, зонд и устройство для гибридной эластографии | |
US20240090875A1 (en) | System and Method for Non-Invasive Determination of Bladder Overactivity Using Ultrasound Vibrometry | |
JP2014030774A (ja) | 関心領域の粘弾性の平均値を測定するための方法および装置 | |
Barboza | Evaluation of Elastographic techniques generated by means of external vibration | |
Yao et al. | Match filter design for improved strain estimation in transient elastography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 16877471 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 3009241 Country of ref document: CA |
|
WWE | Wipo information: entry into national phase |
Ref document number: MX/A/2018/007600 Country of ref document: MX |
|
ENP | Entry into the national phase |
Ref document number: 2018532567 Country of ref document: JP Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
REG | Reference to national code |
Ref country code: BR Ref legal event code: B01A Ref document number: 112018012589 Country of ref document: BR |
|
ENP | Entry into the national phase |
Ref document number: 2016378228 Country of ref document: AU Date of ref document: 20161027 Kind code of ref document: A |
|
ENP | Entry into the national phase |
Ref document number: 20187021155 Country of ref document: KR Kind code of ref document: A |
|
WWE | Wipo information: entry into national phase |
Ref document number: 1020187021155 Country of ref document: KR |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2016877471 Country of ref document: EP Ref document number: 2018126777 Country of ref document: RU |
|
ENP | Entry into the national phase |
Ref document number: 2016877471 Country of ref document: EP Effective date: 20180724 |
|
REG | Reference to national code |
Ref country code: BR Ref legal event code: B01E Ref document number: 112018012589 Country of ref document: BR Free format text: APRESENTAR O RELATORIO DESCRITIVO. |
|
REG | Reference to national code |
Ref country code: BR Ref legal event code: B01Y Ref document number: 112018012589 Country of ref document: BR Kind code of ref document: A2 Free format text: ANULADA A PUBLICACAO CODIGO 1.5 NA RPI NO 2500 DE 04/12/2018 POR TER SIDO INDEVIDA. |
|
ENP | Entry into the national phase |
Ref document number: 112018012589 Country of ref document: BR Kind code of ref document: A2 Effective date: 20180620 |