Michaël Lemaire
4a6d9bb001
git-svn-id: https://subversion.assembla.com/svn/thunderk/paysages@356 b1fd45b6-86a6-48da-8261-f70d1f35bdcc
462 lines
15 KiB
C
462 lines
15 KiB
C
#include "sky.h"
|
|
|
|
#include <stdlib.h>
|
|
#include <math.h>
|
|
|
|
#include "shared/types.h"
|
|
#include "color.h"
|
|
#include "clouds.h"
|
|
#include "euclid.h"
|
|
#include "lighting.h"
|
|
#include "render.h"
|
|
#include "tools.h"
|
|
|
|
#define SPHERE_SIZE 1000.0
|
|
|
|
/******************************** Configuration ********************************/
|
|
void skyInit()
|
|
{
|
|
}
|
|
|
|
void skyQuit()
|
|
{
|
|
}
|
|
|
|
void skySave(PackStream* stream, SkyDefinition* definition)
|
|
{
|
|
packWriteInt(stream, (int*)&definition->model);
|
|
packWriteDouble(stream, &definition->daytime);
|
|
colorSave(stream, &definition->sun_color);
|
|
packWriteDouble(stream, &definition->sun_radius);
|
|
packWriteDouble(stream, &definition->sun_halo_size);
|
|
curveSave(stream, definition->sun_halo_profile);
|
|
|
|
packWriteInt(stream, &definition->model_custom.auto_from_daytime);
|
|
colorSave(stream, &definition->model_custom.zenith_color);
|
|
colorSave(stream, &definition->model_custom.haze_color);
|
|
packWriteDouble(stream, &definition->model_custom.haze_height);
|
|
packWriteDouble(stream, &definition->model_custom.haze_smoothing);
|
|
|
|
packWriteDouble(stream, &definition->model_preetham.turbidity);
|
|
}
|
|
|
|
void skyLoad(PackStream* stream, SkyDefinition* definition)
|
|
{
|
|
packReadInt(stream, (int*)&definition->model);
|
|
packReadDouble(stream, &definition->daytime);
|
|
colorLoad(stream, &definition->sun_color);
|
|
packReadDouble(stream, &definition->sun_radius);
|
|
packReadDouble(stream, &definition->sun_halo_size);
|
|
curveLoad(stream, definition->sun_halo_profile);
|
|
|
|
packReadInt(stream, &definition->model_custom.auto_from_daytime);
|
|
colorLoad(stream, &definition->model_custom.zenith_color);
|
|
colorLoad(stream, &definition->model_custom.haze_color);
|
|
packReadDouble(stream, &definition->model_custom.haze_height);
|
|
packReadDouble(stream, &definition->model_custom.haze_smoothing);
|
|
|
|
packReadDouble(stream, &definition->model_preetham.turbidity);
|
|
|
|
skyValidateDefinition(definition);
|
|
}
|
|
|
|
SkyDefinition skyCreateDefinition()
|
|
{
|
|
SkyDefinition def;
|
|
|
|
def.model = SKY_MODEL_CUSTOM;
|
|
def.daytime = 0.0;
|
|
def.sun_color = COLOR_BLACK;
|
|
def.sun_radius = 1.0;
|
|
def.sun_halo_size = 0.0;
|
|
def.sun_halo_profile = curveCreate();
|
|
def.model_custom.auto_from_daytime = 0;
|
|
def.model_custom.zenith_color = COLOR_BLACK;
|
|
def.model_custom.haze_color = COLOR_BLACK;
|
|
def.model_custom.haze_height = 0.0;
|
|
def.model_custom.haze_smoothing = 0.0;
|
|
def.model_custom._sky_gradation = colorGradationCreate();
|
|
def.model_preetham.turbidity = 0.0;
|
|
|
|
skyValidateDefinition(&def);
|
|
|
|
return def;
|
|
}
|
|
|
|
void skyDeleteDefinition(SkyDefinition* definition)
|
|
{
|
|
curveDelete(definition->sun_halo_profile);
|
|
colorGradationDelete(definition->model_custom._sky_gradation);
|
|
}
|
|
|
|
void skyCopyDefinition(SkyDefinition* source, SkyDefinition* destination)
|
|
{
|
|
destination->model = source->model;
|
|
destination->daytime = source->daytime;
|
|
destination->sun_color = source->sun_color;
|
|
destination->sun_radius = source->sun_radius;
|
|
destination->sun_halo_size = source->sun_halo_size;
|
|
destination->model_custom.auto_from_daytime = source->model_custom.auto_from_daytime;
|
|
destination->model_custom.zenith_color = source->model_custom.zenith_color;
|
|
destination->model_custom.haze_color = source->model_custom.haze_color;
|
|
destination->model_custom.haze_height = source->model_custom.haze_height;
|
|
destination->model_custom.haze_smoothing = source->model_custom.haze_smoothing;
|
|
destination->model_preetham.turbidity = source->model_preetham.turbidity;
|
|
|
|
curveCopy(source->sun_halo_profile, destination->sun_halo_profile);
|
|
|
|
skyValidateDefinition(destination);
|
|
}
|
|
|
|
static void _setAutoCustomModel(SkyDefinition* definition)
|
|
{
|
|
ColorGradation* zenith_gradation;
|
|
ColorGradation* haze_gradation;
|
|
|
|
zenith_gradation = colorGradationCreate();
|
|
haze_gradation = colorGradationCreate();
|
|
|
|
colorGradationQuickAddRgb(zenith_gradation, 0.2, 0.03, 0.03, 0.05);
|
|
colorGradationQuickAddRgb(zenith_gradation, 0.25, 0.25, 0.33, 0.37);
|
|
colorGradationQuickAddRgb(zenith_gradation, 0.35, 0.52, 0.63, 0.8);
|
|
colorGradationQuickAddRgb(zenith_gradation, 0.65, 0.52, 0.63, 0.8);
|
|
colorGradationQuickAddRgb(zenith_gradation, 0.75, 0.25, 0.33, 0.37);
|
|
colorGradationQuickAddRgb(zenith_gradation, 0.8, 0.03, 0.03, 0.05);
|
|
colorGradationQuickAddRgb(haze_gradation, 0.2, 0.05, 0.05, 0.08);
|
|
colorGradationQuickAddRgb(haze_gradation, 0.25, 0.55, 0.42, 0.42);
|
|
colorGradationQuickAddRgb(haze_gradation, 0.3, 0.6, 0.6, 0.6);
|
|
colorGradationQuickAddRgb(haze_gradation, 0.4, 0.92, 0.93, 1.0);
|
|
colorGradationQuickAddRgb(haze_gradation, 0.6, 0.92, 0.93, 1.0);
|
|
colorGradationQuickAddRgb(haze_gradation, 0.7, 0.6, 0.6, 0.8);
|
|
colorGradationQuickAddRgb(haze_gradation, 0.75, 0.62, 0.50, 0.42);
|
|
colorGradationQuickAddRgb(haze_gradation, 0.8, 0.05, 0.05, 0.08);
|
|
|
|
definition->model_custom.zenith_color = colorGradationGet(zenith_gradation, definition->daytime);
|
|
definition->model_custom.haze_color = colorGradationGet(haze_gradation, definition->daytime);
|
|
|
|
colorGradationDelete(zenith_gradation);
|
|
colorGradationDelete(haze_gradation);
|
|
}
|
|
|
|
void skyValidateDefinition(SkyDefinition* definition)
|
|
{
|
|
if (definition->model == SKY_MODEL_CUSTOM)
|
|
{
|
|
if (definition->model_custom.auto_from_daytime)
|
|
{
|
|
_setAutoCustomModel(definition);
|
|
}
|
|
colorGradationClear(definition->model_custom._sky_gradation);
|
|
colorGradationQuickAdd(definition->model_custom._sky_gradation, 0.0, &definition->model_custom.haze_color);
|
|
colorGradationQuickAdd(definition->model_custom._sky_gradation, definition->model_custom.haze_height - definition->model_custom.haze_smoothing, &definition->model_custom.haze_color);
|
|
colorGradationQuickAdd(definition->model_custom._sky_gradation, definition->model_custom.haze_height, &definition->model_custom.zenith_color);
|
|
colorGradationQuickAdd(definition->model_custom._sky_gradation, 1.0, &definition->model_custom.zenith_color);
|
|
}
|
|
}
|
|
|
|
int skyGetLights(SkyDefinition* sky, LightDefinition* lights, int max_lights)
|
|
{
|
|
double sun_angle;
|
|
Vector3 sun_direction;
|
|
int nblights = 0;
|
|
|
|
sun_angle = (sky->daytime + 0.75) * M_PI * 2.0;
|
|
sun_direction.x = cos(sun_angle);
|
|
sun_direction.y = sin(sun_angle);
|
|
sun_direction.z = 0.0;
|
|
|
|
/* TODO Night lights */
|
|
|
|
if (max_lights > 0)
|
|
{
|
|
/* Light from the sun */
|
|
lights[0].direction = v3Scale(sun_direction, -1.0);
|
|
lights[0].color = sky->sun_color;
|
|
lights[0].reflection = 1.0;
|
|
lights[0].filtered = 1;
|
|
lights[0].masked = 1;
|
|
lights[0].amplitude = 0.0;
|
|
nblights = 1;
|
|
if (max_lights > 1)
|
|
{
|
|
/* Skydome lighting */
|
|
lights[1].direction.x = 0.0;
|
|
lights[1].direction.y = -1.0;
|
|
lights[1].direction.z = 0.0;
|
|
lights[1].color = sky->model_custom.zenith_color;
|
|
lights[1].color.r *= 0.6;
|
|
lights[1].color.g *= 0.6;
|
|
lights[1].color.b *= 0.6;
|
|
lights[1].reflection = 0.0;
|
|
lights[1].filtered = 1;
|
|
lights[1].masked = 0;
|
|
lights[1].amplitude = M_PI / 2.0;
|
|
nblights = 2;
|
|
}
|
|
}
|
|
|
|
return nblights;
|
|
}
|
|
|
|
/******************************** Preetham Model ********************************/
|
|
static inline double _angleBetween(double thetav, double phiv, double theta, double phi)
|
|
{
|
|
double cospsi = sin(thetav) * sin(theta) * cos(phi-phiv) + cos(thetav) * cos(theta);
|
|
if (cospsi > 1.0) return 0.0;
|
|
if (cospsi < -1.0) return M_PI;
|
|
return acos(cospsi);
|
|
}
|
|
|
|
static inline Color _xyYToRGB(double x, double y, double Y)
|
|
{
|
|
double fX, fY, fZ;
|
|
Color result;
|
|
|
|
fY = Y;
|
|
fX = x / y * Y;
|
|
fZ = ((1.0 - x - y) / y) * Y;
|
|
|
|
double r, g, b;
|
|
r = 3.2404 * fX - 1.5371 * fY - 0.4985 * fZ;
|
|
g = -0.9692 * fX + 1.8759 * fY + 0.0415 * fZ;
|
|
b = 0.0556 * fX - 0.2040 * fY + 1.0573 * fZ;
|
|
|
|
double expo = -(1.0 / 10000.0);
|
|
r = 1.0 - exp(expo * r);
|
|
g = 1.0 - exp(expo * g);
|
|
b = 1.0 - exp(expo * b);
|
|
|
|
if (r < 0.0) r = 0.0;
|
|
if (g < 0.0) g = 0.0;
|
|
if (b < 0.0) b = 0.0;
|
|
|
|
result.r = r;
|
|
result.g = g;
|
|
result.b = b;
|
|
result.a = 1.0;
|
|
|
|
colorNormalize(&result);
|
|
|
|
return result;
|
|
}
|
|
|
|
static Color _preethamApproximate(SkyDefinition* definition, double theta, double phi)
|
|
{
|
|
double thetaSun;
|
|
double phiSun;
|
|
double gamma;
|
|
double turbidity = definition->model_preetham.turbidity;
|
|
|
|
/* Handle angles */
|
|
if (theta > M_PI / 2.0)
|
|
{
|
|
theta = M_PI / 2.0;
|
|
}
|
|
if (definition->daytime <= 0.5)
|
|
{
|
|
thetaSun = M_PI - definition->daytime * 2.0 * M_PI;
|
|
phiSun = 0.0;
|
|
}
|
|
else
|
|
{
|
|
thetaSun = (definition->daytime - 0.5) * 2.0 * M_PI;
|
|
phiSun = M_PI;
|
|
}
|
|
gamma = _angleBetween(theta, phi, thetaSun, phiSun);
|
|
|
|
double cosTheta;
|
|
if (theta > M_PI / 2.0)
|
|
{
|
|
cosTheta = 0.0000001;
|
|
}
|
|
else
|
|
{
|
|
cosTheta = cos(theta);
|
|
}
|
|
|
|
double T = turbidity;
|
|
double T2 = T * T;
|
|
double suntheta = thetaSun;
|
|
double suntheta2 = thetaSun * thetaSun;
|
|
double suntheta3 = thetaSun * suntheta2;
|
|
|
|
double Ax = -0.01925 * T - 0.25922;
|
|
double Bx = -0.06651 * T + 0.00081;
|
|
double Cx = -0.00041 * T + 0.21247;
|
|
double Dx = -0.06409 * T - 0.89887;
|
|
double Ex = -0.00325 * T + 0.04517;
|
|
|
|
double Ay = -0.01669 * T - 0.26078;
|
|
double By = -0.09495 * T + 0.00921;
|
|
double Cy = -0.00792 * T + 0.21023;
|
|
double Dy = -0.04405 * T - 1.65369;
|
|
double Ey = -0.01092 * T + 0.05291;
|
|
|
|
double AY = 0.17872 * T - 1.46303;
|
|
double BY = -0.35540 * T + 0.42749;
|
|
double CY = -0.02266 * T + 5.32505;
|
|
double DY = 0.12064 * T - 2.57705;
|
|
double EY = -0.06696 * T + 0.37027;
|
|
|
|
double cosGamma = cos(gamma);
|
|
cosGamma = cosGamma < 0.0 ? 0.0 : cosGamma;
|
|
double cosSTheta = fabs(cos(thetaSun));
|
|
|
|
double xz = ( 0.00165 * suntheta3 - 0.00375 * suntheta2 + 0.00209 * suntheta + 0.00000) * T2 +
|
|
(-0.02903 * suntheta3 + 0.06377 * suntheta2 - 0.03202 * suntheta + 0.00394) * T +
|
|
( 0.11693 * suntheta3 - 0.21196 * suntheta2 + 0.06052 * suntheta + 0.25886);
|
|
|
|
double yz = ( 0.00275 * suntheta3 - 0.00610 * suntheta2 + 0.00317 * suntheta + 0.00000) * T2 +
|
|
(-0.04214 * suntheta3 + 0.08970 * suntheta2 - 0.04153 * suntheta + 0.00516) * T +
|
|
( 0.15346 * suntheta3 - 0.26756 * suntheta2 + 0.06670 * suntheta + 0.26688);
|
|
|
|
double X = (4.0 / 9.0 - T / 120.0) * (M_PI - 2.0 * suntheta);
|
|
double Yz = ((4.0453 * T - 4.9710) * tan(X) - 0.2155 * T + 2.4192) * 1000.0;
|
|
|
|
double val1, val2;
|
|
|
|
val1 = (1.0 + Ax * exp(Bx / cosTheta)) * (1.0 + Cx * exp(Dx * gamma) + Ex * sqrt(cosGamma));
|
|
val2 = (1.0 + Ax * exp(Bx)) * (1.0 + Cx * exp(Dx * suntheta) + Ex * sqrt(cosSTheta));
|
|
double x = xz * val1 / val2;
|
|
|
|
val1 = (1.0 + Ay * exp(By / cosTheta)) * (1.0 + Cy * exp(Dy * gamma) + Ey * sqrt(cosGamma));
|
|
val2 = (1.0 + Ay * exp(By)) * (1.0 + Cy * exp(Dy * suntheta) + Ey * sqrt(cosSTheta));
|
|
double y = yz * val1 / val2;
|
|
|
|
val1 = (1.0 + AY * exp(BY / cosTheta)) * (1.0 + CY * exp(DY * gamma) + EY * sqrt(cosGamma));
|
|
val2 = (1.0 + AY * exp(BY)) * (1.0 + CY * exp(DY * suntheta) + EY * sqrt(cosSTheta));
|
|
double Y = Yz * val1 / val2;
|
|
|
|
return _xyYToRGB(x, y, Y);
|
|
}
|
|
|
|
/******************************** Rendering ********************************/
|
|
Color skyGetColor(SkyDefinition* definition, Renderer* renderer, Vector3 eye, Vector3 look)
|
|
{
|
|
double dist;
|
|
Vector3 sun_position;
|
|
Color sun_color, sky_color;
|
|
|
|
sun_position = skyGetSunDirection(definition);
|
|
|
|
look = v3Normalize(look);
|
|
dist = v3Norm(v3Sub(look, sun_position));
|
|
|
|
if (definition->model == SKY_MODEL_PREETHAM)
|
|
{
|
|
sky_color = _preethamApproximate(definition, M_PI/2.0 - asin(look.y), atan2(-look.z, look.x));
|
|
}
|
|
else
|
|
{
|
|
sky_color = colorGradationGet(definition->model_custom._sky_gradation, look.y * 0.5 + 0.5);
|
|
}
|
|
|
|
if (dist < definition->sun_radius + definition->sun_halo_size)
|
|
{
|
|
sun_color = definition->sun_color;
|
|
if (dist <= definition->sun_radius)
|
|
{
|
|
return sun_color;
|
|
}
|
|
else
|
|
{
|
|
dist = (dist - definition->sun_radius) / definition->sun_halo_size;
|
|
sun_color.a = curveGetValue(definition->sun_halo_profile, dist);
|
|
colorMask(&sky_color, &sun_color);
|
|
return sky_color;
|
|
}
|
|
}
|
|
else
|
|
{
|
|
return sky_color;
|
|
}
|
|
}
|
|
|
|
static Color _postProcessFragment(Renderer* renderer, Vector3 location, void* data)
|
|
{
|
|
Vector3 direction;
|
|
Color result;
|
|
SkyDefinition* definition;
|
|
|
|
definition = (SkyDefinition*)data;
|
|
|
|
direction = v3Sub(location, renderer->camera_location);
|
|
|
|
result = skyGetColor(definition, renderer, renderer->camera_location, v3Normalize(direction));
|
|
result = renderer->applyClouds(renderer, result, renderer->camera_location, v3Add(renderer->camera_location, v3Scale(direction, 10.0)));
|
|
|
|
return result;
|
|
}
|
|
|
|
void skyRender(SkyDefinition* definition, Renderer* renderer)
|
|
{
|
|
int res_i, res_j;
|
|
int i, j;
|
|
double step_i, step_j;
|
|
double current_i, current_j;
|
|
Vector3 vertex1, vertex2, vertex3, vertex4;
|
|
Vector3 direction;
|
|
|
|
res_i = renderer->render_quality * 40;
|
|
res_j = renderer->render_quality * 20;
|
|
step_i = M_PI * 2.0 / (double)res_i;
|
|
step_j = M_PI / (double)res_j;
|
|
|
|
for (j = 0; j < res_j; j++)
|
|
{
|
|
if (!renderer->addRenderProgress(renderer, 0.0))
|
|
{
|
|
return;
|
|
}
|
|
|
|
current_j = (double)(j - res_j / 2) * step_j;
|
|
|
|
for (i = 0; i < res_i; i++)
|
|
{
|
|
current_i = (double)i * step_i;
|
|
|
|
direction.x = SPHERE_SIZE * cos(current_i) * cos(current_j);
|
|
direction.y = SPHERE_SIZE * sin(current_j);
|
|
direction.z = SPHERE_SIZE * sin(current_i) * cos(current_j);
|
|
vertex1 = v3Add(renderer->camera_location, direction);
|
|
|
|
direction.x = SPHERE_SIZE * cos(current_i + step_i) * cos(current_j);
|
|
direction.y = SPHERE_SIZE * sin(current_j);
|
|
direction.z = SPHERE_SIZE * sin(current_i + step_i) * cos(current_j);
|
|
vertex2 = v3Add(renderer->camera_location, direction);
|
|
|
|
direction.x = SPHERE_SIZE * cos(current_i + step_i) * cos(current_j + step_j);
|
|
direction.y = SPHERE_SIZE * sin(current_j + step_j);
|
|
direction.z = SPHERE_SIZE * sin(current_i + step_i) * cos(current_j + step_j);
|
|
vertex3 = v3Add(renderer->camera_location, direction);
|
|
|
|
direction.x = SPHERE_SIZE * cos(current_i) * cos(current_j + step_j);
|
|
direction.y = SPHERE_SIZE * sin(current_j + step_j);
|
|
direction.z = SPHERE_SIZE * sin(current_i) * cos(current_j + step_j);
|
|
vertex4 = v3Add(renderer->camera_location, direction);
|
|
|
|
/* TODO Triangles at poles */
|
|
renderer->pushQuad(renderer, vertex1, vertex4, vertex3, vertex2, _postProcessFragment, definition);
|
|
}
|
|
}
|
|
}
|
|
|
|
Vector3 skyGetSunDirection(SkyDefinition* definition)
|
|
{
|
|
Vector3 result;
|
|
double sun_angle = (definition->daytime + 0.75) * M_PI * 2.0;
|
|
result.x = cos(sun_angle);
|
|
result.y = sin(sun_angle);
|
|
result.z = 0.0;
|
|
return result;
|
|
}
|
|
|
|
Color skyGetSunColor(SkyDefinition* definition)
|
|
{
|
|
return definition->sun_color;
|
|
}
|
|
|
|
Color skyGetZenithColor(SkyDefinition* definition)
|
|
{
|
|
return definition->model_custom.zenith_color;
|
|
}
|