Validation of 3D models of the outer and inner surfaces of an ovine femur

The validation of Computed Tomography (CT) based 3D models takes an integral part in studies involving 3D models of bones. This is of particular importance when such models are used for Finite Element studies. The validation of 3D models typically involves the generation of a reference model representing the bones outer surface. Several different devices have been utilised for digitising a bone’s outer surface such as mechanical 3D digitising arms, mechanical 3D contact scanners, electro-magnetic tracking devices and 3D laser scanners. However, none of these devices is capable of digitising a bone’s internal surfaces, such as the medullary canal of a long bone. Therefore, this study investigated the use of a 3D contact scanner, in conjunction with a microCT scanner, for generating a reference standard for validating the internal and external surfaces of a CT based 3D model of an ovine femur. One fresh ovine limb was scanned using a clinical CT scanner (Phillips, Brilliance 64) with a pixel size of 0.4 mm2 and slice spacing of 0.5 mm. Then the limb was dissected to obtain the soft tissue free bone while care was taken to protect the bone’s surface. A desktop mechanical 3D contact scanner (Roland DG Corporation, MDX 20, Japan) was used to digitise the surface of the denuded bone. The scanner was used with the resolution of 0.3 × 0.3 × 0.025 mm. The digitised surfaces were reconstructed into a 3D model using reverse engineering techniques in Rapidform (Inus Technology, Korea). After digitisation, the distal and proximal parts of the bone were removed such that the shaft could be scanned with a microCT (µCT40, Scanco Medical, Switzerland) scanner. The shaft, with the bone marrow removed, was immersed in water and scanned with a voxel size of 0.03 mm3. The bone contours were extracted from the image data utilising the Canny edge filter in Matlab (The Mathswork).. The extracted bone contours were reconstructed into 3D models using Amira 5.1 (Visage Imaging, Germany). The 3D models of the bone’s outer surface reconstructed from CT and microCT data were compared against the 3D model generated using the contact scanner. The 3D model of the inner canal reconstructed from the microCT data was compared against the 3D models reconstructed from the clinical CT scanner data. The disparity between the surface geometries of two models was calculated in Rapidform and recorded as average distance with standard deviation. The comparison of the 3D model of the whole bone generated from the clinical CT data with the reference model generated a mean error of 0.19±0.16 mm while the shaft was more accurate(0.08±0.06 mm) than the proximal (0.26±0.18 mm) and distal (0.22±0.16 mm) parts. The comparison between the outer 3D model generated from the microCT data and the contact scanner model generated a mean error of 0.10±0.03 mm indicating that the microCT generated models are sufficiently accurate for validation of 3D models generated from other methods. The comparison of the inner models generated from microCT data with that of clinical CT data generated an error of 0.09±0.07 mm Utilising a mechanical contact scanner in conjunction with a microCT scanner enabled to validate the outer surface of a CT based 3D model of an ovine femur as well as the surface of the model’s medullary canal.

Modeling corneal surfaces with three-dimensional basis functions

Purpose: To ascertain the effectiveness of object-centered three-dimensional representations for the modeling of corneal surfaces. Methods: Three-dimensional (3D) surface decomposition into series of basis functions including: (i) spherical harmonics, (ii) hemispherical harmonics, and (iii) 3D Zernike polynomials were considered and compared to the traditional viewer-centered representation of two-dimensional (2D) Zernike polynomial expansion for a range of retrospective videokeratoscopic height data from three clinical groups. The data were collected using the Medmont E300 videokeratoscope. The groups included 10 normal corneas with corneal astigmatism less than −0.75 D, 10 astigmatic corneas with corneal astigmatism between −1.07 D and 3.34 D (Mean = −1.83 D, SD = ±0.75 D), and 10 keratoconic corneas. Only data from the right eyes of the subjects were considered. Results: All object-centered decompositions led to significantly better fits to corneal surfaces (in terms of the RMS error values) than the corresponding 2D Zernike polynomial expansions with the same number of coefficients, for all considered corneal surfaces, corneal diameters (2, 4, 6, and 8 mm), and model orders (4th to 10th radial orders) The best results (smallest RMS fit error) were obtained with spherical harmonics decomposition which lead to about 22% reduction in the RMS fit error, as compared to the traditional 2D Zernike polynomials. Hemispherical harmonics and the 3D Zernike polynomials reduced the RMS fit error by about 15% and 12%, respectively. Larger reduction in RMS fit error was achieved for smaller corneral diameters and lower order fits. Conclusions: Object-centered 3D decompositions provide viable alternatives to traditional viewer-centered 2D Zernike polynomial expansion of a corneal surface. They achieve better fits to videokeratoscopic height data and could be particularly suited to the analysis of multiple corneal measurements, where there can be slight variations in the position of the cornea from one map acquisition to the next.

Asphalt concrete surfaces macrotexture determination from still images

