CN110313936A - Ultrasonic signal processing method, apparatus, equipment and storage medium - Google Patents

Ultrasonic signal processing method, apparatus, equipment and storage medium Download PDF

Info

Publication number
CN110313936A
CN110313936A CN201910706440.1A CN201910706440A CN110313936A CN 110313936 A CN110313936 A CN 110313936A CN 201910706440 A CN201910706440 A CN 201910706440A CN 110313936 A CN110313936 A CN 110313936A
Authority
CN
China
Prior art keywords
signal
ultrasonic signal
destination organization
frequency
section
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
CN201910706440.1A
Other languages
Chinese (zh)
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.)
Wuxi Hisky Medical Technologies Co Ltd
Original Assignee
Wuxi Hisky Medical Technologies Co Ltd
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 Wuxi Hisky Medical Technologies Co Ltd filed Critical Wuxi Hisky Medical Technologies Co Ltd
Priority to CN201910706440.1A priority Critical patent/CN110313936A/en
Publication of CN110313936A publication Critical patent/CN110313936A/en
Priority to PCT/CN2020/105010 priority patent/WO2021018103A1/en
Pending legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic 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
    • A61B8/0858Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
    • 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
    • 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/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data

Landscapes

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

Abstract

The application provides a kind of ultrasonic signal processing method, apparatus, equipment and storage medium, this method comprises: obtaining the target ultrasonic signal of the corresponding dynamic wideband of destination organization;According to the target ultrasonic signal, the decay characteristics of the destination organization are determined;According to the decay characteristics of the destination organization, evaluation processing is carried out to the destination organization.Pass through the target ultrasonic signal according to dynamic wideband, determine the decay characteristics of destination organization, and evaluation processing is carried out to destination organization according to the decay characteristics of destination organization, due to combining the ultrasonic signal of destination organization at multiple frequencies, so that signal component is rich and varied, the accuracy to destination organization evaluation is improved.

Description

