Imaging method, imaging device, electronic device and storage medium

Abstract

This application discloses an imaging method, imaging device, electronic device, and computer storage medium. The imaging method includes: for each preset frequency within a preset frequency range; transmitting radio frequency pulses of the preset frequencies to a target area, the radio frequency pulses being used to excite water and fat within the target area to a saturated state; after the radio frequency pulse transmission ends, acquiring saturation signals corresponding to the preset frequencies; comparing all saturation signals corresponding to the preset frequencies with pre-acquired initial signals to determine a mixed signal of water and fat within the preset frequency range; calculating a first signal of water; determining a second signal based on the mixed signal and the first signal; and generating a fat peak map of the target area based on the spectral distribution of the second signal, where the second signal is the fat signal. This method can directly perform imaging based on the spectral distribution of the fat signal, thereby improving the spatial resolution of the image.

Description

Technical Field

[0001] The present application belongs to the field of medical imaging technology, and in particular, relates to an imaging method, an imaging device, an electronic device, and a computer-readable storage medium. Background Art

[0002] Fat is an important component of the human body and is involved in many life activities of the human body. The fat metabolism in different parts of the body can reflect the health status of different aspects of the human body. Therefore, fat imaging technology can provide effective auxiliary information for doctors' diagnosis.

[0003] Existing fat imaging technologies mainly include magnetic resonance spectroscopy imaging methods. However, magnetic resonance spectroscopy imaging methods have limited spatial resolution and can only obtain overall information of centimeter-level voxels, that is, the imaging spatial resolution is low. Summary of the invention

[0004] The present application provides an imaging method, an imaging device, an electronic device and a computer-readable storage medium, which can solve the problem of low imaging spatial resolution in existing magnetic resonance spectroscopy imaging methods.

[0005] In a first aspect, the present application provides an imaging method, comprising:

[0006] For each preset frequency within the preset frequency range: emitting a radio frequency pulse of the preset frequency to the target area, the radio frequency pulse being used to stimulate water and fat in the target area so that the water and fat in the target area are stimulated to a saturated state;

[0007] After the transmission of the radio frequency pulse is completed, a saturation signal corresponding to the preset frequency is collected;

[0008] Comparing the saturated signals corresponding to all the preset frequencies with the pre-collected initial signals to determine the mixed signal of the water and the fat within the preset frequency range;

[0009] Calculating the first signal of the water;

[0010] A second signal is determined based on the mixed signal and the first signal, and a fat peak map of the target area is generated based on the frequency spectrum distribution of the second signal, wherein the second signal is a signal of the fat.

[0011] In a second aspect, the present application provides an imaging device, comprising:

[0012] an excitation module, for each preset frequency within the preset frequency range: emitting a radio frequency pulse of the preset frequency to a target area, wherein the radio frequency pulse is used to excite water and fat in the target area, so that the water and fat in the target area are excited to a saturated state;

[0013] An acquisition module, used for acquiring a saturation signal corresponding to the preset frequency after the transmission of the radio frequency pulse is completed;

[0014] A determination module, used for comparing the saturated signals corresponding to all the preset frequencies with the pre-collected initial signals to determine a mixed signal of the water and the fat within the preset frequency range;

[0015] A calculation module, used for calculating the first signal of water based on the spectrum distribution characteristics of water molecules;

[0016] The imaging module is used to determine a second signal based on the mixed signal and the first signal, and to generate a fat peak map of the target area based on the spectrum distribution of the second signal, wherein the second signal is a signal of the fat.

[0017] In a third aspect, the present application provides an electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of the method of the first aspect when executing the computer program.

[0018] In a fourth aspect, the present application provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the method of the first aspect are implemented.

[0019] In a fifth aspect, the present application provides a computer program product, which includes a computer program. When the computer program is executed by one or more processors, it implements the steps of the method of the first aspect.

[0020] Compared with the prior art, the present application has the following beneficial effects: by exciting and saturating the water and fat in the target area with radio frequency pulses, a mixed signal of water and fat can be obtained; then the first signal of water can be calculated based on the water molecule saturation model; after obtaining the mixed signal and the first signal, the spectrum distribution of the first signal can be separated from the spectrum distribution of the mixed signal, thereby obtaining the spectrum distribution of the fat signal, that is, the spectrum distribution of the second signal, and finally the spectrum distribution of the second signal is used to image the target area to generate a corresponding fat peak map, so as to provide auxiliary information for doctors to diagnose patients. In this method, a conventional spin echo type magnetic resonance signal acquisition method can be used to obtain a mixed signal, and by processing the mixed signal, a millimeter-level image can be generated, thereby improving the spatial resolution of the spectrum image. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0022] Figure 1 It is a schematic diagram of a flow chart of an imaging method provided in an embodiment of the present application;

