Processing the image gradient field using a topographic primal sketch approach.

The spatial derivatives of the image intensity provide topographic information that may be used to identify and segment objects. The accurate computation of the derivatives is often hampered in medical images by the presence of noise and a limited resolution. This paper focuses on accurate computation of spatial derivatives and their subsequent use to process an image gradient field directly, from which an image with improved characteristics can be reconstructed. The improvements include noise reduction, contrast enhancement, thinning object contours and the preservation of edges. Processing the gradient field directly instead of the image is shown to have numerous benefits. The approach is developed such that the steps are modular, allowing the overall method to be improved and possibly tailored to different applications. As presented, the approach relies on a topographic representation and primal sketch of an image. Comparisons with existing image processing methods on a synthetic image and different medical images show improved results and accuracy in segmentation. Here, the focus is on objects with low spatial resolution, which is often the case in medical images. The methods developed show the importance of improved accuracy in derivative calculation and the potential in processing the image gradient field directly. Copyright © 2015 John Wiley & Sons, Ltd.

Dynamics of airflow in a short inhalation.

During a rapid inhalation, such as a sniff, the flow in the airways accelerates and decays quickly. The consequences for flow development and convective transport of an inhaled gas were investigated in a subject geometry extending from the nose to the bronchi. The progress of flow transition and the advance of an inhaled non-absorbed gas were determined using highly resolved simulations of a sniff 0.5 s long, 1 l s⁻¹ peak flow, 364 ml inhaled volume. In the nose, the distribution of airflow evolved through three phases: (i) an initial transient of about 50 ms, roughly the filling time for a nasal volume, (ii) quasi-equilibrium over the majority of the inhalation, and (iii) a terminating phase. Flow transition commenced in the supraglottic region within 20 ms, resulting in large-amplitude fluctuations persisting throughout the inhalation; in the nose, fluctuations that arose nearer peak flow were of much reduced intensity and diminished in the flow decay phase. Measures of gas concentration showed non-uniform build-up and wash-out of the inhaled gas in the nose. At the carina, the form of the temporal concentration profile reflected both shear dispersion and airway filling defects owing to recirculation regions.

Sensitivity to outflow boundary conditions and level of geometry description for a cerebral aneurysm.

Mathematical models, namely the flow boundary conditions, as well as the detail of the bounding geometry, can highly influence the computed flow field. In this work, an anatomically realistic portion of cerebral vasculature with a saccular aneurysm, and its geometric idealisation, are considered. The importance of the geometric description, namely including the side branches or modelling them as holes in the main vessel, is studied. Several approaches to prescribe the outflow boundary conditions at the side branches are analysed, including the traction-free condition, zero velocity (hence neglecting the side-branch), and the coupling with simple zero-dimensional and one-dimensional models. Results of the effects of outflow boundary modelling choice on computed haemodynamic parameters are used to identify appropriateness of the models based on the physical interpretation. Estimated range of error-bars associated to outflow boundary model choice and the level of geometric details are presented for patient-specific computational haemodynamics, and can serve as invitation for future studies. The zero-dimensional and one-dimensional models are shown to provide good representations of the side branches in the case of the clipped geometry.

Decomposition and description of the nasal cavity form.

Patient-specific studies of physiological flows rely on anatomically realistic or idealized models. Objective comparison of datasets or the relation of specific to idealized geometries has largely been performed in an ad hoc manner. Here, two rational procedures (based respectively on Fourier descriptors and medial axis (MA) transforms) are presented; each provides a compact representation of a complex anatomical region, specifically the nasal airways. The techniques are extended to furnish average geometries. These retain a sensible anatomical form, facilitating the identification of a specific anatomy as a set of weighted perturbations about the average. Both representations enable a rapid translation of the surface description into a virtual model for computation of airflow, enabling future work to comprehensively investigate the relation between anatomic form and flow-associated function, for the airways or for other complex biological conduits. The methodology based on MA transforms is shown to allow flexible geometric modeling, as illustrated by a local alteration in airway patency. Computational simulations of steady inspiratory flow are used to explore the relation between the flow in individual vs. averaged anatomical geometries. Results show characteristic flow measures of the averaged geometries to be within the range obtained from the original three subjects, irrespective of averaging procedure. However the effective regularization of anatomic form resulting from the shape averaging was found to significantly reduce trans-nasal pressure loss and the mean shear stress in the cavity. It is suggested that this may have implications in attempts to relate model geometries and flow patterns that are broadly representative.

Nasal architecture: form and flow.

Current approaches to model nasal airflow are reviewed in this study, and new findings presented. These new results make use of improvements to computational and experimental techniques and resources, which now allow key dynamical features to be investigated, and offer rational procedures to relate variations in anatomical form. Specifically, both replica and simplified airways of a single subject were investigated and compared with the replica airways of two other individuals with overtly differing geometries. Procedures to characterize and compare complex nasal airway geometry are first outlined. It is then shown that coupled computational and experimental studies, capable of obtaining highly resolved data, reveal internal flow structures in both intrinsically steady and unsteady situations. The results presented demonstrate that the intimate relation between nasal form and flow can be explored in greater detail than hitherto possible. By outlining means to compare complex airway geometries and demonstrating the effects of rational geometric simplification on the flow structure, this work offers a fresh approach to studies of how natural conduits guide and control flow. The concepts and tools address issues that are thus generic to flow studies in other physiological systems.