Ultrasonic signal processing method, apparatus, equipment and storage medium
Technical field
This application involves ultrasonic technique fields more particularly to a kind of ultrasonic signal processing method, apparatus, equipment and storage to be situated between Matter.
Background technique
With the development of science and technology ultrasonic imaging technique is applied extensively in each field, and because its have in real time, it is cheap, The advantages that Noninvasive and Non-ionizing radiation and be widely applied and clinical diagnosis, quantitative ultrasound can provide very for clinician Intuitive quantitative assessment, such as elasticity, blood flow etc..
But quantitative ultrasound to ultrasound itself characteristics of signals it is related, and be easy organized in or surrounding other believe Number interference.Therefore, how accurately to judge that the state of destination organization, the technology for becoming urgent need to resolve are asked using ultrasonic technique Topic.
Summary of the invention
The application provides a kind of ultrasonic signal processing method, apparatus, equipment and storage medium, to solve the prior art to mesh The state for marking tissue judges the defects of inaccurate.
The application first aspect provides a kind of ultrasonic signal processing method, comprising:
Obtain the target ultrasonic signal of the corresponding dynamic wideband of destination organization;
According to the target ultrasonic signal, the decay characteristics of the destination organization are determined;
According to the decay characteristics of the destination organization, evaluation processing is carried out to the destination organization.
The application the second aspect provides a kind of ultrasonic signal processing device, comprising:
Module is obtained, for obtaining the target ultrasonic signal of the corresponding dynamic wideband of destination organization;
Determining module, for determining the decay characteristics of the destination organization according to the target ultrasonic signal;
Processing module carries out evaluation processing to the destination organization for the decay characteristics according to the destination organization.
A kind of computer equipment is provided in terms of the application third, comprising: at least one processor and memory;
The memory stores computer program;At least one described processor executes the computer of the memory storage Program, the method to realize first aspect offer.
The 4th aspect of the application provides a kind of computer readable storage medium, stores in the computer readable storage medium There is computer program, the computer program is performed the method for realizing that first aspect provides.
Ultrasonic signal processing method, apparatus, equipment and storage medium provided by the present application, pass through the mesh according to dynamic wideband Ultrasonic signal is marked, determines the decay characteristics of destination organization, and evaluate destination organization according to the decay characteristics of destination organization Processing, improves so that signal component is rich and varied to mesh due to combining the ultrasonic signal of destination organization at multiple frequencies Mark the accuracy of tissue-estimating.
Detailed description of the invention
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is this Shen Some embodiments please for those of ordinary skill in the art without any creative labor, can be with It obtains other drawings based on these drawings.
Fig. 1 is the flow diagram for the ultrasonic signal processing method that one embodiment of the application provides;
Fig. 2 is the flow diagram for the ultrasonic signal processing method that another embodiment of the application provides;
Fig. 3 is the structural schematic diagram for the ultrasonic signal processing device that one embodiment of the application provides;
Fig. 4 is the structural schematic diagram for the computer equipment that one embodiment of the application provides.
Through the above attached drawings, it has been shown that the specific embodiment of the application will be hereinafter described in more detail.These attached drawings It is not intended to limit the scope of this disclosure concept by any means with verbal description, but is by referring to specific embodiments Those skilled in the art illustrate the concept of the application.
Specific embodiment
To keep the purposes, technical schemes and advantages of the embodiment of the present application clearer, below in conjunction with the embodiment of the present application In attached drawing, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that described embodiment is Some embodiments of the present application, instead of all the embodiments.Based on the embodiment in the application, those of ordinary skill in the art Every other embodiment obtained without creative efforts, shall fall in the protection scope of this application.
Ultrasonic signal processing method provided by the embodiments of the present application is suitable for using dynamic wideband probe detection target group It knits, obtains the ultrasonic signal of the dynamic wideband of destination organization, and the ultrasonic signal of the dynamic wideband based on destination organization, judge mesh Mark the scene of the state of tissue.Destination organization can be liver organization, musculature, adipose tissue, breast tissue, thyroid gland group It knits and its hetero-organization etc., specifically without limitation.The state of destination organization can be normal condition or pathological condition, if disease Change state can also include lesion type and corresponding lesion grade.For example liver organization, pathological condition include fatty liver, tumour Deng fatty liver may include mild fatty liver, moderate fatty liver, severe fatty liver etc..Specific lesion grade can be according to reality The setting of border demand, the embodiment of the present application is without limitation.
In addition, term " first ", " second " etc. are used for description purposes only, it is not understood to indicate or imply relatively important Property or implicitly indicate the quantity of indicated technical characteristic.In the description of following embodiment, the meaning of " plurality " is two More than a, unless otherwise specifically defined.
These specific embodiments can be combined with each other below, may be at certain for the same or similar concept or process It is repeated no more in a little embodiments.Below in conjunction with attached drawing, the embodiment of the present invention is described.
Embodiment one
The present embodiment provides a kind of ultrasonic signal processing methods, handle for the ultrasonic signal to dynamic wideband, with Judge the state of destination organization.The executing subject of the present embodiment is ultrasonic signal processing device, which, which can be set, is calculating In machine equipment.
As shown in Figure 1, be the flow diagram of ultrasonic signal processing method provided in this embodiment, this method comprises:
Step 101, the target ultrasonic signal of the corresponding dynamic wideband of destination organization is obtained.
Specifically, can pop one's head in by dynamic broad-band ultrasonic to detect destination organization, the raw ultrasound of dynamic wideband is obtained Signal, and certain processing is carried out to raw ultrasound signal, obtain the target ultrasonic signal of the corresponding dynamic wideband of destination organization.
Wherein, the signal bandwidth range that dynamic wideband refers to that ultrasonic probe can work is wider, and the centre frequency of probe is dynamic State is adjustable.Specifically, the ratio that the wideband typically refers to the signal frequency that ultrasonic probe can work and centre frequency is greater than etc. In 60%.Correspondingly, the target ultrasonic signal of dynamic wideband refers to using the probe of wideband under different centre frequency excitations Detect the ultrasonic signal that destination organization obtains.
Illustratively, in the medical field, the ratio of ultrasonic probe bandwidth and centre frequency is greater than 60%, is related to signal frequency Rate is 0.1MHz-100MHz, certainly, is not limited to the ultrasonic probe of this frequency range, and the ultrasound for being also possible to other ranges is visited Head, for example for aerial or geology aspect, then frequency includes 20Hz-0.5MHz.For example the frequency range of ultrasonic probe is 1MHz- 10MHz, then the bandwidth of ultrasonic probe is 9MHz, if centre frequency is 5MHz, the ratio of the ultrasonic probe bandwidth and centre frequency Value is 9/5*100%=180%.
Destination organization can for the human bodies such as liver organization, musculature, adipose tissue, breast tissue or parathyroid tissue or Animal tissue is also possible to the tissue of aerial or geology etc., and the present embodiment does not limit.
Ultrasonic probe is connect with supersonic imaging device, realizes the ultrasound echo signal or transmitted wave letter of acquisition destination organization Number, the embodiment of the present application is known as raw ultrasound signal.Raw ultrasound signal can be it is one-dimensional, two-dimentional, three-dimensional etc., with specific reference to reality The setting of border demand.Supersonic imaging device may include emitter, reception device, imaging processing device etc..Because being that dynamic is wide Frequency ultrasonic probe, the raw ultrasound signal for acquiring acquisition includes the signal of at least one frequency, it can wideband Signal.
It is alternatively possible to carry out data calibration to raw ultrasound signal, and extraction of ingredient interested etc. is handled, is come Obtain target ultrasonic signal.Data calibration be there is directive property due to ultrasonic probe itself, the focusing of probe configures or sets up, visits The influence to the raw ultrasound signal of acquisition such as head sensitivity, system gain and other signal processing means, needs to original super Acoustical signal is calibrated, and for restoring original ultrasonic signal, specific data calibration treatment process can be according to ultrasonic device Hardware, software and algorithm process process be configured.
It extracts ingredient interested to specifically refer to raw ultrasound signal extraction interested information, such as to liver group It knits when being detected, the raw ultrasound signal of acquisition may include the ultrasonic signal of subcutaneous tissue and the ultrasound letter of liver organization Number, it needs to extract the ultrasonic signal of liver organization from raw ultrasound signal.For another example, raw ultrasound signal, including The frequency range of 1-20MHz, based on practical experience, frequency too high arrival destination organization when decayed serious, therefore do not had Effect effect, can extract the ultrasonic signal of 1-5MHz range, as target ultrasonic signal, etc. from raw ultrasound signal. It can specifically handle according to the actual situation, the present embodiment does not limit.
After acquisition obtains raw ultrasound signal, the corresponding dynamic wideband of destination organization can be obtained by certain processing Target ultrasonic signal.
Optionally, after obtaining raw ultrasound signal, original can also be obtained using modes such as frequency filtering or wavelet decompositions The signal of the different frequency of beginning ultrasonic signal.
Step 102, according to target ultrasonic signal, the decay characteristics of destination organization are determined.
Specifically, after the target ultrasonic signal for getting the corresponding dynamic wideband of destination organization, then it can be according to target Ultrasonic signal determines the decay characteristics of destination organization.
It is alternatively possible to be segmented to target ultrasonic signal along depth, by calculate at multistage signal energy (for For every section, for example, each frequency signal amplitude decibel value corresponding with the quadratic sum of amplitude, amplitude and etc.), come Determine the decay characteristics of destination organization.Wherein, the corresponding decibel value of amplitude, which refers to, carries out what log-compressed processing obtained to amplitude Decibel DB value.
Illustratively, it is divided into n sections along depth, n >=2, central depths d1, d2 ..., dn, every segment length is △ l, between depth It is divided into △ d.Every section all comprising the signal of multi-frequency.
After obtaining every section of signal energy, the signal energy of multistage is subjected to linear fit, and oblique after being fitted Rate, by slope divided by centre frequency, the i.e. corresponding decay characteristics of acquisition destination organization.Wherein, centre frequency can pass through frequency Analysis obtains, such as the modes such as Fourier transformation, wavelet transformation.The centre frequency of acquisition be the corresponding amplitude of each frequency signal most The frequency of big signal.
Illustratively, it can be fitted using one-dimensional linear for one-dimensional, least square fitting can be used for two dimension, it is right It can be the prior art using Cubic Spline Fitting etc., specific fit approach in three-dimensional, details are not described herein.
Alternatively it is also possible to calculate separately decay characteristics of the destination organization under every kind of frequency.For unifrequency situation, in Frequency of heart is the frequency, for example only calculates decay characteristics of the destination organization at 3MHz, then centre frequency is 3MHz;For Multi-frequency situation, centre frequency be amplitude maximum signal frequency, for example, frequency include 1MHz, 2MHz, 3MHz, 4MHz, Five kinds of frequencies of 5MHz, wherein the signal frequency of amplitude maximum is 4MHz, then centre frequency is 4MHz.
Step 103, according to the decay characteristics of destination organization, evaluation processing is carried out to destination organization.
Specifically, after the decay characteristics of destination organization have been determined, then it can be according to the decay characteristics of destination organization to mesh Mark tissue carries out evaluation processing.
Illustratively, the state of destination organization can be judged according to the decay characteristics of destination organization.The state may include Normal condition and abnormality.
Optionally, wherein abnormality may include Exception Type and corresponding exception level.
Illustratively, the corresponding ginseng of liver organization normal condition can be obtained ahead of time by taking liver organization as an example in destination organization Examine decay characteristics or corresponding with reference to decay characteristics or ginseng with reference to different lesions grade under decay characteristics range, different lesions type Decay characteristics range is examined, the decay characteristics of destination organization are compared with preset with reference to decay characteristics range, destination organization Decay characteristics be in which kind of situation it is corresponding with reference to decay characteristics range then determine destination organization which kind of situation, such as mesh belonged to It is corresponding with reference to decay characteristics range that mark liver organization decay characteristics belong to normal liver, it is determined that destination organization is normal shape State, if the mild fatty liver pathological changes grade belonged under fatty live lesions type is corresponding to refer to decay characteristics range, it is determined that mesh Marking tissue is mild fatty liver, etc..It here is exemplary illustration, specific state and grade can be according to different tissues Actual conditions are configured, and the present embodiment does not limit.Wherein, exception level embodies the degree of the lesion of tissue.Lesion class Type indicates which class exception tissue produces, for example the Exception Type of liver organization has fatty liver, tumour etc., and fatty liver is corresponding different Normal grade includes mild fatty liver, moderate fatty liver, severe fatty liver etc..Specific Exception Type and exception level can bases The actual conditions of different tissues are arranged, and the present embodiment does not limit.
Ultrasonic signal processing method provided in this embodiment determines mesh by the target ultrasonic signal according to dynamic wideband The decay characteristics of tissue are marked, and evaluation processing is carried out to destination organization according to the decay characteristics of destination organization, due to combining mesh The ultrasonic signal of tissue at multiple frequencies is marked, so that signal component is rich and varied, is improved to the accurate of destination organization evaluation Property.
Embodiment two
The method that the present embodiment provides embodiment one does further supplementary explanation.
As shown in Fig. 2, being the flow diagram of ultrasonic signal processing method provided in this embodiment.
As a kind of enforceable mode, on the basis of the above embodiment 1, optionally, target ultrasonic signal includes extremely A kind of signal of few frequency.Step 102 specifically includes:
Step 1021, target ultrasonic signal is subjected to segment processing along depth direction, obtains at least two sections of ultrasonic signals.
Step 1022, the signal energy for obtaining each section of ultrasonic signal is calculated.
Step 1023, according to the signal energy of each section of ultrasonic signal, the decay characteristics of destination organization are determined.
Specifically, after the target ultrasonic signal for getting the corresponding dynamic wideband of destination organization, it can be to target ultrasound Signal is segmented along depth, by calculate multistage at signal energy (for every section, for example, each frequency signal width Value the decibel value corresponding with quadratic sum of amplitude, amplitude and waiting), to determine the decay characteristics of destination organization.
Illustratively, it is divided into n sections along depth, n >=2, central depths d1, d2 ..., dn, every segment length is △ l, between depth It is divided into △ d.Every section all comprising the signal of multi-frequency.
After obtaining every section of signal energy, the signal energy of multistage is subjected to linear fit, and oblique after being fitted Rate, by slope divided by centre frequency, the i.e. corresponding decay characteristics of acquisition destination organization.Wherein, centre frequency can pass through frequency Analysis obtains, such as the modes such as Fourier transformation, wavelet transformation.The centre frequency of acquisition be the corresponding amplitude of each frequency signal most The frequency of big signal.
Optionally, for one-dimensional signal, it can be and be segmented according to depth line;For two dimension, can be divided into more Bar line, is segmented, for example be classified as a line according to each of pixel by line, constitutes two-dimensional ultrasonic signal by a plurality of line, It is also possible to be segmented by the fritter of ultrasonic signal figure;For three-dimensional, it can be divided into a plurality of line, be segmented by line, It is also possible to be segmented according to the volume of ultrasonic signal figure, can be specifically arranged according to actual needs, the present embodiment does not limit It is fixed.
Alternatively it is also possible to calculate separately decay characteristics of the destination organization under every kind of frequency.For unifrequency situation, in Frequency of heart is the frequency, for example only calculates decay characteristics of the destination organization at 3MHz, then centre frequency is 3MHz;For Multi-frequency situation, centre frequency be amplitude maximum signal frequency, for example, frequency include 1MHz, 2MHz, 3MHz, 4MHz, Five kinds of frequencies of 5MHz, wherein the signal frequency of amplitude maximum is 4MHz, then centre frequency is 4MHz.
Illustratively, n=5, the signal energy at d1=1cm is 500, at d2=2cm be 400, be at d3=3cm 300, it is 100 for 200, at d5=5cm at d4=4cm, then carries out linear fit, acquisitions slope is 100/cm, the center of acquisition Frequency is 4MHz, then slope is 25/cm/MHz divided by the decay characteristics that centre frequency obtains.It here is exemplary illustration, and It is non-that it is limited.
Optionally, the signal energy for obtaining each section of ultrasonic signal is calculated, comprising:
Step 2011, for every section of ultrasonic signal, the subsignal energy of the wherein signal of each frequency is calculated.
Step 2012, the signal energy by the sum of each subsignal energy as this section of ultrasonic signal.
Optionally, the signal energy for obtaining each section of ultrasonic signal is calculated, comprising:
Step 2021, the signal energy of each section of ultrasonic signal is obtained using time domain mode.
Optionally, the signal energy for obtaining each section of ultrasonic signal is calculated, comprising:
Step 2031, the signal energy of each section of ultrasonic signal is obtained using frequency domain mode.
Specifically, time domain mode and frequency domain mode obtain the signal energy of each section of ultrasonic signal overall process and above-mentioned mistake Journey is consistent, the difference is that, under frequency domain mode, need to carry out target ultrasonic signal Short Time Fourier Transform, it is more to calculate Signal energy at section, other parts are consistent, and details are not described herein.
Optionally, according to the signal energy of each section of ultrasonic signal, the decay characteristics of destination organization are determined, comprising:
Step 2041, linear fit processing is carried out to the signal energy of each section of ultrasonic signal, the slope after being fitted.
Step 2042, the slope after fitting is obtained into the decaying of destination organization divided by the centre frequency of target ultrasonic signal Feature.
Wherein, centre frequency is the frequency of the signal of amplitude maximum in the signal of each frequency.
Specifically, the signal energy of multistage is carried out linear fit, and be fitted after obtaining every section of signal energy Slope afterwards, by slope divided by centre frequency, the i.e. corresponding decay characteristics of acquisition destination organization.Wherein, centre frequency can lead to Overfrequency analysis obtains, such as the modes such as Fourier transformation, wavelet transformation.The centre frequency of acquisition is that each frequency signal is corresponding The frequency of the signal of amplitude maximum.
Alternatively it is also possible to calculate separately decay characteristics of the destination organization under every kind of frequency.For unifrequency situation, in Frequency of heart is the frequency, for example only calculates decay characteristics of the destination organization at 3MHz, then centre frequency is 3MHz;For Multi-frequency situation, centre frequency be amplitude maximum signal frequency, for example, frequency include 1MHz, 2MHz, 3MHz, 4MHz, Five kinds of frequencies of 5MHz, wherein the signal frequency of amplitude maximum is 4MHz, then centre frequency is 4MHz.
Illustratively, n=5, the signal energy at d1=1cm is 500, at d2=2cm be 400, be at d3=3cm 300, it is 100 for 200, at d5=5cm at d4=4cm, then carries out linear fit, acquisitions slope is 100/cm, the center of acquisition Frequency is 4MHz, then slope is 25/cm/MHz divided by the decay characteristics that centre frequency obtains.It here is exemplary illustration, and It is non-that it is limited.
As another enforceable mode, on the basis of the above embodiment 1, optionally, it is corresponding to obtain destination organization Dynamic wideband target ultrasonic signal, can specifically include:
Step 1011, the raw ultrasound signal of dynamic wideband is obtained.
Step 1012, data calibration processing is carried out to raw ultrasound signal, the signal after being calibrated.
Step 1013, according to the signal after calibration, ingredient interested is extracted, as target ultrasonic signal.
Specifically, can pop one's head in by dynamic broad-band ultrasonic to detect destination organization, the raw ultrasound of dynamic wideband is obtained Signal, and certain processing is carried out to raw ultrasound signal, obtain the target ultrasonic signal of the corresponding dynamic wideband of destination organization.
Illustratively, ultrasonic probe bandwidth is 60%-200%, and being related to signal frequency is 0.5MHz-50MHz, certainly, no Be limited to the ultrasonic probe of this frequency range, be also possible to the ultrasonic probe of other ranges, for example, for aerial or geology in terms of, Then comment includes 20Hz-0.5MHz.
Ultrasonic probe is connect with supersonic imaging device, realizes the ultrasound echo signal or transmitted wave letter of acquisition destination organization Number, the embodiment of the present application is known as raw ultrasound signal.Raw ultrasound signal can be it is one-dimensional, two-dimentional, three-dimensional etc., with specific reference to reality The setting of border demand.Supersonic imaging device may include emitter, reception device, imaging processing device etc..Because being that dynamic is wide Frequency ultrasonic probe, the raw ultrasound signal for acquiring acquisition includes the signal of at least two frequencies.
It is alternatively possible to carry out data calibration to raw ultrasound signal, and extraction of ingredient interested etc. is handled, is come Obtain target ultrasonic signal.Data calibration be there is directive property due to ultrasonic probe itself, the focusing of probe configures or sets up, visits The influence to the raw ultrasound signal of acquisition such as head sensitivity, system gain and other signal processing means, needs to original super Acoustical signal is calibrated, and for restoring original ultrasonic signal, specific data calibration treatment process can be according to ultrasonic device Hardware, software and algorithm process process be configured.
It extracts ingredient interested to specifically refer to raw ultrasound signal extraction interested information, such as to liver group It knits when being detected, the raw ultrasound signal of acquisition may include the ultrasonic signal of subcutaneous tissue and the ultrasound letter of liver organization Number, it needs to extract the ultrasonic signal of liver organization from raw ultrasound signal.For another example, raw ultrasound signal, including The frequency range of 1-20MHz, based on practical experience, frequency too high arrival destination organization when decayed serious, therefore do not had Effect effect, can extract the ultrasonic signal of 1-5MHz range, as target ultrasonic signal, etc. from raw ultrasound signal. It can specifically handle according to the actual situation, the present embodiment does not limit.
After acquisition obtains raw ultrasound signal, the corresponding dynamic wideband of destination organization can be obtained by certain processing Target ultrasonic signal.
Optionally, after obtaining raw ultrasound signal, original can also be obtained using modes such as frequency filtering or wavelet decompositions The signal of the different frequency of beginning ultrasonic signal.
Optionally, according to the decay characteristics of destination organization, evaluation processing is carried out to destination organization, can specifically include:
Step 1031, according to the decay characteristics of destination organization, the state of destination organization, state normal condition and exception are determined State.
Specifically, after the decay characteristics of destination organization have been determined, then it can be according to the decay characteristics of destination organization to mesh Mark tissue carries out evaluation processing.The state of destination organization can be judged according to the decay characteristics of destination organization.The state can wrap Include normal condition and abnormality.
Optionally, wherein abnormality may include Exception Type and corresponding exception level.
Illustratively, illustratively, liver organization normal condition can be obtained ahead of time by taking liver organization as an example in destination organization It is corresponding to decay with reference to decay characteristics or with reference to the corresponding reference of different lesions grade under decay characteristics range, different lesions type Feature refers to decay characteristics range, and the decay characteristics of destination organization are compared with preset with reference to decay characteristics range, The decay characteristics of destination organization are in the corresponding reference decay characteristics range of which kind of situation and then determine which kind of feelings destination organization belongs to Condition, such as target liver organization decay characteristics belong to the corresponding reference decay characteristics range of normal liver, it is determined that destination organization For normal condition, if the mild fatty liver pathological changes grade belonged under fatty live lesions type is corresponding to refer to decay characteristics range, Then determine that destination organization is mild fatty liver, etc..It here is exemplary illustration, specific state and grade can be according to not Actual conditions with tissue are configured, and the present embodiment does not limit.Wherein, exception level embodies the journey of the lesion of tissue Degree.Lesion type indicates which class exception tissue produces, for example the Exception Type of liver organization has fatty liver, tumour etc., fat The corresponding exception level of liver includes mild fatty liver, moderate fatty liver, severe fatty liver etc..Specific Exception Type and exception etc. Grade can be arranged according to the actual conditions of different tissues, and the present embodiment does not limit.
Optionally, after the decay characteristics for obtaining destination organization, display processing can be carried out, specifically can be display numerical value Or numerical value is shown with pcolor, it can specifically be arranged according to actual needs, the present embodiment does not limit.
It should be noted that each enforceable mode can individually be implemented in the present embodiment, it can also be in the feelings not conflicted It is combined in any combination under condition and implements the application without limitation.
Ultrasonic signal processing method provided in this embodiment determines mesh by the target ultrasonic signal according to dynamic wideband The decay characteristics of tissue are marked, and evaluation processing is carried out to destination organization according to the decay characteristics of destination organization, due to combining mesh The ultrasonic signal of tissue at multiple frequencies is marked, so that signal component is rich and varied, is improved to the accurate of destination organization evaluation Property.
Embodiment three
The present embodiment provides a kind of ultrasonic signal processing devices, the method for executing above-described embodiment one.
As shown in figure 3, being the structural schematic diagram of ultrasonic signal processing device provided in this embodiment.The ultrasonic signal processing Device 30 includes obtaining module 31, determining module 32 and processing module 33.
Wherein, the target ultrasonic signal that module 31 is used to obtain the corresponding dynamic wideband of destination organization is obtained;Determining module 32 for determining the decay characteristics of destination organization according to target ultrasonic signal;Processing module 33 is used for declining according to destination organization Subtract feature, evaluation processing is carried out to destination organization.
Device in this present embodiment is closed, wherein modules execute the concrete mode of operation in related this method It is described in detail in embodiment, no detailed explanation will be given here.
According to ultrasonic signal processing device provided in this embodiment, by the target ultrasonic signal according to dynamic wideband, really Set the goal the decay characteristics of tissue, and carries out evaluation processing to destination organization according to the decay characteristics of destination organization, due to synthesis The ultrasonic signal of destination organization at multiple frequencies improves so that signal component is rich and varied to destination organization evaluation Accuracy.
Example IV
The device that the present embodiment provides above-described embodiment three does further supplementary explanation, is mentioned with executing above-described embodiment two The method of confession.
As a kind of enforceable mode, on the basis of above-described embodiment three, optionally, target ultrasonic signal includes extremely A kind of signal of few frequency;
Determining module is specifically used for:
Target ultrasonic signal is subjected to segment processing along depth direction, obtains at least two sections of ultrasonic signals;
Calculate the signal energy for obtaining each section of ultrasonic signal;
According to the signal energy of each section of ultrasonic signal, the decay characteristics of destination organization are determined.
Optionally it is determined that module, is specifically used for:
For every section of ultrasonic signal, the subsignal energy of the wherein signal of each frequency is calculated;
Using the sum of each subsignal energy as the signal energy of this section of ultrasonic signal.
Optionally it is determined that module, is specifically used for:
The signal energy of each section of ultrasonic signal is obtained using time domain mode;Alternatively,
The signal energy of each section of ultrasonic signal is obtained using frequency domain mode.
Optionally it is determined that module, is specifically used for:
Linear fit processing is carried out to the signal energy of each section of ultrasonic signal, the slope after being fitted;
By the slope after fitting divided by the centre frequency of target ultrasonic signal, the decay characteristics of destination organization, center are obtained Frequency is the frequency of the signal of amplitude maximum in the signal of each frequency.
As another enforceable mode, on the basis of above-described embodiment three, optionally, module is obtained, it is specific to use In:
Obtain the raw ultrasound signal of dynamic wideband;
Data calibration processing is carried out to raw ultrasound signal, the signal after being calibrated;
According to the signal after calibration, ingredient interested is extracted, obtains target ultrasonic signal.
As another enforceable mode, on the basis of above-described embodiment three, optionally, processing module is specific to use In:
According to the decay characteristics of destination organization, determine that the state of destination organization, state include normal condition and abnormality.
Device in this present embodiment is closed, wherein modules execute the concrete mode of operation in related this method It is described in detail in embodiment, no detailed explanation will be given here.
It should be noted that each enforceable mode can individually be implemented in the present embodiment, it can also be in the feelings not conflicted It is combined in any combination under condition and implements the application without limitation.
Mesh is determined by the target ultrasonic signal according to dynamic wideband according to the ultrasonic signal processing device of the present embodiment The decay characteristics of tissue are marked, and evaluation processing is carried out to destination organization according to the decay characteristics of destination organization, due to combining mesh The ultrasonic signal of tissue at multiple frequencies is marked, so that signal component is rich and varied, is improved to the accurate of destination organization evaluation Property.
Embodiment five
The present embodiment provides a kind of computer equipments, for executing method provided by the above embodiment.
As shown in figure 4, being the structural schematic diagram of computer equipment provided in this embodiment.The computer equipment 50 includes: At least one processor 51 and memory 52;
Memory stores computer program;At least one processor executes the computer program of memory storage, to realize Method provided by the above embodiment.
Destination organization is determined by the target ultrasonic signal according to dynamic wideband according to the computer equipment of the present embodiment Decay characteristics, and evaluation processing is carried out to destination organization according to the decay characteristics of destination organization, due to combining destination organization Ultrasonic signal at multiple frequencies improves the accuracy to destination organization evaluation so that signal component is rich and varied.
Embodiment six
The present embodiment provides a kind of computer readable storage medium, computer is stored in the computer readable storage medium Program, computer program are performed the method for realizing that any of the above-described embodiment provides.
It is determined according to the computer readable storage medium of the present embodiment by the target ultrasonic signal according to dynamic wideband The decay characteristics of destination organization, and evaluation processing is carried out to destination organization according to the decay characteristics of destination organization, due to combining The ultrasonic signal of destination organization at multiple frequencies improves the standard to destination organization evaluation so that signal component is rich and varied True property.
In several embodiments provided herein, it should be understood that disclosed device and method can pass through it Its mode is realized.For example, the apparatus embodiments described above are merely exemplary, for example, the division of the unit, only Only a kind of logical function partition, there may be another division manner in actual implementation, such as multiple units or components can be tied Another system is closed or is desirably integrated into, or some features can be ignored or not executed.Another point, it is shown or discussed Mutual coupling, direct-coupling or communication connection can be through some interfaces, the INDIRECT COUPLING or logical of device or unit Letter connection can be electrical property, mechanical or other forms.
The unit as illustrated by the separation member may or may not be physically separated, aobvious as unit The component shown may or may not be physical unit, it can and it is in one place, or may be distributed over multiple In network unit.It can select some or all of unit therein according to the actual needs to realize the mesh of this embodiment scheme 's.
It, can also be in addition, each functional unit in each embodiment of the application can integrate in one processing unit It is that each unit physically exists alone, can also be integrated in one unit with two or more units.Above-mentioned integrated list Member both can take the form of hardware realization, can also realize in the form of hardware adds SFU software functional unit.
The above-mentioned integrated unit being realized in the form of SFU software functional unit can store and computer-readable deposit at one In storage media.Above-mentioned SFU software functional unit is stored in a storage medium, including some instructions are used so that a computer It is each that equipment (can be personal computer, server or the network equipment etc.) or processor (processor) execute the application The part steps of embodiment the method.And storage medium above-mentioned includes: USB flash disk, mobile hard disk, read-only memory (Read- Only Memory, ROM), random access memory (Random Access Memory, RAM), magnetic or disk etc. it is various It can store the medium of program code.
Those skilled in the art can be understood that, for convenience and simplicity of description, only with above-mentioned each functional module Division progress for example, in practical application, can according to need and above-mentioned function distribution is complete by different functional modules At the internal structure of device being divided into different functional modules, to complete all or part of the functions described above.On The specific work process for stating the device of description, can refer to corresponding processes in the foregoing method embodiment, and details are not described herein.
Finally, it should be noted that the above various embodiments is only to illustrate the technical solution of the application, rather than its limitations;To the greatest extent Pipe is described in detail the application referring to foregoing embodiments, those skilled in the art should understand that: its according to So be possible to modify the technical solutions described in the foregoing embodiments, or to some or all of the technical features into Row equivalent replacement;And these are modified or replaceed, each embodiment technology of the application that it does not separate the essence of the corresponding technical solution The range of scheme.