[0023] Figure 2 is a schematic diagram of obtaining the frequency spectrum distribution of a fat signal provided in an embodiment of the present application;

[0024] Figure 3 It is a fat peak spectrum of fat in the cone region of an abdominal cross section at a chemical shift of -2.6 ppm provided in an embodiment of the present application;

[0025] Figure 4 The present invention provides a fat peak spectrum of the fat in the cone region of the abdominal cross section at a chemical shift of -3.4 ppm;

[0026] Figure 5 is a schematic structural diagram of an imaging device provided in an embodiment of the present application;

[0027] Figure 6 It is a schematic diagram of the structure of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0028] In the following description, specific details such as specific system structures, technologies, etc. are provided for the purpose of illustration rather than limitation, so as to provide a thorough understanding of the embodiments of the present application. However, it should be clear to those skilled in the art that the present application may also be implemented in other embodiments without these specific details. In other cases, detailed descriptions of well-known systems, devices, circuits, and methods are omitted to prevent unnecessary details from obstructing the description of the present application.

[0029] In the prior art, magnetic resonance spectroscopy imaging can only obtain the overall information of large voxels at the centimeter level, resulting in a low spatial resolution of imaging. Based on this, the present application proposes an imaging method, which can remove the spectral distribution of the first signal (water signal) from the spectral distribution of the mixed signal (water and fat signals) to obtain the spectral distribution of the second signal (fat signal). Based on the spectral distribution of the second signal, the target area can be further imaged to obtain a fat peak map with a high spatial resolution.

[0030] The imaging method provided in the embodiment of the present application can be applied to a magnetic resonance device, which can have an embedded processor to implement the imaging method of the present application, or an external processor and an electronic device such as a computer to implement the imaging method of the present application. Specifically, the embodiment of the present application does not impose any restrictions on the specific type of electronic device.

[0031] In order to illustrate the technical solution proposed in this application, a specific embodiment is provided below for illustration.

[0032] Figure 1 A schematic flow chart of the imaging method provided by the present application is shown, and the imaging method comprises:

[0033] Step 110: for each preset frequency within the preset frequency range: transmitting a radio frequency pulse of the preset frequency to the target area, so that the water and fat in the target area are stimulated to a saturated state.

[0034] The target area is generally the cross-sectional area of ​​the human body to be imaged, such as the cross-sectional area of ​​the human abdomen; of course, it can also be the area formed by a certain part within the cross-sectional area, such as the area where the cone is located within the cross-sectional area of ​​the human abdomen. To obtain the mixed signal of water and fat in the target area, radio frequency pulses can be used within a specific frequency range to excite the water and fat in the target area so that the water and fat are excited to a saturated state.

[0035] The specific frequency range is the preset frequency range, which is a frequency range covering the resonance frequency of water molecules and the resonance frequency of fat. Within the preset frequency range, the magnetic resonance device can be divided into multiple preset frequencies, and at each preset frequency, a radio frequency pulse can be emitted to excite the water and fat in the target area to a saturated state, so as to collect the relevant signals of water and fat at the preset frequency. The number of times the radio frequency pulse is emitted is proportional to the number of preset frequencies, that is, the more the number of preset frequencies, the more the number of times the radio frequency pulse is emitted, and the richer and more accurate the relevant signals obtained are.

[0036] In the process of drawing the spectrum diagram, the chemical shift is used as the horizontal coordinate, so the above-mentioned preset frequency range and preset frequency can be converted. In the embodiment of the present application, the resonance frequency of water is used as the reference frequency, that is, the chemical shift refers to the chemical shift of water and / or fat relative to the resonance frequency of water molecules. For the convenience of subsequent description, the chemical shift of the resonance frequency of water relative to water molecules can be recorded as the first chemical shift of water, and the chemical shift of fat relative to the resonance frequency of water molecules can be recorded as the second chemical shift. Among them, the preset frequency range is converted into chemical shift and can be recorded as a preset chemical shift range. Correspondingly, the preset chemical shift range covers the chemical shift of the second chemical shift and the first chemical shift, for example, between -5ppm and 5ppm. Within the preset chemical shift range, as water and fat gradually approach saturation, the intensity of the mixed signal of water and fat will gradually decrease, that is, the degree of reduction of the mixed signal reflects the change in the intensity of the mixed signal within the preset chemical shift range.

