From Wakapon
Jump to: navigation, search
Line 94: Line 94:
 
You don't change a winning team : ray-marching a volume texture is the key to realistic 3D clouds.
 
You don't change a winning team : ray-marching a volume texture is the key to realistic 3D clouds.
  
== First Tests ==
+
== Review ==
I first tried to ray-march a cube texture of 128x32x128. I had a 2-tiered process :
+
First, let's review the basics of light behavior within a participating medium. I could advise you more to read the great [http://graphics.ucsd.edu/~henrik/papers/book/ "Realistic Image Synthesis using Photon Mapping"] by Henrik Wann Jensen, to which I refer to as soon as I hit an obstacle or forget something (which is, I confess, most of the time).
 +
 
 +
 
 +
 
 +
== Light Diffusion Method ==
 +
I first tried to ray-march a cube texture of 128x32x128. I used a 2-tiered process :
 
* the first process was a "reaction-diffusion" of light in a RGBA16F texture
 
* the first process was a "reaction-diffusion" of light in a RGBA16F texture
 
* the second process was displaying the diffusion texture through ray-marching
 
* the second process was displaying the diffusion texture through ray-marching
  
=== Lighting ===
 
  
 
What is this light diffusion thing, you might ask ?
 
What is this light diffusion thing, you might ask ?
Line 122: Line 126:
 
** The new current cell's directional light component is equal to :
 
** The new current cell's directional light component is equal to :
 
  <math>RGB\prime_c = RGB_d + e^{\tau(\omega_d)}</math>
 
  <math>RGB\prime_c = RGB_d + e^{\tau(\omega_d)}</math>
where:
+
where:
 
  <math>\tau(\omega_d)=\sigma_t\frac{(\rho_d+\rho_c)}{2}\Delta_x</math>
 
  <math>\tau(\omega_d)=\sigma_t\frac{(\rho_d+\rho_c)}{2}\Delta_x</math>
 +
<math>\rho_d = RGB_d\cdot z</math> = cloud density one cell away in light's direction
 +
<math>\rho_c = RGB_c\cdot z</math> = cloud density of current cell
 +
<math>\Delta_x</math> = cell size (that we chose equal to 1 for simplicity)
 +
 +
This corresponds to the directional energy of one cell away toward the light, attenuated by extinction. Extinction being the combination of absorption (almost 0 in a cloud, as seen earlier) and scattering.
  
 +
<math></math>
  
 
[[Image:LightDiffusion.png|400px|thumb|center|The diffusion process within a bounded volume]]
 
[[Image:LightDiffusion.png|400px|thumb|center|The diffusion process within a bounded volume]]

Revision as of 15:47, 30 March 2011

This page is about the complex task of rendering clouds in real-time.

I have loved and studied clouds for a pretty long time now. I was taking pictures of sunsets from my bedroom's window when I was 15, I'm now 35 and nothing has changed : I'm still amazed by the variety and colors and shapes that cloud bring to an otherwise "sad" blue sky.

Cumulo-Nimbus

My favorites are storm clouds : the humongous cumulo-nimbus.

When it comes to rendering, you have 2 difficulties to overcome : lighting and shape. Both are difficult and very important.

Cloud Types

You may underestimate the importance of the shape of a cloud, especially because they look random so you might think a random shape is okay. But then you are plain wrong, because clouds are all but random.

Shape

If you watch a time-lapse video of a cloud, you can see how clouds move, how they are created and destroyed. One main observation to make is that accelerated clouds motion make them look like regular smoke. This is an important fractal pattern : clouds are like smoke, but bigger. And as they are orders of magnitude larger, the time scale is orders of magnitude smaller. It will take several minutes for the huge mass of a cloud to accomplish the same motion as smoke. A bit like the time for humans seems like slow motion as compared to insects.

You can see in the example video above how the cloud is "pulsed" from hot air pockets below. You can also see how there are "sources" and "sinks" that make the cloud appear and disappear.

All these obey the complex laws of thermodynamics : hot humid air rises up and evaporates while cold air falls down and liquefies, all these occurring at precise temperature and pressure, depending on the altitude.

Clouds are composed mainly of water in all its states : vapor, liquid (i.e. small droplets) and solid (i.e. snow and ice crystals). Clouds also need "nuclei" particles to form, like dust or aerosol molecules : despite large humidity quantities in the air, if there is no nuclei (i.e. the air is too clean), there won't be any cloud.

Various cloud types and their respective altitudes

Every Altitude its Cloud

