# US 10,117,597 B2 · Apparatus, methods and articles for four dimensional (4D) flow magnetic resonance imaging using coherency identification for magnetic resonance imaging flow data > Agents: this data is available as JSON from the patentgrep API and MCP server — see https://patentgrep.com/llms.txt - Number: US 10,117,597 B2 - Document ID: US10117597B2 - Status: Granted - Kind: B2 - Application number: 15/112,130 - Filed: 2015-01-16 - Published: 2016-11-24 - Granted: 2018-11-06 - Earliest claimed: 2014-01-17 - Assignees: Arterys Inc. - Applicants: Arterys Inc. - Inventors: Fabien Beckers; Albert Hsiao; John Axerio-Cilies; Torin Arni Taerum; Daniel Marc Raymond Beauchamp - Web page: https://patentgrep.com/patent/US10117597B2 ## CPC - A61B 5/055: involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging - A61B 5/0044: {for the heart} - A61B 5/02007: {Evaluating blood vessel condition, e.g. elasticity, compliance} - A61B 5/02014: {Determining aneurysm} - A61B 5/021: Measuring pressure in heart or blood vessels - A61B 5/0263: {using NMR} - A61B 5/7257: {using Fourier transforms} - G01R 33/56308: {Characterization of motion or flow; Dynamic imaging} - G01R 33/56545: {caused by finite or discrete sampling, e.g. Gibbs ringing, truncation artefacts, phase aliasing artefacts} - G06T 7/12: Edge-based segmentation - G06T 7/168: involving transform domain methods - G06T 7/174: involving the use of two or more images - G16H 10/60: for patient-specific data, e.g. for electronic patient records - G16H 30/20: for handling medical images, e.g. DICOM, HL7 or PACS - G16H 30/40: for processing medical images, e.g. editing - G16Z 99/00: Subject matter not provided for in other main groups of this subclass - A61B 5/08: Measuring devices for evaluating the respiratory organs - A61B 2576/023: for the heart - G01R 33/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 (image data processing or generation, in general G06T)} - G06T 2207/10088: Magnetic resonance imaging [MRI] - G06T 2207/20056: Discrete and fast Fourier transform, [DFT, FFT] - G06T 2207/20076: Probabilistic image processing - G06T 2207/30048: Heart; Cardiac - G06T 2207/30104: Vascular flow; Blood flow; Perfusion - A61B: DIAGNOSIS; SURGERY; IDENTIFICATION (usefulness limited to only animals A61D) - G01R: MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES (indicating correct tuning of resonant circuits H03J3/12) - G06T: IMAGE DATA PROCESSING OR GENERATION, IN GENERAL - G16H: HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA - G16Z: INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR ## Abstract An MRI image processing and analysis system may identify instances of structure in MRI flow data, e.g., coherency, derive contours and/or clinical markers based on the identified structures. The system may be remotely located from one or more MRI acquisition systems, and perform: perform error detection and/or correction on MRI data sets (e.g., phase error correction, phase aliasing, signal unwrapping, and/or on other artifacts); segmentation; visualization of flow (e.g., velocity, arterial versus venous flow, shunts) superimposed on anatomical structure, quantification; verification; and/or generation of patient specific 4-D flow protocols. An asynchronous command and imaging pipeline allows remote image processing and analysis in a timely and secure manner even with complicated or large 4-D flow MRI data sets. ## Claims ### Claim 1 A method of operation for use with magnetic resonance imaging (MRI) based medical imaging systems, the method comprising: receiving a set of MRI data by at least one processor-based device, the set of MRI data comprising respective anatomical structure and blood flow information for each of a plurality of voxels; identifying one or more instances of structure in the set of MRI flow data by at least one processor-based device, wherein identifying one or more instances of structure in the set of MRI flow data comprises identifying one or more instances of coherency in the set of MRI flow data; and deriving contours in the set of MRI flow data based on the identified one or more instances of structure in the set of MRI flow data by at least one processor-based device. ### Claim 2 (depends on 1) The method of claim 1 wherein identifying one or more instances of coherency in the set of MRI flow data comprises identifying one or more instances of directional coherency in the set of MRI flow data. ### Claim 3 (depends on 1) The method of claim 1 wherein identifying one or more instances of coherency in the set of MRI flow data comprises identifying one or more instances of directional pathline or structural coherency in the set of MRI flow data. ### Claim 4 (depends on 1) The method of claim 1 wherein identifying one or more instances of coherency in the set of MRI flow data comprises identifying one or more instances of Discrete Fourier Transform (DFT) component coherency in the set of MRI flow data. ### Claim 5 (depends on 1) The method of claim 1 wherein identifying one or more instances of coherency in the set of MRI flow data comprises identifying one or more instances of acceleration coherency in the set of MRI flow data. ### Claim 6 (depends on 1) The method of claim 1, further comprising: identifying one or more instances of clinical markers in the set of MRI flow data, by the at least one processor, based on the identified one or more instances of structure in the set of MRI flow data. ### Claim 7 (depends on 6) The method of claim 6 wherein identifying one or more instances of clinical markers in the set of MRI flow data comprises identifying one or more instances of anatomical markers and/or temporal markers in the set of MRI flow data. ### Claim 8 (depends on 6) The method of claim 6 wherein identifying one or more instances of clinical markers in the set of MRI flow data comprises identifying one or more instances of aneurysms, stenosis, or plaque in the set of MRI flow data. ### Claim 9 (depends on 6) The method of claim 6 wherein identifying one or more instances of clinical markers in the set of MRI flow data comprises identifying one or more pressure gradients in the set of MRI flow data. ### Claim 10 (depends on 6) The method of claim 6 wherein identifying one or more instances of clinical markers in the set of MRI flow data comprises identifying one or more instances of anatomical landmarks of a heart in the set of MRI flow data. ### Claim 11 (depends on 1) The method of claim 1 wherein deriving contours in the set of MRI flow data based on the identified one or more instances of structure in the set of MRI flow data comprises deriving contours in the set of MRI flow data that represent various bodily tissues. ### Claim 12 (depends on 11) The method of claim 11, further comprising: autonomously segmenting blood bodily tissue from non-blood bodily tissues by at least one processor-based device; and autonomously segmenting air from bodily tissues by at least one processor-based device. ### Claim 13 (depends on 1) The method of claim 1, further comprising: applying a Discrete Fourier Transform (DFT) to each of a plurality of voxels in the set of MRI data over a plurality of available time-points, by the at least one processor-based device; and examining a number of components of the DFT at each of the voxels over the available time-points by the at least one processor-based device. ### Claim 14 (depends on 13) The method of claim 13 wherein applying a DFT to each of a plurality of voxels in the set of MRI data over a plurality of available time-points comprises applying a DFT to three velocity components (x, y, z) of the blood flow information. ### Claim 15 (depends on 13) The method of claim 13 wherein applying a DFT to each of a plurality of voxels in the set of MRI data over a plurality of available time-points comprises applying a DFT to a resulting velocity magnitude of the blood flow information. ### Claim 16 (depends on 13) The method of claim 13, further comprising: segmenting a blood pool from static tissue based at least part on the DFT components, by the at least one processor, wherein segmenting a blood pool from static tissue based at least part on the DFT components comprises segmenting the blood pool from the static tissue based at least part on a set of low order, non-DC DFT components, autonomously by the at least one processor. ### Claim 17 (depends on 13) The method of claim 13, further comprising: combining a number of DFT components together without regard for relative magnitude or phase to produce a general mask to locate all blood flow within a chest scan, by the at least one processor. ### Claim 18 (depends on 13) The method of claim 13, further comprising: combining a number of DFT components together taking into account relative magnitude or phase of the DFT components to produce a refined mask to identify a particular region of a blood pool in the body, by the at least one processor. ### Claim 19 (depends on 18) The method of claim 18, further comprising: comparing a phase of the DFT components to a time-point of peak systole, and assigning a probability to each voxel based on an amount of deviation of the phase from an expected value, by the at least one processor. ### Claim 20 (depends on 19) The method of claim 19, further comprising: distinguishing between a blood flow in the aorta and a blood flow in the pulmonary arteries based at least in part on the refined mask, by the at least one processor. ### Claim 21 (depends on 19) The method of claim 19, further comprising: identifying a probability cutoff value based at least in part on a histogram of the resulting probability values, autonomously by the at least one processor. ### Claim 22 (depends on 19) The method of claim 19, further comprising: identifying a probability cutoff value for an arterial specific mask based at least in part on a histogram of the resulting probability values, autonomously by the at least one processor; and determining a probability cutoff values for at least one other mask based at least in part on the arterial specific mask, wherein determining a probability cutoff values for at least one other mask based at least in part on the arterial specific mask comprises performing a flood fill of a general blood mask to remove extraneous non-connected pieces. ### Claim 23 (depends on 22) The method of claim 22, further comprising: separating the arterial specific mask into two main portions based at least in part on a number of flow directions and a number or gradients and/or a number of path-lines along with the resulting probability values, by the at least one processor. ### Claim 24 (depends on 23) The method of claim 23, further comprising: distinguishing the aorta and the pulmonary artery from one another based at least in part at least one of an average direction of flow or a relative position of the two main portions in space, by the at least one processor. ### Claim 25 (depends on 19) The method of claim 19, further comprising: producing a probability mask for a wall of the heart, autonomously by the at least one processor. ### Claim 26 (depends on 25) The method of claim 25, further comprising: combining the probability mask for the wall of the heart with a blood flow mask, by the at least one processor. ### Claim 27 (depends on 25) The method of claim 25, further comprising: employing the probability mask for the wall of the heart in performing eddy current correction, by the at least one processor. ### Claim 28 (depends on 25) The method of claim 25, further comprising: employing the probability mask for the wall of the heart to provide at least one of a location and/or a size of the heart in an image, by the at least one processor. ### Claim 29 A processor-based device, comprising: at least one processor; and at least one nontransitory processor-readable medium communicatively coupled to the at least one processor, and which stores at least one of processor-executable instructions or date which when executed causes the at least one processor to: receive a set of MRI data by at least one processor-based device, the set of MRI data comprising respective anatomical structure and blood flow information for each of a plurality of voxels; identify one or more instances of structure in the set of MRI flow data by at least one processor-based device, wherein identifying one or more instances of structure in the set of MRI flow data comprises identifying one or more instances of coherency in the set of MRI flow data; and derive contours in the set of MRI flow data based on the identified one or more instances of structure in the set of MRI flow data by at least one processor-based device. ### Claim 30 (depends on 29) The processor-based device of claim 29 wherein the at least one of processor-executable instructions or date, when executed, causes the at least one processor to identify one or more instances of directional coherency in the set of MRI flow data. ### Claim 31 (depends on 29) The processor-based device of claim 29 wherein the at least one of processor-executable instructions or date, when executed, causes the at least one processor to identify one or more instances of directional pathline or structural coherency in the set of MRI flow data. ### Claim 32 (depends on 29) The processor-based device of claim 29 wherein the at least one of processor-executable instructions or date, when executed, causes the at least one processor to identify one or more instances of Discrete Fourier Transform (DFT) component coherency in the set of MRI flow data. ### Claim 33 (depends on 29) The processor-based device of claim 29 wherein the at least one of processor-executable instructions or date, when executed, causes the at least one processor to identifying one or more instances of acceleration coherency in the set of MRI flow data. ### Claim 34 (depends on 29) The processor-based device of claim 29 wherein the at least one of processor-executable instructions or date, when executed, causes the at least one processor further to identify one or more instances of clinical markers in the set of MRI flow data, by the at least one processor, based on the identified one or more instances of structure in the set of MRI flow data. ### Claim 35 (depends on 34) The processor-based device of claim 34 wherein the at least one of processor-executable instructions or date, when executed, causes the at least one processor to identify one or more instances of anatomical markers and/or temporal markers in the set of MRI flow data in order to identify one or more instances of clinical markers in the set of MRI flow data. ### Claim 36 (depends on 34) The processor-based device of claim 34 wherein the at least one of processor-executable instructions or date, when executed, causes the at least one processor to identify one or more instances of aneurysms, stenosis, or plaque in the set of MRI flow data in order to identify one or more instances of clinical markers in the set of MRI flow data. ### Claim 37 (depends on 34) The processor-based device of claim 34 wherein the at least one of processor-executable instructions or date, when executed, causes the at least one processor to identify one or more pressure gradients in the set of MRI flow data in order to identify one or more instances of clinical markers in the set of MRI flow data. ### Claim 38 (depends on 34) The processor-based device of claim 34 wherein the at least one of processor-executable instructions or date, when executed, causes the at least one processor to identify one or more instances of anatomical landmarks of a heart in the set of MRI flow data identify one or more instances of clinical markers in the set of MRI flow data. ### Claim 39 (depends on 29) The processor-based device of claim 29 wherein the at least one of processor-executable instructions or date, when executed, causes the at least one processor to derive contours in the set of MRI flow data that represent various bodily tissues in order to derive contours in the set of MRI flow data based on the identified one or more instances of structure in the set of MRI flow data. ### Claim 40 (depends on 29) The processor-based device of claim 29 wherein the at least one of processor-executable instructions or date, when executed, causes the at least one processor further to: autonomously segment blood bodily tissue from non-blood bodily tissues by at least one processor-based device; and autonomously segment air from bodily tissues by at least one processor-based device.