[0037] In some embodiments, in order to stimulate the water and fat in the target area to a saturated state, a radio frequency pulse with an energy of 0.5 μT to 1.5 μT, preferably 0.9 μT to 1.1 μT, and more preferably 1 μT can be selected; the emission duration of the radio frequency pulse can be between 50 ms and 500 ms, preferably between 100 ms and 200 ms, and more preferably 100 ms. The energy value and the emission duration value of the radio frequency pulse can be freely combined, for example, the energy of the radio frequency pulse is 1 μT and the emission duration is 100 ms, or the energy of the radio frequency pulse is 0.9 μT and the emission duration is 100 ms.

[0038] During the excitation process, although there are other endogenous substances such as amino groups, magnetization transfer, nuclear Ostwald effect, etc. in the target area, even if these endogenous substances have chemical exchange saturation transfer effect, only by continuously exciting the target area with a radio frequency pulse with an emission duration of more than 500ms will these endogenous substances produce obvious interference signals. In other words, by exciting the target area with the radio frequency pulse in this embodiment, the interference of the above endogenous substances on the mixed signal can be reduced.

[0039] Step 120: After the transmission of the radio frequency pulse is completed, a saturation signal corresponding to the preset frequency is collected.

[0040] The emission of radio frequency pulses will change the strength of the mixed signal of water and fat. In order to detect the changed mixed signal, the corresponding saturation signal can be collected at each preset frequency. After collecting the saturation signals corresponding to all preset frequencies, the saturation signals equal to the number of preset frequencies can be obtained to generate subsequent mixed signals of water and fat.

[0041] In some embodiments, the signal acquisition method is a method for acquiring spin echo magnetic resonance data. The magnetic resonance data acquisition method is adopted because the signal data acquired by the method can be processed to obtain millimeter-level imaging. That is, the present application can process the acquired mixed signal to obtain imaging with a higher spatial resolution than magnetic resonance spectroscopy imaging (with a spatial resolution of centimeters).

[0042] Step 130: compare the saturated signals corresponding to all preset frequencies with the pre-collected initial signals to determine the mixed signals of water and fat within the preset frequency range;

[0043] In order to detect the change of the mixed signal of water and fat, the initial signal of water and fat can be collected first, and the saturated signal corresponding to each preset frequency can be compared with the initial signal to determine the mixed signal of water and fat within the preset frequency range. For example, the signal obtained by water and fat in a saturated state can be subtracted from the initial signal of water and fat collected when no radio frequency pulse is emitted, so as to obtain the change of the signal at each preset frequency. Based on this, the overall mixed signal of water and fat can be obtained. It can be considered that the final mixed signal is a signal collection, which contains the mixed signal corresponding to each preset frequency. It can be understood that the more preset frequencies there are in the preset frequency range, the higher the acquisition frequency of the saturated signal, and correspondingly, the more accurate the obtained mixed signal is.

[0044] Step 140: Calculate a first signal of water.

[0045] In order to separate the fat signal from the mixed signal, the first signal of water can be calculated first. For example, the first signal can be fitted by a single-cell Lorentz model, or the first signal can be calculated by a water molecule signal saturation model based on relaxation information and property information of the radio frequency pulse. The specific calculation method of the first signal is not limited in this embodiment.

[0046] In some embodiments, in order to improve the accuracy of the first signal calculation, the above step 140 specifically includes:

[0047] The first signal of water is fitted by the single-cell Lorentz model, where the fitting formula of the first signal is:

[0048]

[0049] Z water (Δω) is the first signal, Δω is the chemical shift of water relative to the resonance frequency of water molecules, A is the amplitude of the first signal, and lw is the bandwidth of the first signal.