Road surface macrotexture is identified as one of the factors contributing to the surface's skid resistance. Existing methods of quantifying the surface macrotexture, such as the sand patch test and the laser profilometer test, are either expensive or intrusive, requiring traffic control. High-resolution cameras have made it possible to acquire good quality images from roads for the automated analysis of texture depth. In this paper, a granulometric method based on image processing is proposed to estimate road surface texture coarseness distribution from their edge profiles. More than 1300 images were acquired from two different sites, extending to a total of 2.96 km. The images were acquired using camera orientations of 60 and 90 degrees. The road surface is modeled as a texture of particles, and the size distribution of these particles is obtained from chord lengths across edge boundaries. The mean size from each distribution is compared with the sensor measured texture depth obtained using a laser profilometer. By tuning the edge detector parameters, a coefficient of determination of up to R2 = 0.94 between the proposed method and the laser profilometer method was obtained. The high correlation is also confirmed by robust calibration parameters that enable the method to be used for unseen data after the method has been calibrated over road surface data with similar surface characteristics and under similar imaging conditions.

Estimation of illuminance on the south facing surfaces for clear skies in Iran

Background: Daylight availability data are essential for designing effectively day lighted buildings. In respect to no available daylight availability data in Iran, illuminance data on the south facing vertical surfaces were estimated using a proper method. Methods: An illuminance measuring set was designed for measuring vertical illuminances for standard times over 15 days at one hour intervals from 9 a.m. to 3 p.m. at three measuring stations (Hamadan, Eshtehard and Kerman). Measuring data were used to confirm predicted by the IESNA method. Results: Measurement of respective illuminances on the south vertical surfaces resulted in minimum values of 10.5 KLx, mean values of 33.59 KLx and maximum values of 79.6 KLx. Conclusion: In this study was developed a regression model between measured and calculated data of south facing vertical illuminance. This model, have a good linear correlation between measured and calculated values (r= 0.892).

Simulating droplet motion on virtual leaf surfaces

A curvilinear thin film model is used to simulate the motion of droplets on a virtual leaf surface, with a view to better understand the retention of agricultural sprays on plants. The governing model, adapted from Roy et al. (2002 J. Fluid Mech. 454, 235–261) with the addition of a disjoining pressure term, describes the gravity- and curvature driven flow of a small droplet on a complex substrate: a cotton leaf reconstructed from digitized scan data. Coalescence is the key mechanism behind spray coating of foliage, and our simulations demonstrate that various experimentally observed coalescence behaviours can be reproduced qualitatively. By varying the contact angle over the domain, we also demonstrate that the presence of a chemical defect can act as an obstacle to the droplet’s path, causing break-up. In simulations on the virtual leaf, it is found that the movement of a typical spray size droplet is driven almost exclusively by substrate curvature gradients. It is not until droplet mass is sufficiently increased via coalescence that gravity becomes the dominating force.

Miniscleral lens wear influences corneal curvature and optics

Purpose Transient changes in corneal topography associated with soft and conventional or reverse geometry rigid contact lens wear have been well documented; however, only a few studies have examined the influence of scleral contact lens wear upon the cornea. Therefore, in this study, we examined the influence of modern miniscleral contact lenses, which land entirely on the sclera and overlying tissues, upon anterior corneal curvature and optics. Methods Anterior corneal topography and elevation data were acquired using Scheimpflug imaging (Pentacam HR, Oculus) immediately prior to and following 8 hours of miniscleral contact lens wear in 15 young healthy adults (mean age 22 ± 3 years, 8 East Asian, 7 Caucasian) with normal corneae. Corneal diurnal variations were accounted for using data collected on a dedicated measurement day without contact lens wear. Corneal clearance was quantified using an optical coherence tomographer (RS-3000, Nidek) following lens insertion and after 8 hours of lens wear. Results Although corneal clearance was maintained throughout the 8 hour lens wear period, significant corneal flattening (up to 0.08 ± 0.04 mm) was observed, primarily in the superior mid-peripheral cornea, which resulted in a slight increase in against-the-rule corneal astigmatism (mean +0.02/-0.15 x 94 for an 8 mm diameter). Higher order aberration terms of horizontal coma, vertical coma and spherical aberration all underwent significant changes for an 8 mm corneal diameter (p ≤ 0.01), which typically resulted in a decrease in RMS error values (mean change in total higher order RMS -0.035 ± 0.046 µm for an 8 mm diameter). There was no association between the magnitude of change in central or mid-peripheral corneal clearance during lens wear and the observed changes in corneal curvature (p > 0.05). However, Asian participants displayed a significantly greater reduction in corneal clearance (p = 0.04) and greater superior-nasal corneal flattening compared to Caucasians (p = 0.048). Conclusions Miniscleral contact lenses that vault the cornea induce significant changes in anterior corneal surface topography and higher order aberrations following 8 hours of lens wear. The region of greatest corneal flattening was observed in the superior-nasal mid-periphery, more so in Asian participants. Practitioners should be aware that corneal measurements obtained following miniscleral lens removal may mask underlying corneal steepening.

Some surface profiles of a strip after plane-strain indentation by rigid bodies with serrated surfaces

This paper is concerned with the surface profiles of a strip after rigid bodies with serrated (saw-teeth) surfaces indent the strip and are subsequently removed. Plane-strain conditions are assumed. This has application in roughness transfer of final metal forming process. The effects of the semi-angle of the teeth, the depth of indentation and the friction on the contact surface on the profile are considered.