This article was originally published in a journal published by Elsevier, and the attached copy is provided by Elsevier for the author’s benefit and for the benefit of the author’s institution, for non-commercial research and educational use including without limitation use in instruction at your institution, sending it to specific colleagues that you know, and providing a copy to your institution’s administrator. All other uses, reproduction and distribution, including without limitation commercial reprints, selling or licensing copies or access, or posting on open internet sites, your personal or institution’s website or repository, are prohibited. For exceptions, permission may be sought for such use through Elsevier’s permissions site at:
http://www.elsevier.com/locate/permissionusematerial ARTICLE IN PRESS
Computers & Graphics 30 (2006) 994–1010 www.elsevier.com/locate/cag Technical Section Simulation of atmospheric phenomena
Diego GutierrezÃ, Francisco J. Seron , Adolfo Munoz , Oscar Anson
Department of Computer Science and Systems Engineering, University of Zaragoza, C/Maria de Luna 1, 50018 Zaragoza, Spain
Received 19 January 2006; accepted 12 May 2006
Abstract
This paper presents a physically based simulation of atmospheric phenomena. It takes into account the physics of non-homogeneous media in which the index of refraction varies continuously, creating curved light paths. As opposed to previous research on this area, we solve the physically based differential equation that describes the trajectory of light. We develop an accurate expression of the index of refraction in the atmosphere as a function of wavelength, based on real measured data. We also describe our atmosphere profile manager, which lets us mimic the initial conditions of real-world scenes for our simulations. The method is validated both visually (by comparing the images with the real pictures) and numerically (with the extensive literature from other areas of research such as optics or meteorology). The phenomena simulated include the inferior and superior mirages, the Fata Morgana, the Novaya–Zemlya, the Viking’s end of the world, the distortions caused by heat waves and the green flash. r 2006 Elsevier Ltd. All rights reserved.
Keywords: Atmospheric phenomena; Nonlinear ray tracing; Inhomogeneous media
1. Introduction apply it here to simulate six different effects: the inferior and superior mirages, the Fata Morgana, the Novaya- Simulating natural phenomena has always been one of the Zemlya, the Viking’s end of the world and the green flash, most challenging goals of computer graphics. In particular, and propose several others, such as the distortions caused the effect of the atmosphere on sunlight, under certain by heat waves or fire. circumstances, causes a number of really spectacular phenomena; some are well-known, like the mirages, whereas 1.1. Overview some of them remain mostly unknown, such as the Fata Morgana or the Novaya-Zemlya. As far as the authors The atmosphere is a non-homogeneous medium, with its know some of the atmospheric phenomena presented in this properties varying continuously from point to point. The paper have never been simulated before, including the last most important property in terms of defining the light two mentioned above. This is due to the difficulty of first paths is the index of refraction, which in turn depends reproducing the exact, specific atmospheric conditions strongly on temperature, density and wavelength (and required to reproduce the effect and second, accurately more subtly on other parameters such as humidity). The calculating the nonlinear trajectories of light as it travels continuous variation of the light paths induced by the non- through the atmosphere under such conditions. homogeneous index of refraction can be described by a We present here our ray tracing method to reproduce differential equation, as will be shown in Section 3. some of these spectacular phenomena. The method is Our simulation is based on three main points: general and therefore can be applied to the simulation of any effect caused by nonlinear trajectories of light. We Fermat’s principle, which yields the physically based description of the light path through a non-homoge- ÃCorresponding author. Tel.: +34 976 762 354; fax: +34 976 761 914. E-mail addresses: [email protected] (D. Gutierrez), [email protected] neous medium (with a continuously varying index of (F.J. Seron), [email protected] (A. Munoz), [email protected] refraction), coupled with a fast, efficient numerical (O. Anson). resolution method.
0097-8493/$ - see front matter r 2006 Elsevier Ltd. All rights reserved. doi:10.1016/j.cag.2006.05.002 ARTICLE IN PRESS D. Gutierrez et al. / Computers & Graphics 30 (2006) 994–1010 995