Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Models, Anatomic”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61Linked to original sources

Localized specific absorption rate calculations in a realistic phantom leg at 1-30 MHz using a finite element method.

Protection standards for radiofrequency electromagnetic radiation are principally intended to avoid detrimental thermal effects. To this end the International Commission on Non-Ionizing Radiation Protection (ICNIRP), and national bodies such as the National Radiological Protection Board (NRPB), recommend limitations on the localized specific energy absorption rate (SAR) in various parts of the body. The role of numerical dosimetry is to estimate the SAR from measurable parameters such as external field strengths and total body currents. In recent years there have been significant advances in the sophistication of the anatomical models available, and in our knowledge of the electrical properties of the body tissues. Several groups, including NRPB, have developed mathematical phantoms from medical imaging data, such as MRI scans. It has been known for some time that under certain circumstances SAR restrictions may be violated in the ankle due to the concentration of current in a small area. In this paper the author presents calculations of the SAR distribution in a human leg in the high-frequency (HF) band. This band contains the human whole-body resonance frequency and therefore gives the strongest coupling of the body to the field. The present study uses a finite element model with variable mesh size, derived from a 2 mm resolution voxel phantom of the whole body. It also uses recently acquired data on the electrical properties of the tissues. The results are discussed in the light of the exposure standards promulgated by national and international bodies such as NRPB and ICNIRP, and it is shown that the basic SAR restrictions in the leg are ensured by a current reference level of 100 mA.

Absorptiometry, Photon↗

An EGS4-ready tomographic computational model of a 14-year-old female torso for calculating organ doses from CT examinations.

Fifty-four consecutive CT scans have been used to construct a tomographic computational model of a 14-year-old female torso suitable for the determination of organ doses from CT. The model, known as ADELAIDE, is in the form of an input file compatible with user codes based on XYZDOS.MOR from the readily available EGS4 Monte Carlo radiation transport code. ADELAIDE's dimensions are close to the Australian averages for her age so the model is representative of a 14-year-old girl. The realistic anatomy in the model differs considerably from that in Cristy's 15-year-old mathematical computational model by having realistically shaped organs that are appropriately located within a real external contour. Average absorbed dose to organs from simulated CT examinations of the chest and abdomen have been calculated for ADELAIDE using EGS4 within a geometry specific to the General Electric Hi-Speed Advantage CT scanner and using an x-ray spectrum calculated using data from the scanner's x-ray tube. The simulations include the scanner's beam shaping filter and patient table. It is suggested that the resulting values have fewer possible sources of uncertainty than organ doses derived from dose coefficients calculated for a MIRD style model with mathematical anatomy and a spectrum that may not match that of the scanner. The organ doses were normalized using the scanner's CTDI measured free-in-air and an EGS4 simulation of the CTDI measurement. Effective dose to the torso from 26-slice chest and 24-slice abdomen examinations (at 120 kV, 200 mAs, 7 mm slices) is 4.6 +/- 0.1 mSv and 4.3 +/- 0.1 mSv respectively.

Adolescent↗

EEG dipole source localization using artificial neural networks.

Localization of focal electrical activity in the brain using dipole source analysis of the electroencephalogram (EEG), is usually performed by iteratively determining the location and orientation of the dipole source, until optimal correspondence is reached between the dipole source and the measured potential distribution on the head. In this paper, we investigate the use of feed-forward layered artificial neural networks (ANNs) to replace the iterative localization procedure, in order to decrease the calculation time. The localization accuracy of the ANN approach is studied within spherical and realistic head models. Additionally, we investigate the robustness of both the iterative and the ANN approach by observing the influence on the localization error of both noise in the scalp potentials and scalp electrode mislocalizations. Finally, after choosing the ANN structure and size that provides a good trade off between low localization errors and short computation times, we compare the calculation times involved with both the iterative and ANN methods. An average localization error of about 3.5 mm is obtained for both spherical and realistic head models. Moreover, the ANN localization approach appears to be robust to noise and electrode mislocations. In comparison with the iterative localization, the ANN provides a major speed-up of dipole source localization. We conclude that an artificial neural network is a very suitable alternative for iterative dipole source localization in applications where large numbers of dipole localizations have to be performed, provided that an increase of the localization errors by a few millimetres is acceptable.

Brain↗

Thermal effects of radiation from cellular telephones.

