The invention relates to a system and method for determining a trigger signal for imaging with an imaging system.
Electrocardiogram (ECG) and contact photoplethysmography (PPG), e.g. measured with a finger oximeter, have been used in standard magnetic resonance imaging (MRI) and computed tomography (CT) systems for cardiac triggering or gating of scans.
MRI typically uses triggering based on ECG measurements for scans which require synchronization to the cardiac pulse and/or heart phase. The ECG measurement holds the so-called R-peak, signaling the electrical activity starting the contraction of the heart muscle, i.e., the start of the systolic phase. Based on the R-peak, an MRI triggering sequence is started consisting of a preparation pulse (pre-pulse) and multiple acquisition windows. Such ECG based triggering may suffer from disadvantages like complex placement of the electrodes as well as multiple wires that can limit patient comfort as well as disturb imaging quality due to magnetic field interference.
Therefore, feasibility of cardiac-triggered MRI with camera-based PPG of the human face has been shown. A fully or semi-automated and contactless camera-PPG solution can significantly simplify the clinical MRI workflow because no contact sensors are required, and it may eliminate the interference of magnetic field to the signal.
However, compared to ECG-based triggering, a possible limitation of camera-PPG based triggering is that a patient-specific delay is present between the R-peak that is commonly used for cardiac triggering, and the trigger derived from PPG signals. This delay is caused by the pulse transit times from the heart to the peripheral skin site and may lead to inaccurate triggering unless compensated for.
The publication US2017055934A1 discloses a method and system for determining a trigger signal. The method for determining a trigger signal for an imaging device is based on receiving a film of a surface of a first body part of a patient. In an embodiment, image values in a region in the images of the film are averaged. The noise of the image values in the region is reduced by way of the averaging. The film is transformed into a temporal signal series on the basis of the averaging. The transformation takes place in such a manner that the temporal series is a measure of the timing pattern of the blood circulation in the area.
Camera-based PPG requires a certain amount processing in order to achieve robust trigger detection. This may lead to costly delays in the workflow before actual triggered scanning can commence.
Hence, there is a need for improved speed and accuracy of triggered imaging.
It is an object of the invention to provide a method and system to improve the speed and accuracy of triggered imaging.
The invention is defined by the independent claims. Advantageous embodiments are defined in the dependent claims.
According to a first aspect of the invention there is provided a method for determining a trigger signal for imaging with an imaging system. The method comprises the steps of:
The patient-specific characteristic is a measured characteristic related to the physiology or physical characteristics of the patient, such as, but not limited to, heart rate, respiration, oxygenation, blood pressure, cardiac cycle, temperature, weight, height, body-mass-index etc. The patient-specific characteristic may be a constant, one or several measurement points, amplitude, frequency, phase, shape, signal-to-noise level etc. of a physiological measurement. This characteristic or combinations of characteristics can be used to model parameters that are specific for the patient in a patient-specific model. Similarly, the triggering waveform is a time-resolved measurement of a physiological parameter of the patient, such as, but not limited to, heart rate, respiration, oxygenation, cardiac cycle, temperature etc. The trigger signal directly or indirectly triggers acquisition of an image of the patient with the imaging system at a moment in time.
The proposed method may be advantageous for speed and accuracy of the imaging triggering workflow since one or multiple patient-specific characteristics can be measured, and generation of a patient-specific model from such characteristics can be initiated, already before the patient is in the second position where the triggering waveform is measured. The second position may also be the position of the patient when the imaging system is triggered to acquire an image of the patient. In this way, the combination of patient-specific information and the measured triggering waveform may form the basis for accurate imaging triggering without causing unnecessary delays to the workflow.
In an embodiment of the invention, in the first position a patient anatomy of interest is located substantially outside of a bore of the imaging system, and in the second position said patient anatomy of interest is located substantially inside the bore of the imaging system.
The patient anatomy of interest is a part of the patient where triggering is determined and/or the body part that is imaged by the imaging system. By measuring patient-specific characteristics when the patient anatomy of interest is outside of the bore of the imaging system, such as e.g. during a patient preparation phase, patient time inside the bore can be shortened. In this way the total workflow time may be shortened, and patient comfort may be increased. The bore of an imaging system may be e.g. the bore of a computed tomography (CT) system or a magnetic resonance imaging (MRI) system or a molecular imaging system, like positron emission tomography (PET) or single photon emission computed tomography (SPECT), or any combination of such systems.
In an embodiment of the invention, the triggering waveform comprises a waveform derived from photoplethysmography (PPG), and the step of measuring the triggering waveform comprises a measurement with a sensor, preferably a camera.
A measurement of the triggering waveform with e.g. camera-based PPG, instead of multiple electrocardiogram (ECG) sensors placed in contact with the patient may have multiple advantages, such as automation or simplification of the workflow, avoidance of magnetic field interference of the ECG signal etc.
The sensor is preferably e.g. a 2D or 3D camera sensor or similar, configured to contactless measure PPG parameters from visible light or infrared light or light of other wavelengths. However, other sensor technologies are also conceivable to measure a waveform derived from PPG, such transmissive or reflective PPG sensors in close proximity to or contact with the patient, contactless measurements with RADAR or LIDAR, etc.
In an embodiment of the invention, the patient-specific characteristic comprises a characteristic derived from PPG, and the step of measuring the patient-specific characteristic comprises a measurement with a sensor, preferably a camera.
Similarly to the triggering waveform, a measurement of patient-specific parameters with e.g. camera-based PPG, alone or in combination with other measurements, may speed up and simplify the imaging and/or patient preparation workflow.
The sensor used to measure a patient-specific characteristic, comprising a characteristic derived from PPG, may be the same sensor as the sensor used to measure the triggering waveform. The sensor may also be a separate sensor of the same, similar, or different type or functionality as the sensor used to measure the triggering waveform. By way of example, both sensors may e.g. be sensors configured to measure camera-based PPG signals.
In an embodiment of the invention, the patient-specific model comprises a patient-specific PPG waveform profile.
Modelling of a patient-specific model comprising a PPG waveform profile, such as a template morphology and/or amplitude of the waveform, is advantageous for fast and accurate imaging triggering. This may be of particular importance when the triggering waveform comprises a waveform derived from PPG.
In an embodiment of the invention, the patient-specific model comprises a patient-specific pulse transit time.
Similarly to the patient specific-waveform profile, information about a patient-specific pulse transit time may be advantageous for fast and accurate imaging triggering.
In an embodiment of the invention, the patient-specific model comprises a patient-specific reliability factor of a camera-based PPG parameter.
The patient-specific reliability factor may provide an indication of the reliability, quality and/or robustness of camera-based PPG for the specific patient and/or circumstances. If e.g. the reliability is low, additional or alternative measurements to camera-based PPG may be required for accurate triggering based on the triggering waveform. Having such information at an early stage may save bore time and/or avoid that sets of image acquisitions need to be repeated.
In an embodiment of the invention, the step of determining the trigger signal comprises calculating multiple PPG markers, wherein the PPG markers are derived from the triggering waveform and/or the patient-specific PPG waveform profile, and wherein the PPG markers are configured to predict an R-peak.
PPG markers are specific PPG features of the PPG waveform, such as inter-beat interval, waveform features, phase etc. Calculating multiple PPG markers may be advantageous, since the PPG markers may improve accuracy and/or robustness when predicting an R-peak from e.g. camera-based PPG.
In an embodiment of the invention, the PPG markers as mentioned above comprise a diastolic phase PPG marker and a flush phase PPG marker.
In an embodiment of the invention, the method further comprises the steps of:
In an embodiment of the invention, the method further comprises an image formation step, wherein the image formation step forms an image of the patient based on processing of the imaging scan data and the timing parameter.
This may advantageous, since using a timing parameter from the triggering waveform and/or the patient-specific model can help to improve the quality and/or speed of image formation.
According to another aspect of the invention, there is provided a computer program element, which, when being executed by at least one processing unit, is adapted to cause the processing unit to perform a for the processing unit any computer-executable method step according to any of the methods described above.
According to another aspect of the invention, there is provided a computer readable medium having stored thereon the above-mentioned computer program element.
According to a second aspect of the invention, there is provided a system for determining a trigger signal for imaging a patient with an imaging system, the system comprising:
In an embodiment of the invention, at least one of the patient-specific characteristic sensor and/or the triggering waveform sensor is a sensor, preferably a camera, and the sensor is configured to measure a sensor-based PPG signal.
According to another aspect of the invention, there is provided an imaging system comprising the system for determining a trigger signal, wherein the imaging system is preferably a computed tomography system or a magnetic resonance imaging system.
These and other aspects of the invention will be apparent from and elucidated with reference to the embodiments described hereinafter.
A schematic of a system 100 for determining a trigger signal, in accordance with embodiments of the invention, is illustrated in
Many imaging system environments already offer a ceiling-mounted camera that allows for monitoring the patient during the preparation phase, before moving the table into the bore. Such camera sensors can therefore be further used for e.g. remote PPG analysis. The patient-specific characteristic sensor 10 and the triggering waveform sensor 20 are connected, via a wired or wireless connection, to a processing unit (not shown in the figure).
In the example illustrated in
In
During the phase illustrated in
In
Several potential advantages can be envisioned from such a system as illustrated in
With knowledge of patient-specific PPG models, such as PPG waveform morphology, PPG amplitude, pulse transit time delays etc., camera-based accurate and fast timing of e.g. the MRI and CT acquisition windows is simplified.
Images of MRI cardiac cine movies are commonly annotated with time stamps relative to the R-peak, which is defined as time zero. This is possible with prior knowledge of the delays.
MRI cardiac cine movies are typically presented for diagnosis and comparison such that they start with the time of the R-peak. This is possible with prior knowledge of the delays.
Functional evaluation of cardiac images as for ejection fraction and wall-motion may be simplified and/or more accurate with knowledge about timing relative to the R-peak.
Triggering with a low delay between R-peak and trigger requires a patient specific PPG template and therefore some initial calibration time to initialize this model. When such calibration can be made in advance by the processing unit, without relying on an in-bore measurement, expensive in-bore time can be saved.
When camera-based PPG signal quality can be evaluated in advance, to calculate a patient-specific reliability factor of a camera-based PPG parameter, this can have significant positive effects on the workflow time. E.g. in an example of triggered MRI, the workflow would be prolonged drastically if the patient is already inside the bore and the in-bore camera 20 detects a patient-specific insufficient PPG signal quality. In that case the patient 30 would have to be taken out again, MR coils removed, ECG-electrodes applied, MR coils applied again, and then the patient 30 would have to be repositioned in the bore.
In the example illustrated in
Current ECG-detection for MRI triggering is typically based on at least four electrodes on the bare chest. At least four electrodes are typically required for robust R-peak detection inside the bore, because here the ECG signal may be corrupted due to interference with MR fields. For synchronous measurement of ECG with PPG outside of the bore for the purpose of measuring patient-specific characteristics and initialization of calculating a patient-specific model, it is possible to apply only two ECG-electrodes. The electrodes may be applied e.g. at the wrists or hands of the patient 30. Such an ECG set-up may simplify the workflow tremendously. Furthermore, in such cases where the electrodes are removed before the patient enters the bore, since in-bore triggering waveform is measured with e.g. camera-PPG, the ECG electrodes, leads and detection circuitry do not have to be made MR-compatible. It can be envisioned that in some cases the ECG electrodes may just be simple handles or conductive fields at both sides of the patient support that may be grabbed for a few seconds by the patient for as long as data is required to calibrate the model.
Similarly.
As shown in
Similarly, as shown in
A trigger signal is generated based on the identified cardiac cycles and R-peaks, such that the imaging system is triggered to acquire images of the patient.
In the example illustrated in
Additional parameters like patient height, age, body-mass index etc. may also be used to further improve the models.
In order to improve the robustness of PPG marker detection from the triggering waveform, multiple PPG markers can be generated in combination with patient-specific PPG waveform profiles. In this way one marker may predict the next and thereby improve the robustness of the model.
Instead of creating one single patient-specific PPG waveform profile, it is also possible to create a multitude of profiles associated with specific cardiac cycle durations or for ranges of cardiac cycle durations. Creating different profiles per ranges of cycle duration may be motivated by the fact that the time curves of cardiac beat and related physiological quantities as arterial blood velocity and PPG signal do not simply scale linearly if the heart rate changes. Instead, some parts of the cardiac cycle are extended in time (as typically the resting phase in diastole) and some parts remain more or less unchanged (as the cardiac contraction phase).
In each waveform profile, one or more PPG markers can be defined. Of particular relevance are the zero-crossings in this signal. The negative-going zero-crossing is associated with the highest blood volume change at the observation point, e.g. forehead with camera-PPG, in response to a contraction of the heart. The positive-going zero-crossing is associated with the diastolic phase and typically precedes the ECG R-peak. A plurality of markers may by generated from a single waveform profile thus enabling tracking the current phase of the PPG signal in the cardiac cycle.
Detection of PPG markers can e.g. be done as follows. The waveform profile can be shifted so that it starts and stops with a zero-crossing. In this way two waveform profiles are created, each with a full cardiac waveform: one beginning and ending with a positive-going zero-crossing and the second one beginning and ending with a negative-going zero-crossing.
Both waveforms are shown in
The patient-specific PPG signal may be decomposed using a patient-specific model. The model describes the waveform as being composed of three parts: the PPG waveform profile, a low frequency (LF) signal component and a noise component. The LF signal component can be modelled as a low-order polynomial. More formally, let the discrete-time PPG signal segment be denoted as s(n), n=[1, . . . , N], the prototype as p(n), the LF signal component as b(n), then the signal is described as
with c(n) the error signal describing noise (e.g., measurement noise) and disturbances, and with a being a scale factor. The signal component b captures the signal drifting in the PPG and PPG variation associated with respiration (and thus of lower frequency than the heart rate). Various ways to arrive at such a decomposition are known, e.g., a least-squares fitting procedure may be used to attain said decomposition, where windowing may be used to stress the importance of certain parts of the segment in the fit.
The decomposition signals can be analyzed for their strength, e.g., in terms of energy. A good fit is attained if the energy associated with the PPG component p is significantly larger than that associated with e. Further evidence supporting the generation of the marker would be that the coefficient a is positive. Further evidence would be that when considering the next segment consisting of one simple shifted input signal, the coefficient a has become lower, or the balance between the energies in (a p) and e has become less favorable. In these or other ways, a decision may be taken that the marker has occurred.
Since the length of the current cardiac cycle is a priori unknown, the system operates at various lengths. Thus, the detection has an extra dimension, namely that along the scale (cardiac cycle length). The various cues for the existence of the marker also have to be seen across this extra dimension. A hard decision may be taken based on thresholding, more statistical-like approaches may be used, or neural nets may be trained to create the final decision mechanism
Examples of multiple generated PPG markers to improve R-peak prediction are illustrated in
In the above description, several embodiments of the invention are elaborated on. An aspect of the invention is also illustrated by
Similarly,
Besides achieving a more robust triggering from the measured triggering waveform, patient-specific models, e.g. including PPG markers, may contain timing parameters that can be further leveraged during the imaging workflow. An image formation process S60, where scan data is combined into an image or video sequence, may use such timing parameters as extra information to control the image formation process. E.g. by weighting or as a compensation factor for the data obtained at different data acquisition windows.
E.g. in MRI, knowledge about timing parameters as the distance between pre-pulse generation and data acquisition timing can be used to improve image formation. Alternatively, the trigger mechanism can be configured to keep stricter requirements on constant pre-pulse to acquisition window distance. In this case the information of the occurrence of each MR sequence event, such as the motion state of the myocardium, can be used during the image formation process.
It should be noted that the above-mentioned embodiments illustrate rather than limit the invention, and that those skilled in the art will be able to design many alternative embodiments without departing from the scope of the appended claims. In the claims, any reference signs placed between parentheses shall not be construed as limiting the claim. The word “comprising” does not exclude the presence of elements or steps other than those listed in a claim. The word “a” or “an” preceding an element does not exclude the presence of a plurality of such elements. The invention may be implemented by means of hardware comprising several distinct elements, and/or by means of a suitably programmed processor. In the device claim enumerating several means, several of these means may be embodied by one and the same item of hardware. Measures recited in mutually different dependent claims may advantageously be used in combination.
Number | Date | Country | Kind |
---|---|---|---|
22183591.1 | Jul 2022 | EP | regional |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2023/068277 | 7/4/2023 | WO |