CN108877902B - Ultrasonic image brightness adjusting method and adjusting system - Google Patents
Ultrasonic image brightness adjusting method and adjusting system Download PDFInfo
- Publication number
- CN108877902B CN108877902B CN201810873293.2A CN201810873293A CN108877902B CN 108877902 B CN108877902 B CN 108877902B CN 201810873293 A CN201810873293 A CN 201810873293A CN 108877902 B CN108877902 B CN 108877902B
- Authority
- CN
- China
- Prior art keywords
- gain compensation
- partition
- compensation coefficient
- amplitude value
- variance
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/20—ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- 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/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
Abstract
The invention provides an ultrasonic image brightness adjusting method and an ultrasonic image brightness adjusting system, wherein the method comprises the following steps: acquiring an echo RF signal and performing gain compensation on the echo RF signal to form an envelope amplitude signal; partitioning the two-dimensional data to form a plurality of partitions, and acquiring an energy spectrum density function and an amplitude value variance corresponding to each partition according to the obtained envelope amplitude signal; judging whether the ratio of the number of the amplitude value variances smaller than the preset amplitude value variance in the amplitude value variance sum in each interval to the amplitude value variance sum is larger than a preset ratio threshold value or not, if so, acquiring a first gain compensation coefficient corresponding to each interval according to an energy spectral density function, and performing brightness compensation on the image by using the first gain compensation coefficient and displaying the image; if not, acquiring a second gain compensation coefficient corresponding to each partition according to the amplitude value variance and the first gain compensation coefficient, and replacing the preset gain compensation coefficient with the second gain compensation coefficient to reprocess the echo RF signal. The invention can achieve the purpose of improving the image quality.
Description
Technical Field
The invention belongs to the technical field of medical ultrasound, and mainly relates to an ultrasonic image brightness adjusting method and an ultrasonic image brightness adjusting system.
Background
Color blood flow imaging of a color ultrasonic diagnostic apparatus (B-ultrasonic machine) becomes one of indispensable auxiliary diagnostic means in modern medicine due to the unique real-time dynamic characteristic of the color blood flow imaging, and becomes a judgment standard of certain diseases in clinical diagnosis.
In the ultrasonic imaging process, an ultrasonic diagnostic apparatus transmits ultrasonic signals to tissues, the tissues are reflected to form echo signals and processed to form radio frequency RF signals, the radio frequency RF signals are subjected to RF filtering, the attenuation of the signals in the time direction is compensated through a time gain compensation module, the signals subjected to gain compensation are sent to an orthogonal demodulation module to be demodulated and processed to form envelope signals, and the envelope signals are further subjected to logarithmic compression, image post-processing and coordinate transformation in sequence and then displayed and output.
In order to make the output image clearer, the brightness of the image is generally required to be adjusted to facilitate recognition; in the traditional image brightness adjustment, gain compensation is adjusted based on gray scale, the basic data of the method is image gray scale, and the precision of the method is limited to 256 gray scale; meanwhile, in the process of adjusting the image brightness by adopting the method, a plurality of factors can influence the gray scale, so that the display effect of finally obtaining the output image is still to be improved.
Disclosure of Invention
The invention aims to provide an ultrasonic image brightness adjusting method and an ultrasonic image brightness adjusting system.
In order to achieve one of the above objects, an ultrasound image brightness adjusting method according to an embodiment of the present invention includes: s1, acquiring echo RF signals, and performing gain compensation on the RF signals in the transverse direction and/or the longitudinal direction respectively through preset gain compensation coefficients to form envelope amplitude signals; the preset gain compensation coefficient comprises: presetting a transverse gain compensation coefficient and/or a longitudinal gain compensation coefficient;
s2, partitioning the two-dimensional data to form a plurality of partitions, wherein the partitions comprise transversely distributed transverse partitions and/or longitudinally distributed longitudinal partitions; acquiring image brightness characteristic parameters corresponding to each partition according to the obtained envelope amplitude signals, wherein the image brightness characteristic parameters comprise: an energy spectral density function, and a variance of amplitude values;
s3, analyzing the variance of the amplitude value corresponding to each partition, judging whether the ratio of the number of the variance of the amplitude value corresponding to each partition, which is smaller than the preset variance of the amplitude value, to the total number of the variance of the amplitude value is larger than the preset ratio threshold value,
if yes, obtaining a first gain compensation coefficient corresponding to each partition according to the energy spectrum density function, wherein the first gain compensation coefficient comprises: a first lateral gain compensation factor and/or a first longitudinal gain compensation factor; performing brightness compensation on the image by using the first gain compensation coefficient corresponding to each subarea, and displaying the image;
if not, acquiring a first gain compensation coefficient corresponding to each subarea according to the energy spectrum density function, and acquiring a second gain compensation coefficient corresponding to each subarea according to the amplitude value variance and the first gain compensation coefficient, wherein the second gain compensation coefficient comprises: a second lateral gain compensation factor and/or a second longitudinal gain compensation factor; and replacing a preset gain compensation coefficient with a second gain compensation coefficient corresponding to the current partition to perform gain compensation on the echo RF signal to form a new envelope amplitude signal until the ratio of the number of the amplitude value variances corresponding to each partition, which is smaller than the preset amplitude value variance, to the total number of the amplitude value variances is larger than a preset ratio threshold.
As a further improvement of an embodiment of the present invention, acquiring the image brightness characteristic parameter corresponding to each partition according to the obtained envelope amplitude signal specifically includes:
the envelope amplitude signal corresponding to any one scanning point is expressed by env (t, sln),
the longitudinal energy spectrum density function psd (i) corresponding to any partition in the longitudinal direction is:
the variance v (i) of the longitudinal amplitude value corresponding to any partition in the longitudinal direction is as follows:
the method comprises the steps of obtaining two-dimensional data, wherein i is 0,1,2.. M in the longitudinal direction, M represents the number of longitudinal partition intervals in the two-dimensional data, t represents a scanning point, ti represents a first scanning point corresponding to the current partition interval, step0 represents the number of sampling points in each partition interval in the longitudinal direction, sln represents a scanning line, and lines represents the number of transverse scanning lines.
As a further improvement of an embodiment of the present invention, acquiring the image brightness characteristic parameter corresponding to each partition according to the obtained envelope amplitude signal specifically includes:
the envelope amplitude signal corresponding to any scanning point is represented by env (t, sln), the transverse energy spectrum density function corresponding to any transverse subarea is represented by psd (i),
then:
the variance v (i) of the longitudinal amplitude value corresponding to any one horizontal partition is as follows:
n, N represents the number of horizontal partition intervals in the two-dimensional data, sln represents a scanning line, li represents a first scanning line corresponding to the current partition interval, step1 represents the number of scanning lines in each partition interval in the horizontal direction, t represents scanning points, and samples represents the number of longitudinal sampling points.
As a further improvement of the embodiment of the present invention, the obtaining the first gain compensation coefficient corresponding to each sub-region according to the energy spectral density function specifically includes:
enabling a first gain compensation coefficient corresponding to any one partition to be equal to a quotient value of a maximum energy spectrum density function in the direction of the current partition and an energy spectrum density function corresponding to the current partition;
namely: the first gain compensation coefficient corresponding to any one of the partitions is represented by C1(i), then,
c1(i) ═ max (psd (i))/psd (i), where max (psd (i)) represents the maximum energy spectral density function in the direction in which the current partition is located.
As a further improvement of the embodiment of the present invention, the obtaining the second gain compensation coefficient corresponding to each partition according to the variance of the amplitude value and the first gain compensation coefficient specifically includes:
the second gain compensation coefficient corresponding to any one of the partitions is represented by C2(i), then,
where max (v (i)) represents the maximum variance of the amplitude value in the direction of the current partition, and min (v (i)) represents the minimum variance of the amplitude value in the direction of the current partition.
In order to achieve the above object, according to another embodiment of the present invention, there is provided an ultrasound image brightness adjustment system, including: the device comprises an acquisition module, a gain compensation module and a control module, wherein the acquisition module is used for acquiring an echo RF signal and respectively carrying out gain compensation on the RF signal in the transverse direction and/or the longitudinal direction through a preset gain compensation coefficient so as to form an envelope amplitude signal; the preset gain compensation coefficient comprises: presetting a transverse gain compensation coefficient and/or a longitudinal gain compensation coefficient;
the processing module is used for partitioning the two-dimensional data to form a plurality of partitions, wherein the partitions comprise transversely distributed transverse partitions and/or longitudinally distributed longitudinal partitions; acquiring image brightness characteristic parameters corresponding to each partition according to the obtained envelope amplitude signals, wherein the image brightness characteristic parameters comprise: an energy spectral density function, and a variance of amplitude values;
the analysis output module is used for analyzing the variance of the amplitude value corresponding to each interval, judging whether the ratio of the number of the variance of the amplitude value corresponding to each interval, which is smaller than the preset variance of the amplitude value, to the total number of the variance of the amplitude value is larger than a preset ratio threshold value or not,
if yes, obtaining a first gain compensation coefficient corresponding to each partition according to the energy spectrum density function, wherein the first gain compensation coefficient comprises: a first lateral gain compensation factor and/or a first longitudinal gain compensation factor; performing brightness compensation on the image by using the first gain compensation coefficient corresponding to each subarea, and displaying the image;
if not, acquiring a first gain compensation coefficient corresponding to each subarea according to the energy spectrum density function, and acquiring a second gain compensation coefficient corresponding to each subarea according to the amplitude value variance and the first gain compensation coefficient, wherein the second gain compensation coefficient comprises: a second lateral gain compensation factor and/or a second longitudinal gain compensation factor; and replacing a preset gain compensation coefficient with a second gain compensation coefficient corresponding to the current partition to perform gain compensation on the echo RF signal to form a new envelope amplitude signal until the ratio of the number of the amplitude value variances corresponding to each partition, which is smaller than the preset amplitude value variance, to the total number of the amplitude value variances is larger than a preset ratio threshold.
As a further improvement of an embodiment of the present invention, the processing module is specifically configured to:
the envelope amplitude signal corresponding to any one scanning point is expressed by env (t, sln),
the longitudinal energy spectrum density function psd (i) corresponding to any partition in the longitudinal direction is:
the variance v (i) of the longitudinal amplitude value corresponding to any partition in the longitudinal direction is as follows:
the method comprises the steps of obtaining two-dimensional data, wherein i is 0,1,2.. M in the longitudinal direction, M represents the number of longitudinal partition intervals in the two-dimensional data, t represents a scanning point, ti represents a first scanning point corresponding to the current partition interval, step0 represents the number of sampling points in each partition interval in the longitudinal direction, sln represents a scanning line, and lines represents the number of transverse scanning lines.
As a further improvement of an embodiment of the present invention, the processing module is specifically configured to:
the envelope amplitude signal corresponding to any scanning point is represented by env (t, sln), the transverse energy spectrum density function corresponding to any transverse subarea is represented by psd (i),
then:
the variance v (i) of the longitudinal amplitude value corresponding to any one horizontal partition is as follows:
n, N represents the number of horizontal partition intervals in the two-dimensional data, sln represents a scanning line, li represents a first scanning line corresponding to the current partition interval, step1 represents the number of scanning lines in each partition interval in the horizontal direction, t represents scanning points, and samples represents the number of longitudinal sampling points.
As a further improvement of an embodiment of the present invention, the analysis output module is specifically configured to:
enabling a first gain compensation coefficient corresponding to any one partition to be equal to a quotient value of a maximum energy spectrum density function in the direction of the current partition and an energy spectrum density function corresponding to the current partition;
namely: the first gain compensation coefficient corresponding to any one of the partitions is represented by C1(i), then,
c1(i) ═ max (psd (i))/psd (i), where max (psd (i)) represents the maximum energy spectral density function in the direction in which the current partition is located.
As a further improvement of an embodiment of the present invention, the analysis output module is further configured to:
the second gain compensation coefficient corresponding to any one of the partitions is represented by C2(i), then,
where max (v (i)) represents the maximum variance of the amplitude value in the direction of the current partition, and min (v (i)) represents the minimum variance of the amplitude value in the direction of the current partition.
Compared with the prior art, the ultrasonic image brightness adjusting method and the ultrasonic image brightness adjusting system acquire the gain compensation coefficient for adjusting the image brightness through the energy spectrum density function, and correct the acquired gain compensation coefficient according to the amplitude value variance, so that the image brightness can be adaptively adjusted according to specific requirements, the problem of image quality reduction caused by individual difference is effectively solved, the image performance is improved, the diagnosis effect is enhanced, and the user satisfaction is improved.
Drawings
FIG. 1 is a flowchart illustrating a method for adjusting brightness of an ultrasound image according to an embodiment of the present invention;
FIG. 2 is a block diagram of an ultrasound image brightness adjustment system according to an embodiment of the present invention.
Detailed Description
The present invention will be described in detail below with reference to embodiments shown in the drawings. These embodiments are not intended to limit the present invention, and structural, methodological, or functional changes made by those skilled in the art according to these embodiments are included in the scope of the present invention.
It should be noted that the present invention is mainly applied to an ultrasound device, and correspondingly, the object to be measured may be a tissue to be measured, which is not described in detail herein.
In the ultrasonic imaging process; transmitting a pulse signal into a tissue through a probe, wherein the pulse signal is reflected in the tissue to form an ultrasonic signal, the ultrasonic signal is converted into an electrical analog signal through different elements of a probe transducer, the electrical analog signal is amplified by a front-end amplification module, and the electrical analog signal is converted into a digital signal by an A/D (analog-to-digital) conversion module; the digital signals of different elements are synthesized into radio frequency signals through a beam synthesis module; after the radio frequency signal is subjected to RF filtering, the attenuation of the signal in the time direction is compensated through a time gain compensation module, the signal after gain compensation is sent to an orthogonal demodulation module for demodulation processing, the I/Q signal of the orthogonal demodulation result is sent to a subsequent imaging processing module, and the generated image is displayed and output.
Referring to fig. 1, the method for adjusting brightness of an ultrasound image according to the present invention includes: s1, acquiring echo RF signals, and performing gain compensation on the RF signals in the transverse direction and/or the longitudinal direction respectively through preset gain compensation coefficients to form envelope amplitude signals; the preset gain compensation coefficient comprises: presetting a transverse gain compensation coefficient and/or a longitudinal gain compensation coefficient;
the preset gain compensation coefficient is a gain compensation coefficient obtained in a conventional manner, for example: the preset gain compensation coefficient obtaining process is obtained through empirical values, is the prior art, and is not described herein any more, and the processing process of the preset gain compensation coefficient obtaining process does not affect the core invention part of the present invention. Correspondingly, the preset gain compensation coefficient can selectively comprise a preset transverse gain compensation coefficient and/or a preset longitudinal gain compensation coefficient; in the data processing process, the more types of the selected preset gain compensation coefficients, the more accurate the result is, and details are not described herein.
It should be noted that the invention uses the envelope amplitude signal as the basic parameter, which not only contains the important organization information of imaging, but also avoids the change of image brightness information caused by different parameter settings in the log compression link and the post-processing link.
S2, partitioning the two-dimensional data to form a plurality of partitions, wherein the partitions comprise transversely distributed transverse partitions and/or longitudinally distributed longitudinal partitions; acquiring image brightness characteristic parameters corresponding to each partition according to the obtained envelope amplitude signals, wherein the image brightness characteristic parameters comprise: an energy spectral density function, and a variance of amplitude values.
In the specific implementation of the invention, in order to make the obtained result more accurate, the brightness compensation is respectively carried out on the image in the transverse direction and the longitudinal direction; in a preferred embodiment of the present invention, the two-dimensional data is partitioned to form a plurality of partitions, and the two-dimensional data is further processed by the partitions, but in other embodiments of the present invention, the same processing may be performed on each two-dimensional data without considering the calculation amount, and details are not described here.
In a specific embodiment of the present invention, an envelope amplitude signal corresponding to any scanning point is represented by env (t, sln), and a longitudinal energy spectral density function psd (i) corresponding to any partition in the longitudinal direction is:
the variance v (i) of the longitudinal amplitude value corresponding to any partition in the longitudinal direction is as follows:
the method comprises the steps of obtaining two-dimensional data, wherein i is 0,1,2.. M in the longitudinal direction, M represents the number of longitudinal partition intervals in the two-dimensional data, t represents a scanning point, ti represents a first scanning point corresponding to the current partition interval, step0 represents the number of sampling points in each partition interval in the longitudinal direction, sln represents a scanning line, and lines represents the number of transverse scanning lines.
And (3) representing the transverse energy spectrum density function corresponding to any one transverse subarea by psd (i), then:
the variance v (i) of the longitudinal amplitude value corresponding to any one horizontal partition is as follows:
n, N represents the number of horizontal partition intervals in the two-dimensional data, sln represents a scanning line, li represents a first scanning line corresponding to the current partition interval, step1 represents the number of scanning lines in each partition interval in the horizontal direction, t represents scanning points, and samples represents the number of longitudinal sampling points.
It can be understood that, the two-dimensional data is divided into a plurality of partitions in the longitudinal direction and/or the transverse direction, and when the two-dimensional data is divided in the longitudinal direction, the number of sampling points in each partition can be the same or different; in a preferred embodiment of the present invention, the number of each inter-partition sampling point is determined according to accuracy and robustness, and can be modified according to image consistency requirements, applications, and the like, which is not described in detail herein.
Further, the method further comprises: s3, analyzing the amplitude value variance corresponding to each partition, and judging whether the ratio of the number of the amplitude value variances corresponding to each partition, which is smaller than the preset amplitude value variance, to the total number of the amplitude value variances is larger than a preset ratio threshold, if so, acquiring a first gain compensation coefficient corresponding to each partition according to an energy spectrum density function, wherein the first gain compensation coefficient comprises: a first lateral gain compensation factor and/or a first longitudinal gain compensation factor; performing brightness compensation on the image by using the first gain compensation coefficient corresponding to each subarea, and displaying the image; if not, acquiring a first gain compensation coefficient corresponding to each subarea according to the energy spectrum density function, and acquiring a second gain compensation coefficient corresponding to each subarea according to the amplitude value variance and the first gain compensation coefficient, wherein the second gain compensation coefficient comprises: a second lateral gain compensation factor and/or a second longitudinal gain compensation factor; and replacing a preset gain compensation coefficient with a second gain compensation coefficient corresponding to the current partition to perform gain compensation on the echo RF signal to form a new envelope amplitude signal until the proportion of the amplitude value variance smaller than the preset amplitude value variance is larger than a preset proportion threshold.
In the preferred embodiment of the present invention, stability of pdeof (probability reliability estimate of the estimated magnitude of the signal) is used to determine whether to correct the obtained first gain compensation coefficient, pdeof stably represents that the image has better brightness uniformity, in the specific embodiment of the present invention, statistical analysis is performed in partitions, whether pdeof is stable is determined by variance of magnitude values corresponding to each partition, specifically, variance of magnitude values in each partition is used as one of factors affecting the gain compensation coefficient, and the preset variance of magnitude values can be set according to an empirical value, or can be an average value or a weighted value of variance of magnitude values corresponding to each partition; the preset proportion threshold value can be adjusted and optimized according to different users, different application scenes and the like, and is specifically specified, in a specific embodiment of the invention, the value range is expressed by percentage and can be any value between (50% and 100%), and in a preferred embodiment of the invention, the preset proportion threshold value is 70%.
Specifically, the obtaining of the first gain compensation coefficient corresponding to each partition according to the energy spectral density function specifically includes: the first gain compensation coefficient corresponding to any one partition is equal to the quotient of the maximum energy spectrum density function in the direction of the current partition and the energy spectrum density function corresponding to the current partition;
namely: the first gain compensation coefficient corresponding to any one of the partitions is represented by C1(i), then,
c1(i) ═ max (psd (i))/psd (i), where max (psd (i)) represents the maximum energy spectral density function in the direction in which the current partition is located.
It is to be understood that when partitioning a vertical partition, i is 0,1,2.. M, M represents the number of partitions in the vertical partition in the two-dimensional data; when the data is divided into the horizontal partitions, i is 0,1,2.
The obtaining of the second gain compensation coefficient corresponding to each partition according to the amplitude value variance and the first gain compensation coefficient specifically includes: the second gain compensation coefficient corresponding to any one of the partitions is represented by C2(i), then,
wherein max (v (i)) represents the maximum variance of the amplitude value in the direction of the current partition, min (v (i)) represents the minimum variance of the amplitude value in the direction of the current partition, and when the partition is vertical, i is 0,1,2. When the data is divided into the horizontal partitions, i is 0,1,2.
For convenience of understanding, a specific example is described to facilitate understanding, the two-dimensional data is divided into 10 partitions in the horizontal direction and the longitudinal direction respectively, the preset proportion threshold is 70%, after calculation, the amplitude variance between 8 partitions in the horizontal direction is larger than the preset amplitude variance, and the amplitude variance between 5 partitions in the longitudinal direction is larger than the preset amplitude variance; therefore, the ratio of the lateral amplitude variance is 8/10 × 100 ═ 80%, and the ratio of the longitudinal amplitude variance is 5/10 × 100 ═ 50%, so that the image is luminance-compensated in the lateral direction by the first lateral gain compensation coefficient corresponding to each division, and displayed; in the longitudinal direction, replacing a preset gain compensation coefficient with a second gain compensation coefficient corresponding to each interval to perform gain compensation on the original echo RF signal to form a new envelope amplitude signal, and performing cyclic calculation until whether the ratio of the number of the amplitude value variances smaller than the preset amplitude value variance in the amplitude value variances corresponding to each interval to the total number of the amplitude value variances is larger than a preset ratio threshold value.
And further, sequentially carrying out brightness compensation on the graph by using the obtained gain compensation coefficients corresponding to the partitions by using methods of logarithmic compression, image post-processing, coordinate conversion and display so as to output and display the graph.
Referring to fig. 2, an ultrasound image brightness adjustment system according to an embodiment of the present invention includes: the device comprises an acquisition module 100, a processing module 200 and an analysis output module 300.
The obtaining module 100 is configured to obtain an echo RF signal, and perform gain compensation on the RF signal in a transverse direction and/or a longitudinal direction through a preset gain compensation coefficient to form an envelope amplitude signal; the preset gain compensation coefficient comprises: presetting a transverse gain compensation coefficient and/or a longitudinal gain compensation coefficient;
the preset gain compensation coefficient is a gain compensation coefficient obtained in a conventional manner, for example: the preset gain compensation coefficient obtaining process is obtained through empirical values, is the prior art, and is not described herein any more, and the processing process of the preset gain compensation coefficient obtaining process does not affect the core invention part of the present invention. Correspondingly, the preset gain compensation coefficient can selectively comprise a preset transverse gain compensation coefficient and/or a preset longitudinal gain compensation coefficient; in the data processing process, the more types of the selected preset gain compensation coefficients, the more accurate the result is, and details are not described herein.
It should be noted that the invention uses the envelope amplitude signal as the basic parameter, which not only contains the important organization information of imaging, but also avoids the change of image brightness information caused by different parameter settings in the log compression link and the post-processing link.
The processing module 200 is configured to partition the two-dimensional data to form a plurality of partitions, where the partitions include laterally distributed partitions and/or longitudinally distributed partitions; acquiring image brightness characteristic parameters corresponding to each partition according to the obtained envelope amplitude signals, wherein the image brightness characteristic parameters comprise: an energy spectral density function, and a variance of amplitude values.
In the specific implementation of the invention, in order to make the obtained result more accurate, the brightness compensation is respectively carried out on the image in the transverse direction and the longitudinal direction; in a preferred embodiment of the present invention, the processing module 200 partitions the two-dimensional data to form a plurality of partitions, and further partitions process the two-dimensional data, and certainly, in other embodiments of the present invention, the same processing may be performed on each two-dimensional data without considering the calculation amount, which is not described in detail herein.
In a specific embodiment of the present invention, the processing module 200 represents the envelope amplitude signal corresponding to any scanning point by env (t, sln), and a longitudinal energy spectral density function psd (i) corresponding to any partition in the longitudinal direction is:
the variance v (i) of the longitudinal amplitude value corresponding to any partition in the longitudinal direction is as follows:
the method comprises the steps of obtaining two-dimensional data, wherein i is 0,1,2.. M in the longitudinal direction, M represents the number of longitudinal partition intervals in the two-dimensional data, t represents a scanning point, ti represents a first scanning point corresponding to the current partition interval, step0 represents the number of sampling points in each partition interval in the longitudinal direction, sln represents a scanning line, and lines represents the number of transverse scanning lines.
The processing module 200 is further configured to represent a lateral energy spectral density function corresponding to any one of the lateral sub-regions in psd (i), and then:
the variance v (i) of the longitudinal amplitude value corresponding to any one horizontal partition is as follows:
n, N represents the number of horizontal partition intervals in the two-dimensional data, sln represents a scanning line, li represents a first scanning line corresponding to the current partition interval, step1 represents the number of scanning lines in each partition interval in the horizontal direction, t represents scanning points, and samples represents the number of longitudinal sampling points.
It can be understood that, the two-dimensional data is divided into a plurality of partitions in the longitudinal direction and/or the transverse direction, and when the two-dimensional data is divided in the longitudinal direction, the number of sampling points in each partition can be the same or different; in a preferred embodiment of the present invention, the number of each inter-partition sampling point is determined according to accuracy and robustness, and can be modified according to image consistency requirements, applications, and the like, which is not described in detail herein.
The analysis output module 300 is configured to analyze the variance of the amplitude value corresponding to each partition, determine whether a ratio of the number of the variance of the amplitude value corresponding to each partition, which is smaller than a preset variance of the amplitude value, to a total number of the variance of the amplitude value is greater than a preset ratio threshold, and if so, obtain a first gain compensation coefficient corresponding to each partition according to an energy spectral density function, where the first gain compensation coefficient includes: a first lateral gain compensation factor and/or a first longitudinal gain compensation factor; performing brightness compensation on the image by using the first gain compensation coefficient corresponding to each subarea, and displaying the image; if not, acquiring a first gain compensation coefficient corresponding to each subarea according to the energy spectrum density function, and acquiring a second gain compensation coefficient corresponding to each subarea according to the amplitude value variance and the first gain compensation coefficient, wherein the second gain compensation coefficient comprises: a second lateral gain compensation factor and/or a second longitudinal gain compensation factor; and replacing a preset gain compensation coefficient with a second gain compensation coefficient corresponding to the current partition to perform gain compensation on the echo RF signal to form a new envelope amplitude signal until the proportion of the amplitude value variance smaller than the preset amplitude value variance is larger than a preset proportion threshold.
In the preferred embodiment of the present invention, stability of pdeof (probability reliability estimate of the estimated magnitude of the signal) is used to determine whether to correct the obtained first gain compensation coefficient, pdeof stably represents that the image has better brightness uniformity, in the specific embodiment of the present invention, statistical analysis is performed in partitions, whether pdeof is stable is determined by variance of magnitude values corresponding to each partition, specifically, variance of magnitude values in each partition is used as one of factors affecting the gain compensation coefficient, and the preset variance of magnitude values can be set according to an empirical value, or can be an average value or a weighted value of variance of magnitude values corresponding to each partition; the preset proportion threshold value can be adjusted and optimized according to different users, different application scenes and the like, and is specifically specified, in a specific embodiment of the invention, the value range is expressed by percentage and can be any value between (50% and 100%), and in a preferred embodiment of the invention, the preset proportion threshold value is 70%.
Specifically, the obtaining, by the analysis output module 300, the first gain compensation coefficient corresponding to each partition according to the energy spectrum density function specifically includes: enabling a first gain compensation coefficient corresponding to any one partition to be equal to a quotient value of a maximum energy spectrum density function in the direction of the current partition and an energy spectrum density function corresponding to the current partition;
namely: the first gain compensation coefficient corresponding to any one of the partitions is represented by C1(i), then,
c1(i) ═ max (psd (i))/psd (i), where max (psd (i)) represents the maximum energy spectral density function in the direction in which the current partition is located.
It is to be understood that when partitioning a vertical partition, i is 0,1,2.. M, M represents the number of partitions in the vertical partition in the two-dimensional data; when the data is divided into the horizontal partitions, i is 0,1,2.
The obtaining, by the analysis output module 300, the second gain compensation coefficient corresponding to each partition according to the variance of the amplitude value and the first gain compensation coefficient specifically includes: the second gain compensation coefficient corresponding to any one of the partitions is represented by C2(i), then,
wherein max (v (i)) represents the maximum variance of the amplitude value in the direction of the current partition, min (v (i)) represents the minimum variance of the amplitude value in the direction of the current partition, and when the partition is vertical, i is 0,1,2. When the data is divided into the horizontal partitions, i is 0,1,2.
Further, the analysis output module 300 is further configured to perform brightness compensation on the graph by sequentially using the obtained gain compensation coefficients corresponding to the respective partitions through log compression, image post-processing, coordinate transformation, and display, so as to perform output display.
In summary, the ultrasound image brightness adjustment method and the ultrasound image brightness adjustment system of the present invention obtain the gain compensation coefficient for adjusting the image brightness through the energy spectrum density function, and correct the obtained gain compensation coefficient according to the variance of the amplitude value, so that the image brightness can be adaptively adjusted according to specific requirements, the problem of image quality reduction caused by individual differences is effectively improved, the image performance is improved, the diagnostic effect is enhanced, and the user satisfaction is improved.
It can be clearly understood by those skilled in the art that, for convenience and brevity of description, the specific working process of the system described above may refer to the corresponding process in the foregoing method embodiment, and is not described herein again.
In the several embodiments provided in this application, it should be understood that the disclosed system, and method may be implemented in other ways. For example, the system embodiments described above are merely illustrative, and for example, the division of the modules is merely a logical division, and in actual implementation, there may be another division, for example, multiple modules or components may be combined or integrated into another system, or some features may be omitted, or not implemented. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, systems or modules, and may be in an electrical, mechanical or other form.
The modules described as separate parts may or may not be physically separate, and parts displayed as modules may or may not be physical modules, may be located in one place, or may be distributed on a plurality of network modules. Some or all of the modules may be selected according to actual needs to achieve the purpose of the embodiment.
In addition, each functional module in the embodiments of the present application may be integrated into one processing module, or each module may exist alone physically, or 2 or more modules may be integrated into one module. The integrated module can be realized in a hardware form, and can also be realized in a form of hardware and a software functional module.
The integrated module implemented in the form of a software functional module may be stored in a computer-readable storage medium. The software functional module is stored in a storage medium and includes several instructions for causing a computer system (which may be a personal computer, a server, or a network system) or a processor (processor) to execute some steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: various media capable of storing program codes, such as a usb disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solutions of the present application, and not for limiting the same; although the present application has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may be modified or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present application.
Claims (10)
1. An ultrasound image brightness adjustment method, characterized in that the method comprises the following steps:
s1, acquiring echo RF signals, and performing gain compensation on the RF signals in the transverse direction and/or the longitudinal direction respectively through preset gain compensation coefficients to form envelope amplitude signals; the preset gain compensation coefficient comprises: presetting a transverse gain compensation coefficient and/or a longitudinal gain compensation coefficient;
s2, partitioning the two-dimensional data to form a plurality of partitions, wherein the partitions comprise transversely distributed transverse partitions and/or longitudinally distributed longitudinal partitions; acquiring image brightness characteristic parameters corresponding to each partition according to the obtained envelope amplitude signals, wherein the image brightness characteristic parameters comprise: an energy spectral density function, and a variance of amplitude values;
s3, analyzing the variance of the amplitude value corresponding to each partition, judging whether the ratio of the number of the variance of the amplitude value corresponding to each partition, which is smaller than the preset variance of the amplitude value, to the total number of the variance of the amplitude value is larger than the preset ratio threshold value,
if yes, obtaining a first gain compensation coefficient corresponding to each partition according to the energy spectrum density function, wherein the first gain compensation coefficient comprises: a first lateral gain compensation factor and/or a first longitudinal gain compensation factor; performing brightness compensation on the image by using the first gain compensation coefficient corresponding to each subarea, and displaying the image;
if not, acquiring a first gain compensation coefficient corresponding to each subarea according to the energy spectrum density function, and acquiring a second gain compensation coefficient corresponding to each subarea according to the amplitude value variance and the first gain compensation coefficient, wherein the second gain compensation coefficient comprises: a second lateral gain compensation factor and/or a second longitudinal gain compensation factor; and replacing a preset gain compensation coefficient with a second gain compensation coefficient corresponding to the current partition to perform gain compensation on the echo RF signal to form a new envelope amplitude signal, and returning to the step S1 until the ratio of the number of the amplitude value variances smaller than the preset amplitude value in the amplitude value variances corresponding to each partition to the total number of the amplitude value variances is greater than a preset ratio threshold.
2. The method for adjusting brightness of an ultrasound image according to claim 1, wherein obtaining the image brightness characteristic parameter corresponding to each partition according to the obtained envelope magnitude signal specifically includes:
the envelope amplitude signal corresponding to any one scanning point is expressed by env (t, sln),
the longitudinal energy spectrum density function psd (i) corresponding to any partition in the longitudinal direction is:
the variance v (i) of the longitudinal amplitude value corresponding to any partition in the longitudinal direction is as follows:
the method comprises the steps of obtaining two-dimensional data, wherein i is 0,1,2.. M in the longitudinal direction, M represents the number of longitudinal partition intervals in the two-dimensional data, t represents a scanning point, ti represents a first scanning point corresponding to the current partition interval, step0 represents the number of sampling points in each partition interval in the longitudinal direction, sln represents a scanning line, and lines represents the number of transverse scanning lines.
3. The method for adjusting brightness of an ultrasound image according to claim 1, wherein obtaining the image brightness characteristic parameter corresponding to each partition according to the obtained envelope magnitude signal specifically includes:
the envelope amplitude signal corresponding to any scanning point is represented by env (t, sln), the transverse energy spectrum density function corresponding to any transverse subarea is represented by psd (i),
then:
the variance v (i) of the longitudinal amplitude value corresponding to any one horizontal partition is as follows:
n, N represents the number of horizontal partition intervals in the two-dimensional data, sln represents a scanning line, li represents a first scanning line corresponding to the current partition interval, step1 represents the number of scanning lines in each partition interval in the horizontal direction, t represents scanning points, and samples represents the number of longitudinal sampling points.
4. The method according to claim 2 or 3, wherein the obtaining of the first gain compensation coefficient corresponding to each sub-region according to the energy spectral density function specifically comprises:
enabling a first gain compensation coefficient corresponding to any one partition to be equal to a quotient value of a maximum energy spectrum density function in the direction of the current partition and an energy spectrum density function corresponding to the current partition;
namely: the first gain compensation coefficient corresponding to any one of the partitions is represented by C1(i), then,
c1(i) ═ max (psd (i))/psd (i), where max (psd (i)) represents the maximum energy spectral density function in the direction in which the current partition is located.
5. The method for adjusting ultrasound image brightness according to claim 4, wherein obtaining the second gain compensation coefficient corresponding to each partition according to the variance of the amplitude value and the first gain compensation coefficient specifically includes:
the second gain compensation coefficient corresponding to any one of the partitions is represented by C2(i), then,
where max (v (i)) represents the maximum variance of the amplitude value in the direction of the current partition, and min (v (i)) represents the minimum variance of the amplitude value in the direction of the current partition.
6. An ultrasound image brightness adjustment system, the system comprising:
the device comprises an acquisition module, a gain compensation module and a control module, wherein the acquisition module is used for acquiring an echo RF signal and respectively carrying out gain compensation on the RF signal in the transverse direction and/or the longitudinal direction through a preset gain compensation coefficient so as to form an envelope amplitude signal; the preset gain compensation coefficient comprises: presetting a transverse gain compensation coefficient and/or a longitudinal gain compensation coefficient;
the processing module is used for partitioning the two-dimensional data to form a plurality of partitions, and the partitions comprise transversely distributed partitions and/or longitudinally distributed partitions; acquiring image brightness characteristic parameters corresponding to each partition according to the obtained envelope amplitude signals, wherein the image brightness characteristic parameters comprise: an energy spectral density function, and a variance of amplitude values;
the analysis output module is used for analyzing the variance of the amplitude value corresponding to each interval, judging whether the ratio of the number of the variance of the amplitude value corresponding to each interval, which is smaller than the preset variance of the amplitude value, to the total number of the variance of the amplitude value is larger than a preset ratio threshold value or not,
if yes, obtaining a first gain compensation coefficient corresponding to each partition according to the energy spectrum density function, wherein the first gain compensation coefficient comprises: a first lateral gain compensation factor and/or a first longitudinal gain compensation factor; performing brightness compensation on the image by using the first gain compensation coefficient corresponding to each subarea, and displaying the image;
if not, acquiring a first gain compensation coefficient corresponding to each subarea according to the energy spectrum density function, and acquiring a second gain compensation coefficient corresponding to each subarea according to the amplitude value variance and the first gain compensation coefficient, wherein the second gain compensation coefficient comprises: a second lateral gain compensation factor and/or a second longitudinal gain compensation factor; replacing a preset gain compensation coefficient with a second gain compensation coefficient corresponding to the current partition to perform gain compensation on the echo RF signal to form a new envelope amplitude signal;
the processing module is also used for circularly executing and acquiring the image brightness characteristic parameters corresponding to each subarea after the analysis output module forms a new envelope amplitude signal; the analysis output module is further used for circularly analyzing the amplitude value variance corresponding to each interval after the processing module circularly executes the image brightness characteristic parameters, until the ratio of the number of the amplitude value variances corresponding to each interval, which is smaller than the preset amplitude value variance, to the total number of the amplitude value variances is larger than the preset ratio threshold.
7. The ultrasound image brightness adjustment system of claim 6,
the processing module is specifically configured to:
the envelope amplitude signal corresponding to any one scanning point is expressed by env (t, sln),
the longitudinal energy spectrum density function psd (i) corresponding to any partition in the longitudinal direction is:
the variance v (i) of the longitudinal amplitude value corresponding to any partition in the longitudinal direction is as follows:
the method comprises the steps of obtaining two-dimensional data, wherein i is 0,1,2.. M in the longitudinal direction, M represents the number of longitudinal partition intervals in the two-dimensional data, t represents a scanning point, ti represents a first scanning point corresponding to the current partition interval, step0 represents the number of sampling points in each partition interval in the longitudinal direction, sln represents a scanning line, and lines represents the number of transverse scanning lines.
8. The ultrasound image brightness adjustment system of claim 6,
the processing module is specifically configured to:
the envelope amplitude signal corresponding to any scanning point is represented by env (t, sln), the transverse energy spectrum density function corresponding to any transverse subarea is represented by psd (i),
then:
the variance v (i) of the longitudinal amplitude value corresponding to any one horizontal partition is as follows:
n, N represents the number of horizontal partition intervals in the two-dimensional data, sln represents a scanning line, li represents a first scanning line corresponding to the current partition interval, step1 represents the number of scanning lines in each partition interval in the horizontal direction, t represents scanning points, and samples represents the number of longitudinal sampling points.
9. The ultrasound image brightness adjustment system according to claim 7 or 8,
the analysis output module is specifically configured to:
enabling a first gain compensation coefficient corresponding to any one partition to be equal to a quotient value of a maximum energy spectrum density function in the direction of the current partition and an energy spectrum density function corresponding to the current partition;
namely: the first gain compensation coefficient corresponding to any one of the partitions is represented by C1(i), then,
c1(i) ═ max (psd (i))/psd (i), where max (psd (i)) represents the maximum energy spectral density function in the direction in which the current partition is located.
10. The ultrasound image brightness adjustment system of claim 9,
the analysis output module is further configured to:
the second gain compensation coefficient corresponding to any one of the partitions is represented by C2(i), then,
where max (v (i)) represents the maximum variance of the amplitude value in the direction of the current partition, and min (v (i)) represents the minimum variance of the amplitude value in the direction of the current partition.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810873293.2A CN108877902B (en) | 2018-08-02 | 2018-08-02 | Ultrasonic image brightness adjusting method and adjusting system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810873293.2A CN108877902B (en) | 2018-08-02 | 2018-08-02 | Ultrasonic image brightness adjusting method and adjusting system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108877902A CN108877902A (en) | 2018-11-23 |
| CN108877902B true CN108877902B (en) | 2020-12-04 |
Family
ID=64307648
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810873293.2A Active CN108877902B (en) | 2018-08-02 | 2018-08-02 | Ultrasonic image brightness adjusting method and adjusting system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108877902B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110148090B (en) * | 2019-01-19 | 2023-04-07 | 深圳蓝影医学科技股份有限公司 | Automatic optimization method and device for B-type image gain |
| CN110327076B (en) * | 2019-07-05 | 2022-08-16 | 深圳开立生物医疗科技股份有限公司 | Blood flow gain adjusting method, device, equipment and readable storage medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6059727A (en) * | 1995-06-15 | 2000-05-09 | The Regents Of The University Of Michigan | Method and apparatus for composition and display of three-dimensional image from two-dimensional ultrasound scan data |
| CN101888876A (en) * | 2007-10-01 | 2010-11-17 | 因赛泰克有限公司 | Motion compensated image-guided focused ultrasound therapy system |
| CN102499711A (en) * | 2011-09-28 | 2012-06-20 | 无锡祥生医学影像有限责任公司 | Three-dimensional or four-dimensional automatic ultrasound image optimization and adjustment method |
| CN105528764A (en) * | 2015-12-01 | 2016-04-27 | 沈阳东软医疗系统有限公司 | Image brightness adjustment method, installation and device |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5762675A (en) * | 1980-10-03 | 1982-04-15 | Matsushita Electric Ind Co Ltd | Magnetic recorder and reproducer |
| US6679844B2 (en) * | 2002-06-20 | 2004-01-20 | Acuson Corporation | Automatic gain compensation for multiple mode or contrast agent imaging |
-
2018
- 2018-08-02 CN CN201810873293.2A patent/CN108877902B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6059727A (en) * | 1995-06-15 | 2000-05-09 | The Regents Of The University Of Michigan | Method and apparatus for composition and display of three-dimensional image from two-dimensional ultrasound scan data |
| CN101888876A (en) * | 2007-10-01 | 2010-11-17 | 因赛泰克有限公司 | Motion compensated image-guided focused ultrasound therapy system |
| CN102499711A (en) * | 2011-09-28 | 2012-06-20 | 无锡祥生医学影像有限责任公司 | Three-dimensional or four-dimensional automatic ultrasound image optimization and adjustment method |
| CN105528764A (en) * | 2015-12-01 | 2016-04-27 | 沈阳东软医疗系统有限公司 | Image brightness adjustment method, installation and device |
Non-Patent Citations (1)
| Title |
|---|
| 医学超声成像中自适应波束形成方法研究;许梦玲;《中国博士学位论文全文数据库》;20150731;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108877902A (en) | 2018-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8403850B2 (en) | Rapid two/three-dimensional sector strain imaging | |
| JP2002526141A (en) | Adaptive Cross Section Calculation Method Using Statistical Discriminant Characteristics | |
| CN106709967B (en) | Endoscopic imaging algorithm and control system | |
| JP5079177B2 (en) | Ultrasonic observation apparatus, operation method of ultrasonic observation apparatus, and operation program of ultrasonic observation apparatus | |
| CN110378982B (en) | Reconstructed image processing method, device, equipment and storage medium | |
| CN108877902B (en) | Ultrasonic image brightness adjusting method and adjusting system | |
| CN109754397B (en) | Blood vessel section center path extraction method and device, terminal equipment and storage medium | |
| CN108648164A (en) | Ultrasound image optimization method, apparatus and computer readable storage medium | |
| US20100331698A1 (en) | Ultrasonic diagnostic apparatus, ultrasonic elasticity information processing method and ultrasonic elasticity information processing program | |
| KR20130006011A (en) | Apparatus for processing an ultrasound muscle image | |
| CN111513754A (en) | Ultrasonic imaging equipment and quality evaluation method of ultrasonic image | |
| EP3639752A1 (en) | Analyzing apparatus and analyzing program | |
| GB2461991A (en) | Positron emission tomography data with time activity derived framing intervals. | |
| CN108366784A (en) | The working procedure of ultrasound observation apparatus, the working method of ultrasound observation apparatus and ultrasound observation apparatus | |
| US20220142615A1 (en) | Ultrasonic signal processing method and apparatus, device and storage medium | |
| CN115299986A (en) | Ultrasonic imaging equipment and ultrasonic inspection method thereof | |
| CN111383323B (en) | Ultrasonic imaging method and system and ultrasonic image processing method and system | |
| KR101627319B1 (en) | medical image processor and method thereof for medical diagnosis | |
| JP5931414B2 (en) | Ultrasonic diagnostic apparatus, image generation method, and image processing apparatus | |
| CN117152181B (en) | Tumor image segmentation method, device, electronic equipment and readable storage medium | |
| CN107736895B (en) | Protocol parameter configuration method, device and terminal | |
| KR102204448B1 (en) | Method for calculating oleic acid content of domestic animals from ultrasonic image and electronic device thereof | |
| CN111882515B (en) | Ultrasonic signal processing method, ultrasonic signal processing apparatus, and storage medium | |
| RU2795207C1 (en) | Method and device for processing ultrasonic signals, device and information media | |
| TW201732288A (en) | Calibration method of the quantified anechoic-area and hyperechoic-foci features for images acquired from different ultrasound systems |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province Patentee after: Feiyinuo Technology (Suzhou) Co.,Ltd. Address before: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province Patentee before: VINNO TECHNOLOGY (SUZHOU) Co.,Ltd. Address after: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province Patentee after: Feiyinuo Technology Co.,Ltd. Address before: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province Patentee before: Feiyinuo Technology (Suzhou) Co.,Ltd. |
|
| CP01 | Change in the name or title of a patent holder |