A finite element thermal model of the head has been developed to calculate temperature rises generated in the brain by radiation from cellular telephones and similar electromagnetic devices. A 1 mm resolution MRI dataset was segmented semiautomatically, assigning each volume element to one of ten tissue types. A finite element mesh was then generated using a fully automatic tetrahedral mesh generator developed at NRPB. There are two sources of heat in the model: firstly the natural metabolic heat production; and secondly the power absorbed from the electromagnetic field. The SAR was derived from a finite difference time domain model of the head, coupled to a model 'mobile phone', namely a quarter-wavelength antenna mounted on a metal box. The steady-state temperature distribution was calculated using the standard Pennes 'bioheat equation'. In the normal cerebral cortex the high blood perfusion rate serves to provide an efficient cooling mechanism. In the case of equipment generally available to the public, the maximum temperature rise found in the brain was about 0.1 degrees C. These results will help in the further development of criteria for exposure guidelines, and the technique developed may be used to assess temperature rises associated with SARs for different types of RF exposure.

Absorption↗

Tissue resistivity estimation in the presence of positional and geometrical uncertainties.

Geometrical uncertainties (organ boundary variation and electrode position uncertainties) are the biggest sources of error in estimating electrical resistivity of tissues from body surface measurements. In this study, in order to decrease estimation errors, the statistically constrained minimum mean squared error estimation algorithm (MiMSEE) is constrained with a priori knowledge of the geometrical uncertainties in addition to the constraints based on geometry, resistivity range, linearization and instrumentation errors. The MiMSEE calculates an optimum inverse matrix, which maps the surface measurements to the unknown resistivity distribution. The required data are obtained from four-electrode impedance measurements, similar to injected-current electrical impedance tomography (EIT). In this study, the surface measurements are simulated by using a numerical thorax model. The data are perturbed with additive instrumentation noise. Simulated surface measurements are then used to estimate the tissue resistivities by using the proposed algorithm. The results are compared with the results of conventional least squares error estimator (LSEE). Depending on the region, the MiMSEE yields an estimation error between 0.42% and 31.3% compared with 7.12% to 2010% for the LSEE. It is shown that the MiMSEE is quite robust even in the case of geometrical uncertainties.

Algorithms↗

A comparison of quantitative measures of arterial tortuosity using sine wave simulations and 3D wire models.

This paper discusses the quantitative measurement of arterial tortuosity. The measures described are a tortuosity coefficient, TC, based on the second derivative of the vessel midline and the distance factor, DF, defined as L/D - 1, where L is the vessel length and D is the straight line distance between the vessel end points. A tortuosity scale derived from sine wave simulations was used to determine the ability of these two measures to quantify tortuosity. The tortuosity coefficient was demonstrated to be an appropriate measure whilst the distance factor was found to be a measure of vessel elongation rather than tortuosity. A tomographic method of deriving three-dimensional measures of TC and DF is described using a wire model. The x, y coordinates along the length of the wire were obtained from CT slices and the three-dimensional values calculated as the geometric mean of the TC and DF values measured in the x and y directions. Rotation of the wire through 90 degrees demonstrated that the three-dimensional derivations of TC and DF were independent of the orientation of the wire.

Arteries↗

Voxel effects within digital images of trabecular bone and their consequences on chord-length distribution measurements.

Chord-length distributions through the trabecular regions of the skeleton have been investigated since the early 1960s. These distributions have become important features for bone marrow dosimetry; as such, current models rely on the accuracy of their measurements. Recent techniques utilize nuclear magnetic resonance (NMR) microscopy to acquire 3D images of trabecular bone that are then used to measure 3D chord-length distributions by Monte Carlo methods. Previous studies have shown that two voxel effects largely affect the acquisition of these distributions within digital images. One is particularly pertinent as it dramatically changes the shape of the distribution and reduces its mean. An attempt was made to reduce this undesirable effect and good results were obtained for a single-sphere model using minimum acceptable chord (MAC) methods (Jokisch et al 2001 Med. Phys. 28 1493-504). The goal of the present work is to extend the study of these methods to more general models in order to better quantify their consequences. First, a mathematical model of a trabecular bone sample was used to test the usefulness of the MAC methods. The results showed that these methods were not efficient for this simulated bone model. These methods were further tested on a single voxelized sphere over a large range of voxel sizes. The results showed that the MAC methods are voxel-size dependent and overestimate the mean chord length for typical resolutions used with NMR microscopy. The study further suggests that bone and marrow chord-length distributions currently utilized in skeletal dosimetry models are most likely affected by voxel effects that yield values of mean chord length lower than their true values.

