WO2018227469A1 - 磁共振波谱交互方法、系统及计算机可读的存储媒介 - Google Patents
磁共振波谱交互方法、系统及计算机可读的存储媒介 Download PDFInfo
- Publication number
- WO2018227469A1 WO2018227469A1 PCT/CN2017/088425 CN2017088425W WO2018227469A1 WO 2018227469 A1 WO2018227469 A1 WO 2018227469A1 CN 2017088425 W CN2017088425 W CN 2017088425W WO 2018227469 A1 WO2018227469 A1 WO 2018227469A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- concentration
- complex
- metabolites
- operators
- composite
- Prior art date
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/742—Details of notification to user or communication with user or patient ; user input means using visual displays
- A61B5/7435—Displaying user selection data, e.g. icons in a graphical user interface
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/46—NMR spectroscopy
- G01R33/4625—Processing of acquired signals, e.g. elimination of phase errors, baseline fitting, chemometric analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
- G01R33/485—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy based on chemical shift information [CSI] or spectroscopic imaging, e.g. to acquire the spatial distributions of metabolites
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/742—Details of notification to user or communication with user or patient ; user input means using visual displays
- A61B5/7425—Displaying combinations of multiple images regardless of image source, e.g. displaying a reference anatomical image with a live image
Definitions
- the present application relates to a method and system for magnetic resonance spectroscopy processing, and in particular to a method and system for magnetic resonance spectroscopy interaction and analysis.
- Magnetic resonance spectroscopy is a method for performing quantitative analysis of specific nuclei and compounds using magnetic resonance phenomena and chemical shift effects. As a non-invasive technique for measuring the metabolism of a compound developed in recent years, magnetic resonance spectroscopy detects abnormality of a certain tissue by measuring the concentration of a certain compound or complex.
- the interaction mode of the calculation complex has the disadvantages that the interface is not friendly enough, the user operation is limited, the editing of the composite is inflexible, etc., which affects the user's diagnosis experience and reduces the user's diagnosis efficiency. Based on this, it is necessary to design an interactive method and system that is interactive, clear, and can preview the edited composite, and can flexibly edit the composite to meet the user's diagnostic needs.
- a magnetic resonance spectroscopy system comprising a computer readable storage medium configured to store executable modules, including an editing module, configured to select two or more based on user input a metabolite and one or more operators to generate a complex; a calculation module configured to calculate a concentration of the complex; and a display module configured to display a concentration calculation result of the complex to a user; and a process
- the processor is capable of executing the executable module of the computer readable storage medium storage.
- Some embodiments of the present application further disclose a magnetic resonance spectroscopy system, wherein the editing module includes a metabolite list unit configured to provide a plurality of selectable metabolites.
- Some embodiments of the present application also disclose a magnetic resonance spectroscopy system, wherein the editing module further includes an operator unit configured to provide a plurality of selectable operators.
- Some embodiments of the present application also disclose a magnetic resonance spectroscopy system, wherein the editing module further includes a preview unit configured to display the one or more compounds and the one or more operations symbol.
- Some embodiments of the present application also disclose a magnetic resonance spectroscopy system, wherein the computing module is configured to accept the composite from an editing module and calculate a concentration of the composite.
- Some embodiments of the present application also disclose a magnetic resonance spectroscopy system, wherein the calculation module parses the complex according to an inverse Polish algorithm to obtain the two or more metabolites and the one or more operations And calculating the concentration of the complex based on a single concentration of the two or more metabolites and the one or more operators.
- Some embodiments of the present application also disclose a magnetic resonance spectroscopy system, wherein the display module is configured to display a spatial distribution of the concentration of the complex within the subject of the test subject.
- Some embodiments of the present application disclose a magnetic resonance spectroscopy method comprising the steps of selecting two or more metabolites and one or more operators to generate a complex based on user input; calculating the concentration of the complex And displaying the concentration calculation result of the complex to the user.
- Some embodiments of the present application disclose a magnetic resonance spectroscopy method, further comprising displaying the two or more metabolites and the one or more operators.
- Some embodiments of the present application disclose a magnetic resonance spectroscopy method, wherein calculating the concentration of the complex further comprises resolving the complex according to an inverse Polish algorithm to obtain the two or more metabolites and the one Or a plurality of operators and calculating the concentration of the complex based on a single concentration of the two or more metabolites and the one or more operators.
- Some embodiments of the present application disclose a magnetic resonance spectroscopy method, wherein displaying a concentration calculation result of the complex further comprises displaying a spatial distribution of the concentration of the complex in a subject to be detected.
- Some embodiments of the present application disclose a computer readable storage medium configured to store an executable program to perform operations of selecting two or more metabolites and one or more operator generation based on user input.
- a complex calculating a concentration of the complex; displaying a concentration calculation result of the complex to a user; and a processor capable of executing the computer readable storage medium Executable executable program.
- Some embodiments of the present application disclose a computer readable storage medium, wherein the executable program is further operative to perform operations of displaying the two or more metabolites and the one or more operators.
- Some embodiments of the present application disclose a computer readable storage medium, wherein the executable program is further operative to: calculate a concentration of the composite comprising analyzing the composite according to an inverse Polish algorithm to obtain the Two or more metabolites and the one or more operators, and calculating the concentration of the complex based on a single concentration of the two or more metabolites and the one or more operators.
- Some embodiments of the present application disclose a computer readable storage medium, wherein the executable program is further executable to: display a concentration calculation result of the complex further comprising displaying a concentration of the complex at a detection object Spatial distribution in the body.
- FIG. 1A is a schematic diagram of a magnetic resonance spectroscopy processing system shown in accordance with some embodiments of the present application;
- FIG. 1B is an architectural diagram of a computer device of a compound editing interactive system shown in accordance with some embodiments of the present application;
- 2A is a diagram showing an example of a time domain spectrum signal collected by a spectrum signal acquisition device
- 2B is a diagram showing an example of a frequency domain spectrum signal processed by Fourier transform
- FIG. 3 is an exemplary schematic diagram of a compound editing interaction system shown in accordance with some embodiments of the present application.
- FIG. 4 is an exemplary schematic diagram of an editing module shown in accordance with some embodiments of the present application.
- FIG. 5 is a flow chart showing a concentration of a generated complex to a user, shown in some embodiments of the present application.
- FIG. 6 is a flow diagram of a method of magnetic resonance spectroscopy analysis and interaction of an editable composite, in accordance with some embodiments of the present application.
- FIG. 7 is a spatial distribution image of a composite concentration obtained according to Flowchart 6 of an exemplary composite, in accordance with some embodiments of the present application.
- the "scan area" represents the actual area where the scan is performed.
- the region of interest of the detected object can be scanned.
- the region of interest may be an entirety or a part of the detection object, for example, a head, a chest, an abdomen, a heart, a liver, an upper limb, a lower limb, a spine, a bone, a blood vessel, a lesion, a tumor site, etc., or the above random combination.
- a scanned image of a plurality of time points is acquired for a selected layer of the region of interest.
- the "element” represents the smallest component in the image matrix
- the "voxel” represents the smallest component in the actual region. Unless the context clearly indicates an exception, the "element” in the image matrix and the “voxel” in the actual region corresponding to the image matrix in this application may mean the same meaning and may be replaced.
- an aspect of the present application relates to a magnetic resonance spectroscopy analysis and exchange system.
- the magnetic resonance analysis and exchange system can include an editing module, a computing module, and a display module.
- Another aspect of the present application relates to a method of displaying a change in concentration distribution or concentration distribution of an edited complex.
- the magnetic resonance spectroscopy interaction method described in the present application includes performing concentration calculation and visual display on a complex that needs to be edited based on user input.
- Embodiments of the present application can be applied to different image processing systems.
- Different image processing systems may include a Magnetic Resonance Imaging (MRI) system, a nuclear magnetic resonance-positron emission tomography hybrid system (MR-PET system), and the like.
- MRI Magnetic Resonance Imaging
- MR-PET nuclear magnetic resonance-positron emission tomography hybrid system
- Magnetic Resonance Spectroscopy as a method for non-invasive detection of biochemical characteristics of living tissue, can be used to study the pathological and physiological changes of human cell metabolism. In many diseases, metabolic changes precede pathological changes, and MRS is potentially sensitive to such metabolic changes, providing information to detect lesions early.
- Magnetic resonance spectroscopy processing system 100 may include one or more spectral signal acquisition devices 110, one or more networks 120, one or more compound editing interaction systems 130, one or more databases 140.
- the spectral signal acquisition device 110 can scan the detection object to obtain scan data, which can be sent to the compound editing interaction system 130 through the network 120 for further processing, or can be directly stored in the database 140 through the network 120. In some embodiments, the spectral signal acquisition device 110 can be directly coupled to the compound editing interaction system 130. The scan data can be sent directly to the compound editing interaction system 130.
- the detection object may include a human body, an animal, or the like.
- the spectral signal acquisition device 110 may include, but is not limited to, a Magnetic Resonance Imaging (MRI) device or a nuclear magnetic resonance-positron emission computed tomography hybrid system (MR-PET system) or the like.
- MRI Magnetic Resonance Imaging
- MR-PET system nuclear magnetic resonance-positron emission computed tomography hybrid system
- the compound editing interaction system 130 can input the data (eg, spectral signals obtained by the spectral signal acquisition device 110 and/or stored in the database 140, metabolite concentration data, user selected Operators, etc.) are processed and analyzed to generate processing results.
- the compound editing interaction system 130 can generate a complex based on the selected plurality of metabolites.
- the compound editing interaction system 130 can visualize the obtained composite concentration data to generate a spatial distribution map of the complex concentration.
- the scanned image may be a two-dimensional image or a three-dimensional image.
- the compound editing interaction system 130 can include a processor and input/output devices (not shown).
- the processor can be a server or a server group.
- a server group can be centralized, such as a data center.
- a server group can also be distributed, such as a distributed system.
- the processor may be one or a combination of a cloud server, a file server, a database server, a File Transfer Protocol (FTP) server, an application server, a proxy server, a mail server, and the like.
- the processor can be local or remote.
- the local processor can include a processor integrated in the compound editing interaction system 130.
- the remote processor can include a processor coupled to the compound editing interaction system 130 over a network (e.g., network 120).
- the processor can access information stored in database 140 (eg, a medical image stored in database 140), information in spectral signal acquisition device 110 (eg, medical imaging taken by spectral signal acquisition device 110) image).
- the input/output device may input data to the processor, or may receive data output by the processor, and output the data in numbers, characters, images, videos, animations, sounds, and the like. The form is expressed.
- the input/output devices may include, but are not limited to, one or a combination of input devices, output devices, and the like.
- the input device may include, but is not limited to, a character input device (eg, a keyboard), an optical reading device (eg, an optical tag reader, an optical character reader), a graphic input device (eg, a mouse, a joystick, a stylus), an image A combination of one or more of an input device (eg, a video camera, a scanner, a fax machine), an analog input device (eg, a language analog to digital conversion recognition system), and the like.
- the output device may include, but is not limited to, a combination of one or more of a display device, a printing device, a plotter, an image output device, a voice output device, a magnetic recording device, and the like.
- the compound editing interaction system 130 can further include a storage device (not shown).
- the storage device can store various information such as programs, data, and the like.
- data and/or processing results eg, anatomical images, frequency domain spectral signals, density spatial distribution images, etc.
- the storage device can also be output through an input/output device.
- Database 140 can be broadly referred to as a device having storage capabilities.
- the database 140 can store the scan data collected from the spectral signal acquisition device 110 and various data generated in the operation of the compound editing interaction system 130.
- Database 140 can be local or remote.
- the local database may include a device having a storage function integrated in the database 140.
- the remote database may include a storage-enabled device connected to the database 140 over a network (eg, network 120).
- Database 140 may include, but is not limited to, a combination of one or more of a hierarchical database, a networked database, and a relational database.
- the database 140 can digitize the information and store it in a storage device that utilizes electrical, magnetic or optical means.
- the database 140 can be used to store various information such as programs and data.
- the database 140 may be a device that stores information by means of electrical energy, such as various memories, random access memory (RAM), read only memory (ROM), and the like.
- the random access memory includes but is not limited to a decimal counter tube, a selection tube, a delay line memory, a Williams tube, a dynamic random access memory (DRAM), a static random access memory (SRAM), a thyristor random access memory (T-RAM), and a zero capacitor.
- DRAM dynamic random access memory
- SRAM static random access memory
- T-RAM thyristor random access memory
- Z-RAM random access memory
- Read-only memory includes, but is not limited to, bubble memory, magnetic button line memory, thin film memory, magnetic plate line memory, magnetic core memory, drum memory, optical disk drive, hard disk, magnetic tape, early non-volatile memory (NVRAM), phase change Memory, magnetoresistive random storage memory, ferroelectric random access memory, nonvolatile SRAM, flash memory, electronic erasable rewritable read only memory, erasable programmable read only memory, programmable read only memory, shielded A combination of one or more of a heap read memory, a floating connection gate random access memory, a nano random access memory, a track memory, a variable resistive memory, a programmable metallization cell, and the like.
- the database 140 may be a device that stores information using magnetic energy, such as a hard disk, a floppy disk, a magnetic tape, a magnetic core memory, a magnetic bubble memory, a USB flash drive, a flash memory, or the like.
- the database 140 may be a device that optically stores information, such as a CD or a DVD.
- the database 140 may be a device that stores information using magneto-optical means, such as a magneto-optical disk or the like.
- the access mode of the database 140 may be one or a combination of random storage, serial access storage, read-only storage, and the like.
- the database 140 can be a non-permanent memory or a permanent memory.
- the above-mentioned storage devices are merely examples, and the database that can be used in the magnetic resonance spectrum analysis 100 is not limited thereto.
- Network 120 can be a single network or a combination of multiple networks.
- Network 120 may include, but is not limited to, a local area network, a wide area network, a public network, a private network, a wireless local area network, a virtual network, an urban A combination of one or more of a metropolitan area network, a public switched telephone network, and the like.
- Network 120 may include a variety of network access points, such as wired or wireless access points, base stations, or network switching points, through which the data sources are connected to network 120 and transmitted over the network.
- database 140 may be a cloud computing platform with data storage capabilities, including but not limited to public clouds, private clouds, community clouds, hybrid clouds, and the like. Variations such as these are within the scope of the present application.
- FIG. 1B is an architectural diagram of a computer device of the compound editing interaction system 130, in accordance with some embodiments of the present application.
- the computer can be a general purpose computer or a computer with a specific purpose.
- Compound editing interaction system 130 can be implemented by the computer device architecture through its hardware devices, software programs, firmware, and combinations thereof. For ease of illustration, only one computer device architecture is depicted in FIG. 1B, but the related computer functions of the compound editing interaction system 110 may be distributed by multiple computer devices.
- the computer device architecture can include a communication port 230 to which a network (eg, network 120 in FIG. 1A) can be implemented for data communication.
- the computer device architecture may also include a central processing system (CPU) unit 240 for executing program instructions comprised of one or more processors.
- the computer device architecture includes an internal communication bus 270, different forms of program storage units and data storage units, such as a hard disk 210, a read only memory (ROM) 250, a random access memory (RAM) 260, which can be configured as a storage device.
- the computer device architecture can also include an input/output component 220 that supports data and/or information interaction between the computer device architecture and an external (eg, user).
- the computer device architecture can also accept programs and data over a communication network.
- Magnetic resonance spectroscopy can provide metabolic information about the tissue.
- the tissue referred to herein may be a human or animal organ, body fluid, nerve, cell, or the like, or a random combination of the above.
- the working principle of magnetic resonance spectroscopy is to use the difference in the precession frequency of the same magnetic nucleus caused by the different structures of different molecules. Specifically, since the magnetic nucleus of the magnetic nucleus is fixed, the magnetic nucleus is affected by the external static magnetic field in addition to the external static magnetic field, and is also affected by the electron cloud around the nucleus and other surrounding electron clouds. The perceived magnetic field strength is slightly lower than the applied static magnetic field, so the precession frequency is also slightly reduced.
- the magnetic shielding effect of the electron cloud on the magnetic nucleus also differs, showing the difference in the precession frequency.
- the phenomenon that the difference in the precession frequency of the same magnetic nucleus due to the different molecular structure is called the chemical shift phenomenon.
- a designed RF pulse is applied to the target region, and the frequency range is required to cover the precession frequency of the proton to be detected in the metabolite.
- the MR signal (which may be an FID signal or an echo signal) of the region is collected, and the MR signal is derived from protons in various metabolites.
- the proton precession frequency in different metabolites is slightly different. Fourier transform can obtain information on the mass spectrum of different substances.
- the line includes a series of relatively narrow peaks.
- the abscissa indicates the precession frequency of protons in different materials, usually expressed in parts per million (PPM) (based on the proton precession frequency of the standard, the proton precession frequency in other metabolites and the proton precession frequency in the standard)
- PPM parts per million
- the area under the peak of a narrow wave is proportional to the content or concentration of a particular metabolite in the target area. From this, the concentration of the metabolite can be calculated.
- 1 H, 31 P, 12 C, 23 Na, and 19 F can be collected for magnetic resonance spectroscopy for brain tumors, cerebral ischemic diseases, prostate cancer, and diffuse liver disease. , renal function analysis, renal transplant rejection, brain changes in metabolic diseases, diagnosis and differential diagnosis of recurrence after treatment of brain tumors.
- a relatively stable chemical may be selected as a reference for the precession frequency of a magnetic nucleus-related metabolite.
- 1 H MRS is often selected from trimethylsilane as a standard, and 31 P MRS may be used. Creatine phosphate was used as a standard.
- the concentration of the metabolite is low and the signal generated is almost one ten-thousandth of the normal water signal, the sensitivity of the MRS is relatively low, and it is necessary to combine water inhibition and fat suppression to obtain a signal.
- the metabolites of brain 1 H magnetic resonance spectroscopy can include (1) NAA (N-acetylaspartate), which is mainly found in neurons and their axons, and can be used as an internal standard for neurons, and its content can reflect nerves. The functional status of the element, the decrease in content indicates that the neuron is damaged; (2) Creatine (Cr), an energy metabolite, is relatively stable in brain tissue, and is generally used as an internal reference for brain tissue 1 H MRS, and other metabolites are commonly used. The ratio to Cr reflects changes in other metabolites. (3) Choline (CHo) is mainly present in the cell membrane, and its content changes reflect changes in cell membrane metabolism, and its content increases when the cell membrane is degraded or synthesized.
- NAA N-acetylaspartate
- Lactic acid the end product of glycolysis.
- 1 H magnetic resonance spectroscopy may not have a significant Lac peak, but in cerebral ischemia or malignancy, the sugar anaerobic glycolysis process is enhanced and the Lac content may be increased.
- the acquired time domain spectral signal is converted to a frequency domain signal by Fourier transform (FFT transform) by the collected metabolite time domain spectral signal, and then subjected to baseline correction, phase correction and curve.
- FFT transform Fourier transform
- the metabolite is fitted to a processing step such as fitting, and the concentration of the metabolite can be calculated.
- FIG. 2A shows a time domain spectrum signal acquired by the spectral signal acquisition device 110. While Figure 2B shows the frequency domain spectral signal, NAA, Cr, Cho and other types of metabolites or compounds are fitted at different PPMs after the post-treatment step.
- the compound editing interaction system can include, but is not limited to, an editing module 310, a computing module 320, and a display module 330.
- the editing module 310, the computing module 320, and the display module 330 can be implemented by the CPU 240 as shown in FIG.
- the editing module 310 can edit the compounds or metabolites listed in a list to generate a complex that the user needs to view.
- the metabolites in this list may be derived from intermediates or compounds generated by magnetic resonance spectroscopy using specific chemicals including, but not limited to, 1 H, 31 P, 12 C, 23 Na, and 19 F.
- the metabolites in the list can be associated with people associated with magnetic resonance spectroscopy.
- the metabolites in the list can be used to display metabolites associated with magnetic resonance spectroscopy only in patients with this brain tumor, depending on the identity information of the brain tumor patient undergoing magnetic resonance spectroscopy. In some embodiments, this list will intelligently suggest a complex of possible metabolites.
- the compound editing interaction system 130 pops up a pop-up window showing the possible combination of metabolite A and other metabolites B1. , B2, ... and so on.
- the editing module can provide various operators for editing the metabolites listed in the list. For example, add, subtract, multiply, divide, left and right parentheses, delete and confirm operations.
- the editing module can also provide a preview function for the composite of the predetermined edit. For example, the editing module can display selected metabolites and operators in real time, and/or display a predetermined generated complex.
- the editing module can display the composites that need to be edited in a planar or stereoscopic manner.
- the editing module can animate the resulting complex, where the atoms and/or molecules in the metabolite are dynamically simulated with spheres of different colors and radii.
- the calculation module 320 may be directed to the input of a nuclear magnetic resonance anatomical image, and/or the time domain spectral signal of the related metabolite involved in the complex transmitted by the spectral signal acquisition device 110.
- Information to calculate the spatial concentration of the complex in the anatomical image are simultaneously obtained information about the detected object. For example, for a single layer of Chemical Shift Imaging (CSI) 18*14, there are 252 voxels. By calculating the complex concentration on each voxel, the calculation module can generate an 18*14 complex concentration matrix. This complex concentration matrix can be demonstrated by certain visualization methods.
- CSI Chemical Shift Imaging
- the display position of the CSI data on the anatomical image can be located first by the spatial positional relationship of the CSI and the anatomical image.
- CSI can be imagined as a cuboid in three-dimensional space, identified by a vertex and a vector of three edges.
- the display of the CSI on a two-dimensional image in three-dimensional space can be represented by calculating the intersection of the CSI and the two-dimensional image.
- the calculation unit can interpolate the 18*14 concentration matrix to the size of the FOV and anatomize
- the image is fused.
- the specific calculation process can be seen in the flow chart 6 section below.
- Display module 330 can be used to visualize the results of the calculation of the composite, showing spatial variations in the concentration of the composite.
- the calculation module can transmit the result of the fusion of the calculated concentration matrix and the anatomical image to the display module.
- the display module can display the spatial distribution of the concentration within the region of interest.
- the temporal variation of the spatial distribution of the displayed concentrations can also be demonstrated by the display module. For example, by inputting a series of anatomical images with time stamps, after the calculation and concentration processing of the concentration matrix in the region of interest by the calculation module, the fused anatomical images of the series can show the spatial distribution of the complex concentration. Change of time.
- the editing module can include, but is not limited to, one metabolite list unit 410, one operator unit 420, and one preview unit 430.
- the metabolite list unit 410, the operator unit 420, and the preview unit 430 can be implemented by a computer in FIG. 2 by an integral or partial component in the CPU 240.
- the metabolite list unit 410 can display a list of metabolites involved in the magnetic resonance spectroscopy analysis.
- the metabolites in this list can be derived from the use of specific chemicals, including but not limited to 1 H, 31 P, 12 C, 23 Na, and 19 F (or similar to the above), in the middle of magnetic resonance spectroscopy. Or compound.
- specific chemicals including but not limited to 1 H, 31 P, 12 C, 23 Na, and 19 F (or similar to the above), in the middle of magnetic resonance spectroscopy. Or compound.
- the metabolites in this list may also include other compounds that may be involved in brain magnetic resonance. Updates to this list can be real-time or non-real-time.
- This list can be a simple list, for example, the listed metabolites are arranged in a row or column.
- This list can also be a complex list.
- this list can be a drop-down list that expands the sublists by clicking on the intermediate nodes in some lists.
- the metabolites in the list can be associated with people associated with magnetic resonance spectroscopy.
- the metabolites in the list can exhibit metabolites that are only relevant to this guest magnetic resonance spectroscopy analysis, depending on the identity information of the object subjected to magnetic resonance spectroscopy.
- the object referred to herein may be a person, an animal, an organ, a tissue, a body fluid, or a combination similar to the above, or a combination of the above.
- Operator unit 420 can provide various operators for editing the metabolites listed in the list. For example, add, subtract, multiply, divide, left and right parentheses, delete and confirm operations.
- the operators mentioned can also include the logarithm, the exponential, the derivation, the integration, the square, the cube, the fourth power, the square root, the cube root, and the fourth root.
- metabolite list unit 310 and operator unit 320 may also be combined into one unit.
- metabolite unit 310 and operator unit 320 can employ an inverse Polish expression.
- a binary operator is always placed between two operands associated with it (eg, 1+1), so this notation is also called an infix representation.
- the Polish logician J. Lukasiewicz proposed another method of expressing expressions in 1929, called inverse Polish notation. In the inverse Polish notation, all operators are placed. Behind the operand, it is also known as the suffix notation. For example, a formula labeled a+b can be written as a, b, + in the inverse Polish expression.
- Inverse Polish expressions can transform complex expressions into expressions that rely on simple operations to get the result of the calculation. Inverse Polish expressions can be used to perform ordinary expression operations on only two simple operations, stacking and popping. The push and pop stacks here are two standard operations on stack operations in computer data structures.
- the preview unit 430 can display the selected metabolites and operators, and/or display the generated complexes. This display can be real-time or non-real-time.
- the editing module can display the composites that need to be edited in a planar or stereoscopic manner. For example, the editing module can animate the resulting complex, where the atoms and/or molecules in the metabolite are dynamically simulated with spheres of different colors and radii.
- preview unit 430 can provide an intelligent information prompt for the generated composite. For example, the preview unit 430 can prompt whether the generated composite is a composite that is actually present in the physical world. As another example, preview unit 430 can display relevant chemical properties of the generated complex, and/or physiological or pathological information of the subject in question, and the like.
- an editable generated composite unit may also be present in the editing module 310, showing all of the composites that can be combined using the metabolite list unit 310.
- the composite concentration display can be implemented by processor 200.
- processor 200 may first select two or more metabolites and one or more operators to generate a composite based on user input in step 510.
- the user input includes two or more metabolites selected, and one or more operators. The user generates an expected generated complex based on the selected two or more metabolites and one or more operators.
- the user input may also include a time domain spectral signal corresponding to the selected two or more metabolites.
- the time domain spectrum signal may be obtained by scanning by MRI, MRI-PET, or the like, or may be input by the spectrum signal collecting device 110. Alternatively, the time domain spectrum signal can also be obtained by other means.
- the user input may further comprise a frequency domain spectral signal generated by Fourier transform according to the time domain spectral signal.
- the spatial concentration distribution of two or more metabolites selected by the user may also be input in step 510. Step 510 can be implemented by editing module 310.
- Step 520 can include calculating a concentration of the complex based on the generated complex.
- the calculated complex concentration may be a spatial distribution of the complex concentration or a change in the spatial distribution of the complex concentration.
- the inverse Polish algorithm is used to resolve the formulas NAA, Cr, and Cho, and the operators "/", “+”, "(", and ")", Then, based on the analytical result, the concentration of the complex is calculated. For example, on a certain voxel, the concentration of NAA is 43, the concentration of Cr is 30, and the concentration of Cho is 47, and the complex NAA/(Cr+Cho) is on the voxel.
- step 520 can be implemented by computing module 320.
- Step 530 can display the concentration calculation result of the complex to the user.
- the calculated complex concentration can be a spatial distribution of the complex concentration or a spatial distribution of the complex concentration over time.
- the results of the calculations can be presented by visualization techniques to show spatial variations in concentration. Taking a single layer of CSI 18*14 as an example, there are 252 voxels.
- a composite concentration matrix of 18*14 can be generated by the calculation module 320. Through the spatial positional relationship between the CSI and the anatomical image, the display position of the CSI data on the anatomical image and the size of the FOV and the VOI can be located, and the 18*14 complex concentration matrix is interpolated into the FOV size and displayed in fusion with the anatomical image. The concentration change is only displayed in the VOI.
- step 530 can be implemented by display module 330.
- FIG. 6 is a flow diagram of a method of magnetic resonance spectroscopy analysis and interaction of an editable composite, in accordance with some embodiments of the present application.
- Step 610 can select two or more metabolites to be edited by metabolite list unit 410 based on user input.
- the user input may also include a time domain spectral signal corresponding to the selected two or more metabolites.
- the time domain spectrum signal may be obtained by scanning by MRI, MRI-PET, or the like, or may be input by the spectrum signal collecting device 110. Alternatively, the time domain spectrum signal can also be obtained by other means.
- the user input may further include a frequency domain spectrum signal generated by Fourier transform according to the time domain spectrum signal.
- Step 620 can select an operator to participate in the operation from operator unit 420. These operators can be selected by the user, or the various options selected by the compound editing interaction system 130 can be selected by the user. Alternatively, the compound editing interaction system 130 can give a default operator to participate in the operation and prompt the user for the current selection as the default option.
- Step 630 determines whether the currently edited composite meets the requirements. If the currently edited composite does not meet the user's requirements, for example, an erroneous input occurs, or the corresponding composite needs to be determined according to the diagnostic requirements of the detected object, or an abnormal situation needs to be processed, the flow proceeds to step 610 or step 620. Reselect the metabolites and/or operators to be edited.
- the anomaly may be that the edited complex lacks physical or chemical significance, such as non-existence, or the concentration of the metabolite as a denominator is zero, resulting in anomalies in mathematical calculations.
- step 640 the edited composite is added to display module 330.
- the edited composite can be stored in a composite collection unit, such as in database 140.
- the concentration of the edited composite is calculated by calculation unit 320.
- the calculated complex concentration may be a spatial distribution of the complex concentration or a change in the spatial distribution of the complex concentration. Taking a single layer of CSI 18*14 as an example, there are 252 voxels.
- a composite concentration matrix of 18*14 can be generated by the calculation module 320. Through the spatial positional relationship between the CSI and the anatomical image, the display position of the CSI data on the anatomical image and the size of the FOV and the VOI can be located, and the 18*14 complex concentration matrix can be interpolated to the FOV size, and the edited generated composite can be obtained. The distribution of the spatial concentration of the object on the FOV within the anatomical image.
- step 660 it is shown that the distribution of the spatial concentration of the complex calculated in step 650 on the FOV within the anatomical image can be displayed in fusion with the anatomical image, and the concentration distribution of the complex can be displayed in the VOI.
- temporal changes in concentration profile can also be displayed. For example, the above separate processing is performed in chronological order by anatomical images taken at different times to obtain a graph of the concentration distribution of the edited complex in the VOI as a function of time.
- Figure 7 is a spatial distribution image of the composite NAA/(Cr+Cho) concentration obtained according to Flowchart 6.
- the distribution of the complex NAA/(Cr+Cho) values on each voxel is simultaneously displayed on this image. It can be seen that the portion with the highest PPM value in the figure is roughly at the lower left part of the image, indicating that this part may be lesionous. One of the places.
- the present application uses specific words to describe embodiments of the present application.
- a "one embodiment,” “an embodiment,” and/or “some embodiments” means a feature, structure, or feature associated with at least one embodiment of the present application. Therefore, it should be emphasized and noted that “an embodiment” or “an embodiment” or “an alternative embodiment” that is referred to in this specification two or more times in different positions does not necessarily refer to the same embodiment. . Furthermore, some of the features, structures, or characteristics of one or more embodiments of the present application can be combined as appropriate.
- aspects of the present application can be illustrated and described by a number of patentable categories or conditions, including any new and useful process, machine, product, or combination of materials, or Any new and useful improvements. Accordingly, various aspects of the present application can be performed entirely by hardware, entirely by software (including firmware, resident software, microcode, etc.) or by a combination of hardware and software.
- the above hardware or software may be referred to as a "data block,” “module,” “engine,” “unit,” “component,” or “system.”
- aspects of the present application may be embodied in a computer product located in one or more computer readable medium(s) including a computer readable program code.
- a computer readable signal medium may contain a propagated data signal containing a computer program code, for example, on a baseband or as part of a carrier.
- the propagated signal may have a variety of manifestations, including electromagnetic forms, optical forms, and the like, or a suitable combination.
- the computer readable signal medium may be any computer readable medium other than a computer readable storage medium that can be communicated, propagated, or transmitted for use by connection to an instruction execution system, apparatus, or device.
- Computer readable signal The qualitative program code can be propagated through any suitable medium, including radio, cable, fiber optic cable, RF, or similar medium, or a combination of any of the foregoing.
- the computer program code required for the operation of various parts of the application can be written in any one or more programming languages, including object oriented programming languages such as Java, Scala, Smalltalk, Eiffel, JADE, Emerald, C++, C#, VB.NET, Python. Etc., regular programming languages such as C, Visual Basic, Fortran 2003, Perl, COBOL 2002, PHP, ABAP, dynamic programming languages such as Python, Ruby and Groovy, or other programming languages.
- the program code can run entirely on the user's computer, or run as a stand-alone software package on the user's computer, or partially on the user's computer, partly on a remote computer, or entirely on a remote computer or server.
- the remote computer can be connected to the user's computer via any network, such as a local area network (LAN) or wide area network (WAN), or connected to an external computer (eg via the Internet), or in a cloud computing environment, or as a service.
- LAN local area network
- WAN wide area network
- an external computer eg via the Internet
- SaaS software as a service
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Optics & Photonics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Signal Processing (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Radiology & Medical Imaging (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Human Computer Interaction (AREA)
- Chemical & Material Sciences (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
Abstract
一种磁共振波谱交互方法,包括:在基于用户输入选择两个或多个代谢物以及一个或多个运算符生成一个复合物(510)后,计算复合物的浓度及空间分布情况(520),并向用户显示复合物的浓度计算结果(530)。
Description
本申请涉及一种磁共振波谱处理的方法及系统,特别地,涉及一种磁共振波谱交互及分析的方法与系统。
磁共振波谱法是利用磁共振现象和化学位移效应进行特定原子核以及化合物定量分析的一种方法。作为近年来发展起来的一种非侵入性测定化合物代谢的技术,磁共振波谱法通过测定某种化合物或复合物的浓度来检测某种组织的异常。
在目前的波谱分析的处理方式中,计算复合物的交互方式具有界面不够友好,用户操作受限,对复合物的编辑不灵活等缺点,影响用户的诊断体验,降低用户的诊断效率。基于此,设计一种交互友好、流程清晰、可预览已编辑复合物,且可对复合物进行灵活编辑以满足用户诊断需求的交互方法和系统是有必要的。
简述
本申请的一部分附加特性可以在下面的描述中进行说明。通过对以下描述和相应附图的检查或者对实施例的生产或操作的了解,本申请的一部分附加特性对于本领域技术人员是明显的。本披露的特性可以通过对以下描述的具体实施例的各种方面的方法、手段和组合的实践或使用得以实现和达到。
本申请的一些实施例披露了一种磁共振波谱分析系统,包含一种计算机可读的存储媒介,被配置为存储可执行模块,包括编辑模块,被配置为基于用户输入选择两个或多个代谢物以及一个或多个运算符生成一个复合物;计算模块,被配置为计算所述复合物的浓度;以及显示模块,被配置为向用户显示所述复合物的浓度计算结果;以及一个处理器,所述处理器能够执行所述计算机可读的存储媒介存储的可执行模块。
本申请的一些实施例进一步披露了一种磁共振波谱分析系统,其中所述编辑模块包含代谢物列表单元,所述代谢物列表单元被配置为提供多个可选代谢物。
本申请的一些实施例还披露了一种磁共振波谱分析系统,其中所述编辑模块进一步包含运算符单元,所述运算符单元被配置为提供多个可选运算符。
本申请的一些实施例还披露了一种磁共振波谱分析系统,其中所述编辑模块进一步包含预览单元,所述预览单元被配置为显示所述一个或多个化合物以及所述一个或多个运算符。
本申请的一些实施例还披露了一种磁共振波谱分析系统,其中所述计算模块被配置为从编辑模块接受所述复合物,并计算所述复合物的浓度。
本申请的一些实施例还披露了一种磁共振波谱分析系统,其中所述计算模块依据逆波兰式算法解析所述复合物获得所述两个或多个代谢物以及所述一个或多个运算符,并依据所述两个或多个代谢物的单一浓度以及所述一个或多个运算符计算所述复合物的浓度。
本申请的一些实施例还披露了一种磁共振波谱分析系统,其中所述显示模块被配置为显示所述复合物的浓度在检测对象体内的空间分布。
本申请的一些实施例披露了一种磁共振波谱分析方法,包含以下步骤:基于用户输入选择两个或多个代谢物以及一个或多个运算符生成一个复合物;计算所述复合物的浓度;以及向所述用户显示所述复合物的浓度计算结果。
本申请的一些实施例披露了一种磁共振波谱分析方法,其中进一步包括显示所述两个或多个代谢物以及所述一个或多个运算符。
本申请的一些实施例披露了一种磁共振波谱分析方法,其中计算所述复合物的浓度进一步包括依据逆波兰式算法解析所述复合物获得所述两个或多个代谢物以及所述一个或多个运算符,并依据所述两个或多个代谢物的单一浓度以及所述一个或多个运算符计算所述复合物的浓度。
本申请的一些实施例披露了一种磁共振波谱分析方法,其中显示所述复合物的浓度计算结果进一步包括显示所述复合物的浓度在检测对象体内的空间分布。
本申请的一些实施例披露了一种计算机可读的存储媒介,被配置为存储可执行程序,以完成如下的操作:基于用户输入选择两个或多个代谢物以及一个或多个运算符生成一个复合物;计算所述复合物的浓度;向用户显示所述复合物的浓度计算结果;以及一个处理器,所述处理器能够执行所述计算机可读的存储媒
介存储的可执行程序。
本申请的一些实施例披露了一种计算机可读的存储媒介,其中所述可执行程序进一步可以执行显示所述两个或多个代谢物以及所述一个或多个运算符的操作。
本申请的一些实施例披露了一种计算机可读的存储媒介,其中所述可执行程序进一步可以执行如下操作:计算所述复合物的浓度包括依据逆波兰式算法解析所述复合物获得所述两个或多个代谢物以及所述一个或多个运算符,并依据所述两个或多个代谢物的单一浓度以及所述一个或多个运算符计算所述复合物的浓度。
本申请的一些实施例披露了一种计算机可读的存储媒介,其中所述可执行程序进一步可以执行如下操作:显示所述复合物的浓度计算结果进一步包括显示所述复合物的浓度在检测对象体内的空间分布。
附图描述
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单的介绍。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图将本发明应用于其它类似情景。除非从语言环境中显而易见或另做说明,图中相同标号代表相同结构和操作。
图1A是根据本申请的一些实施例所示的一个磁共振波谱处理系统的示意图;
图1B是根据本申请的一些实施例所示的化合物编辑交互系统的一种计算机设备的架构图;
图2A是波谱信号采集设备采集到的时域波谱信号的示例图;
图2B是经过傅里叶变换处理后的频域波谱信号的示例图;
图3是根据本申请的一些实施例所示的化合物编辑交互系统的一种示例性示意图;
图4是根据本申请的一些实施例所示的编辑模块的一种示例性示意图;
图5是根据本申请的一些实施例所示的向用户展示一个所生成复合物浓度的流程图;
图6是根据本申请的一些实施例所示的对可编辑复合物进行磁共振波谱分析与交互的方法的流程图;
图7是根据本申请的一些实施例所示的一种示例复合物按照流程图6得到的复合物浓度的空间分布图像。
具体描述
为了更清楚地说明本申请的实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单的介绍。显而易见地,下面描述中的附图仅仅是本申请的一些示例或实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图将本申请应用于其他类似情景。除非从语言环境中显而易见或另做说明,图中相同标号代表相同结构或操作。
如本申请和权利要求书中所示,除非上下文明确提示例外情形,“一”、“一个”、“一种”和/或“该”等词并非特指单数,也可包括复数。一般说来,术语“包括”与“包含”仅提示包括已明确标识的步骤和元素,而这些步骤和元素不构成一个排它性的罗列,方法或者设备也可能包含其他的步骤或元素。
虽然本申请对根据本申请的实施例的系统中的某些模块做出了各种引用,然而,任何数量的不同模块可以被使用并运行在客户端和/或服务器上。所述模块仅是说明性的,并且所述系统和方法的不同方面可以使用不同模块。
本申请中使用了流程图用来说明根据本申请的实施例的系统所执行的操作。应当理解的是,前面或下面操作不一定按照顺序来精确地执行。相反,可以按照倒序或同时处理各种步骤。同时,也可以将其他操作添加到这些过程中,或从这些过程移除某一步或数步操作。
“扫描区域”代表了进行扫描的实际区域。在一些实施例中,可以对检测对象的感兴趣区域进行扫描。所述感兴趣区域可以是检测对象的整体或其中的一部分,例如,头部、胸部、腹部、心脏、肝脏、上肢、下肢、脊椎、骨骼、血管、病变部位、肿瘤部位等,或者上述部位的任意组合。在一些实施例中,可以
对感兴趣区域的选定层面采集多个时间点的扫描图像。
“元素”代表了图像矩阵中的最小的成分,“体素”代表了实际区域中最小的成分。除非上下文明确提示例外情形,在本申请中图像矩阵中的“元素”和与图像矩阵相对应的实际区域中的“体素”可以表示相同的意思并可以进行替换。
在一些实施例中,本申请一方面涉及一种磁共振波谱分析与交换系统。该磁共振分析与交换系统可以包括编辑模块、计算模块和显示模块。本申请另一方面涉及一种可以显示被编辑复合物的浓度分布或者浓度分布的变化的方法。本申请所述的磁共振波谱交互方法包括基于用户输入对需要编辑生成的复合物进行浓度计算和可视化展示。
本申请的实施例可以应用于不同的图像处理系统。不同的图像处理系统可以包括磁共振成像(Magnetic Resonance Imaging,MRI)系统、核磁共振-正电子发射计算机断层显像混合系统(MR-PET系统)等。
根据本申请的一些实施例,图1A所示的是一个磁共振波谱处理系统的示意图。磁共振波谱分析(Magnetic Resonance Spectroscopy,简称为MRS)作为一种能无创伤地探测活体组织化学特性的方法,可以用于研究人体细胞代谢的病理及生理改变。在许多疾病中,代谢改变先于病理形态改变,而MRS对这种代谢改变的潜在敏感性很高,故能提供信息以早期检测病变。磁共振波谱处理系统100可以包括一个或多个波谱信号采集设备110、一个或多个网络120、一个或多个化合物编辑交互系统130、一个或多个数据库140。
波谱信号采集设备110可以对检测对象进行扫描,得到扫描数据,所述扫描数据可以通过网络120发送到化合物编辑交互系统130进行进一步处理,也可以通过网络120直接存储到数据库140。在一些实施例中,波谱信号采集设备110可以与化合物编辑交互系统130直接连接。所述扫描数据可以直接发送至化合物编辑交互系统130。所述检测对象可以包括人体、动物等。波谱信号采集设备110可以包括但不限于磁共振成像(Magnetic Resonance Imaging,MRI)设备或核磁共振-正电子发射计算机断层显像混合系统(MR-PET系统)等。
化合物编辑交互系统130可以对输入的数据(例如,波谱信号采集设备110得到的和/或数据库140中存储的波谱信号、代谢物浓度数据、用户所选择的
运算符等)进行处理、分析以生成处理结果。例如,化合物编辑交互系统130可以根据所选择的多个代谢物生成复合物。又例如,化合物编辑交互系统130可以对得到的复合物浓度数据进行可视化处理,生成复合物浓度的空间分布图。所述扫描图像可以是二维图像,也可以是三维图像。化合物编辑交互系统130可以包括处理器和输入/输出装置(图中未画出)。在一些实施例中,所述处理器可以是一个服务器,也可以是一个服务器群组。一个服务器群组可以是集中式的,例如数据中心。一个服务器群组也可以是分布式的,例如一个分布式系统。所述处理器可以是云服务器、文件服务器、数据库服务器、File Transfer Protocol(FTP)服务器、应用程序服务器、代理服务器、邮件服务器等中的一种或几种的组合。所述处理器可以是本地的,也可以是远程的。本地处理器可以包括集成在化合物编辑交互系统130中的处理器。远程处理器可以包括通过网络(例如,网络120)与化合物编辑交互系统130连接的处理器。在一些实施例中,所述处理器可以访问数据库140中存储的信息(例如,存储在数据库140中的医学图像)、波谱信号采集设备110中的信息(例如,波谱信号采集设备110拍摄的医学图像)。在一些实施例中,所述输入/输出装置可以向所述处理器输入数据,也可以接收所述处理器输出的数据,并将输出的数据以数字、字符、图像、录像,动画,声音等形式表示出来。在一些实施例中,所述输入/输出装置可以包括但不限于输入装置、输出装置等中的一种或几种的组合。所述输入装置可以包括但不限于字符输入装置(例如,键盘)、光学阅读装置(例如,光学标记阅读机、光学字符阅读机)、图形输入装置(例如,鼠标器、操作杆、光笔)、图像输入装置(例如,摄像机、扫描仪、传真机)、模拟输入装置(例如,语言模数转换识别系统)等中的一种或几种的组合。所述输出装置可以包括但不限于显示装置、打印装置、绘图仪、影像输出装置、语音输出装置、磁记录装置等中的一种或几种的组合。在一些实施例中,化合物编辑交互系统130可以进一步包括存储装置(图中未画出)。所述存储装置可以存储各种信息,例如,程序和数据等。在一些实施例中,化合物编辑交互系统130产生的数据和/或处理结果(例如,解剖图像、频域波谱信号、浓度空间分布图像等)可以存储在数据库140和/或化合物编辑交互系统130的存储装置中,也可以通过输入/输出装置输出。
数据库140可以泛指具有存储功能的设备。数据库140可以存储从波谱信号采集设备110收集的扫描数据和化合物编辑交互系统130工作中产生的各种数据。数据库140可以是本地的,也可以是远程的。本地数据库可以包括集成在数据库140中的具有存储功能的设备。远程数据库可以包括通过网络(例如,网络120)与数据库140连接的具有存储功能的设备。数据库140可以包括但不限于层次式数据库、网络式数据库和关系式数据库等其中的一种或几种的组合。数据库140可以将信息数字化后再以利用电、磁或光学等方式的存储设备加以存储。数据库140可以用来存放各种信息,例如程序和数据等。数据库140可以是利用电能方式存储信息的设备,例如各种存储器、随机存取存储器(Random Access Memory(RAM))、只读存储器(Read Only Memory(ROM))等。其中随机存储器包括但不限于十进计数管、选数管、延迟线存储器、威廉姆斯管、动态随机存储器(DRAM)、静态随机存储器(SRAM)、晶闸管随机存储器(T-RAM)、零电容随机存储器(Z-RAM)等中的一种或几种的组合。只读存储器包括但不限于磁泡存储器、磁钮线存储器、薄膜存储器、磁镀线存储器、磁芯内存、磁鼓存储器、光盘驱动器、硬盘、磁带、早期非易失存储器(NVRAM)、相变化内存、磁阻式随机存储式内存、铁电随机存储内存、非易失SRAM、闪存、电子抹除式可复写只读存储器、可擦除可编程只读存储器、可编程只读存储器、屏蔽式堆读内存、浮动连接门随机存取存储器、纳米随机存储器、赛道内存、可变电阻式内存、可编程金属化单元等中的一种或几种的组合。数据库140可以是利用磁能方式存储信息的设备,例如硬盘、软盘、磁带、磁芯存储器、磁泡存储器、U盘、闪存等。数据库140可以是利用光学方式存储信息的设备,例如CD或DVD等。数据库140可以是利用磁光方式存储信息的设备,例如磁光盘等。数据库140的存取方式可以是随机存储、串行访问存储、只读存储等中的一种或几种的组合。数据库140可以是非永久记忆存储器,也可以是永久记忆存储器。以上提及的存储设备只是列举了一些例子,在磁共振波谱分析100中可以使用的数据库并不局限于此。
网络120可以是单一网络,也可以是多种网络的组合。网络120可以包括但不限于局域网、广域网、公用网络、专用网络、无线局域网、虚拟网络、都市
城域网、公用开关电话网络等中的一种或几种的组合。网络120可以包括多种网络接入点,如有线或无线接入点、基站或网络交换点,通过以上接入点使数据源连接网络120并通过网络发送信息。
需要注意的是,以上对于磁共振波谱分析的描述,仅为示例,并不能把本申请限制在所举实施例范围之内。可以理解的是,对于本领域的技术人员来说,在了解该系统的原理后,可能在不背离这一原理的情况下,对各个模块进行任意组合,或者构成子系统与其他模块连接,对实施上述方法和系统的应用领域形式和细节上的各种修正和改变。例如,在一些实施例中,数据库140可以是具有数据存储功能的云计算平台,包括但不限于公用云、私有云、社区云和混合云等。诸如此类的变形,均在本申请的保护范围之内。
根据本申请的一些实施例,图1B是化合物编辑交互系统130的一种计算机设备的架构图。所述计算机可以是一个通用目的的计算机,或是一个有特定目的的计算机。化合物编辑交互系统130能够被所述计算机设备架构通过其硬件设备、软件程序、固件以及它们的组合所实现。为了示例方便,图1B中只绘制了一台计算机设备的架构,但是化合物编辑交互系统110的相关计算机功能可以由多个计算机设备分布式实施。
所述计算机设备架构可以包括通信端口230,与之相连的可以是实现数据通信的网络(例如,图1A中网络120)。所述计算机设备架构还可以包括一个中央处理系统(CPU)单元240用于执行程序指令,由一个或多个处理器组成。所述计算机设备架构包括一个内部通信总线270,不同形式的程序储存单元以及数据储存单元,例如硬盘210,只读存储器(ROM)250,随机存取存储器(RAM)260,能够被配置为存储所述计算机设备架构处理和/或通信使用的各种数据文件,以及CPU单元240所执行的可能的程序指令。所述计算机设备架构还可以包括一个输入/输出组件220,支持所述计算机设备架构与外部(例如,用户)之间的数据和/或信息交互。所述计算机设备架构也可以通过通信网络接受程序及数据。
磁共振波谱分析可以提供组织的代谢信息。这里所说的组织可以是人或者动物的器官,体液,神经,细胞,或者类似结构,或者以上所举诸例的随意组合。磁共振波谱分析的工作原理是利用不同分子的结构不同所造成的同一磁性原子核
进动频率的差异。具体来说,由于磁性原子核的磁旋比是固定不变的,磁性原子核在外磁场强度中除了受外加静磁场影响外,还受原子核周围的电子云和周围其他原子电子云的影响,使磁性原子核所感受的磁场强度略低于外加静磁场的强度,因而其进动频率也略有降低。同一种磁性原子核如果处于不同的分子中,由于分子化学结构的不同,电子云对磁性原子核的磁屏蔽作用也存在差别,表现出其进动频率的差别。由于所处的分子结构不同造成同一磁性原子核进动频率差异的现象称为化学位移现象。以H质子为例,对目标区域施加经过设计的射频脉冲,其频率范围要求涵盖所要检测代谢物中质子的进动频率。然后采集该区域的MR信号(可以是FID信号或回波信号),该MR信号来源于多种代谢物中的质子,由于化学位移效应,不同的代谢物中质子进动频率有轻微差别,通过傅里叶转换可得到不同物质谱谱线的信息。谱线包括一系列相对比较窄的波峰。其横坐标表示不同物质中质子的进动频率,通常用parts per million(PPM)表示(以标准物的质子进动频率为基准,其他代谢物中质子进动频率与标准物中质子进动频率的差别,以百万分几(parts per million,简写为PPM)来表示)。某一窄波的波峰下面积与目标区域内某特定代谢物的含量,或者浓度,成正比关系。由此可以计算得出代谢物的浓度。(具体描述可以看下文相关段落)例如,可以采集1H、31P、12C、23Na及19F来进行磁共振波谱分析,用于脑肿瘤、脑缺血疾病、前列腺癌、弥漫性肝病、肾脏功能分析、肾移植排斥反应、代谢性疾病的脑改变、脑肿瘤治疗后复发的诊断和鉴别诊断等。在一些实施例中,可以选择一种比较稳定的化学物质作为某种磁性原子核相关代谢物的进动频率参照的标准物,如1H MRS常选用三甲基硅烷作为标准物,31P MRS可采用磷酸肌酸作为标准物。在一些实施例中,由于代谢物的浓度较低,产生的信号几乎是正常水信号的万分之一,所以MRS的敏感性相对较低,需要结合水抑制与脂肪抑制才能得到信号。
大脑1H磁共振波谱分析的代谢物可以包括(1)NAA(N-乙酰门冬氨酸),主要存在于神经元及其轴突,可作为神经元的内标物,其含量可反映神经元的功能状态,含量降低表示神经元受损;(2)肌酸(Cr),能量代谢产物,在脑组织中其浓度相对稳定,一般作为脑组织1H MRS的内参物,常用其他代谢产物与Cr的比值反映其他代谢产物的变化。(3)胆碱(CHo),主要存在于细胞膜,其含
量变化反映细胞膜代谢变化,在细胞膜降解或合成旺盛时其含量增加。在脑肿瘤时,常有Cho升高和NAA降低,因此Cho/NAA升高,尤以恶性肿瘤更为明显。多发硬化等脱髓鞘病变如果Cho升高,往往提示病变活动。(4)乳酸(Lac),为糖酵解的终产物。在一些实施例中,1H磁共振波谱分析可能没有明显的Lac峰,但在脑缺血或恶性肿瘤时,糖无氧酵解过程加强,Lac含量可能会增高。
在一些实施例中,通过将收集到的代谢物时域波谱信号,经过傅里叶变换(FFT变换)将采集到的时域波谱信号转换为频域信号,然后经过基线校正,相位校正和曲线拟合等处理步骤拟合该代谢物,可以计算出该代谢物的浓度。例如,图2A所示为波谱信号采集设备110采集到的时域波谱信号。而图2B所示为频域波谱信号,经过后处理步骤之后在不同PPM处拟合NAA,Cr,Cho及其他类型的代谢物或者化合物。
根据本申请的一些实施例,图3所示的是化合物编辑交互系统130的一种示例性示意图。化合物编辑交互系统可以包括但不限于一个编辑模块310、一个计算模块320、一个显示模块330。编辑模块310、计算模块320和显示模块330可以被如图2中的计算机通过CPU240所实现。
编辑模块310可以对一个列表中所列出的化合物或者代谢物进行编辑,以生成用户需要查看的复合物。这个列表中的代谢物可以来源于使用特定化学物质,包括但不限于1H、31P、12C、23Na及19F,进行磁共振波谱分析生成的中间物或者化合物。在一些实施例中,列表中的代谢物可以与磁共振波谱分析相关的人员有所联系。例如,列表中的代谢物可以依据接受磁共振波谱分析的脑肿瘤病人身份信息的不同,展示出只与此脑肿瘤病人磁共振波谱分析中相关的代谢物。在一些实施例中,这个列表会智能提示代谢物可能组成的复合物。比如,对应列表中的代谢物A,当用户点击或高亮代谢物A所在的位置时,化合物编辑交互系统130会弹出一个弹窗,显示代谢物A与其他代谢物可能组成的复合物组B1,B2,…等等。在一些实施例中,编辑模块可以提供对列表中所列代谢物进行编辑的各种运算符。例如,加,减,乘,除,左右括号,删除和确认操作。编辑模块也可以提供对预定编辑的复合物的预览功能。比如,编辑模块可以实时显示所选择的代谢物和运算符,和/或显示预定生成的复合物。在一些实施例中,编辑模块可以
以平面或者立体的方式展示需要编辑的复合物。例如,编辑模块可以以动画的方式展示所生成的复合物,其中的代谢物中的原子和/或分子以不同颜色和半径的球状物进行动态模拟。
对于编辑模块所生成的复合物,计算模块320可以针对所输入的一幅核磁共振解剖图像,和/或波谱信号采集设备110所传送来的该复合物所涉及的相关代谢物的时域波谱信号信息,来计算出该复合物在该解剖图像中的空间浓度。其中核磁共振解剖图像与时域波谱信号信息为同时获得的关于检测对象的信息。例如,对于一个单层的多体素(Chemical Shift Imaging,CSI)18*14来说,共有252个体素。通过计算每一个体素上的复合物浓度,计算模块可以生成一个18*14的复合物浓度矩阵。这个复合物浓度矩阵可以通过一定的可视化手段进行展示。例如,可以首先通过CSI和解剖图像的空间位置关系,定位CSI数据在解剖图像上的显示位置。比如,CSI在三维空间中可以想象成一个长方体,通过顶点和三条边的向量来标识。CSI在三维空间中的一个二维图像上的显示可以通过计算CSI与该二维图像对应的相交面来表示。
在给定了视域(Field of View,简称为FOV)以及感兴趣区域(Volume of interest,简称为VOI)的大小后,计算单元可以将18*14的浓度矩阵插值为FOV的大小并与解剖图像进行融合。与解剖图像的融合可以指显示出来的值是根据参考图像与浓度彩图根据一定的比例计算。比如,在某点上图像像素值为ImageValue,浓度彩图的值为MetaboliteValue,融合比例为factor,则实际显示出来的值为ActualValue=ImageValue*(1-factor)+MetaboliteValue*factor。具体计算过程可见下文中的流程图6部分。
显示模块330可以用于将该复合物的计算结果可视化,显示该复合物浓度的空间变化。例如,计算模块可以将计算出的浓度矩阵与解剖图像融合后的结果传送给显示模块。显示模块可以在感兴趣区域内显示浓度的空间分布。在一些实施例中,所显示浓度的空间分布的随时间变化也可以通过显示模块进行展示。例如,通过输入带有时间戳的一系列解剖图像,在经过计算模块针对感兴趣区域内的浓度矩阵计算和融合处理后,融合后的该系列解剖图像即可显示复合物浓度的空间分布的随时间变化。
根据本申请的一些实施例,图4所示的是编辑模块310的一种示例性示意图。编辑模块可以包括但不限于一个代谢物列表单元410、一个运算符单元420、一个预览单元430。代谢物列表单元410、运算符单元420和预览单元430可以被如图2中的计算机通过CPU240中的整体或部分组件所实现。
代谢物列表单元410可以将磁共振波谱分析中涉及到的代谢物进行列表展示。这个列表中的代谢物可以来源于使用特定化学物质,包括但不限于1H、31P、12C、23Na及19F(或与以上所举物质类似),进行磁共振波谱分析生成的中间物或者化合物。例如,NAA,Cho,GLu,Cr,ml,Lac,或者类似的代谢物。更进一步地,这个列表中的代谢物也可以包括脑部磁共振可能涉及到的其他化合物。这个列表的更新可以是实时的,也可以是非实时的。这个列表可以是一个简单列表,例如,所列的代谢物以横行或者竖列的方式进行排列。这个列表也可以是一个复杂列表,例如,这个列表可以是一个下拉式列表,通过点击一些列表中的中间节点来展开其中的子列表。在一些实施例中,列表中的代谢物可以与磁共振波谱分析相关的人员有所联系。例如,列表中的代谢物可以依据接受磁共振波谱分析的客体身份信息的不同,展示出只与此客体磁共振波谱分析中相关的代谢物。这里所说的客体,可以是一个人,一种动物,一个器官,一个组织,一种体液,或与以上所举类似,或以上所举情况的一种组合。
运算符单元420可以提供对列表中所列代谢物进行编辑的各种运算符。例如,加,减,乘,除,左右括号,删除和确认操作。所举运算符也可以包括求对数,求指数,求导,求积分,求平方,求立方,求四次方,求平方根,求立方根,求四次方根。更进一步地,运算符单元可以具有编程功能,可以按照用户设定的方式编辑一个将运算符进行组合形成的一个函数。例如,f(x)=ax*x+bx+c就是一个对复合物浓度进行计算的函数。
在一些实施例中,代谢物列表单元310与运算符单元320也可以结合为一个单元。例如,代谢物单元310与运算符单元320可以采用逆波兰式表达法。在通常的表达式中,二元运算符总是置于与之相关的两个运算对象之间(如:1+1),所以这种表示法也称为中缀表示。波兰逻辑学家J.Lukasiewicz于1929年提出了另一种表示表达式的方法,称为逆波兰记法,在逆波兰记法中,所有操作符置于
操作数的后面,因此也被称为后缀表示法。例如,一个标示为a+b的计算式在逆波兰式表达法中可以记为a,b,+。逆波兰表达式可以将复杂表达式转换为依靠简单的操作即可得到计算结果的表达式。逆波兰式表达式可以只用两种简单操作,入栈和出栈就可以达到普通表达式的运算。这里的入栈和出栈为计算机数据结构中关于堆栈操作的两种标准操作。
预览单元430可以显示所选择的代谢物和运算符,和/或显示所生成的复合物。这种显示可以是实时的,也可以是非实时的。在一些实施例中,编辑模块可以以平面或者立体的方式展示需要编辑的复合物。例如,编辑模块可以以动画的方式展示所生成的复合物,其中的代谢物中的原子和/或分子以不同颜色和半径的球状物进行动态模拟。在一些实施例中,预览单元430可以提供对所生成复合物的智能信息提示。比如,预览单元430可以提示所生成的复合物是否为物理世界中真实存在的复合物。又比如,预览单元430可以显示所生成的复合物的相关化学属性,和/或所涉及客体的生理或病理信息,等等。
以上的描述仅仅是本发明的具体实施例,不应被视为是唯一的实施例。显然,对于本领域的专业人员来说,在了解本发明内容和原理后,都可能在不背离本发明原理、结构的情况下,进行形式和细节上的各种修正和改变,但是这些修正和改变仍在本发明的权利要求保护范围之内。例如,可选的,一个可编辑生成的复合物单元也可以出现在编辑模块310中,显示利用代谢物列表单元310所能组合生成的所有复合物。
图5是根据本申请的一些实施例所示的向用户展示一个所生成复合物浓度的流程图。在一些实施例中,所述复合物浓度展示可以由处理器200实现。如图5所示,处理器200可以首先在步骤510中基于用户输入选择两个或多个代谢物以及一个或多个运算符生成一个复合物。在一些实施例中,所述用户输入包括所选择的两个或多个代谢物,以及一个或多个运算符。用户基于所选择的两个或多个代谢物以及一个或多个运算符,来生成一个预期生成的复合物。在一些实施例中,用户输入还可以包括所选择的两个或多个代谢物所对应的时域波谱信号。所述时域波谱信号可以可以通过MRI、MRI-PET等扫描获得,也可以由波谱信号采集设备110输入。可选地,所述时域波谱信号也可以通过其他方式获得。进一
步地,用户输入还可以包括根据时域波谱信号进行傅里叶变换所生成的频域波谱信号。在其他一些实施例中,用户所选择的两个或多个代谢物的空间浓度分布也可以在步骤510中进行输入。步骤510可以由编辑模块310来实现。
步骤520可以包括根据所生成的复合物来计算所述复合物的浓度。所计算的复合物浓度可以是复合物浓度的空间分布,也可以是复合物浓度的空间分布的变化。例如,以复合物NAA/(Cr+Cho)为例,通过逆波兰式算法解析到公式中包含化合物NAA、Cr及Cho,运算符“/”、“+”、“(”及“)”,然后根据解析结果计算复合物的浓度,如在某一个体素上,NAA的浓度为43,Cr浓度为30,Cho浓度为47,则复合物NAA/(Cr+Cho)在此体素上的浓度为43/(30+47),约等于0.56。如果NAA,Cr,Cho的浓度随时间发生了变化,则变化后的复合物NAA/(Cr+Cho)浓度随时间的变化仍可以在步骤520中计算出来。在一些实施例中,步骤520可以通过计算模块320实现。
步骤530可以向用户显示所述复合物的浓度计算结果。如上所描述,所计算的复合物浓度可以是复合物浓度的空间分布,也可以是复合物浓度的空间分布随时间的变化。在一些实施例中,计算结果可以通过可视化技术进行展现,以显示浓度的空间变化。以单层CSI 18*14为例,共有252个体素。通过计算模块320可以生成18*14的复合物浓度矩阵。通过CSI和解剖图像的空间位置关系可以定位CSI数据在解剖图像上的显示位置及FOV,VOI的大小,将18*14的复合物浓度矩阵插值为FOV的大小并与解剖图像进行融合显示,同时只在VOI内显示浓度变化。在一些实施例中,步骤530可以通过显示模块330实现。
图6是根据本申请的一些实施例所示的对可编辑复合物进行磁共振波谱分析与交互的方法的流程图。
步骤610可以基于用户输入,通过代谢物列表单元410选择要编辑的两个或多个代谢物。可选地,用户输入还可以包括所选择的两个或多个代谢物所对应的时域波谱信号。所述时域波谱信号可以可以通过MRI、MRI-PET等扫描获得,也可以由波谱信号采集设备110输入。可选地,所述时域波谱信号也可以通过其他方式获得。进一步地,用户输入还可以包括根据时域波谱信号进行傅里叶变换所生成的频域波谱信号。
步骤620可以从运算符单元420选择要参与运算的运算符。这些运算符可以由用户进行选择,也可以由化合物编辑交互系统130给出选择的各种选项,再由用户进行选择。可选地,化合物编辑交互系统130可以给出缺省的参与运算的运算符,并提示用户目前的选择为缺省选项。
步骤630对目前所编辑的复合物是否符合要求进行判断。如果目前所编辑的复合物不符合用户要求,例如出现误输入,或者需要根据检测对象的诊断要求来决定对应的复合物,或者出现异常情况需要处理,则流程会转至步骤610或者步骤620进行重新选择要参与编辑的代谢物和/或运算符。这里的异常情况可以是所编辑的复合物缺乏物理或化学意义,比如说不存在,或者作为分母的代谢物的浓度为零,从而数学计算上出现异常。
如果目前所编辑的复合物符合用户要求,则流程会转至步骤640,在显示模块330中添加已编辑完成的复合物。已编辑完成的复合物可以存储在一个复合物收集单元,比如放在数据库140中。
在步骤650中,所编辑生成的复合物的浓度由计算单元320进行计算得出。所计算的复合物浓度可以是复合物浓度的空间分布,也可以是复合物浓度的空间分布的变化。以单层CSI 18*14为例,共有252个体素。通过计算模块320可以生成18*14的复合物浓度矩阵。通过CSI和解剖图像的空间位置关系可以定位CSI数据在解剖图像上的显示位置及FOV,VOI的大小,将18*14的复合物浓度矩阵插值为FOV的大小,则可得到所编辑生成的复合物在解剖图像内FOV上的空间浓度的分布情况。
在步骤660中,显示,步骤650所计算的到的复合物在解剖图像内FOV上的空间浓度的分布情况可以与解剖图像进行融合显示,则可以在VOI内显示复合物的浓度分布情况。在一些实施例中,浓度分布的时间变化也可以被显示出来。例如,通过将不同时间拍摄的解剖图像按时间次序进行以上的分别处理以得到所编辑复合物在VOI内的浓度分布随时间变化图。
图7是按照流程图6得到的复合物NAA/(Cr+Cho)浓度的空间分布图像。在此图像上同时显示每个体素上的复合物NAA/(Cr+Cho)值分布情况。可以看到图中PPM值最高的部分大致在图像的左下部位,标示这一部位可能是有病灶的
地方之一。
以上的描述仅仅是本发明的具体实施例,不应被视为是唯一的实施例。显然,对于本领域的专业人员来说,在了解本发明内容和原理后,都可能在不背离本发明原理、结构的情况下,进行形式和细节上的各种修正和改变,但是这些修正和改变仍在本发明的权利要求保护范围之内。
上文已对基本概念做了描述,显然,对于本领域技术人员来说,上述发明披露仅仅作为示例,而并不构成对本申请的限定。虽然此处并没有明确说明,本领域技术人员可能会对本申请进行各种修改、改进和修正。该类修改、改进和修正在本申请中被建议,所以该类修改、改进、修正仍属于本申请示范实施例的精神和范围。
同时,本申请使用了特定词语来描述本申请的实施例。如“一个实施例”、“一实施例”、和/或“一些实施例”意指与本申请至少一个实施例相关的某一特征、结构或特点。因此,应强调并注意的是,本说明书中在不同位置两次或多次提及的“一实施例”或“一个实施例”或“一替代性实施例”并不一定是指同一实施例。此外,本申请的一个或多个实施例中的某些特征、结构或特点可以进行适当的组合。
此外,本领域技术人员可以理解,本申请的各方面可以通过若干具有可专利性的种类或情况进行说明和描述,包括任何新的和有用的工序、机器、产品或物质的组合,或对他们的任何新的和有用的改进。相应地,本申请的各个方面可以完全由硬件执行、可以完全由软件(包括固件、常驻软件、微码等)执行、也可以由硬件和软件组合执行。以上硬件或软件均可被称为“数据块”、“模块”、“引擎”、“单元”、“组件”或“系统”。此外,本申请的各方面可能表现为位于一个或多个计算机可读介质中的计算机产品,该产品包括计算机可读程序编码。
计算机可读信号介质可能包含一个内含有计算机程序编码的传播数据信号,例如在基带上或作为载波的一部分。该传播信号可能有多种表现形式,包括电磁形式、光形式等等、或合适的组合形式。计算机可读信号介质可以是除计算机可读存储介质之外的任何计算机可读介质,该介质可以通过连接至一个指令执行系统、装置或设备以实现通讯、传播或传输供使用的程序。位于计算机可读信号介
质上的程序编码可以通过任何合适的介质进行传播,包括无线电、电缆、光纤电缆、RF、或类似介质、或任何上述介质的组合。
本申请各部分操作所需的计算机程序编码可以用任意一种或多种程序语言编写,包括面向对象编程语言如Java、Scala、Smalltalk、Eiffel、JADE、Emerald、C++、C#、VB.NET、Python等,常规程序化编程语言如C语言、Visual Basic、Fortran 2003、Perl、COBOL 2002、PHP、ABAP,动态编程语言如Python、Ruby和Groovy,或其他编程语言等。该程序编码可以完全在用户计算机上运行、或作为独立的软件包在用户计算机上运行、或部分在用户计算机上运行部分在远程计算机运行、或完全在远程计算机或服务器上运行。在后种情况下,远程计算机可以通过任何网络形式与用户计算机连接,比如局域网(LAN)或广域网(WAN),或连接至外部计算机(例如通过因特网),或在云计算环境中,或作为服务使用如软件即服务(SaaS)。
此外,除非权利要求中明确说明,本申请所述处理元素和序列的顺序、数字字母的使用、或其他名称的使用,并非用于限定本申请流程和方法的顺序。尽管上述披露中通过各种示例讨论了一些目前认为有用的发明实施例,但应当理解的是,该类细节仅起到说明的目的,附加的权利要求并不仅限于披露的实施例,相反,权利要求旨在覆盖所有符合本申请实施例实质和范围的修正和等价组合。例如,虽然以上所描述的系统组件可以通过硬件设备实现,但是也可以只通过软件的解决方案得以实现,如在现有的服务器或移动设备上安装所描述的系统。
同理,应当注意的是,为了简化本申请披露的表述,从而帮助对一个或多个发明实施例的理解,前文对本申请实施例的描述中,有时会将多种特征归并至一个实施例、附图或对其的描述中。但是,这种披露方法并不意味着本申请对象所需要的特征比权利要求中提及的特征多。实际上,实施例的特征要少于上述披露的单个实施例的全部特征。
一些实施例中使用了描述成分、属性数量的数字,应当理解的是,此类用于实施例描述的数字,在一些示例中使用了修饰词“大约”、“近似”或“大体上”来修饰。除非另外说明,“大约”、“近似”或“大体上”表明所述数字允许有±20%的变化。相应地,在一些实施例中,说明书和权利要求中使用的数值参数均为近
似值,该近似值根据个别实施例所需特点可以发生改变。在一些实施例中,数值参数应考虑规定的有效数位并采用一般位数保留的方法。尽管本申请一些实施例中用于确认其范围广度的数值域和参数为近似值,在具体实施例中,此类数值的设定在可行范围内尽可能精确。
针对本申请引用的每个专利、专利申请、专利申请公开物和其他材料,如文章、书籍、说明书、出版物、文档、物件等,特此将其全部内容并入本申请作为参考。与本申请内容不一致或产生冲突的申请历史文件除外,对本申请权利要求最广范围有限制的文件(当前或之后附加于本申请中的)也除外。需要说明的是,如果本申请附属材料中的描述、定义、和/或术语的使用与本申请所述内容有不一致或冲突的地方,以本申请的描述、定义和/或术语的使用为准。
最后,应当理解的是,本申请中所述实施例仅用以说明本申请实施例的原则。其他的变形也可能属于本申请的范围。因此,作为示例而非限制,本申请实施例的替代配置可视为与本申请的教导一致。相应地,本申请的实施例不仅限于本申请明确介绍和描述的实施例。
Claims (20)
- 一种磁共振波谱分析系统,包含:一种计算机可读的存储媒介,被配置为存储可执行模块,包括:编辑模块,被配置为基于用户输入选择两两个或多个代谢物以及一个或多个运算符生成一个复合物;计算模块,被配置为计算所述复合物的浓度;以及显示模块,被配置为向用户显示所述复合物的浓度计算结果;以及一个处理器,所述处理器能够执行所述计算机可读的存储媒介存储的可执行模块。
- 依据权利要求1的一个系统,其中所述编辑模块包含代谢物列表单元,所述代谢物列表单元被配置为提供多个可选代谢物。
- 依据权利要求2的一个系统,其中所述编辑模块进一步包含运算符单元,所述运算符单元被配置为提供多个可选运算符。
- 依据权利要求3的一个系统,其中所述编辑模块进一步包含预览单元,所述预览单元被配置为显示所述一个或多个化合物以及所述一个或多个运算符。
- 依据权利要求1的一个系统,其中所述计算模块被配置为从编辑模块接收所述复合物,并计算所述复合物的浓度。
- 依据权利要求5的一个系统,其中所述计算模块依据逆波兰式算法解析所述复合物获得所述两个或多个代谢物以及所述一个或多个运算符,并依据所述两个两个或多个代谢物的单一浓度以及所述一个或多个运算符计算所述复合物的浓度。
- 依据权利要求1的一个系统,其中所述显示模块被配置为显示所述复合物的浓度在检测对象体内的空间分布。
- 一种磁共振波谱分析方法,包含以下步骤:基于用户输入选择两个或多个代谢物以及一个或多个运算符生成一个复合物;计算所述复合物的浓度;以及显示所述复合物的浓度计算结果。
- 依据权利要求8的一种方法,进一步包括显示所述两两个或多个代谢物以及所述一个或多个运算符。
- 依据权利要求8的一种方法,其中计算所述复合物的浓度进一步包括依据逆波兰式算法解析所述复合物获得所述两两个或多个代谢物以及所述一个或多个运算符,并依据所述两两个或多个代谢物的单一浓度以及所述一个或多个运算符计算所述复合物的浓度。
- 依据权利要求10的一种方法,其中显示所述复合物的浓度计算结果进一步包括显示所述复合物的浓度在检测对象体内的空间分布。
- 依据权利要求11的一种方法,其中所述用户输入包括一幅或者多幅核磁共振解剖图像。
- 依据权利要求12的一种方法,其中所述用户输入包括所述复合物涉及的相关代谢物的时域波谱信号信息。
- 依据权利要求13的一种方法,其中所述核磁共振解剖图像与所述时域波谱信号信息为同时获得的关于所述检测对象的信息。
- 一种计算机可读的存储媒介,被配置为存储可执行程序,以完成如下的操作:基于用户输入选择两两个或多个代谢物以及一个或多个运算符生成一个复合物; 计算所述复合物的浓度;以及向用户显示所述复合物的浓度计算结果;以及一个处理器,所述处理器能够执行所述计算机可读的存储媒介存储的可执行程序。
- 依据权利要求15的一种存储媒介,其中所述可执行程序进一步可以执行显示所述两两个或多个代谢物以及所述一个或多个运算符的操作。
- 依据权利要求15的一种存储媒介,其中所述可执行程序进一步可以执行如下操作:计算所述复合物的浓度包括依据逆波兰式算法解析所述复合物获得所述两两个或多个代谢物以及所述一个或多个运算符,并依据所述两两个或多个代谢物的单一浓度以及所述一个或多个运算符计算所述复合物的浓度。
- 依据权利要求17的一种存储媒介,其中所述可执行程序进一步可以执行如下操作:显示所述复合物的浓度计算结果进一步包括显示所述复合物的浓度在检测对象体内的空间分布。
- 依据权利要求18的一种存储媒介,其中所述可执行程序进一步可以执行如下操作:所述用户输入包括一幅或者多幅核磁共振解剖图像。
- 依据权利要求19的一种存储媒介,其中所述可执行程序进一步可以执行如下操作:所述用户输入包括所述复合物涉及的相关代谢物的时域波谱信号信息。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2017/088425 WO2018227469A1 (zh) | 2017-06-15 | 2017-06-15 | 磁共振波谱交互方法、系统及计算机可读的存储媒介 |
EP17913892.0A EP3637094A4 (en) | 2017-06-15 | 2017-06-15 | MAGNETIC RESONANCE SPECTROSCOPY INTERACTION METHOD AND SYSTEM, AND COMPUTER-READABLE STORAGE MEDIUM |
US16/710,055 US20200110142A1 (en) | 2017-06-15 | 2019-12-11 | Methods, systems, and computer-readable storage media for interaction in magnetic resonance spectroscopy |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2017/088425 WO2018227469A1 (zh) | 2017-06-15 | 2017-06-15 | 磁共振波谱交互方法、系统及计算机可读的存储媒介 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US16/710,055 Continuation US20200110142A1 (en) | 2017-06-15 | 2019-12-11 | Methods, systems, and computer-readable storage media for interaction in magnetic resonance spectroscopy |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2018227469A1 true WO2018227469A1 (zh) | 2018-12-20 |
Family
ID=64658766
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2017/088425 WO2018227469A1 (zh) | 2017-06-15 | 2017-06-15 | 磁共振波谱交互方法、系统及计算机可读的存储媒介 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20200110142A1 (zh) |
EP (1) | EP3637094A4 (zh) |
WO (1) | WO2018227469A1 (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130149388A1 (en) * | 2010-01-04 | 2013-06-13 | Joseph Louie MCCLAY | Metabolomics-based biomarkers for lung function |
US20140156573A1 (en) * | 2011-07-27 | 2014-06-05 | The Research Foundation Of State University Of New York | Methods for generating predictive models for epithelial ovarian cancer and methods for identifying eoc |
CN106153739A (zh) * | 2015-03-25 | 2016-11-23 | 上海博云生物科技有限公司 | 一种基于代谢组学技术诊断肝癌的模型 |
US20170023575A1 (en) * | 2013-08-21 | 2017-01-26 | Purdue Research Foundation | Identification of blood based metabolite biomarkers of pancreatic cancer |
EP3136100A2 (en) * | 2015-08-31 | 2017-03-01 | The University of Hong Kong | Pleural fluid markers for malignant pleural effusions |
CN106600579A (zh) * | 2016-11-30 | 2017-04-26 | 上海联影医疗科技有限公司 | 磁共振图像波谱分析方法及装置 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6963769B1 (en) * | 1999-04-21 | 2005-11-08 | The United States Of America As Represented By The Secretary Of The Department Of Health And Human Services | Method for enhancing contrast produced by MRI |
US7432926B1 (en) * | 2004-11-23 | 2008-10-07 | Cherkas Barry M | Complete function graphing system and method |
JP5608650B2 (ja) * | 2008-08-15 | 2014-10-15 | コーニンクレッカ フィリップス エヌ ヴェ | 磁気共鳴分光画像データを用いたpet又はspect放射線撮像システムの減衰補正 |
DE102011076776A1 (de) * | 2011-05-31 | 2012-12-06 | Siemens Aktiengesellschaft | Verfahren zu einer Bilddatenaufnahme |
DE102013226244A1 (de) * | 2013-12-17 | 2015-06-18 | Siemens Aktiengesellschaft | Medizintechnische Steuerung |
DE102015204628B4 (de) * | 2015-03-13 | 2023-09-07 | Siemens Healthcare Gmbh | Verfahren zur Planung einer medizinischen Bildgebungsuntersuchung |
-
2017
- 2017-06-15 EP EP17913892.0A patent/EP3637094A4/en active Pending
- 2017-06-15 WO PCT/CN2017/088425 patent/WO2018227469A1/zh unknown
-
2019
- 2019-12-11 US US16/710,055 patent/US20200110142A1/en not_active Abandoned
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130149388A1 (en) * | 2010-01-04 | 2013-06-13 | Joseph Louie MCCLAY | Metabolomics-based biomarkers for lung function |
US20140156573A1 (en) * | 2011-07-27 | 2014-06-05 | The Research Foundation Of State University Of New York | Methods for generating predictive models for epithelial ovarian cancer and methods for identifying eoc |
US20170023575A1 (en) * | 2013-08-21 | 2017-01-26 | Purdue Research Foundation | Identification of blood based metabolite biomarkers of pancreatic cancer |
CN106153739A (zh) * | 2015-03-25 | 2016-11-23 | 上海博云生物科技有限公司 | 一种基于代谢组学技术诊断肝癌的模型 |
EP3136100A2 (en) * | 2015-08-31 | 2017-03-01 | The University of Hong Kong | Pleural fluid markers for malignant pleural effusions |
CN106600579A (zh) * | 2016-11-30 | 2017-04-26 | 上海联影医疗科技有限公司 | 磁共振图像波谱分析方法及装置 |
Non-Patent Citations (1)
Title |
---|
See also references of EP3637094A4 * |
Also Published As
Publication number | Publication date |
---|---|
EP3637094A1 (en) | 2020-04-15 |
US20200110142A1 (en) | 2020-04-09 |
EP3637094A4 (en) | 2020-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Glatard et al. | A virtual imaging platform for multi-modality medical image simulation | |
Filograna et al. | Identification of the most significant magnetic resonance imaging (MRI) radiomic features in oncological patients with vertebral bone marrow metastatic disease: a feasibility study | |
Li et al. | Blooming artifact reduction in coronary artery calcification by a new de-blooming algorithm: initial study | |
US11696701B2 (en) | Systems and methods for estimating histological features from medical images using a trained model | |
Gan et al. | Extracting three-dimensional orientation and tractography of myofibers using optical coherence tomography | |
He et al. | Automatic segmentation and quantification of epicardial adipose tissue from coronary computed tomography angiography | |
US7876940B2 (en) | Universal image processing | |
Tiago et al. | A data augmentation pipeline to generate synthetic labeled datasets of 3D echocardiography images using a GAN | |
Braig et al. | Preclinical 4D-flow magnetic resonance phase contrast imaging of the murine aortic arch | |
Schaefferkoetter et al. | Deep learning for whole-body medical image generation | |
US10964074B2 (en) | System for harmonizing medical image presentation | |
Weng et al. | Deep learning-based segmentation of the lung in MR-images acquired by a stack-of-spirals trajectory at ultra-short echo-times | |
Lim et al. | CardiSort: a convolutional neural network for cross vendor automated sorting of cardiac MR images | |
Sendra-Balcells et al. | Domain generalization in deep learning for contrast-enhanced imaging | |
Benameur et al. | Interpretation of cardiac wall motion from cine-MRI combined with parametric imaging based on the Hilbert transform | |
Warner et al. | Supervised segmentation of polycystic kidneys: a new application for stereology data | |
Bones et al. | Workflow for automatic renal perfusion quantification using ASL‐MRI and machine learning | |
Wang et al. | Fusion of color Doppler and magnetic resonance images of the heart | |
CN107389720B (zh) | 磁共振波谱交互方法及系统 | |
Gulamhussene et al. | Transfer-learning is a key ingredient to fast deep learning-based 4D liver MRI reconstruction | |
WO2018227469A1 (zh) | 磁共振波谱交互方法、系统及计算机可读的存储媒介 | |
Sun et al. | Left atrium wall-mapping application for wall thickness visualisation | |
Nicol et al. | Echocardiographic evaluation of left ventricular function using an automated analysis algorithm is feasible for beginners and experts: comparison with invasive and non-invasive methods | |
Ravi et al. | Accelerated MRI using intelligent protocolling and subject-specific denoising applied to Alzheimer's disease imaging | |
Maillot et al. | Automated inversion time selection for black-blood late gadolinium enhancement cardiac imaging in clinical practice |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 17913892 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2017913892 Country of ref document: EP Effective date: 20191209 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |