scholarly journals Hybrid Scene Structuring for Accelerating 3D Radiative Transfer Simulations

2019 ◽  
Vol 11 (22) ◽  
pp. 2637 ◽  
Author(s):  
Jianbo Qi ◽  
Tiangang Yin ◽  
Donghui Xie ◽  
Jean-Philippe Gastellu-Etchegorry

Three-dimensional (3D) radiative transfer models are the most accurate remote sensing models. However, presently the application of 3D models to heterogeneous Earth scenes is a computationally intensive task. A common approach to reduce computation time is abstracting the landscape elements into simpler geometries (e.g., ellipsoid), which, however, may introduce biases. Here, a hybrid scene structuring approach is proposed to accelerate the radiative transfer simulations while keeping the scene as realistic as possible. In a first step, a 3D description of the Earth landscape with equal-sized voxels is optimized to keep only non-empty voxels (i.e., voxels that contain triangles) and managed using a bounding volume hierarchy (BVH). For any voxel that contains triangles, within-voxel BVHs are created to accelerate the ray–triangle intersection tests. The hybrid scheme is implemented in the Discrete Anisotropic Radiative Transfer (DART) model by integrating the Embree ray-tracing kernels developed at Intel. In this paper, the performance of the hybrid algorithm is compared with the original uniform grid approach implemented in DART for a 3D city scene and a forest scene. Results show that the removal of empty voxels can accelerate urban simulation by 1.4×~3.7×, and that the within-voxel BVH can accelerate forest simulations by up to 258.5×.

2021 ◽  
Author(s):  
Dominik Göldner ◽  
Fotios Alexandros Karakostis ◽  
Armando Falcucci