Bone Marrow↗

Dynamics and interaction of filaments in a computational model of re-entrant ventricular fibrillation.

Ventricular fibrillation (VF) is a lethal cardiac arrhythmia. Re-entry, in which action potential wavefronts rotate around filaments, is believed to sustain VF. In this study we used a computational model of multiple wavelet fibrillation in the thin-walled right ventricle (10 mm thick) and the thicker walled left ventricle (16 mm thick) to investigate the effect of tissue thickness and initiation protocol on re-entry, and to examine whether filament dynamics and interaction in the model could explain why re-entry is both rarely observed and short-lived in experimental studies that map electrical activation on the heart surface. We found (i) that the density of filaments, the proportion of transmural filaments and the proportion of filaments visible on the model surface were all higher in the 10 mm simulation, (ii) that the initiation protocol influences the rate of filament breakdown but not the number of filaments present after 1 s, and (iii) that although many filaments are visible on the surface of the model, the majority are visible for less than one rotation. This study shows that tissue thickness, geometry and initiation protocol influence electrical activation during VF, and that the rapid motion and interaction of filaments result in transient appearance of surface re-entry.

Anisotropy↗

Modified ionic models of cardiac tissue for efficient large scale computations.

Recirculation of excitation, or re-entry, is one of the most important mechanisms of life-threatening cardiac arrhythmias and fibrillation. Modelling these phenomena requires large scale computations in two and three-dimensional slabs of cardiac tissue. Because of computational constraints, most of the studies use simplified (non-ionic) models of cardiac tissue, which are electrophysiologically less accurate than the detailed ionic models. In this paper, we propose a method to modify ionic models of cardiac tissue into an intermediate class of models, which are almost as efficient for computations as simplified models, and retain most of the properties of the original ionic models, such as the shape of the action potential, the restitution of action potential duration and of the conduction velocity, as well as unchanged description of most of the ionic currents.

Algorithms↗

In vivo metabolism and partitioning of 6-[18F]fluoro-L-meta-tyrosine in whole blood: a unified compartment model.

Physiological quantification of dynamic PET data requires the determination of an input function, preferably from plasma. A compartmental model relating a parent radiotracer, its radiolabelled metabolites and their exchange between plasma and erythrocytes is presented. This model allows for the time course of radioactivity measured in whole blood to be transformed into the time course of the radiotracer in plasma. The utility of this approach is illustrated with blood data collected on 30 human subjects injected with 6-[18F]fluoro-L-meta-tyrosine (FmT), a pre-synaptic dopaminergic radiotracer. A three-compartment four-parameter model is shown to yield significantly better fits to the blood data than related lower and higher order models. This model is found to be robust to measurement noise, and yet sensitive to metabolic changes induced by pretreatment with carbidopa. For FmT, the between-subject variations are shown to be small enough to warrant the use of a population-based correction; tissue time-activity curves were simulated to verify that this correction does not significantly affect the precision and accuracy of the derived rate constants. The unified blood model can be adapted for radiotracers other than FmT as long as the blood partition ratio of the parent radiotracer differs from that of its metabolites and/or the rate at which they equilibrate between plasma and erythrocytes is different.

Adult↗

Precise modelling of the eye for proton therapy of intra-ocular tumours.

A new method is described that allows precise modelling of organs at risk and target volume for radiation therapy of intra-ocular tumours. The aim is to optimize the dose distribution and thus to reduce normal tissue complication probability. A geometrical 3D model based on elliptic shapes was developed that can be used for multimodal model-based segmentation of 3D patient data. The tumour volume cannot be clearly identified in CT and MR data, whereas the tumour outline can be discriminated very precisely in fundus photographs. Therefore, a multimodal 2D fundus diagram was developed, which allows us to correlate and display simultaneously information extracted from the eye model, 3D data and the fundus photograph. Thus, the connection of fundus diagram and 3D data is well-defined and the 3D volume can be calculated directly from the tumour outline drawn onto the fundus photograph and the tumour height measured by ultrasound. The method allows the calculation of a precise 3D eye model of the patient, including the different structures of the eye as well as the tumour volume. The method was developed as part of the new 3D treatment planning system OCTOPUS for proton therapy of ocular tumours within a national research project together with the Hahn-Meitner-Institut Berlin.

Eye Neoplasms↗

The effect of volume currents due to myocardial anisotropy on body surface potentials.