[0050] In the spectrum distribution of the fat signal, the fat peak is mainly concentrated in the chemical shift range of -3.8ppm to -2ppm, and the chemical shift range of -2ppm and above is the dominant range of the spectrum distribution of the first signal, that is, the water peak is mainly concentrated in the chemical shift range of -2ppm and above. In order to calculate the first signal, the specific resonance frequency range can be determined according to the first chemical shift, so as to accurately calculate the first signal.

[0051] In the process of determining the chemical shift calculation range, it should be ensured that the chemical shift calculation range includes the first chemical shift, and the determined chemical shift calculation range is the dominant range of the spectral distribution of the first signal. For example, it can be a chemical shift of -2ppm and above, expressed as Δω>-2ppm. The chemical shift calculation range here is determined from the preset chemical shift range in step 110. After the chemical shift calculation range is determined, the first signal can be further determined according to the above fitting formula. It can be understood that the chemical shift calculation range can also be divided into multiple chemical shift points, and for each chemical shift point, a corresponding first signal can be calculated. Similar to the mixed signal, the final first signal is also regarded as a signal collection, which includes the first signal corresponding to each chemical shift point.

[0052] In some embodiments, in addition to obtaining the first signal by fitting the single-cell Lorentz model, the first signal can also be obtained by the following method:

[0053] Step 141: Collect longitudinal relaxation data and transverse relaxation data of water.

[0054] Step 142: Calculate the longitudinal relaxation rate R of water based on the longitudinal relaxation data. 1w .

[0055] Step 143: Calculate the transverse relaxation rate R of water based on the transverse relaxation data. 2w .

[0056] Step 144: Calculate a first signal based on the longitudinal relaxation rate, the transverse relaxation rate, and the attribute information of the radio frequency pulse.

[0057] By additionally acquiring two types of relaxation data of water, calculating the relaxation rates of the two types of relaxation, and combining the two relaxation rates with the property information of the radio frequency pulse, the first signal is estimated, and the formula is:

[0058]

[0059] Among them, R 1ρ =R 1w cos 2 θ+R 2w sin 2 θ, θ = tan-1 (ω1 / Δω), ω1 is the saturation pulse intensity, TR is the scan repetition time, T s is the RF pulse saturation time.

[0060] Among them, the two types of relaxation data are longitudinal relaxation data and transverse relaxation data, specifically, the longitudinal relaxation data is T1 mapping, and the transverse relaxation data is T2 mapping. Simply put, T1 mapping is the data obtained by weighting the signal with the T1 relaxation time, and T2 mapping is the data obtained by weighting the signal with the T2 relaxation time. After obtaining T1mapping and T2 mapping, the longitudinal relaxation rate and transverse relaxation rate of water can be calculated based on the obtained relevant data; the two relaxation rates are combined with the attribute information of the radio frequency pulse to further estimate the first signal. Among them, the attribute information of the radio frequency pulse includes the energy of the radio frequency pulse and the emission duration of the radio frequency pulse.

[0061] Step 150: determine a second signal based on the mixed signal and the first signal, and generate a fat peak map of the target area based on the frequency spectrum distribution of the second signal.

[0062] After obtaining the mixed signal and the first signal, the spectrum distribution of the fat signal can be determined based on the spectrum distributions corresponding to the two signals. That is, the spectrum distribution of the second signal is determined so that the magnetic resonance device can image the target area based on the spectrum distribution of the second signal and obtain a fat peak spectrum of the target area. The fat peak spectrum can characterize the distribution of fat and provide doctors with relevant information about fat metabolism in the target area, which helps doctors further diagnose the patient's condition.

[0063] In some embodiments, the above step 150 specifically includes: subtracting the spectral distribution of the mixed signal from the spectral distribution of the first signal to obtain the spectral distribution of the second signal, and generating a fat peak map of the target area based on the spectral distribution of the second signal.

[0064] The process of determining the spectrum distribution of the second signal can be: firstly, the mixed signal and the first signal are subtracted to obtain the second signal; if the first signal is denoted as Z water (Δω), the mixed signal is denoted as ΔZ(Δω), and the second signal is denoted as Z fat (Δω), then Z fat (Δω)=ΔZ(Δω)-Z water (Δω), then we get Z fat (Δω), it is equivalent to obtaining the spectrum distribution of the second signal.

