Develop blood oxygen level dependent signal by metabolic/hemodynamic model using numerical methods

Background and objective: The metabolic/hemodynamic (MH) model describes the blood flow mechanisms as well as the coupling between the hemodynamic responses and the metabolic activities in a blood vessel in the human brain. In the existing MH model, the blood flow out from a blood vessel is formulat...

Full description

Bibliographic Details
Main Authors: Md. Roni Islam, Sheikh Md. Rabiul Islam
Format: Article
Language:English
Published: Elsevier 2023-01-01
Series:Informatics in Medicine Unlocked
Subjects:
Online Access:http://www.sciencedirect.com/science/article/pii/S235291482300059X
_version_ 1797845075604013056
author Md. Roni Islam
Sheikh Md. Rabiul Islam
author_facet Md. Roni Islam
Sheikh Md. Rabiul Islam
author_sort Md. Roni Islam
collection DOAJ
description Background and objective: The metabolic/hemodynamic (MH) model describes the blood flow mechanisms as well as the coupling between the hemodynamic responses and the metabolic activities in a blood vessel in the human brain. In the existing MH model, the blood flow out from a blood vessel is formulated as dependent only on the capillary volume. In fact, the blood flow out from a blood vessel depends not only on the capillary volume but also on the blood flow into the capillary bed. For this reason, the blood flow out formula of the existing model has been modified. In addition to implementing existing model modification to obtain better accuracy, we have used new methods to solve the model instead of conventional methods. Method: The MH model describes physical phenomena of a blood vessel by eight processes equations (PEs). These PEs are often solved by using a local linearization (LL) scheme and the Taylor series method. In addition to the previously used Taylor series method, we have also used the Euler method and the Runge–Kutta (RK) method to solve the model instead of a LL scheme for estimating dynamical variables (DVs). By using these DVs, a Blood Oxygen Level Dependent (BOLD) signal is generated through a well-defined observation equation (OE). There are two OEs, called Obata and Friston. The Friston OE produces a BOLD signal from the cerebral blood volume and deoxy-hemoglobin content with their nonlinear properties; conversely, the Obata OE produces a BOLD signal without considering nonlinear properties. For this reason, we have used the Friston OE instead of the Obata OE to estimate the BOLD signal perfectly. Results: At 20% resting oxygen extraction fraction (ROEF), the BOLD signals of the modified and the existing model are identical, but when the ROEF increased up to 50% at its standard value, the modified model accuracy is increased by 16.12%–23.07% more than that of the existing model. The Euler and RK methods generate a BOLD signal 6.95% more accurately than that of Taylor series method from the modified model. Conclusion: In the model inversion process, this research will be helpful to estimate the model parameters and hidden states accurately.
first_indexed 2024-04-09T17:32:39Z
format Article
id doaj.art-96a4add3ba3f4a99a8dc90960d56521d
institution Directory Open Access Journal
issn 2352-9148
language English
last_indexed 2024-04-09T17:32:39Z
publishDate 2023-01-01
publisher Elsevier
record_format Article
series Informatics in Medicine Unlocked
spelling doaj.art-96a4add3ba3f4a99a8dc90960d56521d2023-04-18T04:08:56ZengElsevierInformatics in Medicine Unlocked2352-91482023-01-0138101217Develop blood oxygen level dependent signal by metabolic/hemodynamic model using numerical methodsMd. Roni Islam0Sheikh Md. Rabiul Islam1Department of Electronics and Communication Engineering, Khulna University of Engineering and Technology, Khulna, 9203, BangladeshCorresponding author.; Department of Electronics and Communication Engineering, Khulna University of Engineering and Technology, Khulna, 9203, BangladeshBackground and objective: The metabolic/hemodynamic (MH) model describes the blood flow mechanisms as well as the coupling between the hemodynamic responses and the metabolic activities in a blood vessel in the human brain. In the existing MH model, the blood flow out from a blood vessel is formulated as dependent only on the capillary volume. In fact, the blood flow out from a blood vessel depends not only on the capillary volume but also on the blood flow into the capillary bed. For this reason, the blood flow out formula of the existing model has been modified. In addition to implementing existing model modification to obtain better accuracy, we have used new methods to solve the model instead of conventional methods. Method: The MH model describes physical phenomena of a blood vessel by eight processes equations (PEs). These PEs are often solved by using a local linearization (LL) scheme and the Taylor series method. In addition to the previously used Taylor series method, we have also used the Euler method and the Runge–Kutta (RK) method to solve the model instead of a LL scheme for estimating dynamical variables (DVs). By using these DVs, a Blood Oxygen Level Dependent (BOLD) signal is generated through a well-defined observation equation (OE). There are two OEs, called Obata and Friston. The Friston OE produces a BOLD signal from the cerebral blood volume and deoxy-hemoglobin content with their nonlinear properties; conversely, the Obata OE produces a BOLD signal without considering nonlinear properties. For this reason, we have used the Friston OE instead of the Obata OE to estimate the BOLD signal perfectly. Results: At 20% resting oxygen extraction fraction (ROEF), the BOLD signals of the modified and the existing model are identical, but when the ROEF increased up to 50% at its standard value, the modified model accuracy is increased by 16.12%–23.07% more than that of the existing model. The Euler and RK methods generate a BOLD signal 6.95% more accurately than that of Taylor series method from the modified model. Conclusion: In the model inversion process, this research will be helpful to estimate the model parameters and hidden states accurately.http://www.sciencedirect.com/science/article/pii/S235291482300059XBOLD signalMH modelTaylor series methodEuler methodRunge–Kutta method
spellingShingle Md. Roni Islam
Sheikh Md. Rabiul Islam
Develop blood oxygen level dependent signal by metabolic/hemodynamic model using numerical methods
Informatics in Medicine Unlocked
BOLD signal
MH model
Taylor series method
Euler method
Runge–Kutta method
title Develop blood oxygen level dependent signal by metabolic/hemodynamic model using numerical methods
title_full Develop blood oxygen level dependent signal by metabolic/hemodynamic model using numerical methods
title_fullStr Develop blood oxygen level dependent signal by metabolic/hemodynamic model using numerical methods
title_full_unstemmed Develop blood oxygen level dependent signal by metabolic/hemodynamic model using numerical methods
title_short Develop blood oxygen level dependent signal by metabolic/hemodynamic model using numerical methods
title_sort develop blood oxygen level dependent signal by metabolic hemodynamic model using numerical methods
topic BOLD signal
MH model
Taylor series method
Euler method
Runge–Kutta method
url http://www.sciencedirect.com/science/article/pii/S235291482300059X
work_keys_str_mv AT mdroniislam developbloodoxygenleveldependentsignalbymetabolichemodynamicmodelusingnumericalmethods
AT sheikhmdrabiulislam developbloodoxygenleveldependentsignalbymetabolichemodynamicmodelusingnumericalmethods