This protocol presents the first detailed step-by-step pipeline for the 3D scanning and post processing of large batches of lithic artefacts using a micro-computed tomography (micro-CT) scanner (i.e., a Phoenix v-tome-x S model by General Electronics MCC, Boston MA) and an Artec Space Spider scanner (Artec Inc., Luxembourg). This protocol was used to scan and analyze ca. 700 lithic artefacts from the Protoaurignacian layers at Fumane Cave in north-eastern Italy (Falcucci et al., in preparation). For this study several costly scanners and proprietary software packages were employed. Although it is not easy to find a low-budget alternative for the scanners, it is possible to use free and open-source software programs, such as 3D-Slicer (https://www.slicer.org/) or MorphoDig (https://morphomuseum.com/morphodig), to process CT data as well as MeshLab (Cignoni et al. 2008) to interact with the 3D models in general. However, if alternative software is used, the steps and their order described in this protocol might diverge significantly. A cost-effective alternative to create 3D models is digital photogrammetry using commercial cameras and freely available software like Meshroom (https://alicevision.org). Although photogrammetry is an affordable technique to create accurate 3D models of objects, this method might not be useful when scanning large batches of artefacts, as it will require a lot of computation time and processing capacity. Likewise, it could be difficult to generate accurate 3D models of very small and/or detailed tool shapes using 3D surface scanners because stone tools are often much smaller than the recommended minimum field of view. Similarly, the resolution of conventional medical CT scanners might not be sufficient to capture minor details of stone tools, such as the outline or dorsal scars. Thus, high-resolution micro-CT technology is the only reliable way to accurately capture the overall morphology of small stone tools. This protocol aims at providing the first detailed procedure dedicated to the scanning of small lithic implements for further three-dimensional analysis. Note that some of the steps must be repeated at different working stages throughout this protocol. In cases where a task must be done in the exact same way as described in a previous step, a reference to that step is provided. When slight changes were made, the step was modified and reported entirely. This protocol contains a few red and green colours (e.g., arrows or within-program colours) which might be perceived differently by people with dyschromatopsia. However, the display of these colours has been kept to a minimum. We recommend the reader to go over the entire protocol carefully, even if only some specific parts are required. A few points are in fact interdependent, and some of them must be applied simultaneously. Content: Part 1 – Styrofoam preparation Part 2 – Micro-CT scanning Part 3 – 3D model extraction of CT scanned stone artifacts using Avizo Part 4 – Cropping extracted surface model to separate Face A and B in Artec Studio Part 5 – Cropping Face A to separate the lines in Artec Studio Part 6 – Cropping each stone artefact from the lines in Artec Studio Part 7 – Virtually control measurements in MeshLab Part 8 – Artec scanning of larger artifacts Part 9 – Export meshes as non-binary ply models for successive analysis in geomorph Three-dimensional example (in ply format) of the effectivity of the StyroStone Protocol: You can download an example of one Styrofoam line in 3D obtained using our protocol to appreciate the result that can be achieved. We have selected a line where objects are characterized by different metric and morphological attributes. Notice the retouching well visible in the last five smaller artifacts (counting from the left when artifact are oriented with the dorsal face in front of the observer and the butt down), as well as the platforms and bulbs of all artifacts. For more information and examples, feel free to contact us!


2009 ◽  
Vol 66 (10) ◽  
pp. 3131-3146 ◽  
Author(s):  
Robert Pincus ◽  
K. Franklin Evans

Abstract This paper examines the tradeoffs between computational cost and accuracy for two new state-of-the-art codes for computing three-dimensional radiative transfer: a community Monte Carlo model and a parallel implementation of the Spherical Harmonics Discrete Ordinate Method (SHDOM). Both codes are described and algorithmic choices are elaborated. Two prototype problems are considered: a domain filled with stratocumulus clouds and another containing scattered shallow cumulus, absorbing aerosols, and molecular scatterers. Calculations are performed for a range of resolutions and the relationships between accuracy and computational cost, measured by memory use and time to solution, are compared. Monte Carlo accuracy depends primarily on the number of trajectories used in the integration. Monte Carlo estimates of intensity are computationally expensive and may be subject to large sampling noise from highly peaked phase functions. This noise can be decreased using a range of variance reduction techniques, but these techniques can compromise the excellent agreement between the true error and estimates obtained from unbiased calculations. SHDOM accuracy is controlled by both spatial and angular resolution; different output fields are sensitive to different aspects of this resolution, so the optimum accuracy parameters depend on which quantities are desired as well as on the characteristics of the problem being solved. The accuracy of SHDOM must be assessed through convergence tests and all results from unconverged solutions may be biased. SHDOM is more efficient (i.e., has lower error for a given computational cost) than Monte Carlo when computing pixel-by-pixel upwelling fluxes in the cumulus scene, whereas Monte Carlo is more efficient in computing flux divergence and downwelling flux in the stratocumulus scene, especially at higher accuracies. The two models are comparable for downwelling flux and flux divergence in cumulus and upwelling flux in stratocumulus. SHDOM is substantially more efficient when computing pixel-by-pixel intensity in multiple directions; the models are comparable when computing domain-average intensities. In some cases memory use, rather than computation time, may limit the resolution of SHDOM calculations.


Author(s):  
Petr Chelishchev ◽  
Knut Sørby

This paper presents algorithms for estimating the minimum volume bounding box based on a three-dimensional point set measured by a coordinate measuring machine. A new algorithm, which calculates the minimum volume with high accuracy and reduced number of computations, is developed. The algorithm is based on the convex hull operation and established theories about a minimum bounding box circumscribing a convex polyhedron. The new algorithm includes a pre-processing operation that removes convex polyhedron faces located near the edges of the measured object. As showed in the paper, the solution of the minimum bonding box is not based on faces located near the edges; therefore, we can save computation time by excluding them from the convex polyhedron data set. The algorithms have been demonstrated on physical objects measured by a coordinate measuring machine, and on theoretical 3D models. The results show that the algorithm can be used when high accuracy is required, for example in calibration of reference standards.


1997 ◽  
Vol 3 (S2) ◽  
pp. 275-276
Author(s):  
Mark H. Ellisman ◽  
Stephen J. Young ◽  
G. Y. Fan ◽  
Guy Perkins ◽  
Steve Lamont ◽  
...  

The intermediate high-voltage electron microscope (IVEM) located at the National Center for Microscopy and Imaging Research at San Diego (NCMIR) can image relatively thick specimens that contain substantial three-dimensional (3-D) structure. Electron tomography is an important tool used at NCMIR for deriving 3-D cellular and subcellular structure from IVEM images. Reconstruction algorithms commonly used in electron tomography include weighted back projection, and iterative algebraic reconstruction techniques such as ART and SIRT. Improvements in reconstruction quality are possible using the iterative algorithms. Because these algorithms are computationally intensive, we have ported them to massively parallel computers at the San Diego Supercomputer Center, reducing the computation time over that required with workstation level machines.The quality of tomographic data for the 3-D reconstruction of biological structures is also being enhanced by NCMIR projects to improve the microscope. We have designed and constructed special electron optics and microscope control systems for the JEOL 4000EX.


Author(s):  
H. David Miranda ◽  
John Orr ◽  
Chris Williams

Numerical implementations of bond-based peridynamics are computationally intensive. We propose a new class of fast interaction functions for constitutive modelling that reduce calculation time when compared to other formulations in the literature. This is achieved by substituting the stretch definition from the original interaction functions with a new stretch measure that we call modified stretch. The resultant interaction functions are proven to approximate the existing formulations, and proven to require equivalent stability and convergence conditions under explicit time integration. Gains of speed greater than 11% were obtained in numerical tests that compared the new functions with those in the literature. The new approach was verified against classical elastic theory using simple examples and shows good agreement. Examples describing three-dimensional quasi-brittle structures are also presented. The proposed fast interaction functions lead to improvements in the ability to calculate the load response of realistic structures, since they usually require fine discretisation and large computation time.


2013 ◽  
Vol 819 ◽  
pp. 299-303
Author(s):  
Zhi Yu Chen ◽  
Jian Zhong Fu ◽  
Hong Yao Shen ◽  
Wen Feng Gan

Amongst the flourishing Delaunay Triangulation methods, growth algorithm has been widely accepted because of its reputation of being simple and elegant. However, the parallelization of growth algorithm has not been fully exploited. In this work, a novel Growth algorithm of Delaunay Triangulation is proposed. The point cloud is first divided into two parts by a suitable curve and the separated areas are calculated by incremental algorithm. Triangles which cross with the curve are generated by a growth algorithm associated with uniform grid. At the process of merging, these grew triangles are used to detect incorrect triangles of the incremental algorithm areas. Method about generating triangles on curve is elaborated and a simple way to detect interferential triangles is also explained. With above method, triangulation calculation can be parallelized. Unlike the traditional divide-and-conquer method, no flip operation is needed in the proposed methodology. Thus, three dimensional applications are also made possible. A comparative research between tradition incremental algorithm and the proposed method has been conducted. Results show, the algorithm has a higher performance with less computation time.


2016 ◽  
Vol 73 (10) ◽  
pp. 4137-4149 ◽  
Author(s):  
Laurent Makké ◽  
Luc Musson-Genon ◽  
Bertrand Carissimo ◽  
Pierre Plion ◽  
Maya Milliez ◽  
...  

Abstract The atmospheric radiation field has seen the development of more accurate and faster methods to take into account absorption. Modeling fog formation, where infrared radiation is involved, requires accurate methods to compute cooling rates. Radiative fog appears under clear-sky conditions owing to a significant cooling during the night where absorption and emission are the dominant processes. Thanks to high-performance computing, high-resolution multispectral approaches to solving the radiative transfer equation are often used. Nevertheless, the coupling of three-dimensional radiative transfer with fluid dynamics is very computationally expensive. Radiation increases the computation time by around 50% over the pure computational fluid dynamics simulation. To reduce the time spent in radiation calculations, a new method using analytical absorption functions fitted by Sasamori on Yamamoto’s radiation chart has been developed to compute an equivalent absorption coefficient (spectrally integrated). Only one solution of the radiative transfer equation is needed against Nband × Ngauss for an Nband model with Ngauss quadrature points on each band. A comparison with simulation data has been done and the new parameterization of radiative properties proposed in this article shows the ability to handle variations of gas concentrations and liquid water.


2008 ◽  
Author(s):  
Hamzah Asyrani Sulaiman ◽  
Abdullah Bade ◽  
Daut Daman ◽  
Mohd Shahrizal Sunar

Sign in / Sign up

Export Citation Format

Share Document