[0065] The embodiment of the present application uses radio frequency pulses to saturate the water and fat in the target area, thereby obtaining a mixed signal of water and fat; then the first signal of water is calculated based on the resonance frequency of water molecules; after obtaining the mixed signal and the first signal, the spectrum distribution of the first signal is separated from the spectrum distribution of the mixed signal, thereby obtaining the spectrum distribution of the fat signal, that is, the spectrum distribution of the second signal; finally, the spectrum distribution of the second signal can be used to image the target area. The mixed signal in this method can be collected based on a conventional magnetic resonance signal acquisition method, and by processing the mixed signal, a millimeter-level image can be generated, thereby improving the spatial resolution of the spectral image.

[0066] In order to further prove the feasibility of this application and the beneficial effects that this application can bring, the following is an explanation with specific experimental data. Among them, the target area is the area where the cone is located in the cross-section of the abdomen. The above-mentioned first signal represents the feedback signal of water in a saturated state, which can be called a water signal; the second signal represents the feedback signal of fat in a water-saturated state, which can be called a fat signal.

[0067] like Figure 2 As shown, Figure 2 The schematic diagram of obtaining the spectrum distribution of the fat signal is shown. Among them, the spectrum distribution curve of the mixed signal is subtracted from the spectrum distribution curve of the water signal (obtained by fitting the single-cell Lorentz function), and the spectrum distribution curve of the fat signal can be obtained. From the spectrum distribution curve of the fat signal, it can be seen that obvious peaks can be seen at chemical shifts of 0.6ppm, -0.5ppm, -2ppm, -2.6ppm, -3.4ppm and -3.8ppm, and after calculating the resonance frequencies corresponding to these chemical shifts, it is found that the calculated resonance frequencies are consistent with the main resonance frequencies of fat, that is, the results verify that the method for determining the spectrum distribution of fat signals proposed in this application is reliable and accurate.

[0068] Figure 3 The fat peak spectrum with a chemical shift of -2.6 ppm is shown. Figure 4 The fat peak spectrum with a chemical shift of -3.4ppm is shown, which is comparable to ordinary magnetic resonance images, and the spatial resolution can reach the millimeter level, which is better than magnetic resonance spectroscopy imaging. In other words, the fat peak spectrum further verifies that the imaging method proposed in this application can improve the spatial resolution of imaging.

[0069] The imaging method provided in the above embodiments of the present application can bring the following beneficial effects:

[0070] (1) A mixed signal is obtained by using a conventional magnetic resonance signal acquisition method, and after processing the mixed signal, millimeter-level imaging can be achieved;

[0071] (2) The complete spectrum distribution of fat signal can be obtained.

[0072] Corresponding to the imaging method of the above embodiment, Figure 5 A structural block diagram of an imaging device 5 provided in an embodiment of the present application is shown. For ease of explanation, only the parts related to the embodiment of the present application are shown.

[0073] Reference Figure 5 , the imaging device 5 comprises:

[0074] The excitation module 51 is used for transmitting a radio frequency pulse of a preset frequency to a target area for each preset frequency within a preset frequency range, wherein the radio frequency pulse is used to excite water and fat in the target area so that the water and fat in the target area are excited to a saturated state;

[0075] The acquisition module 52 acquires a saturation signal corresponding to a preset frequency after the transmission of the radio frequency pulse ends;

[0076] A determination module 53, used to compare the saturated signals corresponding to all preset frequencies with the pre-collected initial signals to determine a mixed signal of water and fat within a preset frequency range;

[0077] A calculation module 54, used for calculating a first signal of water;

[0078] The imaging module 55 is used to determine the second signal based on the mixed signal and the first signal, and generate a fat peak map of the target area based on the frequency spectrum distribution of the second signal, where the second signal is a fat signal.

[0079] Optionally, the calculation module 54 may include:

[0080] The first computing unit is used to

[0081] The first signal of water is fitted by the single-cell Lorentz model, where the fitting formula of the first signal is:

[0082]

[0083] Z water (Δω) is the first signal, Δω is the chemical shift of water relative to the resonance frequency of water molecules, A is the amplitude of the first signal, and lw is the bandwidth of the first signal.

[0084] Optionally, the calculation module 54 may include:

[0085] A collection unit, used for collecting longitudinal relaxation data and transverse relaxation data of water;

[0086] The second calculation unit is used to calculate the longitudinal relaxation rate R of water based on the longitudinal relaxation data. 1w ;

