# -*- coding: utf-8 -*-
################################ Begin license #################################
# Copyright (C) Laboratory of Imaging technologies,
# Faculty of Electrical Engineering,
# University of Ljubljana.
#
# This file is part of PyXOpto.
#
# PyXOpto is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# PyXOpto is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with PyXOpto. If not, see <https://www.gnu.org/licenses/>.
################################# End license ##################################
import os.path
import numpy as np
from .base import Absorption, Interpolator
from xopto import DATA_PATH
[docs]class DeOxyHem(Absorption):
material = 'deoxygenated hemoglobin'
def __init__(self):
super().__init__()
filename = os.path.join(
DATA_PATH, 'materials', 'absorption', 'blood_deoxy_absorption.npy')
self._interpolator = Interpolator.fromfile(filename)
def __call__(self, wavelength: float or np.ndarray, t: float) \
-> np.ndarray or float:
'''
Computes the absorption coefficient of deoxygenated blood.
Parameters
----------
wavelength: float, np.ndarray
Wavelength of light (m) at which to compute the absorption
coefficient of oxygenated blood.
t: float
Temperature not supported by this model.
Returns
-------
mua: float, np.ndarray
Absorption coefficient (1/m) at the specified wavelength (m).
'''
self._interpolator.check_range(wavelength)
return self._interpolator(wavelength)
[docs] def plot(self, wavelength: np.ndarray = None, show: bool = True):
'''
Plot the absorption coefficient at the given wavelengths.
Parameters
----------
wavelength: np.ndarray
Wavelengths of light. If None, use the reference values.
show: bool
Show the plot window if True.
'''
self._interpolator.plot(wavelength, label=self.material, show=show)
default = DeOxyHem()