Photoacoustic mammary gland image reconstruction method and system
Technical Field
The invention relates to the technical field of breast imaging, in particular to a photoacoustic breast image reconstruction method and a photoacoustic breast image reconstruction system.
Background
The breast diseases are diseases derived from related tissues of the breast such as breast glands, fat, lymph, blood vessels, nipples and the like, and comprise inflammatory diseases of the breast, benign lesions of the breast, malignant tumors of the breast, congenital dysplasia, male breast development and the like. There are many methods for examining breast diseases, including molybdenum target X-ray imaging, nuclear magnetic resonance imaging, ultrasonic imaging, infrared imaging, photoacoustic imaging, etc.
The existing photoacoustic breast formation image reconstruction method adopts image acquisition, image cutting and image analysis methods, and although the method can solve the defects of other imaging technology in accuracy and time, the method directly performs acquisition, cutting and analysis on data, and has the following problems that a large amount of interference or useless information exists in the acquired image, and the image information needs to be processed firstly. Although the noise can be reduced by adopting the existing filtering method, the tumor signal is easily filtered, and the imaging accuracy is influenced.
Disclosure of Invention
The invention provides a photoacoustic breast image reconstruction method and system for solving the technical problems.
The invention is realized by the following technical scheme:
an image reconstruction method of a photoacoustic breast, comprising the steps of:
A. acquiring multi-channel ultrasonic information;
B. respectively carrying out time-frequency transformation on each path of ultrasonic information, and transforming the ultrasonic information to a time-frequency domain;
C. analyzing the time-frequency diagram and locking the position of the focus;
D. attenuating positions except the focus position in the time-frequency diagram;
E. inversely transforming the time-frequency domain information into a time domain to obtain noise-reduced photoacoustic information;
F. and cutting and analyzing the photoacoustic information to obtain a mammary gland imaging image.
The tumor photoacoustic signal is characterized more obviously in the time-frequency domain than in the surrounding normal tissue, i.e. the frequency component and the amplitude are higher in a short time in the time-frequency domain. According to the scheme, based on a time-frequency analysis method, the ultrasonic information is subjected to time-frequency transformation, the ultrasonic information is converted into a time-frequency domain, the focus tumor position is locked, other positions are attenuated, and then the focus tumor position is converted into the time domain, so that the signal-to-noise ratio of a signal with a reconstruction signal can be reduced, the signal-to-noise ratio of a breast tumor signal is improved, useful information of a tumor is not lost while noise is reduced, and the accuracy of breast image reconstruction is improved.
An image reconstruction system of a photoacoustic breast, comprising:
the data acquisition module is used for acquiring human breast ultrasonic information;
a data processing module for processing the image by adopting the method of the steps A-E;
and the data cutting module is used for imaging the photoacoustic information, judging the position of the focus and cutting data of the position of the focus and the surrounding area of the focus.
Compared with the prior art, the invention has the following advantages and beneficial effects:
1. the method is based on a time-frequency analysis method, firstly carries out time-frequency transformation on ultrasonic information, converts the ultrasonic information into a time-frequency domain, locks the position of the focus tumor, attenuates other positions and then transforms the focus tumor into the time domain, can reduce the signal-to-noise ratio of signals with reconstruction, improves the signal-to-noise ratio of breast tumor signals, does not lose useful information of the tumor while reducing noise, and improves the accuracy of breast image reconstruction.
Drawings
The accompanying drawings, which are included to provide a further understanding of the embodiments of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the principles of the invention.
Fig. 1 is a flow chart of a photoacoustic breast image reconstruction method;
FIG. 2 is a result of image reconstruction in a prior art method;
FIG. 3 shows the image reconstruction result of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to examples and accompanying drawings, and the exemplary embodiments and descriptions thereof are only used for explaining the present invention and are not meant to limit the present invention.
Example 1
The embodiment discloses an image reconstruction method of a photoacoustic breast, as shown in fig. 1, including the following steps:
A. acquiring multi-channel ultrasonic information;
B. respectively carrying out time-frequency transformation on each path of ultrasonic information, and transforming the ultrasonic information to a time-frequency domain;
C. analyzing the time-frequency diagram, and locking the focus position to keep effective information in the image;
D. attenuating positions except the focus position in the time-frequency diagram, and improving the contrast between the focus position and other positions;
E. inversely transforming the time-frequency domain information into a time domain to obtain noise-reduced photoacoustic information;
F. and imaging and cutting the photoacoustic information to obtain a mammary gland imaging image.
The method is based on a time-frequency analysis method, a focus area is locked on a time-frequency domain, other areas are attenuated, focus information is locked on the time-frequency domain, and accuracy of focus information analysis is improved; the method has the advantages that other regions are attenuated, the effect of focus region images cannot be damaged, the image characteristic information of the focus regions is kept, and the actual effect of the focus region images is not damaged.
Example 2
Based on the above method, the present embodiment discloses a specific implementation manner, that is, the method includes the following steps:
A. in the step, the array sensor is adopted to realize data acquisition, so that the data acquisition amount is increased, and the comprehensiveness of data information acquisition is enhanced;
B. respectively performing time-frequency transformation on each path of image information by adopting methods such as short-time fractional Fourier transform, S transformation and the like, and transforming the image information to a time-frequency domain;
C. the time is taken as a horizontal axis and the frequency is taken as a vertical axis in a time-frequency domain, the tumor photoacoustic signal is more obvious in the time-frequency domain, and a region which has high frequency components and high amplitude and is in an aggregation state is marked as a focus position;
D. attenuating positions except the focus position in the time-frequency diagram, specifically, attenuating the time-frequency intensity values except the focus position in the time-frequency diagram by multiplying a weight coefficient;
E. inversely transforming the time-frequency domain information into a time domain to obtain noise-reduced photoacoustic information;
F. and imaging the photoacoustic information by adopting a delay superposition or filtering back projection algorithm, cutting the acquired data according to the position of the suspected focus observed by the imaging graph and the position of the focus, and cutting the data to obtain the data of the suspected focus and the peripheral area of the suspected focus so as to obtain a mammary gland imaging graph.
Based on this, the system for realizing the method comprises a data acquisition module, a data processing module and a data cutting module. The data acquisition module adopts the existing array detector; the data processing module performs time domain-time frequency domain-time domain transformation on the image according to the method of the steps A-E in the embodiment 1 or 2; the data cutting module images the photoacoustic information, judges the focus position and cuts out the focus position and the data of the surrounding area of the focus position. The array detector is used as a part of the existing photoacoustic breast imager, and the photoacoustic breast imager, the array structure and the data acquisition process and method thereof are not described in detail at this time. To facilitate data analysis, the system may further include an image analysis system that quantitatively images the cropped data to obtain the hemoglobin concentration and oxygen saturation of the lesion for diagnosis by a physician.
The photoacoustic image reconstruction is carried out on the mammary tissue of the same patient, the image obtained by the existing method is shown in fig. 2, the image obtained by the method and the system of the invention is shown in fig. 3, and the imaging accuracy is higher compared with the part in the circle of fig. 2 and the corresponding part of fig. 3.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are merely exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.