[0087] The third calculation unit is used to calculate the transverse relaxation rate R of water based on the transverse relaxation data. 2w ;

[0088] A fourth computing unit, for

[0089] Based on the longitudinal relaxation rate R 1w , transverse relaxation rate R 2w The first signal is calculated with the attribute information of the RF pulse, and the formula is:

[0090]

[0091] Among them, R 1ρ =R 1w cos 2 θ+R 2w sin 2 θ, θ = tan -1 (ω1 / Δω), Z water (Δω) is the first signal, ω1 is the saturation pulse intensity, TR is the scanning repetition time, T s is the RF pulse saturation time.

[0092] Optionally, the excitation module 51 is specifically used for: transmitting radio frequency pulses to the target area so that the water and fat in the target area are excited to a saturated state, wherein the energy of the radio frequency pulses is 0.5 μT to 1.5 μT, and the emission duration is 50 ms to 500 ms.

[0093] Optionally, the excitation module 51 is specifically used for: transmitting radio frequency pulses to the target area so that the water and fat in the target area are excited to a saturated state, wherein the energy of the radio frequency pulses is 0.9 μT to 1.1 μT, and the emission duration is 100 ms to 200 ms.

[0094] Optionally, the excitation module 51 is specifically used for: transmitting radio frequency pulses to the target area so that the water and fat in the target area are excited to a saturated state, wherein the energy of the radio frequency pulses is 1 μT and the emission duration is 100 ms.

[0095] Optionally, the imaging module 55 is specifically used to: perform a subtraction between the mixed signal and the first signal to obtain a second signal.

[0096] It should be noted that the information interaction and execution process between the above-mentioned devices / units are based on the same concept as the method embodiment of the present application. Their specific functions and technical effects can be found in the method embodiment part and will not be repeated here.

[0097] Figure 6This is a schematic diagram of the physical structure of an electronic device provided by an embodiment of the present application. Figure 6 As shown, the electronic device 6 of this embodiment includes: at least one processor 60 ( Figure 6 Only one is shown in the figure) a processor, a memory 61, and a computer program 62 stored in the memory 61 and executable on at least one processor 60, when the processor 60 executes the computer program 62, the steps in any of the above-mentioned imaging method embodiments are implemented, for example Figure 1 Steps 110-150 are shown.

[0098] The processor 60 may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc.

[0099] In some embodiments, the memory 61 may be an internal storage unit of the electronic device 6, such as a hard disk or memory of the electronic device 6. In other embodiments, the memory 61 may also be an external storage device of the electronic device 6, such as a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, etc., equipped on the electronic device 6.

[0100] Furthermore, the memory 61 may also include both an internal storage unit of the terminal device 6 and an external storage device. The memory 61 is used to store operating means, application programs, boot loaders, data and other programs, such as program codes of computer programs, etc. The memory 61 may also be used to temporarily store data that has been output or is to be output.

[0101] The technicians in the relevant field can clearly understand that for the convenience and simplicity of description, only the division of the above-mentioned functional units and modules is used as an example for illustration. In practical applications, the above-mentioned function allocation can be completed by different functional units and modules as needed, that is, the internal structure of the above-mentioned device can be divided into different functional units or modules to complete all or part of the functions described above. The functional units and modules in the embodiment can be integrated in a processing unit, or each unit can exist physically separately, or two or more units can be integrated in one unit. The above-mentioned integrated unit can be implemented in the form of hardware or in the form of software functional units. In addition, the specific names of the functional units and modules are only for the convenience of distinguishing each other, and are not used to limit the scope of protection of this application. The specific working process of the units and modules in the above-mentioned system can refer to the corresponding process in the aforementioned method embodiment, which will not be repeated here.

[0102] An embodiment of the present application further provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the steps in the above-mentioned method embodiments can be implemented.

[0103] An embodiment of the present application provides a computer program product. When the computer program product runs on a mobile terminal, the mobile terminal can implement the steps in the above-mentioned method embodiments when executing the computer program product.

[0104] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the present application implements all or part of the processes in the above-mentioned embodiment method, which can be completed by instructing the relevant hardware through a computer program. The above-mentioned computer program can be stored in a computer-readable storage medium, and the computer program can implement the steps of the above-mentioned various method embodiments when executed by the processor. Among them, the above-mentioned computer program includes computer program code, and the above-mentioned computer program code can be in source code form, object code form, executable file or some intermediate form. The above-mentioned computer-readable medium may at least include: any entity or device that can carry the computer program code to the camera / electronic device, recording medium, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electrical carrier signal, telecommunication signal and software distribution medium. For example, a USB flash drive, a mobile hard disk, a disk or an optical disk.