Changes in anterior and posterior body surface potential maps (BSPMs) due to myocardial anisotropy were examined using a highly heterogeneous finite element model of an adult male subject constructed from segmented magnetic resonance images. A total of 23 different tissue types were identified in the whole torso. The myocardial fibre orientations in the human heart wall were mapped from the fibre orientations of a canine heart which are available in the literature using deformable mapping techniques. The current and potential distributions in the whole torso were computed using dipolar sources in the septum, apical area, left ventricular wall or right ventricular wall. For each dipole x, y, z orientations were studied. An adaptive finite element solver was used to compute currents and potential distributions in the whole torso with an element size of 0.78 x 0.78 x 3 mm in the myocardium and larger elements in other parts of the torso. For each dipole position two cases were studied. In one case the myocardium was isotropic and in the other it was anisotropic. It was found that BSPMs showed a very notable difference between the isotropic and the anisotropic myocardium for all dipole positions with the largest difference for the apical dipoles. The correlation coefficients for the BSPMs between the isotropic and anisotropic cases ranged from 0.83 for an apical dipole to 0.99 for an RV wall dipole. These results suggest that myocardial fibre anisotropy plays an important role in determining the body surface potentials.

Adult↗

Optimizing the assessment of age-related changes in trabecular bone.

The goal of this study was to develop an optimal procedure to determine age-related changes in trabecular bone. The investigations were based on two-dimensional images of the human vertebral trabecular bone specimens. The following indices of trabecular structure were considered: bone volume/total volume, star volume of the marrow cavity, Euler number and the probability of disconnection (straightforwardly connected with the number of separated parts of the network). To follow precisely the changes in the trabecular structure with age, a computer simulation model was used. Up to 35 years of physiological remodelling were simulated. The validation of the model calculations was based on a quantitative comparison with the data measured for older individuals. The simulations confirmed that the description of the age-related changes in the trabecular bone by means of the architectural parameter (star volume) constitutes a promising tool for subjects older than approximately 50 years. For individuals younger than approximately 50 years bone mineral density (bone volume/total volume) seems to be the best suited descriptor. The results suggest that the optimal diagnostic procedure is age-dependent and should not be limited to the bone mineral density measurement. The clinical usefulness of the procedure has been validated by examination of the CT images.

Adult↗

Effect of a hands-free wire on specific absorption rate for a waist-mounted 1.8 GHz cellular telephone handset.

A common feature of cellular telephony is the use of a 'hands-free' audio extension lead connected to a waist-worn handset. Interaction between the transmitting antenna, the wire and the user's body can occur, with detrimental effects including polar pattern degradation, reduced efficiency and localized increases in specific absorption rate (SAR). Using a realistic full-body model of an adult male, finite difference time domain analysis was employed to investigate the coupling between a hip-mounted 1.8 GHz handset fitted with a monopole antenna and a 1 m long wire representing a hands-free wire. Conduction current densities were computed for three identifiable coupling modes: magnetic-only, conductive-only and combined conductive-and-magnetic. Magnetic-only coupling was dominant. Without the lead, placing the handset at waist height led to a 42.8% increase in the total energy deposited in the body, compared to use at the head. Introducing the lead further increased the body loss, with a reduction in system radiation efficiency from 52% to 43.7%. Without the hands-free lead, the peak 1 g and 10 g SARs were 0.450 W kg(-1) and 0.265 W kg(-1), respectively, for 125 mW transmit power. With the hands-free lead connected, these values increased to 1.14 W kg(-1) and 0.430 W kg(-1), respectively.

Adult↗

A semi-automatic method for developing an anthropomorphic numerical model of dielectric anatomy by MRI.

Complex permittivity values have a dominant role in the overall consideration of interaction between radiofrequency electromagnetic fields and living matter, and in related applications such as electromagnetic dosimetry. There are still some concerns about the accuracy of published data and about their variability due to the heterogeneous nature of biological tissues. The aim of this study is to provide an alternative semi-automatic method by which numerical dielectric human models for dosimetric studies can be obtained. Magnetic resonance imaging (MRI) tomography was used to acquire images. A new technique was employed to correct nonuniformities in the images and frequency-dependent transfer functions to correlate image intensity with complex permittivity were used. The proposed method provides frequency-dependent models in which permittivity and conductivity vary with continuity--even in the same tissue--reflecting the intrinsic realistic spatial dispersion of such parameters. The human model is tested with an FDTD (finite difference time domain) algorithm at different frequencies; the results of layer-averaged and whole-body-averaged SAR (specific absorption rate) are compared with published work, and reasonable agreement has been found. Due to the short time needed to obtain a whole body model, this semi-automatic method may be suitable for efficient study of various conditions that can determine large differences in the SAR distribution, such as body shape, posture, fat-to-muscle ratio, height and weight.

