CN117607770B - Magnetic resonance image reconstruction method, device, electronic device and storage medium - Google Patents
Magnetic resonance image reconstruction method, device, electronic device and storage medium Download PDFInfo
- Publication number
- CN117607770B CN117607770B CN202410085787.XA CN202410085787A CN117607770B CN 117607770 B CN117607770 B CN 117607770B CN 202410085787 A CN202410085787 A CN 202410085787A CN 117607770 B CN117607770 B CN 117607770B
- Authority
- CN
- China
- Prior art keywords
- information
- dimensional contour
- space information
- view
- specific field
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 230000035945 sensitivity Effects 0.000 claims abstract description 92
- 238000012545 processing Methods 0.000 claims abstract description 37
- 239000011159 matrix material Substances 0.000 claims description 108
- 230000009466 transformation Effects 0.000 claims description 31
- 238000006243 chemical reaction Methods 0.000 claims description 21
- 239000003550 marker Substances 0.000 claims description 14
- 230000015654 memory Effects 0.000 claims description 11
- 238000001514 detection method Methods 0.000 claims description 9
- 238000007405 data analysis Methods 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 19
- 238000004590 computer program Methods 0.000 description 13
- 230000006870 function Effects 0.000 description 13
- 238000002595 magnetic resonance imaging Methods 0.000 description 12
- 230000008569 process Effects 0.000 description 9
- 238000004891 communication Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 4
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 3
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 3
- 101000827703 Homo sapiens Polyphosphoinositide phosphatase Proteins 0.000 description 2
- 102100023591 Polyphosphoinositide phosphatase Human genes 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000004806 packaging method and process Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 101100012902 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) FIG2 gene Proteins 0.000 description 1
- 101100233916 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) KAR5 gene Proteins 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013475 authorization Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 229910000365 copper sulfate Inorganic materials 0.000 description 1
- ARUVKPQLZAKDPS-UHFFFAOYSA-L copper(II) sulfate Chemical compound [Cu+2].[O-][S+2]([O-])([O-])[O-] ARUVKPQLZAKDPS-UHFFFAOYSA-L 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- 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
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
-
- 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
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明提供了一种磁共振图像重建方法、装置、电子设备及存储介质。该方法包括:获取扫描对象在第一特定视野范围内扫描切片的第一三维轮廓空间信息和磁共振扫描数据;对第一三维轮廓空间信息进行处理,得到处理后的第一三维轮廓空间信息;将处理后的第一三维轮廓空间信息输入至已构建的线圈灵敏度估计模型,得到与扫描对象在第一特定视野范围内扫描切片对应的多通道线圈灵敏度信息;根据多通道线圈灵敏度信息和磁共振扫描数据,对扫描对象进行磁共振图像重建。
The present invention provides a magnetic resonance image reconstruction method, device, electronic device and storage medium. The method comprises: obtaining first three-dimensional contour space information and magnetic resonance scanning data of a scanned slice of a scanned object within a first specific field of view; processing the first three-dimensional contour space information to obtain processed first three-dimensional contour space information; inputting the processed first three-dimensional contour space information into a constructed coil sensitivity estimation model to obtain multi-channel coil sensitivity information corresponding to the scanned slice of the scanned object within the first specific field of view; and reconstructing a magnetic resonance image of the scanned object according to the multi-channel coil sensitivity information and the magnetic resonance scanning data.
Description
技术领域Technical Field
本发明涉及磁共振图像重建领域,具体地,涉及一种磁共振图像重建方法、装置、电子设备及存储介质。The present invention relates to the field of magnetic resonance image reconstruction, and in particular to a magnetic resonance image reconstruction method, device, electronic equipment and storage medium.
背景技术Background Art
磁共振成像技术是目前比较成熟的成像技术,特别是磁共振并行成像技术,使用多个具有唯一灵敏度图的接收线圈采集信息,加速数据采集,减少成像时间。Magnetic resonance imaging technology is a relatively mature imaging technology at present, especially magnetic resonance parallel imaging technology, which uses multiple receiving coils with unique sensitivity maps to collect information, accelerate data acquisition, and reduce imaging time.
磁共振并行成像技术一般可以基于图像域、K空间域或混合域等方法实现图像重建。基于图像域的方法中,需要额外的扫描或预扫描以获取精准线圈灵敏度图,再基于精准线圈灵敏度图,灵活结合图像先验(例如,图像稀疏性)进行图像重建,但这会使得在实际操作中纳入额外的扫描时间;基于混合域的方法,需要K空间域中心校准数据,利用特征值分解从截断提取的子空间中衍生线圈灵敏度图,并进行掩码操作以在空间上定义图像背景区域,但该方法需要手动调整阈值来截断子空间以及选择特征值来掩盖特征向量图,手动调整阈值不同,以致于衍生线圈灵敏度图不同,重建图像效果不同;基于K空间的方法,也需要K空间中心校准数据,限制了加速因子,同时不利于结合图像稀疏性的先验知识。Magnetic resonance parallel imaging technology can generally achieve image reconstruction based on methods such as image domain, K-space domain or hybrid domain. In the image domain-based method, additional scanning or pre-scanning is required to obtain accurate coil sensitivity maps, and then based on the accurate coil sensitivity maps, image reconstruction is flexibly combined with image priors (for example, image sparsity), but this will cause additional scanning time to be included in actual operations; the hybrid domain-based method requires K-space domain center calibration data, uses eigenvalue decomposition to derive coil sensitivity maps from truncated extracted subspaces, and performs mask operations to spatially define the image background area, but this method requires manual adjustment of thresholds to truncate subspaces and select eigenvalues to mask eigenvector maps. Different manually adjusted thresholds result in different derived coil sensitivity maps and different reconstructed image effects; the K-space-based method also requires K-space center calibration data, which limits the acceleration factor and is not conducive to combining prior knowledge of image sparsity.
发明内容Summary of the invention
有鉴于此,本发明提供了一种磁共振图像重建方法、装置、电子设备及存储介质。In view of this, the present invention provides a magnetic resonance image reconstruction method, device, electronic device and storage medium.
本发明的一方面提供了一种磁共振图像重建方法,包括:获取扫描对象在第一特定视野范围内扫描切片的第一三维轮廓空间信息和磁共振扫描数据;对第一三维轮廓空间信息进行处理,得到处理后的第一三维轮廓空间信息;将处理后的第一三维轮廓空间信息输入至已构建的线圈灵敏度估计模型,得到与扫描对象在第一特定视野范围内扫描切片对应的多通道线圈灵敏度信息;根据多通道线圈灵敏度信息和磁共振扫描数据,对扫描对象进行磁共振图像重建。One aspect of the present invention provides a magnetic resonance image reconstruction method, comprising: obtaining first three-dimensional contour space information and magnetic resonance scanning data of a scanned slice of a scanned object within a first specific field of view; processing the first three-dimensional contour space information to obtain processed first three-dimensional contour space information; inputting the processed first three-dimensional contour space information into a constructed coil sensitivity estimation model to obtain multi-channel coil sensitivity information corresponding to the scanned slice of the scanned object within the first specific field of view; and reconstructing a magnetic resonance image of the scanned object based on the multi-channel coil sensitivity information and the magnetic resonance scanning data.
可选地,在获取扫描对象在第一特定视野范围内扫描切片的第一三维轮廓空间信息和磁共振扫描数据之前,方法还包括:确定第二特定视野范围的坐标系与第一特定视野范围的坐标系之间的第一转换矩阵,其中,第一转换矩阵应用于获取扫描对象在第一特定视野范围内扫描切片的第一三维轮廓空间信息。Optionally, before acquiring first three-dimensional contour space information and magnetic resonance scanning data of scanned slices of the scanned object within the first specific field of view, the method also includes: determining a first transformation matrix between a coordinate system of a second specific field of view and a coordinate system of the first specific field of view, wherein the first transformation matrix is applied to acquiring first three-dimensional contour space information of scanned slices of the scanned object within the first specific field of view.
可选地,确定第二特定视野范围的坐标系与第一特定视野范围的坐标系之间的第一转换矩阵,包括:获取配准装置中目标对象在第二特定视野范围内的第一位置信息,目标对象为可被磁共振扫描识别的标志物;利用预先设置的梯度扫描方向,对配准装置中目标对象在第一特定视野范围内进行扫描,获取与目标对象对应的扫描图像;对扫描图像进行连通域检测处理,得到与扫描图像对应的第二位置信息;根据第一位置信息和第二位置信息,确定第二特定视野范围的坐标系与第一特定视野范围的坐标系之间的第一转换矩阵。Optionally, determining a first transformation matrix between a coordinate system of a second specific field of view and a coordinate system of a first specific field of view includes: obtaining first position information of a target object in a registration device within the second specific field of view, where the target object is a marker that can be identified by magnetic resonance scanning; scanning the target object in the registration device within the first specific field of view using a preset gradient scanning direction to obtain a scanned image corresponding to the target object; performing connected domain detection processing on the scanned image to obtain second position information corresponding to the scanned image; and determining a first transformation matrix between a coordinate system of the second specific field of view and a coordinate system of the first specific field of view based on the first position information and the second position information.
可选地,获取扫描对象在第一特定视野范围内扫描切片的第一三维轮廓空间信息,包括:获取扫描对象在第二特定视野范围内的三维轮廓空间参考信息;确定扫描对象在第一特定视野范围内扫描切片的第一定位信息;根据第一定位信息和三维轮廓空间参考信息,确定扫描对象在第二特定视野范围内扫描切片的第二三维轮廓空间信息;对第二三维轮廓空间信息进行转换处理,得到扫描对象在第一特定视野范围内扫描切片的第一三维轮廓空间信息。Optionally, obtaining first three-dimensional contour space information of scanned slices of the scanned object within a first specific field of view includes: obtaining three-dimensional contour space reference information of the scanned object within a second specific field of view; determining first positioning information of the scanned slices of the scanned object within the first specific field of view; determining second three-dimensional contour space information of the scanned slices of the scanned object within the second specific field of view based on the first positioning information and the three-dimensional contour space reference information; and converting the second three-dimensional contour space information to obtain first three-dimensional contour space information of the scanned slices of the scanned object within the first specific field of view.
可选地,根据第一定位信息和三维轮廓空间参考信息,确定扫描对象在第二特定视野范围内扫描切片的第二三维轮廓空间信息,包括:根据第一定位信息和第一转换矩阵,从三维轮廓空间参考信息中提取与第一定位信息对应的第二定位信息,其中,第一转换矩阵是通过第二特定视野范围的坐标系与第一特定视野范围的坐标系确定的,第二定位信息为扫描对象在第二特定视野范围内扫描切片的信息;根据第二定位信息,确定扫描对象在第二特定视野范围内扫描切片的第二三维轮廓空间信息。Optionally, based on the first positioning information and the three-dimensional contour space reference information, second three-dimensional contour space information of the scanned slices of the scanned object within the second specific field of view is determined, including: based on the first positioning information and the first transformation matrix, second positioning information corresponding to the first positioning information is extracted from the three-dimensional contour space reference information, wherein the first transformation matrix is determined by the coordinate system of the second specific field of view and the coordinate system of the first specific field of view, and the second positioning information is information of the scanned slices of the scanned object within the second specific field of view; based on the second positioning information, second three-dimensional contour space information of the scanned slices of the scanned object within the second specific field of view is determined.
可选地,对第一三维轮廓空间信息进行处理,得到处理后的第一三维轮廓空间信息,包括:确定扫描切片的第一三维轮廓空间信息与目标扫描范围的空间信息之间的第二转化矩阵;根据第二转化矩阵和第一三维轮廓空间信息,得到与第一三维轮廓空间信息对应的第三三维轮廓空间信息,其中,第三三维轮廓空间信息表征扫描对象在第一特定视野范围内扫描切片的三维轮廓空间信息的唯一性;对第三三维轮廓空间信息进行数据分析处理,得到处理后的第一三维轮廓空间信息。Optionally, the first three-dimensional contour space information is processed to obtain the processed first three-dimensional contour space information, including: determining a second conversion matrix between the first three-dimensional contour space information of the scanned slice and the spatial information of the target scanning range; obtaining third three-dimensional contour space information corresponding to the first three-dimensional contour space information based on the second conversion matrix and the first three-dimensional contour space information, wherein the third three-dimensional contour space information characterizes the uniqueness of the three-dimensional contour space information of the scanned slice of the scanned object within the first specific field of view; and performing data analysis processing on the third three-dimensional contour space information to obtain the processed first three-dimensional contour space information.
可选地,将处理后的第一三维轮廓空间信息输入至已构建的线圈灵敏度估计模型,得到与扫描对象在特定视野范围内扫描切片对应的多通道线圈灵敏度信息,包括:根据第一特征矩阵信息、第二特征矩阵信息和第一三维轮廓空间信息,得到与扫描对象在特定视野范围内扫描切片对应的多通道线圈灵敏度信息,其中,第一特征矩阵是通过历史多通道灵敏度信息确定,第二特征矩阵是通过历史第一三维轮廓空间信息确定。Optionally, the processed first three-dimensional contour space information is input into the constructed coil sensitivity estimation model to obtain multi-channel coil sensitivity information corresponding to the scanned slices of the scanned object within a specific field of view, including: obtaining the multi-channel coil sensitivity information corresponding to the scanned slices of the scanned object within a specific field of view based on the first characteristic matrix information, the second characteristic matrix information and the first three-dimensional contour space information, wherein the first characteristic matrix is determined by historical multi-channel sensitivity information, and the second characteristic matrix is determined by historical first three-dimensional contour space information.
本发明的另一方面提供了一种磁共振图像重建装置,包括:获取模块,用于获取扫描对象在第一特定视野范围内扫描切片的第一三维轮廓空间信息和磁共振扫描数据;处理模块,用于对第一三维轮廓空间信息进行处理,得到处理后的第一三维轮廓空间信息;输入模块,用于将处理后的第一三维轮廓空间信息输入至已构建的线圈灵敏度估计模型,得到与扫描对象在第一特定视野范围内扫描切片对应的多通道线圈灵敏度信息;重建模块,用于根据多通道线圈灵敏度信息和磁共振扫描数据,对扫描对象进行磁共振图像重建。Another aspect of the present invention provides a magnetic resonance image reconstruction device, comprising: an acquisition module, used to acquire first three-dimensional contour space information and magnetic resonance scanning data of a scanned slice of a scanned object within a first specific field of view; a processing module, used to process the first three-dimensional contour space information to obtain processed first three-dimensional contour space information; an input module, used to input the processed first three-dimensional contour space information into a constructed coil sensitivity estimation model to obtain multi-channel coil sensitivity information corresponding to the scanned slice of the scanned object within the first specific field of view; and a reconstruction module, used to reconstruct a magnetic resonance image of the scanned object based on the multi-channel coil sensitivity information and the magnetic resonance scanning data.
本发明的另一方面提供了一种电子设备,包括:一个或多个处理器;存储器,用于存储一个或多个程序,其中,当一个或多个程序被一个或多个处理器执行时,使得一个或多个处理器执行上述的方法。Another aspect of the present invention provides an electronic device, comprising: one or more processors; and a memory for storing one or more programs, wherein when the one or more programs are executed by the one or more processors, the one or more processors execute the above method.
本发明的另一方面提供了一种计算机可读存储介质,存储有计算机可执行指令,指令在被执行时用于实现上述的方法。Another aspect of the present invention provides a computer-readable storage medium storing computer-executable instructions, which are used to implement the above method when executed.
本发明的另一方面提供了一种计算机程序产品,计算机程序产品包括计算机可执行指令,指令在被执行时用于实现上述的方法。Another aspect of the present invention provides a computer program product, which includes computer executable instructions, and the instructions are used to implement the above method when executed.
通过将得到的第一三维轮廓空间信息输入至已构建的线圈灵敏度估计模型,得到与扫描对象在第一特定视野范围内扫描切片对应的多通道线圈灵敏度信息;根据多通道线圈灵敏度信息和磁共振扫描数据,对扫描对象进行磁共振图像重建。实现了在不需要中心校准区域的情况下,能够利用已构建的线圈灵敏度估计模型快速得到在第一特定视野范围内扫描切片对应的多通道线圈灵敏度信息,并基于得到的多通道线圈灵敏度信息实现对扫描对象图像的快速重建,减小在实际操作中的扫描时间,提高了无校准磁共振图像重建速度。By inputting the obtained first three-dimensional contour space information into the constructed coil sensitivity estimation model, the multi-channel coil sensitivity information corresponding to the scan slice of the scanned object within the first specific field of view is obtained; and the magnetic resonance image of the scanned object is reconstructed according to the multi-channel coil sensitivity information and the magnetic resonance scanning data. It is achieved that the multi-channel coil sensitivity information corresponding to the scan slice within the first specific field of view can be quickly obtained by using the constructed coil sensitivity estimation model without the need for a central calibration area, and the image of the scanned object can be quickly reconstructed based on the obtained multi-channel coil sensitivity information, thereby reducing the scanning time in actual operation and improving the reconstruction speed of the non-calibrated magnetic resonance image.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
通过以下参照附图对本发明实施例的描述,本发明的上述以及其他目的、特征和优点将更为清楚,在附图中:The above and other objects, features and advantages of the present invention will become more apparent through the following description of the embodiments of the present invention with reference to the accompanying drawings, in which:
图1示出了根据本发明实施例的可以应用磁共振图像重建方法和装置的示例性系统架构;FIG1 shows an exemplary system architecture to which a magnetic resonance image reconstruction method and apparatus may be applied according to an embodiment of the present invention;
图2示出了根据本发明实施例的磁共振图像重建方法的流程图;FIG2 shows a flow chart of a magnetic resonance image reconstruction method according to an embodiment of the present invention;
图3示出了根据本发明实施例的扫描切片轮廓三维示意图;FIG3 shows a three-dimensional schematic diagram of a scan slice profile according to an embodiment of the present invention;
图4(a)示出了根据本发明实施例的扫描对象的VTK坐标系的示意图;FIG4( a ) shows a schematic diagram of a VTK coordinate system of a scanned object according to an embodiment of the present invention;
图4(b)示出了根据本发明实施例的扫描对象的平面轮廓的示意图;FIG4( b ) is a schematic diagram showing a plane profile of a scanned object according to an embodiment of the present invention;
图5示出了根据本发明实施例的配准装置的示意图;FIG5 shows a schematic diagram of a registration device according to an embodiment of the present invention;
图6示出了根据本发明实施例的确定第一转换矩阵的流程示意图;FIG6 shows a schematic diagram of a process of determining a first conversion matrix according to an embodiment of the present invention;
图7示出了根据本发明实施例的获取第一三维轮廓空间信息的方法流程图;FIG7 shows a flow chart of a method for obtaining first three-dimensional contour space information according to an embodiment of the present invention;
图8示出了根据本发明的实施例的磁共振图像重建装置的框图;以及FIG8 shows a block diagram of a magnetic resonance image reconstruction apparatus according to an embodiment of the present invention; and
图9示出了根据本发明实施例的适于实现磁共振图像重建方法的电子设备的框图。FIG. 9 shows a block diagram of an electronic device suitable for implementing a magnetic resonance image reconstruction method according to an embodiment of the present invention.
具体实施方式DETAILED DESCRIPTION
以下,将参照附图来描述本发明的实施例。但是应该理解,这些描述只是示例性的,而并非要限制本发明的范围。在下面的详细描述中,为便于解释,阐述了许多具体的细节以提供对本发明实施例的全面理解。然而,明显地,一个或多个实施例在没有这些具体细节的情况下也可以被实施。此外,在以下说明中,省略了对公知结构和技术的描述,以避免不必要地混淆本发明的概念。Below, embodiments of the present invention will be described with reference to the accompanying drawings. However, it should be understood that these descriptions are exemplary only and are not intended to limit the scope of the present invention. In the following detailed description, for ease of explanation, many specific details are set forth to provide a comprehensive understanding of embodiments of the present invention. However, it is apparent that one or more embodiments may also be implemented without these specific details. In addition, in the following description, descriptions of known structures and technologies are omitted to avoid unnecessary confusion of concepts of the present invention.
在此使用的术语仅仅是为了描述具体实施例,而并非意在限制本发明。在此使用的术语“包括”、“包含”等表明了特征、步骤、操作和/或部件的存在,但是并不排除存在或添加一个或多个其他特征、步骤、操作或部件。The terms used herein are only for describing specific embodiments and are not intended to limit the present invention. The terms "comprise", "include", etc. used herein indicate the existence of features, steps, operations and/or components, but do not exclude the existence or addition of one or more other features, steps, operations or components.
在此使用的所有术语(包括技术和科学术语)具有本领域技术人员通常所理解的含义,除非另外定义。应注意,这里使用的术语应解释为具有与本说明书的上下文相一致的含义,而不应以理想化或过于刻板的方式来解释。All terms (including technical and scientific terms) used herein have the meanings commonly understood by those skilled in the art unless otherwise defined. It should be noted that the terms used herein should be interpreted as having a meaning consistent with the context of this specification and should not be interpreted in an idealized or overly rigid manner.
在使用类似于“A、B和C等中至少一个”这样的表述的情况下,一般来说应该按照本领域技术人员通常理解该表述的含义来予以解释(例如,“具有A、B和C中至少一个的系统”应包括但不限于单独具有A、单独具有B、单独具有C、具有A和B、具有A和C、具有B和C、和/或具有A、B、C的系统等)。When using expressions such as "at least one of A, B, and C, etc.", they should generally be interpreted according to the meaning of the expression commonly understood by those skilled in the art (for example, "a system having at least one of A, B, and C" should include but is not limited to a system having A alone, B alone, C alone, A and B, A and C, B and C, and/or A, B, C, etc.).
在本发明的实施例中,所涉及的数据(例如,包括但不限于用户个人信息)的收集、更新、分析、处理、使用、传输、提供、发明、存储等方面,均符合相关法律法规的规定,被用于合法的用途,且不违背公序良俗。特别地,对用户个人信息采取了必要措施,防止对用户个人信息数据的非法访问,维护用户个人信息安全和网络安全。In the embodiments of the present invention, the collection, updating, analysis, processing, use, transmission, provision, invention, storage and other aspects of the data involved (for example, including but not limited to user personal information) are in compliance with the provisions of relevant laws and regulations, are used for legitimate purposes, and do not violate public order and good morals. In particular, necessary measures are taken for user personal information to prevent illegal access to user personal information data and maintain the security of user personal information and network security.
在本发明的实施例中,在获取或采集用户个人信息之前,均获取了用户的授权或同意。In the embodiment of the present invention, before obtaining or collecting the user's personal information, the user's authorization or consent is obtained.
本发明的实施例提供了一种磁共振图像重建方法,包括:获取扫描对象在第一特定视野范围内扫描切片的第一三维轮廓空间信息和磁共振扫描数据;对第一三维轮廓空间信息进行处理,得到处理后的第一三维轮廓空间信息;将处理后的第一三维轮廓空间信息输入至已构建的线圈灵敏度估计模型,得到与扫描对象在第一特定视野范围内扫描切片对应的多通道线圈灵敏度信息;根据多通道线圈灵敏度信息和磁共振扫描数据,对扫描对象进行磁共振图像重建。An embodiment of the present invention provides a magnetic resonance image reconstruction method, comprising: acquiring first three-dimensional contour space information and magnetic resonance scanning data of a scanned slice of a scanned object within a first specific field of view; processing the first three-dimensional contour space information to obtain processed first three-dimensional contour space information; inputting the processed first three-dimensional contour space information into a constructed coil sensitivity estimation model to obtain multi-channel coil sensitivity information corresponding to the scanned slice of the scanned object within the first specific field of view; and reconstructing a magnetic resonance image of the scanned object according to the multi-channel coil sensitivity information and the magnetic resonance scanning data.
图1示出了根据本发明实施例的可以应用磁共振图像重建方法和装置的示例性系统架构。需要注意的是,图1所示仅为可以应用本发明实施例的系统架构的示例,以帮助本领域技术人员理解本发明的技术内容,但并不意味着本发明实施例不可以用于其他设备、系统、环境或场景。FIG1 shows an exemplary system architecture to which the magnetic resonance image reconstruction method and apparatus can be applied according to an embodiment of the present invention. It should be noted that FIG1 is only an example of a system architecture to which the embodiment of the present invention can be applied, to help those skilled in the art understand the technical content of the present invention, but does not mean that the embodiment of the present invention cannot be used in other devices, systems, environments or scenarios.
如图1所示,根据该实施例的系统架构可以包括磁共振成像设备101、第一终端设备102、第二终端设备103,网络104和服务器105。网络104用以在磁共振成像设备101、第一终端设备102、第二终端设备103和服务器105之间提供通信链路的介质。网络104可以包括各种连接类型,例如有线和/或无线通信链路等等。As shown in Fig. 1, the system architecture according to this embodiment may include a magnetic resonance imaging device 101, a first terminal device 102, a second terminal device 103, a network 104 and a server 105. The network 104 is used to provide a medium for communication links between the magnetic resonance imaging device 101, the first terminal device 102, the second terminal device 103 and the server 105. The network 104 may include various connection types, such as wired and/or wireless communication links, etc.
磁共振图像设备101产生检测的特定信号,并通过磁共振图像设备中的接收线圈(图中未示出)将特定信号传输给第一终端设备102、第二终端设备103,第一终端设备102、第二终端设备103通过网络104与服务器105交互,以接收或发送消息等。第一终端设备102、第二终端设备103上可以安装有各种通讯客户端应用,例如购物类应用、网页浏览器应用、搜索类应用、即时通信工具、邮箱客户端和/或社交平台软件等(仅为示例)。The magnetic resonance imaging device 101 generates a specific signal for detection, and transmits the specific signal to the first terminal device 102 and the second terminal device 103 through a receiving coil (not shown in the figure) in the magnetic resonance imaging device. The first terminal device 102 and the second terminal device 103 interact with the server 105 through the network 104 to receive or send messages, etc. Various communication client applications may be installed on the first terminal device 102 and the second terminal device 103, such as shopping applications, web browser applications, search applications, instant messaging tools, email clients and/or social platform software, etc. (only as examples).
第一终端设备102、第二终端设备103可以是具有显示屏并且支持网页浏览的各种电子设备,包括但不限于智能手机、平板电脑、膝上型便携计算机和台式计算机等等。The first terminal device 102 and the second terminal device 103 may be various electronic devices having display screens and supporting web browsing, including but not limited to smart phones, tablet computers, laptop computers, desktop computers, and the like.
服务器105可以是提供各种服务的服务器,例如对用户利用第一终端设备102、第二终端设备103所浏览的网站提供支持的后台管理服务器(仅为示例)。后台管理服务器可以对接收到的用户请求等数据进行分析等处理,并将处理结果(例如根据用户请求获取或生成的网页、信息、或数据等)反馈给终端设备。The server 105 may be a server that provides various services, such as a background management server (only as an example) that provides support for websites browsed by users using the first terminal device 102 and the second terminal device 103. The background management server may analyze and process received data such as user requests, and feed back processing results (such as web pages, information, or data obtained or generated according to user requests) to the terminal device.
需要说明的是,本发明实施例所提供的磁共振图像重建方法一般可以由服务器105执行。相应地,本发明实施例所提供的磁共振图像重建装置一般可以设置于服务器105中。本发明实施例所提供的磁共振图像重建方法也可以由不同于服务器105且能够与第一终端设备102、第二终端设备103和/或服务器105通信的服务器或服务器集群执行。相应地,本发明实施例所提供的磁共振图像重建装置也可以设置于不同于服务器105且能够与第一终端设备102、第二终端设备103和/或服务器105通信的服务器或服务器集群中。或者,本发明实施例所提供的磁共振图像重建方法也可以由第一终端设备102或第二终端设备103执行,或者也可以由不同于第一终端设备102或第二终端设备103的其他终端设备执行。相应地,本发明实施例所提供的磁共振图像重建装置也可以设置于第一终端设备102或第二终端设备103中,或设置于不同于第一终端设备102或第二终端设备103的其他终端设备中。It should be noted that the magnetic resonance image reconstruction method provided in the embodiment of the present invention can generally be executed by the server 105. Accordingly, the magnetic resonance image reconstruction device provided in the embodiment of the present invention can generally be set in the server 105. The magnetic resonance image reconstruction method provided in the embodiment of the present invention can also be executed by a server or server cluster that is different from the server 105 and can communicate with the first terminal device 102, the second terminal device 103 and/or the server 105. Accordingly, the magnetic resonance image reconstruction device provided in the embodiment of the present invention can also be set in a server or server cluster that is different from the server 105 and can communicate with the first terminal device 102, the second terminal device 103 and/or the server 105. Alternatively, the magnetic resonance image reconstruction method provided in the embodiment of the present invention can also be executed by the first terminal device 102 or the second terminal device 103, or can also be executed by other terminal devices different from the first terminal device 102 or the second terminal device 103. Accordingly, the magnetic resonance image reconstruction device provided in the embodiment of the present invention can also be set in the first terminal device 102 or the second terminal device 103, or can be set in other terminal devices different from the first terminal device 102 or the second terminal device 103.
应该理解,图1中的终端设备、网络和服务器的数目仅仅是示意性的。根据实现需要,可以具有任意数目的终端设备、网络和服务器。It should be understood that the number of terminal devices, networks and servers in Figure 1 is only illustrative. Any number of terminal devices, networks and servers may be provided according to implementation requirements.
图2示出了根据本发明实施例的磁共振图像重建方法的流程图。FIG. 2 shows a flow chart of a magnetic resonance image reconstruction method according to an embodiment of the present invention.
如图2所示,该方法包括操作S210~S240。As shown in FIG. 2 , the method includes operations S210 to S240 .
在操作S210,获取扫描对象在第一特定视野范围内扫描切片的第一三维轮廓空间信息和磁共振扫描数据。In operation S210, first three-dimensional contour space information and magnetic resonance scanning data of a scanned slice of a scanned object within a first specific field of view are acquired.
根据本发明的实施例,扫描对象可以是人体,可以是动物,也可以是指定区域,如身体的特定部分,例如,头部、胸部、腹部,或者是其任一组合部位,在本发明中不做限定。According to an embodiment of the present invention, the scanning object can be a human body, an animal, or a specified area, such as a specific part of the body, for example, the head, chest, abdomen, or any combination thereof, which is not limited in the present invention.
根据本发明的实施例,扫描切片可以为通过高速扫描将扫描对象对应的图像进行切片处理得到的。According to an embodiment of the present invention, the scan slice may be obtained by slicing an image corresponding to the scan object through high-speed scanning.
图3示出了根据本发明实施例的扫描切片轮廓三维示意图。FIG. 3 shows a three-dimensional schematic diagram of a scan slice contour according to an embodiment of the present invention.
如图3,以扫描对象头部为示例,得出的针对头部轮廓的扫描切片的三维数据示意图。在不同方向,分别可以获取多个扫描切片的轮廓三维数据。一个扫描对象可以对象多个扫描切片。As shown in Figure 3, taking the head of a scanned object as an example, a schematic diagram of the three-dimensional data of the scan slice for the head contour is obtained. In different directions, the contour three-dimensional data of multiple scan slices can be obtained respectively. A scanned object can have multiple scan slices.
根据本发明的实施例,第一特定视野范围可以为操作人员选定的需要对扫描对象进行扫描的扫描区域,可以认为是操作人员选定的FOV(field of view)范围。According to an embodiment of the present invention, the first specific field of view range may be a scanning area selected by an operator where the scanning object needs to be scanned, and may be considered as a FOV (field of view) range selected by the operator.
根据本发明的实施例,第一三维轮廓空间信息可以为在第一特定视野范围内获得的扫描对象的扫描切片在VTK(The Visualization Toolkit)坐标下的轮廓空间信息。According to an embodiment of the present invention, the first three-dimensional contour space information may be contour space information of a scan slice of the scan object obtained within the first specific field of view in VTK (The Visualization Toolkit) coordinates.
例如,图4(a)示出了根据本发明实施例的扫描对象的VTK坐标系的示意图;图4(b)示出了根据本发明实施例的扫描对象的平面轮廓的示意图。For example, FIG. 4( a ) shows a schematic diagram of a VTK coordinate system of a scanned object according to an embodiment of the present invention; and FIG. 4( b ) shows a schematic diagram of a plane contour of a scanned object according to an embodiment of the present invention.
结合图4(a)及4(b)可知,以扫描对象中心作为原点,扫描对象左侧指向的方向作为X轴正向,扫描对象背后的方向作为Y轴的正向,扫描对象的头部的方向为Z轴的正向。在VTK坐标系下,图4(a)中X-Y平面构成对应图4(b)中L-P平面,图4(a)中X-Z平面构成对应图4(b)中L-S平面,图4(a)中Y-Z平面构成对应图4(b)中P-S平面。Combining Figures 4 (a) and 4 (b), we can see that the center of the scanned object is taken as the origin, the direction of the left side of the scanned object is taken as the positive direction of the X axis, the direction behind the scanned object is taken as the positive direction of the Y axis, and the direction of the head of the scanned object is taken as the positive direction of the Z axis. In the VTK coordinate system, the X-Y plane in Figure 4 (a) corresponds to the L-P plane in Figure 4 (b), the X-Z plane in Figure 4 (a) corresponds to the L-S plane in Figure 4 (b), and the Y-Z plane in Figure 4 (a) corresponds to the P-S plane in Figure 4 (b).
根据本发明的实施例,磁共振扫描数据可以为扫描对象在进行扫描时的采集的欠采样数据以及欠采样方式数据。According to an embodiment of the present invention, the magnetic resonance scanning data may be under-sampling data and under-sampling mode data collected when the scanning object is scanned.
在操作S220,对第一三维轮廓空间信息进行处理,得到处理后的第一三维轮廓空间信息。In operation S220, the first three-dimensional contour space information is processed to obtain processed first three-dimensional contour space information.
在操作S230,将处理后的第一三维轮廓空间信息输入至已构建的线圈灵敏度估计模型,得到与扫描对象在第一特定视野范围内扫描切片对应的多通道线圈灵敏度信息。In operation S230, the processed first three-dimensional contour space information is input into the constructed coil sensitivity estimation model to obtain multi-channel coil sensitivity information corresponding to the scan slice of the scanned object within the first specific field of view.
根据本发明的实施例,可以将第一三维轮廓空间信息转变到最大扫描范围坐标系下,以使得任何扫描对象的扫描切片的轮廓信息都能在该最大扫描范围坐标系下予以表示第一三维轮廓空间信息。According to an embodiment of the present invention, the first three-dimensional contour space information can be transformed into a maximum scanning range coordinate system, so that contour information of a scanned slice of any scanned object can be represented by the first three-dimensional contour space information in the maximum scanning range coordinate system.
根据本发明的实施例,将转变到最大扫描范围坐标系下的第一三维轮廓空间信息转换成已构建的线圈灵敏度估计模型可识别的输入信号,从而得到处理后的第一三维轮廓空间信息。According to an embodiment of the present invention, the first three-dimensional contour space information converted to the maximum scanning range coordinate system is converted into an input signal recognizable by the constructed coil sensitivity estimation model, thereby obtaining the processed first three-dimensional contour space information.
根据本发明的实施例,线圈灵敏度估计模型可以通过对输入的第一三维轮廓空间信息进行处理,从而快速得到多通道线圈灵敏度信息。According to an embodiment of the present invention, the coil sensitivity estimation model can quickly obtain multi-channel coil sensitivity information by processing the input first three-dimensional contour space information.
根据本发明的实施例,针对扫描对象的每个扫描切片,通过线圈灵敏度估计模型均可以获得一个多通道线圈灵敏度信息。多通道线圈灵敏度信息可以表征接收线圈对输入信号的响应程度。其值越高,表示检测微弱信号的能力越强。According to an embodiment of the present invention, for each scan slice of the scanned object, a multi-channel coil sensitivity information can be obtained through the coil sensitivity estimation model. The multi-channel coil sensitivity information can characterize the response degree of the receiving coil to the input signal. The higher the value, the stronger the ability to detect weak signals.
在操作S240,根据多通道线圈灵敏度信息和磁共振扫描数据,对扫描对象进行磁共振图像重建。In operation S240, a magnetic resonance image is reconstructed for the scanned object according to the multi-channel coil sensitivity information and the magnetic resonance scan data.
根据本发明的实施例,利用已构建的线圈灵敏度估计模型快速获取扫描切片的多通道线圈灵敏度信息,并结合磁共振扫描数据,可以进行图像重建,得到最终重建图像。According to an embodiment of the present invention, the constructed coil sensitivity estimation model is used to quickly obtain multi-channel coil sensitivity information of a scan slice, and image reconstruction can be performed in combination with magnetic resonance scan data to obtain a final reconstructed image.
根据本发明的实施例,磁共振扫描数据可以包括欠采样数据y,欠采样方式数据P,多通道线圈灵敏度信息可以表示为C,可以利用下式确定重建图像,即,式(1):According to an embodiment of the present invention, the magnetic resonance scanning data may include undersampling data y , undersampling mode data P , and the multi-channel coil sensitivity information may be expressed as C. The reconstructed image may be determined using the following formula, that is, formula (1):
(1); (1);
其中,F为傅里叶变换,Ψ为稀疏变换,m为所求的重建图像。Where F is the Fourier transform, Ψ is the sparse transform, and m is the desired reconstructed image.
根据本发明的实施例,通过将得到的第一三维轮廓空间信息输入至已构建的线圈灵敏度估计模型,得到与扫描对象在第一特定视野范围内扫描切片对应的多通道线圈灵敏度信息;根据多通道线圈灵敏度信息和磁共振扫描数据,对扫描对象进行磁共振图像重建。实现了在不需要中心校准区域的情况下,能够利用已构建的线圈灵敏度估计模型快速得到在第一特定视野范围内扫描切片对应的多通道线圈灵敏度信息,并基于得到的多通道线圈灵敏度信息实现对扫描对象图像的快速重建,减小在实际操作中的扫描时间,提高了无校准磁共振图像重建速度。According to an embodiment of the present invention, the obtained first three-dimensional contour space information is input into the constructed coil sensitivity estimation model to obtain the multi-channel coil sensitivity information corresponding to the scan slice of the scanned object within the first specific field of view; and the magnetic resonance image of the scanned object is reconstructed according to the multi-channel coil sensitivity information and the magnetic resonance scanning data. It is achieved that the multi-channel coil sensitivity information corresponding to the scan slice within the first specific field of view can be quickly obtained by using the constructed coil sensitivity estimation model without the need for a central calibration area, and the image of the scanned object can be quickly reconstructed based on the obtained multi-channel coil sensitivity information, thereby reducing the scanning time in actual operation and improving the reconstruction speed of the non-calibrated magnetic resonance image.
根据本发明的实施例,在获取扫描对象在第一特定视野范围内扫描切片的第一三维轮廓空间信息和磁共振扫描数据之前,该方法还包括:确定第二特定视野范围的坐标系与第一特定视野范围的坐标系之间的第一转换矩阵,其中,第一转换矩阵应用于获取扫描对象在第一特定视野范围内扫描切片的第一三维轮廓空间信息。According to an embodiment of the present invention, before acquiring first three-dimensional contour space information and magnetic resonance scanning data of scanned slices of the scanned object within a first specific field of view, the method also includes: determining a first transformation matrix between a coordinate system of a second specific field of view and a coordinate system of the first specific field of view, wherein the first transformation matrix is applied to acquiring first three-dimensional contour space information of scanned slices of the scanned object within the first specific field of view.
根据本发明的实施例,确定第二特定视野范围的坐标系与第一特定视野范围的坐标系之间的第一转换矩阵,包括:获取配准装置中目标对象在第二特定视野范围内的第一位置信息,目标对象为可被磁共振扫描识别的标志物;利用预先设置的梯度扫描方向,对配准装置中目标对象在第一特定视野范围内进行扫描,获取与目标对象对应的扫描图像;对扫描图像进行连通域检测处理,得到与扫描图像对应的第二位置信息;根据第一位置信息和第二位置信息,确定第二特定视野范围的坐标系与第一特定视野范围的坐标系之间的第一转换矩阵。According to an embodiment of the present invention, determining a first transformation matrix between a coordinate system of a second specific field of view and a coordinate system of a first specific field of view includes: acquiring first position information of a target object in a registration device within the second specific field of view, where the target object is a marker that can be identified by magnetic resonance scanning; scanning the target object in the registration device within the first specific field of view using a preset gradient scanning direction to acquire a scanned image corresponding to the target object; performing connected domain detection processing on the scanned image to obtain second position information corresponding to the scanned image; and determining a first transformation matrix between the coordinate system of the second specific field of view and the coordinate system of the first specific field of view based on the first position information and the second position information.
根据本发明的实施例,第二特定视野范围可以为扫描对象在深度摄像设备下的扫描区域。第二特定视野范围的坐标系可以为扫描对象在深度摄像设备下的坐标系,以下可以简称为“深度摄像设备坐标系”。According to an embodiment of the present invention, the second specific field of view may be a scanning area of the scanned object under the depth camera device. The coordinate system of the second specific field of view may be a coordinate system of the scanned object under the depth camera device, which may be referred to as the "depth camera device coordinate system" hereinafter.
根据本发明的实施例,深度摄像设备可以为磁兼容深度摄像设备,例如,可以为深度摄像机、深度照相机等,本发明在此不做限定。According to an embodiment of the present invention, the depth camera device may be a magnetically compatible depth camera device, for example, a depth camera, a depth still camera, etc., and the present invention is not limited thereto.
根据本发明的实施例,第一转换矩阵可以为从深度摄像设备坐标系转换到VTK坐标系的矩阵,以实现深度摄像设备坐标系与VTK坐标系的配准。According to an embodiment of the present invention, the first conversion matrix may be a matrix for converting from the depth camera device coordinate system to the VTK coordinate system, so as to achieve registration between the depth camera device coordinate system and the VTK coordinate system.
根据本发明的实施例,该操作可以基于预先设置在磁共振成像设备上的配准装置来确定第一转换矩阵。According to an embodiment of the present invention, the operation may determine the first transformation matrix based on a registration device pre-set on the magnetic resonance imaging device.
例如,图5示出了根据本发明实施例的配准装置的示意图。For example, FIG5 shows a schematic diagram of a registration device according to an embodiment of the present invention.
如图5所示,该配准装置500是由配准物501、深度摄像设备502和标志物503构成。将深度摄像设备502镶嵌在特殊的配准物501上,同时配准物501上还镶嵌多个可被磁共振成像设备504扫描识别的标志物503,其中,以配准物501位参考对象,深度摄像设备502和多个标志物503之间的相对位置是固定的。配准装置500可以固定于磁共振成像设备504的扫描腔内。As shown in FIG5 , the registration device 500 is composed of a registration object 501, a depth camera device 502 and a marker 503. The depth camera device 502 is embedded in a special registration object 501, and a plurality of markers 503 that can be scanned and identified by a magnetic resonance imaging device 504 are also embedded in the registration object 501, wherein the registration object 501 is used as a reference object, and the relative positions between the depth camera device 502 and the plurality of markers 503 are fixed. The registration device 500 can be fixed in the scanning cavity of the magnetic resonance imaging device 504.
根据本发明的实施例,配准物可以为具有特定形状的装置,例如,可以为长方体或正方体的装置,标志物可以为硫酸铜溶液。According to an embodiment of the present invention, the registration object may be a device having a specific shape, for example, a device in the shape of a cuboid or a cube, and the marker may be a copper sulfate solution.
根据本发明的实施例,目标对象可以为该配准装置中可被磁共振成像设备扫描识别的标志物。将标志物镶嵌在配准物上,与配准物上的深度摄像设备具有相对固定位置,可以先前获得目标对象在深度摄像设备的视野范围内的坐标,也就是第二特定视野范围内的第一位置信息,可以表示为MarksInCamera。According to an embodiment of the present invention, the target object may be a marker in the registration device that can be scanned and identified by a magnetic resonance imaging device. The marker is embedded on the registration object and has a relatively fixed position with the depth camera on the registration object. The coordinates of the target object within the field of view of the depth camera can be obtained in advance, that is, the first position information within the second specific field of view, which can be expressed as MarksInCamera .
根据本发明的实施例,可以将配准装置固定在磁共振成像设备扫描腔体内,利用收发一体线圈,通过预先设置的梯度扫描方向(G x 、G y 、G z ),对配准装置中的目标对象在第一特定视野范围内进行扫描,可以获取与目标对象对应的扫描图像。According to an embodiment of the present invention, the registration device can be fixed in the scanning cavity of the magnetic resonance imaging device, and the target object in the registration device can be scanned within a first specific field of view by using a transceiver-in-one coil through a preset gradient scanning direction ( Gx , Gy , Gz ), so that a scanned image corresponding to the target object can be obtained.
根据本发明的实施例,与目标对象对应的扫描图像可以为目标对象为亮点的背景噪声图像。According to an embodiment of the present invention, the scanned image corresponding to the target object may be a background noise image in which the target object is a bright spot.
根据本发明的实施例,连通域检测处理可以为Blob检测,可以通过对背景噪声图像进行Blob检测,得到目标对象在VTK坐标系下的坐标,也就是第一特定视野范围内的第二位置信息,可以表示为MarksInVTK。According to an embodiment of the present invention, the connected domain detection process may be Blob detection, and the coordinates of the target object in the VTK coordinate system may be obtained by performing Blob detection on the background noise image, that is, the second position information within the first specific field of view, which may be expressed as MarksInVTK .
根据本发明的实施例,可以基于目标对象的第一位置信息MarksInCamera和目标对象的第二位置信息MarksInVTK,确定第二特定视野范围的坐标系与第一特定视野范围的坐标系之间的第一转换矩阵,第一转换矩阵可以表示为CameraToVTK。具体地,可以由下式(2)得到,即:According to an embodiment of the present invention, a first transformation matrix between a coordinate system of a second specific field of view and a coordinate system of a first specific field of view can be determined based on the first position information MarksInCamera of the target object and the second position information MarksInVTK of the target object. The first transformation matrix can be expressed as CameraToVTK . Specifically, it can be obtained by the following formula (2), that is:
(2); (2);
其中,vtkClassTrf表示求解第二特定视野范围的坐标系与第一特定视野范围的坐标系之间的第一转换矩阵的函数。Here, vtkClassTrf represents a function for solving a first transformation matrix between a coordinate system of the second specific field of view and a coordinate system of the first specific field of view.
图6示出了根据本发明实施例的确定第一转换矩阵的流程示意图。FIG. 6 shows a schematic diagram of a process of determining a first conversion matrix according to an embodiment of the present invention.
结合图5和图6所示,在配准装置中,可以以配准物601为参考对象,基于配准物上标志物602与深度摄像设备603的相对固定位置,可以先前获得标志物在深度摄像设备的视野范围内的坐标,也就是第二特定视野范围内的第一位置信息604,可以表示为Marks InCamera。可将配准装置固定在磁共振成像设备扫描腔体内,利用收发一体线圈,设置特定的梯度扫描方向,对标志物602进行FOV扫描,获取标志物为亮点的背景噪声图像,也即与标志物对应的扫描图像605,对扫描图像605进行连通域检测处理,得到标志物在VTK坐标系下的坐标,也就是第一特定视野范围内的第二位置信息606,可以表示为MarksInVTK。根据第一位置信息604和第二位置信息606,确定第二特定视野范围的坐标系与第一特定视野范围的坐标系之间的第一转换矩阵607。As shown in combination with FIG. 5 and FIG. 6 , in the registration device, the registration object 601 can be used as a reference object. Based on the relative fixed position of the marker 602 on the registration object and the depth camera 603, the coordinates of the marker within the field of view of the depth camera can be obtained in advance, that is, the first position information 604 within the second specific field of view, which can be expressed as Marks InCamera . The registration device can be fixed in the scanning cavity of the magnetic resonance imaging device, and a specific gradient scanning direction is set by using the transceiver integrated coil to perform FOV scanning on the marker 602, and obtain a background noise image with the marker as a bright spot, that is, a scanned image 605 corresponding to the marker. The scanned image 605 is processed by connected domain detection to obtain the coordinates of the marker in the VTK coordinate system, that is, the second position information 606 within the first specific field of view, which can be expressed as Marks InVTK . According to the first position information 604 and the second position information 606, the first conversion matrix 607 between the coordinate system of the second specific field of view and the coordinate system of the first specific field of view is determined.
图7示出了根据本发明实施例的获取第一三维轮廓空间信息的方法流程图。FIG. 7 shows a flow chart of a method for acquiring first three-dimensional contour space information according to an embodiment of the present invention.
如图7所示,该方法包括操作S710~操作S740。As shown in FIG. 7 , the method includes operations S710 to S740 .
在操作S710,获取扫描对象在第二特定视野范围内的三维轮廓空间参考信息。In operation S710, three-dimensional contour space reference information of a scan object within a second specific field of view is acquired.
在操作S720,确定扫描对象在第一特定视野范围内扫描切片的第一定位信息。In operation S720, first positioning information of a scan slice of the scan object within a first specific field of view is determined.
在操作S730,根据第一定位信息和三维轮廓空间参考信息,确定扫描对象在第二特定视野范围内扫描切片的第二三维轮廓空间信息。In operation S730, second three-dimensional contour space information of a scanned slice of the scanned object within a second specific field of view is determined according to the first positioning information and the three-dimensional contour space reference information.
在操作S740,对第二三维轮廓空间信息进行转换处理,得到扫描对象在第一特定视野范围内扫描切片的第一三维轮廓空间信息。In operation S740, the second three-dimensional contour space information is converted to obtain first three-dimensional contour space information of a scanned slice of the scanned object within a first specific field of view.
根据本发明的实施例,第二特定视野范围可以为深度摄像设备下可扫描区域。三维轮廓空间参考信息可以为扫描对象的三维表面轮廓空间信息。According to an embodiment of the present invention, the second specific field of view may be a scannable area under the depth camera device. The three-dimensional contour space reference information may be three-dimensional surface contour space information of the scanned object.
根据本发明的实施例,三维轮廓空间参考信息可以通过磁共振成像设备中安装的配准的深度摄像设备捕获扫描对象获取的。According to an embodiment of the present invention, the three-dimensional contour spatial reference information can be obtained by capturing the scanned object through a registered depth camera installed in a magnetic resonance imaging device.
根据本发明的实施例,第一特定视野范围可以为在VTK坐标下选择的可扫描区域。可以根据扫描对象检测部位确定在VTK坐标下扫描切片的扫描区域,也即,第一定位信息。According to an embodiment of the present invention, the first specific field of view may be a scannable area selected under VTK coordinates. The scanning area for scanning slices under VTK coordinates may be determined according to the detected part of the scanned object, that is, the first positioning information.
根据本发明的实施例,可以根据扫描切片在VTK坐标下的第一定位信息,从获取扫描对象在深度摄像设备坐标系下扫描切片第二三维轮廓空间信息,可以表示为SpatialContourInCamera。第二三维轮廓空间信息表征该扫描切片在深度摄像设备内的截面空间信息。According to an embodiment of the present invention, the second three-dimensional contour space information of the scanned slice in the depth camera device coordinate system can be obtained based on the first positioning information of the scanned slice in the VTK coordinate system, which can be expressed as SpatialContourInCamera . The second three-dimensional contour space information represents the cross-sectional space information of the scanned slice in the depth camera device.
根据本发明的实施例,可以基于第一转换矩阵,对第二三维轮廓空间信息进行转换处理,得到扫描对象在第一特定视野范围内扫描切片的第一三维轮廓信息。第一三轮轮廓信息可以表示为SpatialContour。具体地,可以由式(3)得到,即:According to an embodiment of the present invention, the second three-dimensional contour space information can be transformed based on the first transformation matrix to obtain the first three-dimensional contour information of the scanned slice of the scanned object within the first specific field of view. The first three-round contour information can be expressed as SpatialContour . Specifically, it can be obtained by formula (3), that is:
(3); (3);
其中,CameraToVTK为第一转换矩阵。Among them, CameraToVTK is the first transformation matrix.
根据本发明的实施例,根据第一定位信息和三维轮廓空间参考信息,确定扫描对象在第二特定视野范围内扫描切片的第二三维轮廓空间信息,包括:根据第一定位信息和第一转换矩阵,从三维轮廓空间参考信息中提取与第一定位信息对应的第二定位信息,其中,第一转换矩阵是通过第二特定视野范围的坐标系与第一特定视野范围的坐标系确定的,第二定位信息为扫描对象在第二特定视野范围内扫描切片的信息;根据第二定位信息,确定扫描对象在第二特定视野范围内扫描切片的第二三维轮廓空间信息。According to an embodiment of the present invention, second three-dimensional contour space information of scanned slices of a scanned object within a second specific field of view is determined based on first positioning information and three-dimensional contour space reference information, including: extracting second positioning information corresponding to the first positioning information from the three-dimensional contour space reference information based on the first positioning information and a first transformation matrix, wherein the first transformation matrix is determined by a coordinate system of the second specific field of view and a coordinate system of the first specific field of view, and the second positioning information is information of scanned slices of the scanned object within the second specific field of view; and determining second three-dimensional contour space information of scanned slices of the scanned object within the second specific field of view based on the second positioning information.
根据本发明的实施例,扫描切片可以为具有规则形状的切片,例如,可以为长方形。以单一扫描切片为示例,扫描切片的第一定位信息可以包括该扫描切片在第一特定视野范围内的定位信息,定位信息可以包括扫描中心点(ScanCenter)、扫描时的频率编码方向(FreqDir)、扫描时的相位编码方向(PhaseDir)和扫描长度(FOVLength),其中,扫描长度(FOVLength)可以包括频率编码方向的扫描长度(FreqLength)和相位编码方向的扫描长度(PhaseLength)。According to an embodiment of the present invention, the scan slice may be a slice having a regular shape, for example, a rectangle. Taking a single scan slice as an example, the first positioning information of the scan slice may include positioning information of the scan slice within a first specific field of view, and the positioning information may include a scan center point ( ScanCenter ), a frequency encoding direction ( FreqDir ) during scanning, a phase encoding direction ( PhaseDir ) during scanning, and a scan length ( FOVLength ), wherein the scan length ( FOVLength ) may include a scan length ( FreqLength ) in a frequency encoding direction and a scan length ( PhaseLength ) in a phase encoding direction.
根据本发明的实施例,可以对第一转换矩阵进行逆处理,得到第一转换矩阵的逆矩阵,可以表示为CameraToVTK -1。According to an embodiment of the present invention, the first conversion matrix may be inversely processed to obtain an inverse matrix of the first conversion matrix, which may be expressed as CameraToVTK −1 .
根据本发明的实施例,第二定位信息可以为获取扫描对象扫描时的扫描切片在第二特定视野范围内的截面的坐标信息。According to an embodiment of the present invention, the second positioning information may be coordinate information of a cross section of a scan slice within a second specific field of view when the scan object is scanned.
根据本发明的实施例,可以利用第一转换矩阵的逆矩阵,确定与第一定位信息对应的在第二特定视野范围内的扫描切片的截面的第二定位信息,其中,第一定位信息可以为扫描对象扫描时扫描切片在VTK坐标系下的信息。具体地,可以通过下式(4)~(7)得到,即:According to an embodiment of the present invention, the inverse matrix of the first transformation matrix can be used to determine the second positioning information of the cross section of the scan slice within the second specific field of view corresponding to the first positioning information, wherein the first positioning information can be the information of the scan slice in the VTK coordinate system when the scan object is scanned. Specifically, it can be obtained by the following formulas (4) to (7), that is:
(4); (4);
(5); (5);
(6); (6);
(7); (7);
其中,ScanCenterInCamera为转换到第二特定视野范围内的扫描中心点;FreqDirInCamera为转换到第二特定视野范围内的频率编码方向;PhaseDirInCamera为转换到第二特定视野范围内的相位编码方向;FOVLength为在第二特定视野范围内的扫描长度。Among them, ScanCenterInCamera is the scanning center point converted to the second specific field of view; FreqDirInCamera is the frequency encoding direction converted to the second specific field of view; PhaseDirInCamera is the phase encoding direction converted to the second specific field of view; FOVLength is the scanning length in the second specific field of view.
根据本发明的实施例,第二定位信息可以包括在第二特定视野范围内的扫描切片截面的四个方位的定位信息,例如,四个方位的定位信息可以包括左上角定位信息、左下角定位信息、右上角定位信息、右下角定位信息。According to an embodiment of the present invention, the second positioning information may include positioning information of four positions of the scanned slice section within the second specific field of view. For example, the positioning information of the four positions may include upper left corner positioning information, lower left corner positioning information, upper right corner positioning information, and lower right corner positioning information.
根据本发明的实施例,可以以单一扫描切片为示例,例如,确定扫描切片对应的在第二特定视野范围内的切面截面的左上角的定位信息,可以由下式(8)~(10)得到,即:According to an embodiment of the present invention, a single scan slice can be taken as an example. For example, the positioning information of the upper left corner of the cut section corresponding to the scan slice within the second specific field of view can be obtained by the following formulas (8) to (10), that is:
(8); (8);
(9); (9);
(10); (10);
其中,TopLeftCorner[0]为左上角的X轴方向坐标信息;TopLeftCorner[1]为左上角的Y轴方向坐标信息;TopLeftCorner[2]为左上角的Z轴方向坐标信息。Among them, TopLeftCorner [0] is the X-axis coordinate information of the upper left corner; TopLeftCorner [1] is the Y-axis coordinate information of the upper left corner; TopLeftCorner [2] is the Z-axis coordinate information of the upper left corner.
根据本发明的实施例,例如,确定扫描切片对应的在第二特定视野范围内的切面截面的左下角的定位信息,可以由下式(11)~(13)得到,即:According to an embodiment of the present invention, for example, the positioning information of the lower left corner of the cut section corresponding to the scan slice within the second specific field of view can be obtained by the following formulas (11) to (13), that is:
(11); (11);
(12); (12);
(13); (13);
其中,BottomLeftCorner[0]为左下角的X轴方向坐标信息;BottomLeftCorner[1]为左下角的Y轴方向坐标信息;BottomLeftCorner[2]为左下角的Z轴方向坐标信息。Among them, BottomLeftCorner [0] is the X-axis coordinate information of the lower left corner; BottomLeftCorner [1] is the Y-axis coordinate information of the lower left corner; BottomLeftCorner [2] is the Z-axis coordinate information of the lower left corner.
根据本发明的实施例,例如,确定扫描切片对应的在第二特定视野范围内的切面截面的右上角的定位信息,可以由下式(14)~(16)得到,即:According to an embodiment of the present invention, for example, the positioning information of the upper right corner of the cut section corresponding to the scan slice within the second specific field of view can be obtained by the following equations (14) to (16), that is:
(14); (14);
(15); (15);
(16); (16);
其中,TopRightCorner[0]为右上角的X轴方向坐标信息;TopRightCorner[1]为右上角的Y轴方向坐标信息;TopRightCorner[2]为右上角的Z轴方向坐标信息。Among them, TopRightCorner [0] is the X-axis coordinate information of the upper right corner; TopRightCorner [1] is the Y-axis coordinate information of the upper right corner; TopRightCorner [2] is the Z-axis coordinate information of the upper right corner.
根据本发明的实施例,例如,确定扫描切片对应的在第二特定视野范围内的切面截面的右下角的定位信息,可以由下式(17)~(19)得到,即:According to an embodiment of the present invention, for example, the positioning information of the lower right corner of the cut section corresponding to the scan slice within the second specific field of view can be obtained by the following equations (17) to (19), that is:
(17); (17);
(18); (18);
(19); (19);
其中,BottomRightCorner[0]为右下角的X轴方向坐标信息;BottomRightCorner[1]为右下角的Y轴方向坐标信息;BottomRightCorner[2]为右下角的Z轴方向坐标信息。Among them, BottomRightCorner [0] is the X-axis coordinate information of the lower right corner; BottomRightCorner [1] is the Y-axis coordinate information of the lower right corner; BottomRightCorner [2] is the Z-axis coordinate information of the lower right corner.
根据本发明的实施例,可以根据上述得到的扫描切片的截面的四个方位的定位信息,从三维轮廓空间参考信息中提取该四个方位的定位信息,得到第二定位信息。According to an embodiment of the present invention, the positioning information of the four orientations of the cross section of the scanned slice obtained above can be extracted from the three-dimensional contour space reference information to obtain the second positioning information.
根据本发明的实施例,可以通过第二定位信息获取扫描对象在第二特定视野范围内的扫描切片的第二三维轮廓空间信息,可以表示为SpatialContourInCamera。According to an embodiment of the present invention, second three-dimensional contour spatial information of a scan slice of the scanned object within a second specific field of view may be acquired through the second positioning information, which may be represented as SpatialContourInCamera .
根据本发明的实施例,对第一三维轮廓空间信息进行处理,得到处理后的第一三维轮廓空间信息,包括:确定扫描切片的第一三维轮廓空间信息与目标扫描范围的空间信息之间的第二转化矩阵;根据第二转化矩阵和第一三维轮廓空间信息,得到与第一三维轮廓空间信息对应的第三三维轮廓空间信息,其中,第三三维轮廓空间信息表征扫描对象在第一特定视野范围内扫描切片的三维轮廓空间信息的唯一性;对第三三维轮廓空间信息进行数据分析处理,得到处理后的第一三维轮廓空间信息。According to an embodiment of the present invention, the first three-dimensional contour space information is processed to obtain the processed first three-dimensional contour space information, including: determining a second transformation matrix between the first three-dimensional contour space information of the scanned slice and the spatial information of the target scanning range; obtaining third three-dimensional contour space information corresponding to the first three-dimensional contour space information based on the second transformation matrix and the first three-dimensional contour space information, wherein the third three-dimensional contour space information represents the uniqueness of the three-dimensional contour space information of the scanned slice of the scanned object within the first specific field of view; performing data analysis processing on the third three-dimensional contour space information to obtain the processed first three-dimensional contour space information.
根据本发明的实施例,目标扫描范围可以为能够将任何扫描对象的第一三维轮廓空间信息进行唯一表示的最大扫描范围。According to an embodiment of the present invention, the target scanning range may be a maximum scanning range that can uniquely represent the first three-dimensional contour space information of any scanned object.
根据本发明的实施例,第二转化矩阵可以为将第一三维轮廓空间信息转变到目标扫描范围的坐标系下的矩阵,可以表示为TrfToLargestRegion。According to an embodiment of the present invention, the second conversion matrix may be a matrix for converting the first three-dimensional contour space information into a coordinate system of the target scanning range, and may be expressed as TrfToLargestRegion .
根据本发明的实施例,可以基于第二转化矩阵,通过平移的方式将第一三维轮廓空间信息转变到最大扫描范围坐标系下,得到与第一三维轮廓空间信息对应的第三三维轮廓空间信息,可以表示为UniqueSpatialContour,第三三维轮廓空间信息可以表征扫描对象在第一特定视野范围内扫描切片的三维轮廓空间信息的唯一性。第三三维轮廓空间信息可以由下式(20)得到,即:According to an embodiment of the present invention, based on the second transformation matrix, the first three-dimensional contour space information can be transformed into the maximum scanning range coordinate system by translation, and the third three-dimensional contour space information corresponding to the first three-dimensional contour space information can be obtained, which can be expressed as UniqueSpatialContour . The third three-dimensional contour space information can characterize the uniqueness of the three-dimensional contour space information of the scanned slice of the scanned object within the first specific field of view. The third three-dimensional contour space information can be obtained by the following formula (20), that is:
(20); (20);
根据本发明的实施例,可以对第三三维轮廓空间信息进行网格化和逻辑化,得到适合于已构建的线圈灵敏度估计模型的输入代替信号。According to the embodiment of the present invention, the third three-dimensional contour space information may be gridded and logicized to obtain an input substitute signal suitable for the constructed coil sensitivity estimation model.
根据本发明的实施例,例如,可以以网格化目标扫描区域(如),并将第三三维轮廓空间信息在VTK坐标系下和与VTK坐标系相同的对比坐标下进行空间信息对比,将不为null的置为1,其余置为0,即可得到经网格化和逻辑化处理后的第一三维轮廓空间信息,可以表示为S。具体地,可以由下式(21)得到,即:According to an embodiment of the present invention, for example, Grid target scanning area (such as ), and compare the third 3D contour space information in the VTK coordinate system with the comparison coordinate system that is the same as the VTK coordinate system, set those that are not null to 1, and the rest to 0, and the first 3D contour space information after gridding and logic processing can be obtained, which can be expressed as S. Specifically, it can be obtained by the following formula (21), that is:
(21); (twenty one);
其中,Grid为网格化;Logical为逻辑化。Among them, Grid means gridding; Logical means logic.
根据本发明的实施例,将处理后的第一三维轮廓空间信息输入至已构建的线圈灵敏度估计模型,得到与扫描对象在特定视野范围内扫描切片对应的多通道线圈灵敏度信息,包括:根据第一特征矩阵信息、第二特征矩阵信息和第一三维轮廓空间信息,得到与扫描对象在特定视野范围内扫描切片对应的多通道线圈灵敏度信息,其中,第一特征矩阵是通过历史多通道灵敏度信息确定,第二特征矩阵是通过历史第一三维轮廓空间信息确定。According to an embodiment of the present invention, the processed first three-dimensional contour space information is input into a constructed coil sensitivity estimation model to obtain multi-channel coil sensitivity information corresponding to a scan slice of the scanned object within a specific field of view, including: obtaining multi-channel coil sensitivity information corresponding to a scan slice of the scanned object within a specific field of view based on first characteristic matrix information, second characteristic matrix information and first three-dimensional contour space information, wherein the first characteristic matrix is determined by historical multi-channel sensitivity information, and the second characteristic matrix is determined by historical first three-dimensional contour space information.
根据本发明的实施例,第一特征矩阵信息可以为基于历史多通道线圈灵敏度,在构建线圈灵敏度估计模型过程中确定得到的,可以表示为E C ;第二特征矩阵信息可以为基于历史第一三维轮廓空间信息,在构建线圈灵敏度估计模型过程中确定得到的,可以表示为E S 。According to an embodiment of the present invention, the first characteristic matrix information may be determined in the process of building a coil sensitivity estimation model based on historical multi-channel coil sensitivities, and may be expressed as EC ; the second characteristic matrix information may be determined in the process of building a coil sensitivity estimation model based on historical first three-dimensional contour space information, and may be expressed as E S.
根据本发明的实施例,可以利用第一特征矩阵信息、第二特征矩阵信息和第一三维轮廓空间信息,得到与扫描对象在特定视野范围内扫描切片对应的多通道线圈灵敏度信息。具体地,可以由下式(22)得到,即:According to an embodiment of the present invention, the first characteristic matrix information, the second characteristic matrix information and the first three-dimensional contour space information can be used to obtain the multi-channel coil sensitivity information corresponding to the scan slice of the scanned object within a specific field of view. Specifically, it can be obtained by the following formula (22), that is:
(22); (twenty two);
其中,S为第一三维轮廓空间信息,C为多通道线圈灵敏度信息,E C 为第一特征矩阵信息,即,为与扫描切片对应的多通道线圈灵敏度信息的特征矩阵信息;E S 为第二特征矩阵信息,即,与扫描切片对应的第一三维轮廓空间信息的特征矩阵信息。Among them, S is the first three-dimensional contour space information, C is the multi-channel coil sensitivity information , EC is the first characteristic matrix information, that is, the characteristic matrix information of the multi-channel coil sensitivity information corresponding to the scan slice; ES is the second characteristic matrix information, that is, the characteristic matrix information of the first three-dimensional contour space information corresponding to the scan slice.
根据本发明的实施例,已构建的线圈灵敏度估计模型可以通过如下方法进行构建,从而确定第一特征矩阵信息和第二特征矩阵信息。According to an embodiment of the present invention, the constructed coil sensitivity estimation model may be constructed by the following method, thereby determining the first characteristic matrix information and the second characteristic matrix information.
根据本发明的实施例,获取扫描对象在第一特定视野范围内的扫描切片的历史多通道线圈灵敏度信息。According to an embodiment of the present invention, historical multi-channel coil sensitivity information of a scan slice of a scanned object within a first specific field of view is acquired.
根据本发明的实施例,获取扫描对象在第一特定视野范围内的扫描切片的历史第一三维轮廓空间信息。According to an embodiment of the present invention, historical first three-dimensional contour space information of a scan slice of a scanned object within a first specific field of view is acquired.
根据本发明的实施例,获取扫描对象在第一特定视野范围内的扫描切片的历史第一三维轮廓空间信息与上述获取扫描对象在第一特定视野范围内扫描切片的历史第一三维轮廓空间信息的方法相同,本发明在此不再赘述。According to an embodiment of the present invention, obtaining the historical first three-dimensional contour space information of the scanned slices of the scanned object within the first specific field of view is the same as the above-mentioned method of obtaining the historical first three-dimensional contour space information of the scanned slices of the scanned object within the first specific field of view, and the present invention will not be repeated here.
根据本发明的实施例,获取扫描对象在第一特定视野范围内的扫描切片的历史第一三维轮廓空间信息,包括:获取扫描对象的扫描切片的历史平面轮廓图像,对历史平面轮廓图像进行形态学闭操作处理,得到扫描切片的历史平面轮廓信息,确定从历史平面轮廓信息转换到VTK坐标下的历史第一三维轮廓空间信息的第三转换矩阵,根据历史第一三维轮廓空间信息的第三转换矩阵和历史平面轮廓信息,得到扫描对象在第一特定视野范围内的扫描切片的历史第一三维轮廓空间信息。According to an embodiment of the present invention, the historical first three-dimensional contour space information of the scanned slices of the scanned object within the first specific field of view is obtained, including: obtaining the historical plane contour image of the scanned slices of the scanned object, performing morphological closing operation on the historical plane contour image to obtain the historical plane contour information of the scanned slices, determining a third conversion matrix from the historical plane contour information to the historical first three-dimensional contour space information under VTK coordinates, and obtaining the historical first three-dimensional contour space information of the scanned slices of the scanned object within the first specific field of view according to the third conversion matrix of the historical first three-dimensional contour space information and the historical plane contour information.
根据本发明的实施例,第三转换矩阵可以根据历史平面轮廓图像中图像标准格式信息中的历史扫描原点(Origin)、历史扫描距离(Spacing)和历史扫描方向(Orientation)。可以利用式(23)得到,即:According to an embodiment of the present invention, the third transformation matrix can be obtained according to the historical scanning origin ( Origin ), historical scanning distance ( Spacing ) and historical scanning direction ( Orientation ) in the image standard format information in the historical plane contour image using formula (23), that is:
(23); (twenty three);
其中,IndexToVTK为从历史平面轮廓信息转换到VTK坐标下的历史第一三维轮廓空间信息的第三转换矩阵。Among them, IndexToVTK is a third conversion matrix for converting the historical plane contour information into the historical first three-dimensional contour space information under VTK coordinates.
根据本发明的实施例,扫描对象在第一特定视野范围内的扫描切片的历史第一三维轮廓空间信息,可以利用下式(24)得到,即:According to an embodiment of the present invention, the historical first three-dimensional contour space information of the scanned slice of the scanned object within the first specific field of view can be obtained using the following formula (24), that is:
(24); (twenty four);
其中,SpatialContour¢ 为历史第一三维轮廓空间信息;Contour为历史平面轮廓信息。Among them, SpatialContour¢ is the historical first three-dimensional contour spatial information; Contour is the historical plane contour information.
根据本发明的实施例,可以对历史第一三维轮廓空间信息转换至目标扫描范围坐标系下,得到历史第三三维轮廓空间信息,并对历史第三三维轮廓空间信息进行网格化和逻辑化处理,得到处理后的历史第一三维轮廓空间信息。According to an embodiment of the present invention, the historical first three-dimensional contour space information can be converted to the target scanning range coordinate system to obtain the historical third three-dimensional contour space information, and the historical third three-dimensional contour space information can be gridded and logically processed to obtain the processed historical first three-dimensional contour space information.
根据本发明的实施例,以第j个扫描切片得到的历史多通道线圈灵敏度信息为示例,历史多通道线圈灵敏度信息可由下式(25)表示,即:According to an embodiment of the present invention, the historical multi-channel coil sensitivity information obtained from the j -th scan slice is For example, historical multi-channel coil sensitivity information It can be expressed by the following formula (25):
(25); (25);
其中,j为第j个扫描切片,L为第j个扫描切片对应的多通道总数;M为每个通道线圈上的体素总数;l为第j个扫描切片对应的第l个通道;m为第l个通道上的第m个体素的数值;为第j个扫描切片对应的多通道线圈灵敏度信息中第l个通道线圈灵敏度信息中的第m个体素的数值;1≤j≤J,J表示扫描切片总数,为大于等于1的正整数;1≤l≤L,1≤m≤M,L、M均为大于等于1的正整数。Wherein, j is the jth scan slice, L is the total number of multi-channels corresponding to the jth scan slice; M is the total number of voxels on each channel coil; l is the lth channel corresponding to the jth scan slice; m is the value of the mth voxel on the lth channel; is the value of the mth voxel in the lth channel coil sensitivity information in the multi-channel coil sensitivity information corresponding to the jth scan slice; 1≤ j ≤ J , J represents the total number of scan slices, which is a positive integer greater than or equal to 1 ; 1≤ l ≤ L , 1≤ m ≤ M , L and M are both positive integers greater than or equal to 1.
根据本发明的实施例,以第j个扫描切片的历史第一三维轮廓空间信息为示例,历史第一三维轮廓空间信息可由下式(26)表示,即:According to an embodiment of the present invention, the historical first three-dimensional contour space information of the j -th scanning slice is used. For example, the first three-dimensional contour space information It can be expressed by the following formula (26):
(26); (26);
其中,N为第j个扫描切片对应的历史第一三维轮廓空间信息中的体素总数;为第j个扫描切片对应的的历史第一三维轮廓空间信息中第n个体素的数值;1≤n≤N;N为大于等于1的正整数。Wherein, N is the total number of voxels in the historical first three-dimensional contour space information corresponding to the j -th scan slice; is the value of the nth voxel in the historical first three-dimensional contour space information corresponding to the jth scan slice; 1≤ n ≤ N; N is a positive integer greater than or equal to 1.
根据本发明的实施例,可以通过上式(25)~(26)构建如式(27)的与第j个扫描切片对应的矩阵,即:According to an embodiment of the present invention, the matrix corresponding to the j -th scan slice as shown in formula (27) can be constructed by using the above formulas (25) to (26): ,Right now:
(27); (27);
根据本发明的实施例,以扫描切片总数为J为示例,J为大于等于1的正整数。根据上述式(27)可以得到构建该线圈灵敏度估计模型的矩阵P,如下式(28):According to an embodiment of the present invention, the total number of scan slices is J as an example, where J is a positive integer greater than or equal to 1. According to the above formula (27), the matrix P for constructing the coil sensitivity estimation model can be obtained as shown in the following formula (28):
(28); (28);
其中,;;为对与第j个扫描切片对应的矩阵进行去中心化;为对与J个扫描切片对应的矩阵的均值。in, ; ; To decentralize the matrix corresponding to the j -th scan slice; is the mean of the matrices corresponding to J scan slices.
根据本发明的实施例,可以结合上式(25)~(28)可知,构建该线圈灵敏度估计模型的矩阵P的大小可以为。According to the embodiment of the present invention, it can be known from the above equations (25) to (28) that the size of the matrix P for constructing the coil sensitivity estimation model can be .
根据本发明的实施例,可以根据矩阵P构建协方差矩阵。协方差矩阵是半正定的,其特征值为非负。由于协方差矩阵的大小为,因此不能直接进行特征值分解;According to an embodiment of the present invention, the covariance matrix can be constructed according to the matrix P. . Covariance matrix is semi-positive definite and its eigenvalues are non-negative. Since the covariance matrix The size is , so the eigenvalue decomposition cannot be performed directly;
根据本发明的实施例,因为,所以至多有J-1 个无关的特征向量。According to an embodiment of the present invention, because ,so There are at most J -1 irrelevant eigenvectors.
根据本发明的实施例,可以假设有特征值为λ的特征向量X,将PX乘以,可以得到,因此PX是的特征向量,可以表示为式(29),即:According to an embodiment of the present invention, it can be assumed that There is an eigenvector X with eigenvalue λ , multiply PX by , we can get , so PX is The characteristic vector of can be expressed as formula (29), that is:
(29); (29);
根据本发明的实施例,特征矩阵E可以表征与扫描切片对应的历史多通道线圈灵敏度信息和历史第一三维轮廓空间信息的特征。According to an embodiment of the present invention, the feature matrix E may characterize features of historical multi-channel coil sensitivity information and historical first three-dimensional contour space information corresponding to the scan slice.
根据本发明的实施例,可以假设与每个扫描切片对应的矩阵,都可以使用的前K个最大特征值对应的特征向量加权近似,可以表示为式(30),即:According to an embodiment of the present invention, it can be assumed that the matrix corresponding to each scan slice is , you can use The weighted approximation of the eigenvectors corresponding to the first K largest eigenvalues can be expressed as formula (30), that is:
(30); (30);
其中,为第k个特征向量对象的加权值,为第k个特征向量,1≤k≤K;K为大于等于1的正整数。in, is the weighted value of the kth eigenvector object, is the kth eigenvector, 1≤ k ≤ K; K is a positive integer greater than or equal to 1.
根据本发明的实施例,再次结合上式(25)~(28),及式(30),可得式(31)~(33):According to an embodiment of the present invention, by combining the above equations (25) to (28) and equation (30), equations (31) to (33) can be obtained:
(31); (31);
(32); (32);
(33); (33);
其中,为矩阵中前行的多通道线圈灵敏度信息;为矩阵中后N行的第一三维轮廓空间信息。in, For the matrix Multi-channel coil sensitivity information in the middle forward row; For the matrix The first three-dimensional contour space information of the last N lines.
根据本发明的实施例,基于上式(31)~(33)可得出如下式(34)~(35),即:According to an embodiment of the present invention, based on the above equations (31) to (33), the following equations (34) to (35) can be obtained, namely:
(34); (34);
(35); (35);
其中,W为特征矩阵E对应的加权值,可表示为;E C 为与扫描切片对应的多通道线圈灵敏度信息的特征矩阵信息;E S 为与扫描切片对应的第一三维轮廓空间信息的特征矩阵信息。Among them, W is the weighted value corresponding to the feature matrix E , which can be expressed as ; EC is the characteristic matrix information of the multi-channel coil sensitivity information corresponding to the scan slice; ES is the characteristic matrix information of the first three-dimensional contour space information corresponding to the scan slice .
根据本发明的实施例,E C 和E S 可以分别来自特征矩阵E的前行和后N行,利用上式(34)~(35)可得构建的线圈灵敏度估计模型,即,如下式(36):According to an embodiment of the present invention, EC and ES can be respectively derived from the front of the feature matrix E The coil sensitivity estimation model constructed by using the above equations (34) to (35) can be obtained, that is, as shown in the following equation (36):
(36); (36);
即,上述式(22):(22);That is, the above formula (22): (twenty two);
其中,E C 为的矩阵;E S 的逆为的矩阵;S为的矩阵。Among them, E C is The inverse of E S is The matrix of The matrix of .
根据本发明的实施例,可以利用该模型,基于扫描对象在第一特定视野范围内扫描切片的第一三维轮廓空间信息,可以快速得到与扫描切片对应的多通道线圈灵敏度信息,从而快速实现图像重建。According to an embodiment of the present invention, the model can be used to quickly obtain multi-channel coil sensitivity information corresponding to the scanned slice based on the first three-dimensional contour space information of the scanned slice of the scanned object within a first specific field of view, thereby quickly realizing image reconstruction.
图8示出了根据本发明的实施例的磁共振图像重建装置的框图。FIG8 shows a block diagram of a magnetic resonance image reconstruction apparatus according to an embodiment of the present invention.
如图8所示,该装置可以包括:获取模块810、处理模块820、输入模块830和重建模块840。As shown in FIG. 8 , the apparatus may include: an acquisition module 810 , a processing module 820 , an input module 830 and a reconstruction module 840 .
获取模块810,用于获取扫描对象在第一特定视野范围内扫描切片的第一三维轮廓空间信息和磁共振扫描数据。The acquisition module 810 is used to acquire first three-dimensional contour space information and magnetic resonance scanning data of a scanned slice of a scanned object within a first specific field of view.
处理模块820,用于对第一三维轮廓空间信息进行处理,得到处理后的第一三维轮廓空间信息。The processing module 820 is used to process the first three-dimensional contour space information to obtain processed first three-dimensional contour space information.
输入模块830,用于将处理后的第一三维轮廓空间信息输入至已构建的线圈灵敏度估计模型,得到与扫描对象在第一特定视野范围内扫描切片对应的多通道线圈灵敏度信息。The input module 830 is used to input the processed first three-dimensional contour space information into the constructed coil sensitivity estimation model to obtain the multi-channel coil sensitivity information corresponding to the scan slice of the scanned object within the first specific field of view.
重建模块840,用于根据多通道线圈灵敏度信息和磁共振扫描数据,对扫描对象进行磁共振图像重建。The reconstruction module 840 is used to reconstruct the magnetic resonance image of the scanned object according to the multi-channel coil sensitivity information and the magnetic resonance scanning data.
根据本发明的实施例,该装置还包括:确定模块。According to an embodiment of the present invention, the device further includes: a determination module.
确定模块,用于确定第二特定视野范围的坐标系与第一特定视野范围的坐标系之间的第一转换矩阵,其中,第一转换矩阵应用于获取扫描对象在第一特定视野范围内扫描切片的第一三维轮廓空间信息。A determination module is used to determine a first transformation matrix between a coordinate system of a second specific field of view and a coordinate system of a first specific field of view, wherein the first transformation matrix is used to obtain first three-dimensional contour space information of a scanned slice of a scanned object within the first specific field of view.
根据本发明的实施例,确定模块可以包括:第一获取子模块、扫描子模块、第一处理子模块和第一确定子模块。According to an embodiment of the present invention, the determination module may include: a first acquisition submodule, a scanning submodule, a first processing submodule and a first determination submodule.
第一获取子模块,用于获取配准装置中目标对象在第二特定视野范围内的第一位置信息,目标对象为可被磁共振扫描识别的标志物。The first acquisition submodule is used to acquire first position information of a target object in the registration device within a second specific field of view, where the target object is a marker that can be identified by magnetic resonance scanning.
扫描子模块,用于利用预先设置的梯度扫描方向,对配准装置中目标对象在第一特定视野范围内进行扫描,获取与目标对象对应的扫描图像。The scanning submodule is used to scan the target object in the registration device within a first specific field of view using a preset gradient scanning direction to obtain a scanned image corresponding to the target object.
第一处理子模块,用于对扫描图像进行连通域检测处理,得到与扫描图像对应的第二位置信息。The first processing submodule is used to perform connected domain detection processing on the scanned image to obtain second position information corresponding to the scanned image.
第一确定子模块,用于根据第一位置信息和第二位置信息,确定第二特定视野范围的坐标系与第一特定视野范围的坐标系之间的第一转换矩阵。The first determining submodule is used to determine a first conversion matrix between a coordinate system of the second specific visual field and a coordinate system of the first specific visual field according to the first position information and the second position information.
根据本发明的实施例,获取模块810可以包括:第二获取子模块、第二确定子模块、第三确定子模块和第二处理子模块。According to an embodiment of the present invention, the acquisition module 810 may include: a second acquisition submodule, a second determination submodule, a third determination submodule and a second processing submodule.
第二获取子模块,用于获取扫描对象在第二特定视野范围内的三维轮廓空间参考信息。The second acquisition submodule is used to acquire the three-dimensional contour spatial reference information of the scanned object within a second specific field of view.
第二确定子模块,用于确定扫描对象在第一特定视野范围内扫描切片的第一定位信息。The second determining submodule is used to determine first positioning information of a scan slice of the scanned object within a first specific field of view.
第三确定子模块,用于根据第一定位信息和三维轮廓空间参考信息,确定扫描对象在第二特定视野范围内扫描切片的第二三维轮廓空间信息。The third determination submodule is used to determine second three-dimensional contour space information of a scanned slice of the scanned object within a second specific field of view according to the first positioning information and the three-dimensional contour space reference information.
第二处理子模块,用于对第二三维轮廓空间信息进行转换处理,得到扫描对象在第一特定视野范围内扫描切片的第一三维轮廓空间信息。The second processing submodule is used to convert the second three-dimensional contour space information to obtain the first three-dimensional contour space information of the scanned slice of the scanned object within the first specific field of view.
根据本发明的实施例,第三确定子模块可以包括:提取单元和确定单元。According to an embodiment of the present invention, the third determining submodule may include: an extracting unit and a determining unit.
提取单元,用于根据第一定位信息和第一转换矩阵,从三维轮廓空间参考信息中提取与第一定位信息对应的第二定位信息,其中,第一转换矩阵是通过第二特定视野范围的坐标系与第一特定视野范围的坐标系确定的,第二定位信息为扫描对象在第二特定视野范围内扫描切片的信息。An extraction unit is used to extract second positioning information corresponding to the first positioning information from the three-dimensional contour space reference information according to the first positioning information and the first transformation matrix, wherein the first transformation matrix is determined by the coordinate system of the second specific field of view and the coordinate system of the first specific field of view, and the second positioning information is information of the scanned slices of the scanned object within the second specific field of view.
确定单元,用于根据第二定位信息,确定扫描对象在第二特定视野范围内扫描切片的第二三维轮廓空间信息。The determination unit is used to determine the second three-dimensional contour space information of the scanned slice of the scanned object within the second specific field of view according to the second positioning information.
根据本发明的实施例,处理模块820可以包括:第四确定子模块、第一获得子模块和第三处理子模块。According to an embodiment of the present invention, the processing module 820 may include: a fourth determining submodule, a first obtaining submodule and a third processing submodule.
第四确定子模块,用于确定扫描切片的第一三维轮廓空间信息与目标扫描范围的空间信息之间的第二转化矩阵。The fourth determination submodule is used to determine a second conversion matrix between the first three-dimensional contour space information of the scan slice and the space information of the target scanning range.
第一获得子模块,用于根据第二转化矩阵和第一三维轮廓空间信息,得到与第一三维轮廓空间信息对应的第三三维轮廓空间信息,其中,第三三维轮廓空间信息表征扫描对象在第一特定视野范围内扫描切片的三维轮廓空间信息的唯一性。The first acquisition submodule is used to obtain third three-dimensional contour space information corresponding to the first three-dimensional contour space information according to the second transformation matrix and the first three-dimensional contour space information, wherein the third three-dimensional contour space information represents the uniqueness of the three-dimensional contour space information of the scanned slices within the first specific field of view of the scanned object.
第三处理子模块,用于对第三三维轮廓空间信息进行数据分析处理,得到处理后的第一三维轮廓空间信息。The third processing submodule is used to perform data analysis and processing on the third three-dimensional contour space information to obtain processed first three-dimensional contour space information.
根据本发明的实施例,输入模块830可以包括:第二获得子模块。According to an embodiment of the present invention, the input module 830 may include: a second obtaining submodule.
第二获得子模块,用于根据第一特征矩阵信息、第二特征矩阵信息和第一三维轮廓空间信息,得到与扫描对象在特定视野范围内扫描切片对应的多通道线圈灵敏度信息,其中,第一特征矩阵是通过历史多通道灵敏度信息确定,第二特征矩阵是通过历史第一三维轮廓空间信息确定。The second acquisition submodule is used to obtain multi-channel coil sensitivity information corresponding to the scanned slice of the scanned object within a specific field of view based on the first characteristic matrix information, the second characteristic matrix information and the first three-dimensional contour space information, wherein the first characteristic matrix is determined by historical multi-channel sensitivity information, and the second characteristic matrix is determined by historical first three-dimensional contour space information.
根据本发明的实施例的模块、子模块、单元中的任意多个、或其中任意多个的至少部分功能可以在一个模块中实现。根据本发明实施例的模块、子模块、单元、中的任意一个或多个可以被拆分成多个模块来实现。根据本发明实施例的模块、子模块、单元、中的任意一个或多个可以至少被部分地实现为硬件电路,例如现场可编程门阵列(FPGA)、可编程逻辑阵列(PLA)、片上系统、基板上的系统、封装上的系统、专用集成电路(ASIC),或可以通过对电路进行集成或封装的任何其他的合理方式的硬件或固件来实现,或以软件、硬件以及固件三种实现方式中任意一种或以其中任意几种的适当组合来实现。或者,根据本发明实施例的模块、子模块、单元、中的一个或多个可以至少被部分地实现为计算机程序模块,当该计算机程序模块被运行时,可以执行相应的功能。According to the embodiments of the present invention, any multiple of the modules, submodules, and units, or at least part of the functions of any multiple of them can be implemented in one module. According to the embodiments of the present invention, any one or more of the modules, submodules, and units can be split into multiple modules for implementation. According to the embodiments of the present invention, any one or more of the modules, submodules, and units can be at least partially implemented as hardware circuits, such as field programmable gate arrays (FPGAs), programmable logic arrays (PLAs), systems on chips, systems on substrates, systems on packages, application specific integrated circuits (ASICs), or can be implemented by hardware or firmware in any other reasonable way of integrating or packaging the circuit, or by any one of the three implementation methods of software, hardware, and firmware, or by a suitable combination of any of them. Alternatively, according to the embodiments of the present invention, one or more of the modules, submodules, and units can be at least partially implemented as computer program modules, and when the computer program modules are run, the corresponding functions can be executed.
例如,获取模块810、处理模块820、输入模块830和重建模块840中的任意多个可以合并在一个模块/单元/中实现,或者其中的任意一个模块/单元/可以被拆分成多个模块/单元/。或者,这些模块/单元/中的一个或多个模块/单元/的至少部分功能可以与其他模块/单元/的至少部分功能相结合,并在一个模块/单元/中实现。根据本发明的实施例,获取模块810、处理模块820、输入模块830和重建模块840中的至少一个可以至少被部分地实现为硬件电路,例如现场可编程门阵列(FPGA)、可编程逻辑阵列(PLA)、片上系统、基板上的系统、封装上的系统、专用集成电路(ASIC),或可以通过对电路进行集成或封装的任何其他的合理方式等硬件或固件来实现,或以软件、硬件以及固件三种实现方式中任意一种或以其中任意几种的适当组合来实现。或者,获取模块810、处理模块820、输入模块830和重建模块840中的至少一个可以至少被部分地实现为计算机程序模块,当该计算机程序模块被运行时,可以执行相应的功能。For example, any multiple of the acquisition module 810, the processing module 820, the input module 830, and the reconstruction module 840 may be combined in one module/unit/for implementation, or any one of the modules/units/may be split into multiple modules/units/. Alternatively, at least part of the functions of one or more of these modules/units/may be combined with at least part of the functions of other modules/units/and implemented in one module/unit/. According to an embodiment of the present invention, at least one of the acquisition module 810, the processing module 820, the input module 830, and the reconstruction module 840 may be at least partially implemented as a hardware circuit, such as a field programmable gate array (FPGA), a programmable logic array (PLA), a system on a chip, a system on a substrate, a system on a package, an application specific integrated circuit (ASIC), or may be implemented by hardware or firmware such as any other reasonable way of integrating or packaging the circuit, or implemented in any one of the three implementation modes of software, hardware, and firmware, or in a suitable combination of any of them. Alternatively, at least one of the acquisition module 810, the processing module 820, the input module 830 and the reconstruction module 840 may be at least partially implemented as a computer program module, and when the computer program module is executed, the corresponding function may be performed.
需要说明的是,本发明的实施例中磁共振图像重建装置部分与本发明的实施例中磁共振图像重建方法部分是相对应的,磁共振图像重建装置部分的描述具体参考磁共振图像重建方法部分,在此不再赘述。It should be noted that the magnetic resonance image reconstruction device part in the embodiment of the present invention corresponds to the magnetic resonance image reconstruction method part in the embodiment of the present invention. The description of the magnetic resonance image reconstruction device part specifically refers to the magnetic resonance image reconstruction method part, which will not be repeated here.
图9示出了根据本发明实施例的适于实现磁共振图像重建方法的电子设备的框图。图9示出的电子设备仅仅是一个示例,不应对本发明实施例的功能和使用范围带来任何限制。Fig. 9 shows a block diagram of an electronic device suitable for implementing a magnetic resonance image reconstruction method according to an embodiment of the present invention. The electronic device shown in Fig. 9 is only an example and should not bring any limitation to the functions and scope of use of the embodiment of the present invention.
如图9所示,根据本发明实施例的电子设备包括处理器901,其可以根据存储在只读存储器(ROM)902中的程序或者从存储部分908加载到随机访问存储器(RAM)903中的程序而执行各种适当的动作和处理。处理器901例如可以包括通用微处理器(例如CPU)、指令集处理器和/或相关芯片组和/或专用微处理器(例如,专用集成电路(ASIC)),等等。处理器901还可以包括用于缓存用途的板载存储器。处理器901可以包括用于执行根据本发明实施例的方法流程的不同动作的单一处理单元或者是多个处理单元。As shown in Figure 9, the electronic device according to an embodiment of the present invention includes a processor 901, which can perform various appropriate actions and processes according to the program stored in the read-only memory (ROM) 902 or the program loaded from the storage part 908 to the random access memory (RAM) 903. The processor 901 may include, for example, a general-purpose microprocessor (such as a CPU), an instruction set processor and/or a related chipset and/or a special-purpose microprocessor (for example, an application-specific integrated circuit (ASIC)), etc. The processor 901 may also include an onboard memory for caching purposes. The processor 901 may include a single processing unit or multiple processing units for performing different actions of the method flow according to an embodiment of the present invention.
在RAM 903中,存储有电子设备操作所需的各种程序和数据。处理器 901、ROM 902以及RAM 903通过总线904彼此相连。处理器901通过执行ROM 902和/或RAM 903中的程序来执行根据本发明实施例的方法流程的各种操作。需要注意,程序也可以存储在除ROM 902和RAM 903以外的一个或多个存储器中。处理器901也可以通过执行存储在一个或多个存储器中的程序来执行根据本发明实施例的方法流程的各种操作。In RAM 903, various programs and data required for the operation of the electronic device are stored. Processor 901, ROM 902, and RAM 903 are connected to each other via bus 904. Processor 901 performs various operations of the method flow according to the embodiment of the present invention by executing the program in ROM 902 and/or RAM 903. It should be noted that the program can also be stored in one or more memories other than ROM 902 and RAM 903. Processor 901 can also perform various operations of the method flow according to the embodiment of the present invention by executing the program stored in one or more memories.
根据本发明的实施例,电子设备还可以包括输入/输出(I/O)接口905,输入/输出(I/O)接口905也连接至总线904。系统还可以包括连接至I/O接口905的以下部件中的一项或多项:包括键盘、鼠标等的输入部分906;包括诸如阴极射线管(CRT)、液晶显示器(LCD)等以及扬声器等的输出部分907;包括硬盘等的存储部分908;以及包括诸如LAN卡、调制解调器等的网络接口卡的通信部分909。通信部分909经由诸如因特网的网络执行通信处理。驱动器910也根据需要连接至I/O接口905。可拆卸介质911,诸如磁盘、光盘、磁光盘、半导体存储器等等,根据需要安装在驱动器910上,以便于从其上读出的计算机程序根据需要被安装入存储部分908。According to an embodiment of the present invention, the electronic device may further include an input/output (I/O) interface 905, which is also connected to the bus 904. The system may further include one or more of the following components connected to the I/O interface 905: an input section 906 including a keyboard, a mouse, etc.; an output section 907 including a cathode ray tube (CRT), a liquid crystal display (LCD), etc., and a speaker, etc.; a storage section 908 including a hard disk, etc.; and a communication section 909 including a network interface card such as a LAN card, a modem, etc. The communication section 909 performs communication processing via a network such as the Internet. A drive 910 is also connected to the I/O interface 905 as needed. A removable medium 911, such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, etc., is installed on the drive 910 as needed, so that a computer program read therefrom is installed into the storage section 908 as needed.
根据本发明的实施例,根据本发明实施例的方法流程可以被实现为计算机软件程序。例如,本发明的实施例包括一种计算机程序产品,其包括承载在计算机可读存储介质上的计算机程序,该计算机程序包含用于执行流程图所示的方法的程序代码。在这样的实施例中,该计算机程序可以通过通信部分909从网络上被下载和安装,和/或从可拆卸介质911被安装。在该计算机程序被处理器901执行时,执行本发明实施例的系统中限定的上述功能。根据本发明的实施例,上文描述的系统、设备、装置、模块、单元等可以通过计算机程序模块来实现。According to an embodiment of the present invention, the method flow according to an embodiment of the present invention can be implemented as a computer software program. For example, an embodiment of the present invention includes a computer program product, which includes a computer program carried on a computer-readable storage medium, and the computer program includes a program code for executing the method shown in the flowchart. In such an embodiment, the computer program can be downloaded and installed from the network through the communication part 909, and/or installed from the removable medium 911. When the computer program is executed by the processor 901, the above-mentioned functions defined in the system of the embodiment of the present invention are executed. According to an embodiment of the present invention, the system, equipment, device, module, unit, etc. described above can be implemented by a computer program module.
本发明还提供了一种计算机可读存储介质,该计算机可读存储介质可以是上述实施例中描述的设备/装置/系统中所包含的;也可以是单独存在,而未装配入该设备/装置/系统中。上述计算机可读存储介质承载有一个或者多个程序,当上述一个或者多个程序被执行时,实现根据本发明实施例的方法。The present invention also provides a computer-readable storage medium, which may be included in the device/apparatus/system described in the above embodiment; or may exist independently without being assembled into the device/apparatus/system. The above computer-readable storage medium carries one or more programs, and when the above one or more programs are executed, the method according to the embodiment of the present invention is implemented.
根据本发明的实施例,计算机可读存储介质可以是非易失性的计算机可读存储介质。例如可以包括但不限于:便携式计算机磁盘、硬盘、随机访问存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本发明中,计算机可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。According to an embodiment of the present invention, the computer-readable storage medium may be a non-volatile computer-readable storage medium. For example, it may include but is not limited to: a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination thereof. In the present invention, the computer-readable storage medium may be any tangible medium containing or storing a program, which may be used by or in combination with an instruction execution system, an apparatus or a device.
例如,根据本发明的实施例,计算机可读存储介质可以包括上文描述的ROM 902和/或RAM 903和/或ROM 902和RAM 903以外的一个或多个存储器。For example, according to an embodiment of the present invention, the computer-readable storage medium may include the ROM 902 and/or the RAM 903 described above and/or one or more memories other than the ROM 902 and the RAM 903 .
附图中的流程图和框图,图示了按照本发明各种实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分,上述模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个接连地表示的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图或流程图中的每个方框、以及框图或流程图中的方框的组合,可以用执行规定的功能或操作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。本领域技术人员可以理解,本发明的各个实施例中记载的特征可以进行多种组合和/或结合,即使这样的组合或结合没有明确记载于本发明中。特别地,在不脱离本发明精神和教导的情况下,本发明的各个实施例中记载的特征可以进行多种组合和/或结合。所有这些组合和/或结合均落入本发明的范围。The flowcharts and block diagrams in the accompanying drawings illustrate the possible implementation architecture, functions and operations of the systems, methods and computer program products according to various embodiments of the present invention. In this regard, each box in the flowchart or block diagram may represent a module, a program segment, or a part of a code, and the above-mentioned module, program segment, or a part of the code contains one or more executable instructions for implementing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the box may also occur in a different order from the order marked in the accompanying drawings. For example, two boxes represented in succession can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram or flowchart, and the combination of boxes in the block diagram or flowchart, can be implemented with a dedicated hardware-based system that performs the specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions. It will be understood by those skilled in the art that the features recorded in the various embodiments of the present invention can be combined and/or combined in various ways, even if such a combination or combination is not explicitly recorded in the present invention. In particular, without departing from the spirit and teachings of the present invention, the features recorded in the various embodiments of the present invention can be combined and/or combined in various ways. All such combinations and/or combinations fall within the scope of the present invention.
以上对本发明的实施例进行了描述。但是,这些实施例仅仅是为了说明的目的,而并非为了限制本发明的范围。尽管在以上分别描述了各实施例,但是这并不意味着各个实施例中的措施不能有利地结合使用。本发明的范围由所附权利要求及其等同物限定。不脱离本发明的范围,本领域技术人员可以做出多种替代和修改,这些替代和修改都应落在本发明的范围之内。The embodiments of the present invention are described above. However, these embodiments are only for the purpose of illustration, and are not intended to limit the scope of the present invention. Although the embodiments are described above, this does not mean that the measures in the various embodiments cannot be used in combination. The scope of the present invention is defined by the appended claims and their equivalents. Without departing from the scope of the present invention, those skilled in the art may make various substitutions and modifications, which should all fall within the scope of the present invention.
Claims (8)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410085787.XA CN117607770B (en) | 2024-01-22 | 2024-01-22 | Magnetic resonance image reconstruction method, device, electronic device and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410085787.XA CN117607770B (en) | 2024-01-22 | 2024-01-22 | Magnetic resonance image reconstruction method, device, electronic device and storage medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117607770A CN117607770A (en) | 2024-02-27 |
CN117607770B true CN117607770B (en) | 2024-04-16 |
Family
ID=89944691
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410085787.XA Active CN117607770B (en) | 2024-01-22 | 2024-01-22 | Magnetic resonance image reconstruction method, device, electronic device and storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117607770B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1910470A (en) * | 2004-01-15 | 2007-02-07 | 皇家飞利浦电子股份有限公司 | Coil sensitivity estimation for parallel imaging |
CN103371819A (en) * | 2012-04-27 | 2013-10-30 | 株式会社东芝 | Magnetic resonance imaging apparatus and image processing apparatus |
WO2020142109A1 (en) * | 2019-01-04 | 2020-07-09 | University Of Cincinnati | A system and method for reconstruction of magnetic resonance images acquired with partial fourier acquisition |
CN114264997A (en) * | 2021-12-14 | 2022-04-01 | 武汉联影生命科学仪器有限公司 | Gradient sensitivity calibration method, device and magnetic resonance equipment |
CN114419176A (en) * | 2021-12-24 | 2022-04-29 | 杭州微影医疗科技有限公司 | A reconstruction system and method for propeller magnetic resonance imaging |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102010043370B4 (en) * | 2010-11-04 | 2014-03-13 | Siemens Aktiengesellschaft | Combination of MR signals with suppression of unwanted signal components |
-
2024
- 2024-01-22 CN CN202410085787.XA patent/CN117607770B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1910470A (en) * | 2004-01-15 | 2007-02-07 | 皇家飞利浦电子股份有限公司 | Coil sensitivity estimation for parallel imaging |
CN103371819A (en) * | 2012-04-27 | 2013-10-30 | 株式会社东芝 | Magnetic resonance imaging apparatus and image processing apparatus |
WO2020142109A1 (en) * | 2019-01-04 | 2020-07-09 | University Of Cincinnati | A system and method for reconstruction of magnetic resonance images acquired with partial fourier acquisition |
CN114264997A (en) * | 2021-12-14 | 2022-04-01 | 武汉联影生命科学仪器有限公司 | Gradient sensitivity calibration method, device and magnetic resonance equipment |
CN114419176A (en) * | 2021-12-24 | 2022-04-29 | 杭州微影医疗科技有限公司 | A reconstruction system and method for propeller magnetic resonance imaging |
Non-Patent Citations (1)
Title |
---|
《基于多通道并行采集的部分可分离函数 高分辨动态磁共振成像》;史彩云 等;《集成技术》;20130930;1-6 * |
Also Published As
Publication number | Publication date |
---|---|
CN117607770A (en) | 2024-02-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20200214619A1 (en) | Automated co-registration of prostate mri data | |
US10274567B2 (en) | Magnetic resonance imaging apparatus and image processing apparatus | |
US10914798B2 (en) | Eigen-vector approach for coil sensitivity maps estimation | |
US10180476B2 (en) | Systems and methods for segmented magnetic resonance fingerprinting dictionary matching | |
EP3304113B1 (en) | Mri with variable density sampling | |
CN102144923A (en) | Method and magnetic resonance system to reduce distortions in diffusion imaging | |
US7782051B2 (en) | Geometry based field prediction method for susceptibility mapping and phase artifact removal | |
CN110335204B (en) | Thermal imaging image enhancement method | |
CN116097305B (en) | Method, device, computer equipment and storage medium for generating point cloud normal vector | |
US12136145B2 (en) | Signal amplitude feature-based method for fast reconstructing a magnetic particle imaging and device | |
CN111157935B (en) | Magnetic resonance imaging method, magnetic resonance imaging device, storage medium and computer equipment | |
CN107290700B (en) | Phase correction method and device and magnetic resonance system | |
US8195004B2 (en) | Device for processing raw images or tensor images | |
Carmichael et al. | Diffusion tensor smoothing through weighted Karcher means | |
Hoge et al. | Identification of translational displacements between N-dimensional data sets using the high-order SVD and phase correlation | |
CN118566993A (en) | Aviation electromagnetic three-dimensional random inversion method based on multi-scale analysis | |
CN117607770B (en) | Magnetic resonance image reconstruction method, device, electronic device and storage medium | |
US9354339B2 (en) | Systems and methods for combining triaxial geophone data for localizing nearby transient seismic sources | |
CN114994670B (en) | Radar imaging method, device, equipment and medium based on composite regularization parameters | |
Pure et al. | A new image fusion method based on integration of wavelet and fast discrete curvelet transform | |
CN114022877B (en) | Three-dimensional SAR image interested target extraction method based on self-adaptive threshold | |
CN118506106B (en) | A method for accurate identification of rock mass structural surfaces based on boundary noise pre-segmentation of point cloud data | |
Liu et al. | Analysis of statistical properties of atmospheric turbulence-induced image dancing based on Hilbert transform and dense optical flow | |
Singh et al. | Smooth Polynomial Approach for Microwave Imaging in Sparse Processing Framework | |
Shilpa et al. | Compressive Sensing Technique on MRI Reconstruction—Methodical Survey |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20250317 Address after: 230088 5089 Wangjiang West Road, Hefei high tech Zone, Anhui Patentee after: INSTITUTE OF ADVANCED TECHNOLOGY University OF SCIENCE AND TECHNOLOGY OF CHINA Country or region after: China Patentee after: University of Science and Technology of China Address before: 230088 5089 Wangjiang West Road, Hefei high tech Zone, Anhui Patentee before: INSTITUTE OF ADVANCED TECHNOLOGY University OF SCIENCE AND TECHNOLOGY OF CHINA Country or region before: China |
|
TR01 | Transfer of patent right |