Most of the clouds are very compact and fit in a "tiny" altitude range, although the storm cumulo-nimbi clouds span the entire range of altitudes, some of them going up to 15 kilometers (!!).

The scales in these clouds vary dramatically. Thin clouds are easily traversed by light while thick clouds are mostly reflective. Light scattering in clouds is really important because of the scales browsed by clouds : a photon doesn't go in a straight line for a long time when pacing through the many hundreds of meters inside a cloud.

Several ways for modeling clouds

Please, read the excellent Ken Perlin's presentation so you get familiar with the different Perlin noises.

Some are great for introducing regular randomness (plain noise), some others are good for clouds (sum 1/f(noise)) while even others are good for terrain modeling (sum 1/f(|noise|)) and are given funny names like "ridged multi fractal" or "fractal brownian motion".

Once you have your basic noise function, it's all a matter of intelligently combining several noises at several amplitudes and frequencies.


Alto/Cirro cumulus & stratus

These are thin vaporous high-altitude clouds. They are the easiest and most commonly simulated clouds as they are very thin and can be easily mapped to a plane, as shown on the image below.

For alto/cirro cumulus forms, there is a clear distinct "cauliflower" pattern that Perlin noise can easily achieve using the Turbulence pattern, as shown below :

Perlin noise : Turbulence variation.

For a great introduction to the complex light scattering within clouds, and the rendering of thin cloud slabs, I invite you to read the (quite technical) paper "Real-time realistic illumination and shading of stratiform clouds" by Bouthors et al.

You're okay sticking to a 2D cloud plane (even faking thickness through parallax mapping), you can even simulate up to the already thick stratocumulus until you need to display distinct clumps, holes or fly through the clouds. From this point on, you need to really render 3D clouds.


Nimbostratus, Cumulus & Cumulo-nimbus

These are the big guys (and my favorite S13.gif) !

I remember being baffled by the foundation paper "A Simple, Efficient Method for Realistic Animation of Clouds" by Dobashi et al. in 2000, where these guys managed to render impressively realistic clouds using particles while the clouds were modeled using a cellular automaton. This paper served as a base for Mark Harris and his well known clouds.

MarkHarris.png

This team of Japanese, when it comes to imitating reality, is really excellent ! They are also responsible for the foundation paper "Display Method of the Sky Color Taking into Account Multiple Scattering" about the sky model later implemented in real time by Sean O'Neil. On the downside though, like me, their english sucks and I remember re-reading 10 times the same sentence to try and understand what they meant.

Many tests through the years

Of course, I experimented with cloud particles myself at that time and I obtained some nice results but these clouds always kinda lacked the fluffiness I was after. This I learned : you can't really have great fluffy clouds with particles. First because of their disjoint nature (you're not using a mesh or tracing through a continuous volume here, merely splatting flat stuff together), then because of the lighting model that failed to account for important scattering events. Basically, with real-time particle models, you can take into account extinction and 0- or single-scattering events so you get the main features of a cloud, but you lack the general "denseness" and solid-feeling of cumuli and cumulo-nimbi clouds.


I tried many things, with more or less satisfying results :

  • Soft particles with precomputed Spherical Harmonics lighting, the lighting was quite okay but the clouds were obviously too "solid looking"
  • Marching cubes and other "mesh clouds" rendering, my intention was to render these in several depth buffers seen from front and back and to "subtract" the depths so I obtained the thickness of the mesh at a particular point in space (I first had this idea in 1995 and managed to have one of the first volumetric lights ever with my software renderer, I never used it in a demo though and the effect has since been re-used countless times). This could have worked but only computing the absolute thickness of an object, without accounting for holes and concavities always yielded poor results where it seemed light "poured through the volume" and although you got interesting god ray effects, the rendering was really not looking anything like a cloud.
  • Mega-particles, I first encountered this "technique" about 6 months ago (like november 2010) and, although I couldn't see how the effect worked (it looked like quite a "I bullshit you" effect), I decided to give it a try. This yielded interesting results but not because of the effect itself, as it's mainly a post-process and rendering a 3D cloud as if it was a 2D plane gives nasty "scrolling artefacts" when you move the camera. No, the main interesting rendering achievement here was because it was the first time I really correctly used phase functions AND a now quite old technique called "Deep Shadow Maps". This is really an interesting read !


For clouds modeling, I wrote a "recursive sphere aggregation" software where you drew your base cloud shape using big spheres, then the software created smaller "level 1" spheres on the surface of the these "root spheres" and simulated auto-spacing, then recursively created smaller "level 2" spheres and simulated again, and again. I managed to obtain nice cloud shapes with that software but I the rendering was awful. I never managed to have a nice rendering AND modeling at the same time...