Claims (16)

1. a kind of ultrasonic signal processing method characterized by comprising
Obtain the target ultrasonic signal of the corresponding dynamic wideband of destination organization;
According to the target ultrasonic signal, the decay characteristics of the destination organization are determined;
According to the decay characteristics of the destination organization, evaluation processing is carried out to the destination organization.
2. the method according to claim 1, wherein the target ultrasonic signal includes the letter of at least one frequency Number;
It is described according to the target ultrasonic signal, determine the decay characteristics of the destination organization, comprising:
The target ultrasonic signal is subjected to segment processing along depth direction, obtains at least two sections of ultrasonic signals;
Calculate the signal energy for obtaining each section of ultrasonic signal;
According to the signal energy of each section of ultrasonic signal, the decay characteristics of the destination organization are determined.
3. according to the method described in claim 2, it is characterized in that, it is described calculate obtain each section of ultrasonic signal signal energy, Include:
For every section of ultrasonic signal, the subsignal energy of the wherein signal of each frequency is calculated;
Signal energy by the sum of each subsignal energy as this section of ultrasonic signal.
4. according to the method described in claim 2, it is characterized in that, it is described calculate obtain each section of ultrasonic signal signal energy, Include:
The signal energy of each section of ultrasonic signal is obtained using time domain mode;Alternatively,
The signal energy of each section of ultrasonic signal is obtained using frequency domain mode.
5. according to the method described in claim 2, it is characterized in that, the signal energy according to each section of ultrasonic signal, determines The decay characteristics of the destination organization, comprising:
Linear fit processing is carried out to the signal energy of each section of ultrasonic signal, the slope after being fitted;
By the slope after the fitting divided by the centre frequency of the target ultrasonic signal, the decaying for obtaining the destination organization is special Sign, the centre frequency are the frequency of the signal of amplitude maximum in the signal of each frequency.
6. the method according to claim 1, wherein the target for obtaining the corresponding dynamic wideband of destination organization Ultrasonic signal, comprising:
Obtain the raw ultrasound signal of dynamic wideband;
Data calibration processing is carried out to the raw ultrasound signal, the signal after being calibrated;
According to the signal after the calibration, ingredient interested is extracted, as the target ultrasonic signal.
7. method according to claim 1-6, which is characterized in that the decaying according to the destination organization is special Sign, carries out evaluation processing to the destination organization, comprising:
According to the decay characteristics of the destination organization, determine the state of the destination organization, the state include normal condition and Abnormality.
8. a kind of ultrasonic signal processing device characterized by comprising
Module is obtained, for obtaining the target ultrasonic signal of the corresponding dynamic wideband of destination organization;
Determining module, for determining the decay characteristics of the destination organization according to the target ultrasonic signal;
Processing module carries out evaluation processing to the destination organization for the decay characteristics according to the destination organization.
9. device according to claim 8, which is characterized in that the target ultrasonic signal includes the letter of at least one frequency Number;
The determining module, is specifically used for:
The target ultrasonic signal is subjected to segment processing along depth direction, obtains at least two sections of ultrasonic signals;
Calculate the signal energy for obtaining each section of ultrasonic signal;
According to the signal energy of each section of ultrasonic signal, the decay characteristics of the destination organization are determined.
10. device according to claim 9, which is characterized in that the determining module is specifically used for:
For every section of ultrasonic signal, the subsignal energy of the wherein signal of each frequency is calculated;
Signal energy by the sum of each subsignal energy as this section of ultrasonic signal.
11. device according to claim 9, which is characterized in that the determining module is specifically used for:
The signal energy of each section of ultrasonic signal is obtained using time domain mode;Alternatively,
The signal energy of each section of ultrasonic signal is obtained using frequency domain mode.
12. device according to claim 9, which is characterized in that the determining module is specifically used for:
Linear fit processing is carried out to the signal energy of each section of ultrasonic signal, the slope after being fitted;
By the slope after the fitting divided by the centre frequency of the target ultrasonic signal, the decaying for obtaining the destination organization is special Sign, the centre frequency are the frequency of the signal of amplitude maximum in the signal of each frequency.
13. device according to claim 8, which is characterized in that the acquisition module is specifically used for:
Obtain the raw ultrasound signal of dynamic wideband;
Data calibration processing is carried out to the raw ultrasound signal, the signal after being calibrated;
According to the signal after the calibration, ingredient interested is extracted, as the target ultrasonic signal.
14. according to the described in any item devices of claim 8-13, which is characterized in that the processing module is specifically used for:
According to the decay characteristics of the destination organization, determine the state of the destination organization, the state include normal condition and Abnormality.
15. a kind of computer equipment characterized by comprising at least one processor and memory;
The memory stores computer program;At least one described processor executes the computer journey of the memory storage Sequence, to realize method of any of claims 1-7.
16. a kind of computer readable storage medium, which is characterized in that be stored with computer journey in the computer readable storage medium Sequence, the computer program, which is performed, realizes method of any of claims 1-7.
CN201910706440.1A 2019-08-01 2019-08-01 Ultrasonic signal processing method, apparatus, equipment and storage medium Pending CN110313936A (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201910706440.1A CN110313936A (en) 2019-08-01 2019-08-01 Ultrasonic signal processing method, apparatus, equipment and storage medium
PCT/CN2020/105010 WO2021018103A1 (en) 2019-08-01 2020-07-28 Ultrasonic signal processing method and apparatus, device, and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910706440.1A CN110313936A (en) 2019-08-01 2019-08-01 Ultrasonic signal processing method, apparatus, equipment and storage medium

Publications (1)

Publication Number Publication Date
CN110313936A true CN110313936A (en) 2019-10-11

Family

ID=68125174

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910706440.1A Pending CN110313936A (en) 2019-08-01 2019-08-01 Ultrasonic signal processing method, apparatus, equipment and storage medium

Country Status (2)

Country Link
CN (1) CN110313936A (en)
WO (1) WO2021018103A1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021018107A1 (en) * 2019-08-01 2021-02-04 无锡海斯凯尔医学技术有限公司 Ultrasonic signal processing method and ultrasonic signal processing apparatus, and device and storage medium
WO2021018103A1 (en) * 2019-08-01 2021-02-04 无锡海斯凯尔医学技术有限公司 Ultrasonic signal processing method and apparatus, device, and storage medium

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5746209A (en) * 1996-01-26 1998-05-05 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Method of and apparatus for histological human tissue characterizationusing ultrasound
CN103637820A (en) * 2013-11-27 2014-03-19 华南理工大学 Hepatic fibrosis degree relative deviation identification method based on ultrasonic radio frequency (RF) time sequence
CN104042237A (en) * 2013-03-15 2014-09-17 美国西门子医疗解决公司 Fat Fraction Estimation Using Ultrasound with Shear Wave Propagation
CN105377145A (en) * 2013-05-24 2016-03-02 森尼布鲁克研究所 System and method for classifying and characterizing tissues using first-order and second-order statistics of quantitative ultrasound parametric maps
CN107647881A (en) * 2016-07-25 2018-02-02 法国爱科森有限公司 The measuring method of human body or animal organ's viscoelastic parameters
CN109044402A (en) * 2017-09-12 2018-12-21 乐普(北京)医疗器械股份有限公司 A kind of Diagnosis of Fatty liver system based on Ultrasonic Elasticity Imaging

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106170253B (en) * 2015-03-23 2019-03-01 奥林巴斯株式会社 The working method of ultrasound observation apparatus and ultrasound observation apparatus
CN104873221B (en) * 2015-06-05 2018-03-13 无锡海斯凯尔医学技术有限公司 Liver fat quantitative approach and system based on ultrasonic wave
WO2017068892A1 (en) * 2015-10-22 2017-04-27 株式会社日立製作所 Ultrasound diagnostic device and attenuation characteristic measurement method
CN105930665B (en) * 2016-04-26 2018-06-29 北京工业大学 A kind of liver fat method for quantitative measuring based on ultrasonic attenuation coefficient
CN110313936A (en) * 2019-08-01 2019-10-11 无锡海斯凯尔医学技术有限公司 Ultrasonic signal processing method, apparatus, equipment and storage medium

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5746209A (en) * 1996-01-26 1998-05-05 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Method of and apparatus for histological human tissue characterizationusing ultrasound
CN104042237A (en) * 2013-03-15 2014-09-17 美国西门子医疗解决公司 Fat Fraction Estimation Using Ultrasound with Shear Wave Propagation
CN105377145A (en) * 2013-05-24 2016-03-02 森尼布鲁克研究所 System and method for classifying and characterizing tissues using first-order and second-order statistics of quantitative ultrasound parametric maps
CN103637820A (en) * 2013-11-27 2014-03-19 华南理工大学 Hepatic fibrosis degree relative deviation identification method based on ultrasonic radio frequency (RF) time sequence
CN107647881A (en) * 2016-07-25 2018-02-02 法国爱科森有限公司 The measuring method of human body or animal organ's viscoelastic parameters
CN109044402A (en) * 2017-09-12 2018-12-21 乐普(北京)医疗器械股份有限公司 A kind of Diagnosis of Fatty liver system based on Ultrasonic Elasticity Imaging

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
侯宇格 等: "Fibro touch受控衰减参数评价肝细胞脂肪变的相关影响因素", 《世界华人消化杂志》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021018107A1 (en) * 2019-08-01 2021-02-04 无锡海斯凯尔医学技术有限公司 Ultrasonic signal processing method and ultrasonic signal processing apparatus, and device and storage medium
WO2021018103A1 (en) * 2019-08-01 2021-02-04 无锡海斯凯尔医学技术有限公司 Ultrasonic signal processing method and apparatus, device, and storage medium

Also Published As

Publication number Publication date
WO2021018103A1 (en) 2021-02-04

Similar Documents

Publication Publication Date Title
Oosterveld et al. Ultrasound attenuation and texture analysis of diffuse liver disease: methods and preliminary results
EP3806721B1 (en) Method and system for determining fractional fat content of tissue
CN105377145B (en) The single order and second-order statistics that are mapped with quantitative ultrasound parameter classifies and characterizes the system and method for tissue
US11083436B2 (en) Ultrasonic image analysis systems and analysis methods thereof
EP3013244B1 (en) System and method for mapping ultrasound shear wave elastography measurements
US20170251988A1 (en) Ultrasound imaging apparatus
DE102019202545A1 (en) Placement of the region of interest for quantitative ultrasound imaging
EP3400878B1 (en) Method for postural independent location of targets in diagnostic imagines acquired by multimodal acquisitions and system for carrying out the said method
CN110313941B (en) Data processing method, device, equipment and storage medium
CN110313936A (en) Ultrasonic signal processing method, apparatus, equipment and storage medium
EP3701877B1 (en) Ultrasound diagnostic system and ultrasound diagnostic method
CN105025807B (en) For the scanning that stationary fluid pond is identified and is delimited
Daoud et al. Needle detection in curvilinear ultrasound images based on the reflection pattern of circular ultrasound waves
EP2433570A1 (en) Ultrasound system for displaying slice of object and method thereof
Kidera et al. TDOA-based microwave imaging algorithm for real-time microwave ablation monitoring
Kretzek et al. GPU-based 3D SAFT reconstruction including attenuation correction
Daoud et al. Accurate needle localization using two-dimensional power doppler and b-mode ultrasound image analyses: A feasibility study
US20220142615A1 (en) Ultrasonic signal processing method and apparatus, device and storage medium
CN102119341A (en) Method of meshing and calculating a volume in an ultrasound imaging system
US20130315462A1 (en) Use of a Two-Dimensional Analytical Signal in Sonography
CN110368023A (en) Signal processing method, device, equipment and storage medium
JP6483659B2 (en) Beamforming technology for detection of microcalcifications by ultrasound
RU2795207C1 (en) Method and device for processing ultrasonic signals, device and information media
Wang et al. Adaptive sound speed compensation for reflection imaging of ultrasound computed tomography
Abeyratne et al. Ultrasound scatter-spacing based diagnosis of focal diseases of the liver

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: 20191011

RJ01 Rejection of invention patent application after publication