The disclosed subject matter relates to diagnosis and risk evaluation of vertebral artery dissection (VAD).
A VAD is a flap-like tear of the inner lining of the vertebral artery, which is located in the neck and supplies blood to the brain. VAD typically results in blood entering the arterial wall, which in turn may form a blood clot, thickening the artery wall and often impeding blood flow.
Although incidents of VAD are relatively few in the general population (affecting about 1.1 per 100,000 as reported in [1] Schievink W I, Roiter V. Epidemiology of cervical artery dissection, Front Neurol Neurosci. 2005; 20: 12-15, incidents have been rising in the past few decades. Further, VAD may be a primary cause of ischemic stroke in young and middle-aged populations. Moreover, a fairly large percentage of patients who suffer a VAD (approximately 67%-85% of all VAD cases) have a resulting permeant condition or suffer a stroke, as disclosed in [2] Kim Y-K, Schulman S. Cervical artery dissection: pathology, epidemiology and management, Thromb Res. 2009; 123: 810-821.
VAD can be classified into two types: spontaneous VAD and traumatic VAD. Spontaneous VAD generally originates from intrinsic factors, specifically a variety of underlying arteriopathies. Traumatic VAD on the other hand is generally caused by extrinsic factors, such as an injury to the neck. However the cause of a VAD incident can be combinative. That is, Patients with weak vertebral arteries can be more susceptible to extrinsic trauma and traumatic VAD.
Studies, such as [3] Debette S, Leys D. Cervical-artery dissections: predisposing factors, diagnosis, and outcome, Lancet Neurol, 2009; 8: 668-678, have shown possible associations between spontaneous VAD and a number of underlying diseases that can cause weaknesses in the vertebral arterial wall. These diseases include connective tissue diseases, inflammation, and certain types of genetic disorders. However, as noted in [4] Debette S, Markus H S, the genetics of cervical artery dissection: a systematic review, Stroke. 2009; 40: e459-66, none of these possible associations have been convincingly proven.
Vascular imaging techniques, such as computed tomography angiography (CTA) and magnetic resonance angiography (MRA), have typically been used for the diagnosis of VAD. However, so far to date, such imaging techniques have not been used or proposed to be used as a predictive tool for assessing vulnerabilities of the vertebral artery and for evaluating risks of developing future VAD in asymptomatic patients.
As noted above, recent studies have identified a number of potential contributors to the occurrence of VAD. Most of these risk factors are associated with underlying arteriopathies that may lead to vulnerabilities in the vertebral arteries. Such arteriopathies include, but are not limited to hereditary connective tissue diseases, fibromuscular dysplasia, and arterial inflammation. These risk factors have been respectively studied in [5] Rouviére S, Michelini R, Sarda P, Pagés M., Spontaneous carotid artery dissection in two siblings with osteogenesis imperfecta, Cerebrovasc Dis. 2004; 17: 270-272, Cervical-artery dissections: predisposing factors, diagnosis, and outcome (cited above), and [2] Cervical artery dissection: pathology, epidemiology and management (cited above).
Recently, ultrasound imaging has been used to assess arterial wall weakness in carotid arteries, which has been associated with the occurrence of carotid artery dissection, as explained in [6] Calvet D, Boutouyrie P, Touze E, Laloux B, Mas J-L, Laurent S. Increased stiffness of the carotid wall material in patients with spontaneous cervical artery dissection, Stroke. 2004; 35: 2078-2082.
Contrast-enhanced magnetic resonance imaging (MRI) has been used as a follow-up imaging tool post occurrence of a VAD incident, as disclosed in [7] Provenzale J M. MRI and MRA for evaluation of dissection of craniocerebral arteries, lessons from the medical literature. Emerg Radiol, 2009; 16: 185-193.
However, despite the use of these known techniques, there lacks an imaging-based predictive tool to evaluate the risks of VAD for prevention thereof.
The state of the art of four-dimensional phase-contrast magnetic resonance imaging (4D PC-MRI) will now be described.
Disturbed blood flow to the brain has been associated with several neurological diseases, from stroke and vascular diseases to Alzheimer's and cognitive decline.
Two dimensional phase-contrast magnetic resonance imaging (2D PC-MRI) has been clinically used to study blood flow in selected transverse planes of the cerebral arteries as disclosed in [8] Lotz J, Döker R, Noeske R, Schüttert M, Felix R, Galanski M, et al, In vitro validation of phase-contrast flow measurements at 3 T in comparison to 1.5 T: precision, accuracy, and signal-to-noise ratios. J Magn Reson Imaging. 2005; 21: 604-610.
In recent years, 4D PC-MRI has been developed and performed primarily in a research capacity to study time-resolved three dimensional (3D) blood flow in the aorta and carotid arteries, as discussed in [9] Stankovic Z, Allen B D, Garcia J, Jarvis K B, Markl M. 4D flow imaging with MRI, Cardiovasc Diagn Ther. 2014; 4: 173-192.
Due to their smaller size, vertebral arteries are not often targeted for 4D flow imaging. Indeed it is traditionally difficult to target these arteries for 4D flow imaging.
Moreover, because of the difficulties in the acquisition and post-processing of 4D PC-MRI, 4D PC-MRI is rarely used clinically.
Further, 4D PC-MRI has not yet been used to evaluate the risks of the occurrence of VAD.
The instant disclosure has been developed in light of the above circumstances.
Certain embodiments provide a predictive tool for evaluating the risk of developing VAD in certain high-risk asymptomatic patients. This and other embodiments may utilize 4D PC-MRI to visualize and quantify blood flow in the vertebral arteries. In certain embodiments, the time-resolved 3D flow velocity field extracted from the 4D PC-MRI may be processed to extract various hemodynamic variables, such as, for example, pulsatile wave velocity and arterial wall shear stress, to assess the healthiness of vertebral arteries and dynamic interactions between the vertebral arterial wall and blood flow.
In certain embodiments, the hemodynamic variables, may include, but are no way limited to the 4D flow velocity field, flow pulsatile velocity, and time-resolved distribution of the wall shear stress. These illusory hemodynamic variables may be concatenated to form a complex hemodynamic profile of a patient's vertebral arteries.
In certain embodiments, this profile may be used to train the aforementioned predictive tool using machine learning.
Also, in certain embodiments, the flow-imaging-based hemodynamic variables may be further integrated with the other imaging/laboratory/genetic test data, for example, to generate a comprehensive predictive tool for VAD. Further, in certain embodiments, such an integration may be performed through statistical analysis and/or machine learning.
That is, certain embodiments of the instant disclosure provide for a risk evaluation tool for VAD based on analysis of the hemodynamics in the vertebral arteries using 4D PC-MRI and the optional integration of the hemodynamic variables with other anatomic/clinical information.
Accordingly, certain embodiments of the instant disclosure provide for a method, apparatus, and storage medium for imaging time-resolved 3D flow fields with high spatial resolution in relatively small vessels, post-processing of 4D PC-MRI with minimal user interaction, non-invasive assessment of the healthiness and vulnerability of the vertebral arteries of an individual using flow imaging, and information fusion of multi-modality and multi-source data for the evaluation of VAD risks.
One or more embodiments provide a sememe prediction method, a computer device, and a storage medium.
According to an aspect of an embodiment, there is provided a method for vertebral artery dissection risk evaluation that includes obtaining four-dimensional phase-contrast magnetic resonance imaging data, performing pre-processing of the four-dimensional phase-contrast magnetic resonance imaging data, obtaining at least one blood hemodynamic marker from the four-dimensional phase-contrast magnetic resonance imaging data, a classifying the at least one blood hemodynamic marker as a hemodynamic predictor of vertebral artery dissection, and creating a comprehensive risk evaluation of vertebral artery dissection using the hemodynamic predictor.
According to an aspect of an embodiment, there is provided an apparatus for vertebral artery dissection risk evaluation that includes at least one memory configured to store computer program code; at least one hardware processor configured to access said computer program code and operate as instructed by said computer program code, said computer program code including: first obtaining code configured to cause said at least one hardware processor to obtain four-dimensional phase-contrast magnetic resonance imaging data, pre-processing code configured to cause said at least one hardware processor to perform pre-processing of the four-dimensional phase-contrast magnetic resonance imaging data, second obtainment code configured to cause said at least one hardware processor to obtain at least one blood hemodynamic marker from the four-dimensional phase-contrast magnetic resonance imaging data, classification code configured to cause said at least one hardware processor to classify the at least one blood hemodynamic marker as a hemodynamic predictor of vertebral artery dissection, and creation code configured to cause said at least one hardware processor to create a comprehensive risk evaluation of vertebral artery dissection using the hemodynamic predictor.
According to an aspect of an embodiment, there is provided a non-transitory computer-readable medium storing instructions for vertebral artery dissection risk evaluation, the instructions comprising: one or more instructions that, when executed by one or more processors of a device, cause the one or more processors to: obtain four-dimensional phase-contrast magnetic resonance imaging data, pre-process the four-dimensional phase-contrast magnetic resonance imaging data, obtain at least one blood hemodynamic marker from the four-dimensional phase-contrast magnetic resonance imaging data, classify the at least one blood hemodynamic marker as a hemodynamic predictor of vertebral artery dissection, and create a comprehensive risk evaluation of vertebral artery dissection using the hemodynamic predictor.
To make the objectives, technical solutions, and advantages of this application be more clear and comprehensible, embodiments will be further described in detail with reference to the accompany drawings. It should be understood that, the specific implementations described herein are only used for interpreting this application, rather than limiting this application.
Embodiments are not limited to the structure shown in
Referring to
S310: Obtain four-dimensional phase contrast magnetic resonance image (4D PC-MRI) data.
Initially, it should be understood that prior to obtainment of the 4D PC-MRI data, a time-of-flight MR angiography (TOF MRA) or a contrast-enhanced MRA can be performed to serve as a localizer scan, in which a 3D region of interest (ROI) where vertebral arteries reside can be selected. Here, a larger ROI area may be selected to include other important cerebral arteries. However, a larger ROI may increase the scan time. The diameter of the vertebral arteries within the ROI can also be assessed using the localizer MRA, from which a minimal arterial diameter may be used to direct a setting of the spatiotemporal resolution of the PC-MRI. In certain embodiments, the transverse luminal area of the artery may cover enough voxels for a reliable quantification of flow velocity. In certain embodiments, the in-plane spatial resolution of the 4D PC-MRI may be set to 0.22×Diameter_min. The spatial resolution in the axial direction may be set to ≤2 mm. The temporal resolution may be set to <40 ms. The velocity encoding parameter (VENC) may be set to <150 cm/sec. In certain embodiments, 4D PC-MRI of the vertebral arteries may be performed with ECG-gating. The scan parameters of 4D PC-MRI may be determined based on considerations of both image quality and total scan time. When Gadolinium-based MRI contrast is used, for example, in certain embodiments, performing 4D flow imaging after the contrast-enhanced studies can improve the blood-to-tissue contrast and the velocity-to-noise ratio in the 4D PC-MRI images. When available, imaging acceleration methods may be used to shorten the acquisition time and improve the image quality.
S320: Pre-process the four-dimensional phase-contrast magnetic resonance imaging data. In this step, preprocessing is carried out on the obtained 4D PC-MRI data.
A number of sources may contribute to flow quantification errors in raw 4D PC-MRI data. While some sources of these errors may be compensated and corrected automatically on an MRI scanner, for example, typically, there are two phase errors that are addressed in the pre-processing S320 step.
First, the background phase offset induced by eddy currents is compensated. In certain embodiments, regions of static tissues in the 4D flow image, obtained from the 4D PC-MRI data, is identified using thresholding methods. Additionally, or in the alternative, a user can estimate the eddy currents-induced background phase offset errors using polynomial fitting, and subsequently remove the phase offset from the 4D flow data.
Second, phase correction is performed, if necessary, for example, when phase aliasing occurs. Certain embodiments may employ one or several phase-unwrapping algorithms.
In addition to correcting background phase offset and performing phase correction, the pre-processing step S320 may also include segmentation and tracking of the target arteries in certain embodiments. For instance, in some embodiments, flow path-line tracing may be performed only within the boundary of the artery lumen. In some embodiments, arteries in the magnitude image of the PC-MRI are segmented and tracked. Also, in some embodiments, automated segmentation of the arteries may be performed by first tracing the arterial centerlines and then performing the lumen segmentation using deformable models with a tubular shape.
In certain embodiments, not necessarily illustrated in
Referring back to
Since 4D PC-MRI provides full volumetric coverage of the ROI, the vertebral artery dissection (VAD) diagnosis and risk evaluation method provides an unique option of retrospective selection of 2D slices or 3D sub-regions in the 3D field of view for 3D flow quantification and analysis. Thus, besides conventional 2D flow parameters, such as, for example, transvalvular gradient and peak flow velocity, a number of advanced 4D blood hemodynamic markers can be harvested from the 4D PC-MRI image data. Some of these advanced markers are discussed below. However, it should be understood that the markers are not limited to those discussed below.
Shear rate (SR) may be calculated as a spatial gradient of the flow velocity field. It may be associated with the blood thrombus process because it is associated with forces experienced by blood components such as platelets and red blood cells.
Wall shear stress (WSS) is the friction force blood flow exerts on the vertebral arterial wall. It can be estimated in certain embodiments by taking the derivative of 4D flow velocity near the vessel wall boundary. WSS is believed to play an important role in the regulation of the functions of the endothelial cell and the extracellular matrix in the vessel wall. For example, low WSS has been associated with the development of atherosclerosis, and high WSS has been associated with vessel dilation and the formation of aneurysms.
Pulse wave velocity (PWV) is the propagation speed of the systolic wave front through the artery. It is a direct indicator of arterial wall stiffness and an important predictor of arteriopathy progression in patients with hypertension and connective tissue diseases. In order to automatically measure PWV in a 4D flow image, in certain embodiments, velocity waveforms can be measured at selected sites along the centerline of the vertebral artery. Then, PWV may be calculated as the ratio of the distance between measurement sites and the transit-time.
Flow eccentricity (FE) may lead to jet impingement on the vertebral artery wall, and may be associated with weakness in the vessel wall and the occurrence of VAD.
As noted above, it should be understood that the above advanced hemodynamic markers are not all of the advanced hemodynamic markers that may be obtained from the 4D PC-MRI image data. Other advanced hemodynamic markers derived from the 4D PC-MRI image data may include, but are in no way limited to turbulence, kinetic energy, energy dissipation, relative pressure fields, and flow displacement.
Accordingly in S330, at least one of these blood hemodynamic markers is obtained from the 4D PC-MRI image data.
It will be understood that the aforementioned methods of obtaining the aforementioned advanced hemodynamic markers are in no way limiting. Indeed, certain embodiments of the disclosure may obtain the aforementioned advanced hemodynamic markers in different manners.
S340: Classify the at least one blood hemodynamic marker as a hemodynamic predictor of vertebral artery dissection.
Here, the obtained advanced blood hemodynamic marker(s) are classified. In certain embodiments, this classification may be performed by deep learning. However, other classification methods may also be used.
Additionally, when other parameters (e.g. not necessarily advanced blood hemodynamic markers) are available, such as, for example, artery geometric measurements derived from the contrast-enhanced CTA or MRA, patient clinical and medical information, laboratory test results, genetic test results, and potential risk level of the extrinsic factors related to VAD, these additional parameters may also be classified. In certain embodiments, these additional parameters may be classified using deep learning. However, other classification methods may also be used.
This classification process is illustrated in more detail in
In the embodiment depicted in
Referring again to
The above described method uses 4D PC-MRI flow imaging to extract hemodynamic information that is closely related to the healthiness of the vertebral arteries. The discussed embodiments achieves the following functions:
Acquisition of high-resolution 4D PC-MRI image data in the vertebral arteries. Embodiments may provide a guideline to achieve high-resolution 4D PC-MRI image data of the cerebral and extracerebral vessels using commercial MRI scanners.
Post-processing of 4D PC-MRI image data with minimal user interactions for the extraction of the time-resolved 3D flow velocity in the vertebral arteries, in certain embodiments.
Extraction of in-vivo hemodynamic variables that are associated with vulnerability of the vertebral arterial wall, in certain embodiments. A variety of hemodynamic parameters related to arteriopathy may be identified and extracted from the 4D PC-MRI image data in certain embodiments.
Study the healthiness of the vertebral arterial wall in vivo using machine learning techniques with the hemodynamic variables as the input features, as discussed above with reference to certain embodiments.
Integration of hemodynamic information and other imaging/clinical/laboratory/genetic testing results for the comprehensive risk evaluation of VAD, as discussed above with reference to certain embodiments.
In certain high-risk asymptomatic populations, which are prone to either intrinsic, extrinsic, or both factors of VAD risk, identifying patients with an underlying vertebral arteriopathy and advising proactive prevention of VAD may be beneficial. For instance, patients with family histories of spontaneous arterial dissection could benefit from a risk-evaluation test for VAD. Also, for athletes in competitive sports, such a screening tool would be much needed by both the athlete community and the sports industry.
Certain embodiments of the instant disclosure provide for the evaluation of other relatively smaller and deeper arteries (e.g., diameter range: 3-5 mm; not easily accessible by ultrasound).
Also, the above-described embodiments may alternative, or combinatively be modified as follows.
The afore-described classifications may be replaced by other machine learning-based or statistics-based methods, that are not necessarily rooted in deep learning. This may be especially true in embodiments, for which very limited training data is available.
Segmentation of the vertebral arteries in the magnitude image of the 4D PC-MRI may be performed by using other segmentation methods, such as the 3D levelset method.
In embodiments where training data is limited, for classifying the parameters, mean values may be used. Another approach would be to treat the missing data as hidden variables, and use an EM algorithm to estimate them.
The VAD diagnosis and risk evaluation apparatus/method corresponds to the other of the VAD diagnosis and risk evaluation apparatus/method, and the specific technical features that correspond are not repeated herein.
A person of ordinary skill in the art may understand that all or some of the modules, units, components and procedures of the foregoing embodiments may be implemented by a computer program instructing relevant hardware. The program may be stored in a non-volatile computer readable storage medium. When the program is executed, the program may control the hardware to execute the procedures of the embodiments of each foregoing method. Any usage of a memory, storage, a database or other media in each embodiment of this application may include non-volatile and/or volatile memories. The non-volatile memory may include a read-only memory (ROM), a programmable ROM (PROM), an electrically programmable ROM (EPROM), an electrically erasable programmable ROM (EEPROM), or a flash memory. The volatile memory may include a random access memory (RAM) or an external cache memory. For description, rather than for limitation, RAM may be in various forms, for example, a static RAM (SRAM), a dynamic RAM (DRAM), a synchronous DRAM (SDRAM), a double data rate SDRAM (DDRSDRAM), an enhanced SDRAM (ESDRAM), a Synchlink DRAM (SLDRAM), a Rambus direct RAM (RDRAM), a directly memory bus dynamic RAM (DRDRAM), and a memory bus dynamic RAM (RDRAM).
Each technical feature in the foregoing embodiments may be combined randomly. For simplified description, not all possible combinations of each technical feature in the foregoing embodiments are described. However, the combinations of the technical features shall be considered to fall within the scope of the specification as long as the combinations are not contradictory. The foregoing embodiments only describe several implementations of this application, and their description is specific and detailed, but cannot therefore be construed as a limitation to the patent scope of the present disclosure. It should be noted that a person of ordinary skill in the art may further make variations and improvements without departing from the conception of this application, and these all fall within the protection scope of this application. Therefore, the patent protection scope of this application should be subject to the appended claims.