Until finally...

It's only recently that I finally came up with the actual way of rendering nice volumetric clouds, some that you can change dynamically and view from anywhere, even from inside !

To do this, there is no other solution than rendering a 3D volume : a volume that tells you for each position in space whether or not there is a density of "cloud material". It's also very recent that it has become quite tractable to use 3D textures and even more recent that shaders allow you to perform ray-marching steps, the first examples of this being parallax mapping that is simply "terrain marching" to find the intersection of a ray with a height field.

My friend David was the first one (it's been a couple years already) to come up with super nice real-time clouds by ray-marching a volume, as seen in the picture below. His technique is a mix of mesh rendered into a volume texture, coupled with pre-computed spherical harmonics for lighting and various real-time 2D noise for shape variations, and that's brilliant !

U2Clouds.jpg


Rendering

You don't change a winning team : ray-marching a volume texture is the key to realistic 3D clouds.

Review

First, let's review the basics of light behavior within a participating medium. I could advise you more to read the great "Realistic Image Synthesis using Photon Mapping" by Henrik Wann Jensen, to which I refer to as soon as I hit an obstacle or forget something (which is, I confess, most of the time).


Light Diffusion Method

I first tried to ray-march a cube texture of 128x32x128. I used a 2-tiered process :

  • the first process was a "reaction-diffusion" of light in a RGBA16F texture
  • the second process was displaying the diffusion texture through ray-marching


What is this light diffusion thing, you might ask ?

Well I decided to split light into 2 parts : directional and isotropic.

If you read the paper by Bouthors et al. that I mentioned earlier about lighting stratiform clouds, you will read interesting facts about light scattering within clouds. The main fact is that clouds have a very high albedo of about 90%, meaning that they absorb almost no light. Most of the light is scattered. And among these scattering events, 50% are in the original light direction. This means there is a very strong forward scattering, or directional component in light going through a cloud.

The remaining 50% of the light is, grossly, isotropically scattered.

This is what I call the Watoo-Watoo phase function (from, an old cartoon I used to watch when I was a kid) :

WatooWatoo.png


You need to create a double-buffered RGBA16F 3D texture of, say, 128x32x128 pixels. Red will contain the directional light component while Green will contain the isotropic light component. Blue will contain the noise density and Alpha is not used. These buffers are initialized with 0 except the blue component that can be initialized with the initial noise density (although you can also use the shader to compute the noise in real-time).

The diffusion process then goes like this :

  • To propagate the directional component, we :
    • Read the current cell's RGB value (R=Directional G=Isotropic B=Density) => We call this RGB value <math>RGB_c</math>
    • Read the cell at P=CurrentPosition - LightDirection => We call that RGB value <math>RGB_d</math>
    • The new current cell's directional light component is equal to :
<math>RGB\prime_c = RGB_d + e^{\tau(\omega_d)}</math>
where:
<math>\tau(\omega_d)=\sigma_t\frac{(\rho_d+\rho_c)}{2}\Delta_x</math>
<math>\rho_d = RGB_d\cdot z</math> = cloud density one cell away in light's direction
<math>\rho_c = RGB_c\cdot z</math> = cloud density of current cell
<math>\Delta_x</math> = cell size (that we chose equal to 1 for simplicity)

This corresponds to the directional energy of one cell away toward the light, attenuated by extinction. Extinction being the combination of absorption (almost 0 in a cloud, as seen earlier) and scattering.

<math></math>
The diffusion process within a bounded volume


Ray-marching

The absolute key to a nice ray-marching is the amount of ray-marching steps. The more the better !

For a nice rendering I needed about 128 steps but you can't write a shader with a for loop of 128 steps for a 1024x768 surface so I exploited another great feature of clouds : their low frequency nature. The fact that clouds are large and only display small changes in shape made me guess about the impact of rendering not a 1024x768 surface but a downscaled version of 256x192 (downscaled by a factor 4). Ray-marching 128 steps in such a small surface doesn't change the impression that you are looking at clouds, and helps you to maintain a steady frame rate. On top of that, the downscale factor gives a very flexible way to scale the time used by the effect. For low end machines you can downscale by 8, for very high end machines you can downscale by 2. It's all a matter of experimenting with various hardware and use a scale factor appropriate for both the machine and the time you want to allocate to render clouds (as for most games, you can't spend your entire time on rendering).