Algorithms↗

Transit time kinetics in ordered and disordered vascular trees.

Imaging modalities exploit tracer-dilution methods to measure bulk haemodynamic parameters such as blood flow and volume at the level of the microcirculation. Here, we ask the question of whether the kinetics of a tracer can reveal morphological information about the vessels through which the tracers flow. The goal is to relate the acquired time-intensity characteristic to details of the vascular structure that lies below the imaging resolution. Two fractal vascular models are developed that represent organized 'kidney-like' and disorganized 'tumour-like' structures. The models are generated using simple rules of branching and fractal geometry in two dimensions. Blood flow and tracer kinetics are simulated using fundamental laws of haemodynamics. The flow conditions are matched in the two models. The fractal box dimensions of the kidney (D(B) = 1.67 +/- 0.01) and the tumour (D(B) = 1.80 +/- 0.01) vasculatures fall in the range given in the literature (D(B) = 1.61 +/- 0.06 and D(B) = 1.84 +/- 0.04, respectively). The tracer kinetic curves of the kidney and the tumour vasculatures have the same initial slope and final asymptote, corresponding to the same flow rate and vascular volume, but have different forms. The difference in the two curves is related to the distribution function of transit times of the vascular models, and is a consequence of the randomness introduced in vessel diameter and length. In principle, the form of the tracer kinetic curve from a contrast imaging study may offer information relating not only to vascular volume and flow rate, but also to the organization of a microvascular network.

Algorithms↗

Physical mechanisms of gas and perfluoron retinopexy and sub-retinal fluid displacement.

Injection of gas into the eye, followed by face-down positioning, is a common protocol for the reseating of the retina in posterior and superior retinal tears and breaks. The physical mechanism by which injected gas helps reattach retinal flaps is often ascribed to the 'buoyancy' force of the injected gas bubble. The various forces at play in this system (surface tension and buoyancy) were calculated and compared. The results are extended to the case in which the retina is intact (pneumatic displacement of blood) and to the use of intraocular perfluoron (n-perfluorooctane). We show that buoyancy forces are applicable only for gas or n-perfluorooctane bubbles that are smaller than the detached retina and that do not invade underneath the retina. For larger bubbles, as is normally used in reattachment protocols, we show that it is the interfacial tension that reattaches the retina. The range of angles within which patients can position, and still maintain a gas-vitreous interface along a tear is calculated as a function of the volume of injected gas and size of the tear. The maximum retinal flap size that can be reattached using surface tension forces is also estimated.

Fluorocarbons↗

Registration of electronic portal images for patient set-up verification.

Images acquired from an electronic portal imaging device are aligned with digitally reconstructed radiographs (DRRs) or other portal images to verify patient positioning during radiation therapy. Most of the currently available computer aided registration methods are based on the manual placement of corresponding landmarks. The purpose of the paper is twofold: (a) the establishment of a methodology for patient set-up verification during radiotherapy based on the registration of electronic portal images, and (b) the evaluation of the proposed methodology in a clinical environment. The estimation of set-up errors, using the proposed methodology, can be accomplished by matching the portal image of the current fraction of the treatment with the portal image of the baseline treatment (reference portal image) using a nearly automated technique. The proposed registration method is tested on a number of phantom data as well as on data from four patients. The phantom data included portal images that corresponded to various positions of the phantom on the treatment couch. For each patient, a set of 30 portal images was used. For the phantom data (for both transverse and lateral portal images), the maximum absolute deviations of the translational shifts were within 1.5 mm, whereas the in-plane rotation angle error was less than 0.5 degrees. The two-way Anova revealed no statistical significant variability both within observer and between-observer measurements (P > 0.05). For the patient data, the mean values obtained with manual and the proposed registration methods were within 0.5 mm. In conclusion, the proposed registration method has been incorporated within a system, called ESTERR-PRO. Its image registration capability achieves high accuracy and both intra- and inter-user reproducibility. The system is fully operational within the Radiotherapy Department of 'HYGEIA' Hospital in Athens and it could be easily installed in any other clinical environment since it requires standardized hardware specifications and minimal human intervention.

Algorithms↗