[0105] In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described or recorded in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0106] Those of ordinary skill in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of this application.

[0107] In the embodiments provided in the present application, it should be understood that the disclosed devices / network equipment and methods can be implemented in other ways. For example, the device / network equipment embodiments described above are merely schematic. For example, the division of the above modules or units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms.

[0108] The units described above as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0109] The above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. These modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present application, and should all be included in the protection scope of the present application.

Claims

1. An imaging method, characterized in that: include: For each preset frequency within a preset frequency range: transmitting a radio frequency pulse of the preset frequency to a target area, wherein the radio frequency pulse is used to stimulate water and fat in the target area so that the water and fat in the target area are stimulated to a saturated state; After the transmission of the radio frequency pulse is completed, a saturation signal corresponding to the preset frequency is collected; Comparing the saturated signals corresponding to all the preset frequencies with the pre-collected initial signals to determine the mixed signal of the water and the fat within the preset frequency range; calculating a first signal of the water; A second signal is determined based on the mixed signal and the first signal, and a fat peak map of the target area is generated based on the frequency spectrum distribution of the second signal, where the second signal is a signal of the fat.

2. The imaging method according to claim 1, characterized in that The calculating the first water signal comprises: The first signal of the water is fitted by a single-cell Lorentz model, wherein the fitting formula of the first signal is: The Z water (Δω) is the first signal, Δω is the chemical shift of water relative to the resonance frequency of water molecules, A is the amplitude of the first signal, and lw is the bandwidth of the first signal.

3. The imaging method according to claim 1, characterized in that: The calculating the first water signal comprises: Collecting longitudinal relaxation data and transverse relaxation data of the water; The longitudinal relaxation rate R of water is calculated based on the longitudinal relaxation data. 1w ; The transverse relaxation rate R of water is calculated based on the transverse relaxation data. 2w ; Based on the longitudinal relaxation rate R 1w , the transverse relaxation rate R 2w The first signal is calculated based on the attribute information of the radio frequency pulse, and the formula is: Among them, R 1ρ =R 1w cos 2 θ+R 2w sin 2 θ, θ = tan -1 (ω1 / Δω), the Z water (Δω) is the first signal, Δω is the chemical shift of water relative to the resonance frequency of water molecules, ω1 is the saturation pulse intensity, TR is the scanning repetition time, and T s is the RF pulse saturation time.

4. The imaging method according to any one of claims 1 to 3, characterized in that: The energy of the radio frequency pulse is 0.5 μT to 1.5 μT, and the emission duration is 50 ms to 500 ms.

5. The imaging method according to any one of claims 1 to 3, characterized in that: The energy of the radio frequency pulse is 0.9 μT to 1.1 μT, and the emission duration is 100 ms to 200 ms.

6. The imaging method according to any one of claims 1 to 3, characterized in that: The energy of the radio frequency pulse is 1 μT, and the emission duration is 100 ms.

7. The imaging method according to claim 1, characterized in that: The determining of the second signal based on the mixed signal and the first signal comprises: The mixed signal is subtracted from the first signal to obtain the second signal.

8. An imaging device, characterized in that: include: an excitation module, for each preset frequency within a preset frequency range: emitting a radio frequency pulse of the preset frequency to a target area, wherein the radio frequency pulse is used to excite water and fat in the target area so that the water and fat in the target area are excited to a saturated state; An acquisition module, used for acquiring a saturation signal corresponding to the preset frequency after the transmission of the radio frequency pulse is completed; A determination module, used for comparing the saturated signals corresponding to all the preset frequencies with the pre-collected initial signals to determine the mixed signal of the water and the fat within the preset frequency range; A calculation module, used for calculating the first signal of water based on the spectrum distribution characteristics of water molecules; An imaging module is used to determine a second signal based on the mixed signal and the first signal, and to generate a fat peak map of the target area based on the frequency spectrum distribution of the second signal, wherein the second signal is a signal of the fat.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the imaging method according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the imaging method according to any one of claims 1 to 7 is implemented.

Images