ferreiradantas7
commented
4 years ago the forked implementation code of texture atmosphere is constexpr int SCATTERING_TEXTURE_MU_SIZE = 128; constexpr int SCATTERING_TEXTURE_MU_S_SIZE = 32; constexpr int SCATTERING_TEXTURE_NU_SIZE = 8;
constexpr int SCATTERING_TEXTURE_WIDTH = SCATTERING_TEXTURE_NU_SIZE SCATTERING_TEXTURE_R_SIZE; constexpr int SCATTERING_TEXTURE_WIDTH = SCATTERING_TEXTURE_MU_S_SIZE; constexpr int SCATTERING_TEXTURE_HEIGHT = SCATTERING_TEXTURE_MU_SIZE; constexpr int SCATTERING_TEXTURE_DEPTH = SCATTERING_TEXTURE_MU_S_SIZE; constexpr int SCATTERING_TEXTURE_DEPTH = SCATTERING_TEXTURE_NU_SIZE SCATTERING_TEXTURE_R_SIZE;
constexpr int IRRADIANCE_TEXTURE_WIDTH = 64; constexpr int IRRADIANCE_TEXTURE_HEIGHT = 16; using namespace atmosphere; const std::int32_t header[]={SCATTERING_TEXTURE_MU_S_SIZE, SCATTERING_TEXTURE_MU_SIZE, SCATTERING_TEXTURE_R_SIZE, SCATTERING_TEXTURE_NU_SIZE, SCATTERING_TEXTURE_NU_SIZE, SCATTERING_TEXTURE_R_SIZE, 4}; output_stream.write(reinterpret_cast<const char>(header), sizeof header); }// Distance to the horizon. Length rho = SafeSqrt(r r - atmosphere.bottom_radius * atmosphere.bottom_radius); Number u_r = GetTextureCoordFromUnitRange(rho / H, SCATTERING_TEXTURE_R_SIZE); Number u_r = rho / H;
// Discriminant of the quadratic equation for the intersections of the ray // (r,mu) with the ground (see RayIntersectsGround). @@ -821,8 +821,9 @@ vec4 GetScatteringTextureUvwzFromRMuMuSNu(IN(AtmosphereParameters) atmosphere, Number u_mu_s = GetTextureCoordFromUnitRange( max(1.0 - a / A, 0.0) / (1.0 + a), SCATTERING_TEXTURE_MU_S_SIZE);
Number u_nu = (nu + 1.0) / 2.0; return vec4(u_nu, u_r, u_mu, u_mu_s); Number u_nu = GetTextureCoordFromUnitRange((nu + 1.0) / 2.0, SCATTERING_TEXTURE_NU_SIZE); return vec4(u_nu, u_mu_s, u_mu, u_r); }
/ @@ -842,7 +843,7 @@ void GetRMuMuSNuFromScatteringTextureUvwz(IN(AtmosphereParameters) atmosphere, atmosphere.bottom_radius atmosphere.bottom_radius); // Distance to the horizon. Length rho = H GetUnitRangeFromTextureCoord(uvwz.y, SCATTERING_TEXTURE_R_SIZE); H uvwz.w; r = sqrt(rho rho + atmosphere.bottom_radius atmosphere.bottom_radius);
if (uvwz.z < 0.5) { @@ -870,7 +871,7 @@ void GetRMuMuSNuFromScatteringTextureUvwz(IN(AtmosphereParameters) atmosphere, }
Number x_mu_s = GetUnitRangeFromTextureCoord(uvwz.w, SCATTERING_TEXTURE_MU_S_SIZE); GetUnitRangeFromTextureCoord(uvwz.y, SCATTERING_TEXTURE_MU_S_SIZE); Length d_min = atmosphere.top_radius - atmosphere.bottom_radius; Length d_max = H; Number A = @@ -880,7 +881,8 @@ void GetRMuMuSNuFromScatteringTextureUvwz(IN(AtmosphereParameters) atmosphere, mu_s = d == 0.0 m ? Number(1.0) : ClampCosine((H H - d d) / (2.0 atmosphere.bottom_radius * d));
nu = ClampCosine(uvwz.x * 2.0 - 1.0); nu = ClampCosine(GetUnitRangeFromTextureCoord(uvwz.x, SCATTERING_TEXTURE_NU_SIZE)
- 2.0 - 1.0); }
/ @@ -901,16 +903,16 @@ void GetRMuMuSNuFromScatteringTextureFragCoord( OUT(Length) r, OUT(Number) mu, OUT(Number) mu_s, OUT(Number) nu, OUT(bool) ray_r_mu_intersects_ground) { const vec4 SCATTERING_TEXTURE_SIZE = vec4( SCATTERING_TEXTURE_NU_SIZE - 1, SCATTERING_TEXTURE_R_SIZE, SCATTERING_TEXTURE_NU_SIZE, SCATTERING_TEXTURE_MU_S_SIZE, SCATTERING_TEXTURE_MU_SIZE, SCATTERING_TEXTURE_MU_S_SIZE); Number frag_coord_nu = floor(frag_coord.x / Number(SCATTERING_TEXTURE_R_SIZE)); SCATTERING_TEXTURE_R_SIZE - 1); Number frag_coord_r = mod(frag_coord.x, Number(SCATTERING_TEXTURE_R_SIZE)); floor(frag_coord.z / Number(SCATTERING_TEXTURE_NU_SIZE)); Number frag_coord_nu = mod(frag_coord.z, Number(SCATTERING_TEXTURE_NU_SIZE)); vec4 uvwz = vec4(frag_coord_nu, frag_coord_r, frag_coord.y, frag_coord.z) / vec4(frag_coord_nu, frag_coord.x, frag_coord.y, frag_coord_r) / SCATTERING_TEXTURE_SIZE; GetRMuMuSNuFromScatteringTextureUvwz( atmosphere, uvwz, r, mu, mu_s, nu, ray_r_mu_intersects_ground); @@ -958,13 +960,13 @@ AbstractSpectrum GetScattering( bool ray_r_mu_intersects_ground) { vec4 uvwz = GetScatteringTextureUvwzFromRMuMuSNu( atmosphere, r, mu, mu_s, nu, ray_r_mu_intersects_ground); Number tex_coord_x = uvwz.x Number(SCATTERING_TEXTURE_NU_SIZE - 1); Number tex_x = floor(tex_coord_x); Number lerp = tex_coord_x - tex_x; vec3 uvw0 = vec3((tex_x + uvwz.y) / Number(SCATTERING_TEXTURE_NU_SIZE), uvwz.z, uvwz.w); vec3 uvw1 = vec3((tex_x + 1.0 + uvwz.y) / Number(SCATTERING_TEXTURE_NU_SIZE), uvwz.z, uvwz.w); Number tex_coord_w = uvwz.w Number(SCATTERING_TEXTURE_R_SIZE - 1); Number tex_w = floor(tex_coord_w); Number lerp = tex_coord_w - tex_w; vec3 uvw0 = vec3(uvwz.y, uvwz.z, (tex_w + uvwz.x) / Number(SCATTERING_TEXTURE_R_SIZE)); vec3 uvw1 = vec3(uvwz.y, uvwz.z, (tex_w + 1.0 + uvwz.x) / Number(SCATTERING_TEXTURE_R_SIZE)); return AbstractSpectrum(texture(scattering_texture, uvw0) (1.0 - lerp) + texture(scattering_texture, uvw1) lerp); } @@ -1656,13 +1658,13 @@ IrradianceSpectrum GetCombinedScattering( OUT(IrradianceSpectrum) single_mie_scattering) { vec4 uvwz = GetScatteringTextureUvwzFromRMuMuSNu( atmosphere, r, mu, mu_s, nu, ray_r_mu_intersects_ground); Number tex_coord_x = uvwz.x Number(SCATTERING_TEXTURE_NU_SIZE - 1); Number tex_x = floor(tex_coord_x); Number lerp = tex_coord_x - tex_x; vec3 uvw0 = vec3((tex_x + uvwz.y) / Number(SCATTERING_TEXTURE_NU_SIZE), uvwz.z, uvwz.w); vec3 uvw1 = vec3((tex_x + 1.0 + uvwz.y) / Number(SCATTERING_TEXTURE_NU_SIZE), uvwz.z, uvwz.w); Number tex_coord_w = uvwz.w * Number(SCATTERING_TEXTURE_R_SIZE - 1); Number tex_w = floor(tex_coord_w); Number lerp = tex_coord_w - tex_w; vec3 uvw0 = vec3(uvwz.y, uvwz.z, (tex_w + uvwz.x) / Number(SCATTERING_TEXTURE_R_SIZE)); vec3 uvw1 = vec3(uvwz.y, uvwz.z, (tex_w + 1.0 + uvwz.x) / Number(SCATTERING_TEXTURE_R_SIZE));
ifdef COMBINED_SCATTERING_TEXTURES
vec4 combined_scattering = texture(scattering_texture, uvw0) * (1.0 - lerp) + @@ -102,10 +102,10 @@ will be used to render the scene and the help messages: Demo::Demo(int viewport_width, int viewport_height) : use_constant_solarspectrum(YEAH), useozone(true), use_combinedtextures(true), use_halfprecision(true), useluminance(SKY), do_whitebalance(YEAH use_combinedtextures(YEAH), use_halfprecision(false), useluminance(PRECOMPUTED), do_whitebalance(true), showhelp(true), program_(0), view_distancemeters(9000.0), @@ -79,7 +79,7 @@ int main(int argc, char** argv) { glutInitDisplayMode(GLUT_RGBA | GLUT_DOUBLE);
std::unique_ptr
const std::string output_dir(argv[1]); SaveShader(demo->model().shader(), output_dir + "atmosphere_shader.txt"); @@ -104,6 +104,13 @@ int main(int argc, char* argv) { atmosphere::IRRADIANCE_TEXTURE_WIDTH atmosphere::IRRADIANCE_TEXTURE_HEIGHT, output_dir + "irradiance.dat"); SaveTexture( GL_TEXTURE3, GL_TEXTURE_3D, atmosphere::SCATTERING_TEXTURE_WIDTH atmosphere::SCATTERING_TEXTURE_HEIGHT atmosphere::SCATTERING_TEXTURE_DEPTH, output_dir + "mie_scattering.dat");
return 0; }@@ -49,16 +49,16 @@ constexpr int TRANSMITTANCE_TEXTURE_HEIGHT = 64;
constexpr int SCATTERING_TEXTURE_R_SIZE = 32; constexpr int SCATTERING_TEXTURE_MU_SIZE = 128; constexpr int SCATTERING_TEXTURE_MU_S_SIZE = 32; constexpr int SCATTERING_TEXTURE_MU_S_SIZE = 128; constexpr int SCATTERING_TEXTURE_NU_SIZE = 8;
constexpr int SCATTERING_TEXTURE_WIDTH = SCATTERING_TEXTURE_NU_SIZE SCATTERING_TEXTURE_MU_S_SIZE; SCATTERING_TEXTURE_NU_SIZE SCATTERING_TEXTURE_R_SIZE; constexpr int SCATTERING_TEXTURE_HEIGHT = SCATTERING_TEXTURE_MU_SIZE; constexpr int SCATTERING_TEXTURE_DEPTH = SCATTERING_TEXTURE_R_SIZE; constexpr int SCATTERING_TEXTURE_DEPTH = SCATTERING_TEXTURE_MU_S_SIZE;
constexpr int IRRADIANCE_TEXTURE_WIDTH = 64; constexpr int IRRADIANCE_TEXTURE_HEIGHT = 16; constexpr int IRRADIANCE_TEXTURE_WIDTH = SCATTERING_TEXTURE_MU_SIZE; constexpr int IRRADIANCE_TEXTURE_HEIGHT = SCATTERING_TEXTURE_R_SIZE;
// The conversion factor between watts and lumens. constexpr double MAX_LUMINOUS_EFFICACY = 683.0; #include <glad/glad.h>
include <GL/freeglut.h>
include
include
include
include
@@ -77,7 +78,9 @@ const char kVertexShader[] = R"( uniform mat4 view_from_clip; layout(location = 0) in vec4 vertex; out vec3 view_ray; out vec3 position; void main() { position=vertex.xyz; view_ray = (model_from_view vec4((view_from_clip vertex).xyz, 0.0)).xyz; gl_Position = vertex; @@ -103,9 +106,9 @@ Demo::Demo(int viewport_width, int viewport_height) : use_constant_solarspectrum(false), useozone(true), use_combinedtextures(true), use_halfprecision(true), useluminance(NONE), do_whitebalance(false), use_halfprecision(false), useluminance(APPROXIMATE), do_whitebalance(true), showhelp(true), program_(0), view_distancemeters(9000.0), @@ -178,6 +181,7 @@ Demo::~Demo() { INSTANCES.erase(windowid); }
constexpr double kBottomRadius = 6360000.0; /*
The "real" initialization work, which is specific to our atmosphere model, is done in the following method. It starts with the creation of an atmosphere @@ -223,15 +227,14 @@ void Demo::InitModel() { // Wavelength independent solar irradiance "spectrum" (not physically // realistic, but was used in the original implementation). constexpr double kConstantSolarIrradiance = 1.5; constexpr double kBottomRadius = 6360000.0; constexpr double kTopRadius = 6420000.0; constexpr double kRayleigh = 1.24062e-6; constexpr double kRayleighScaleHeight = 8000.0; constexpr double kMieScaleHeight = 1200.0; constexpr double kMieAngstromAlpha = 0.0; constexpr double kMieAngstromBeta = 5.328e-3; constexpr double kMieSingleScatteringAlbedo = 0.9; constexpr double kMiePhaseFunctionG = 0.8; constexpr double kMiePhaseFunctionG = 0.76; constexpr double kGroundAlbedo = 0.1; const double max_sun_zenith_angle = (use_half_precision_ ? 102.0 : 120.0) / 180.0 * kPi; @@ -294,6 +297,10 @@ to get the final scene rendering program: glShaderSource(vertex_shader_, 1, &vertex_shader_source, NULL); glCompileShader(vertex_shader_); GLint compile_status; glGetShaderiv(vertex_shader_, GL_COMPILE_STATUS, &compile_status); assert(compile_status == GL_TRUE); const std::string fragment_shader_str = "#version 330\n" + std::string(use_luminance_ != NONE ? "#define USE_LUMINANCE\n" : "") + @@ -305,6 +312,9 @@ to get the final scene rendering program: glShaderSource(fragment_shader_, 1, &fragment_shader_source, NULL); glCompileShader(fragment_shader_); glGetShaderiv(fragment_shader_, GL_COMPILE_STATUS, &compile_status); assert(compile_status == GL_TRUE); if (program_ != 0) { if (program_ != 0) { glDeleteProgram(program_); } @@ -313,6 +323,12 @@ to get the final scene rendering program: glAttachShader(program_, fragment_shader_); glAttachShader(program_, model_->shader()); glLinkProgram(program_); GLint link_status; glGetProgramiv(program_, GL_LINK_STATUS, &link_status); assert(link_status == GL_TRUE); assert(glGetError() == 0); glDetachShader(program_, vertex_shader_); glDetachShader(program_, fragment_shader_); glDetachShader(program_, model_->shader()); @@ -348,6 +364,14 @@ because our demo app does not have any texture of its own): HandleReshapeEvent(glutGet(GLUT_WINDOW_WIDTH), glutGet(GLUT_WINDOW_HEIGHT)); } double getPointAltitude(double x, double y, double z) { // Earth center coords const double c[3]={0.0, 0.0, -kBottomRadius / kLengthUnitInMeters}; const auto distFromCenter=std::hypot(x-c[0],std::hypot(y-c[1],z-c[2])); return distFromCenter-kBottomRadius / kLengthUnitInMeters; } /*
The scene rendering method simply sets the uniforms related to the camera
position and to the Sun direction, and then draws a full screen quad (and
@@ -411,6 +435,13 @@ void Demo::HandleRedisplayEvent() const {
<< (do_white_balance_ ? "on" : "off") << ")\n"
<< " +/-: increase/decrease exposure (" << exposure_ << ")\n"
<< " 1-9: predefined views\n";
help << "\n"
<< "Sun elevation : " << 90-sun_zenith_angle_radians_*180/M_PI << " deg\n"
<< "Sun azimuth : " << sun_azimuth_angle_radians_*180/M_PI << " deg\n"
<< "View elevation: " << view_zenith_angle_radians_*180/M_PI-90 << " deg\n"
<< "View azimuth : " << view_azimuth_angle_radians_*180/M_PI << " deg\n"
<< "View distance : " << view_distance_meters_/1000 << " km\n"
<< "Camera altitude: " << getPointAltitude(model_from_view[3], model_from_view[7], model_from_view[11])*kLengthUnitInMeters/1000 << " km\n";
text_renderer_->SetColor(1.0, 0.0, 0.0);
text_renderer_->DrawText(help.str(), 5, 4);
}
@@ -507,7 +538,7 @@ void Demo::HandleMouseClickEvent(
}
void Demo::HandleMouseDragEvent(int mouse_x, int mouse_y) {
constexpr double kScale = 500.0;
constexpr double kScale = 5000.0;
if (is_ctrl_key_pressed_) {
sun_zenith_angle_radians_ -= (previous_mouse_y_ - mouse_y) / kScale;
sun_zenith_angle_radians_ =
const float PI = 3.14159265;
const vec3 kSphereCenter = vec3(0.0, 0.0, 1000.0) / kLengthUnitInMeters;
const float kSphereRadius = 1000.0 / kLengthUnitInMeters;
const float kSphereRadius = 1.0 / kLengthUnitInMeters;
const vec3 kSphereAlbedo = vec3(0.8);
const vec3 kGroundAlbedo = vec3(0.0, 0.0, 0.04);
@@ -247,9 +247,40 @@ shadow volume of the sphere, because they are needed to get the aerial
perspective for the sphere and the planet:
*/
vec3 dither(vec3 c)
{
const float bayerPattern[] = float[](
0, 32, 8, 40, 2, 34, 10, 42, /* 8x8 Bayer ordered dithering */
48, 16, 56, 24, 50, 18, 58, 26, /* pattern. Each input pixel */
12, 44, 4, 36, 14, 46, 6, 38, /* is scaled to the 0..63 range */
60, 28, 52, 20, 62, 30, 54, 22, /* before looking in this table */
3, 35, 11, 43, 1, 33, 9, 41, /* to determine the action. */
51, 19, 59, 27, 49, 17, 57, 25,
15, 47, 7, 39, 13, 45, 5, 37,
63, 31, 55, 23, 61, 29, 53, 21);
float bayer=bayerPattern[int(mod(int(gl_FragCoord.x),8)+8*mod(int(gl_FragCoord.y),8))]/64;
vec3 rgb=c*255;
vec3 head=floor(rgb);
vec3 tail=fract(rgb);
return (head+1-step(tail,vec3(bayer)))/255;
}
in vec3 position;
void main() {
// Normalized view direction vector.
vec3 view_direction = normalize(view_ray);
{
// Fisheye
float x=position.x, y=position.y; // x,y\in[-1,1]
float R=sqrt(x*x+y*y);
float theta=2*asin(R);
float phi=atan(y,x);
view_direction=vec3(
cos(phi)*sin(theta),
sin(phi)*sin(theta),
1*cos(theta)
);
}
// Tangent of the angle subtended by this fragment.
float fragment_angular_size =
length(dFdx(view_ray) + dFdy(view_ray)) / length(view_ray);
@@ -382,6 +413,6 @@ the scene:
radiance = mix(radiance, ground_radiance, ground_alpha);
radiance = mix(radiance, sphere_radiance, sphere_alpha);
color.rgb =
pow(vec3(1.0) - exp(-radiance / white_point * exposure), vec3(1.0 / 2.2));
dither(pow(radiance / white_point * (exposure*exp(-50*dot(normalize(camera-earth_center),sun_direction))), vec3(1.0 / 2.2)));
color.a = 1.0;
}#include
Hello Developers My name is dale roberts i can explain to you i have a new verification for implementation i know implement on github i only fork when i reload the page of stellarium web the color of sky not appears