@cesium/engine 6.1.1 → 7.0.0

package/Source/Shaders/Model/AtmosphereStageFS.js

@@ -0,0 +1,113 @@
+//This file is automatically rebuilt by the Cesium build process.
+export default "// robust iterative solution without trig functions\n\
+// https://github.com/0xfaded/ellipse_demo/issues/1\n\
+// https://stackoverflow.com/questions/22959698/distance-from-given-point-to-given-ellipse\n\
+//\n\
+// This version uses only a single iteration for best performance. For fog\n\
+// rendering, the difference is negligible.\n\
+vec2 nearestPointOnEllipseFast(vec2 pos, vec2 radii) {\n\
+    vec2 p = abs(pos);\n\
+    vec2 inverseRadii = 1.0 / radii;\n\
+    vec2 evoluteScale = (radii.x * radii.x - radii.y * radii.y) * vec2(1.0, -1.0) * inverseRadii;\n\
+\n\
+    // We describe the ellipse parametrically: v = radii * vec2(cos(t), sin(t))\n\
+    // but store the cos and sin of t in a vec2 for efficiency.\n\
+    // Initial guess: t = cos(pi/4)\n\
+    vec2 tTrigs = vec2(0.70710678118);\n\
+    vec2 v = radii * tTrigs;\n\
+\n\
+    // Find the evolute of the ellipse (center of curvature) at v.\n\
+    vec2 evolute = evoluteScale * tTrigs * tTrigs * tTrigs;\n\
+    // Find the (approximate) intersection of p - evolute with the ellipsoid.\n\
+    vec2 q = normalize(p - evolute) * length(v - evolute);\n\
+    // Update the estimate of t.\n\
+    tTrigs = (q + evolute) * inverseRadii;\n\
+    tTrigs = normalize(clamp(tTrigs, 0.0, 1.0));\n\
+    v = radii * tTrigs;\n\
+\n\
+    return v * sign(pos);\n\
+}\n\
+\n\
+vec3 computeEllipsoidPositionWC(vec3 positionMC) {\n\
+    // Get the world-space position and project onto a meridian plane of\n\
+    // the ellipsoid\n\
+    vec3 positionWC = (czm_model * vec4(positionMC, 1.0)).xyz;\n\
+\n\
+    vec2 positionEllipse = vec2(length(positionWC.xy), positionWC.z);\n\
+    vec2 nearestPoint = nearestPointOnEllipseFast(positionEllipse, czm_ellipsoidRadii.xz);\n\
+\n\
+    // Reconstruct a 3D point in world space\n\
+    return vec3(nearestPoint.x * normalize(positionWC.xy), nearestPoint.y);\n\
+}\n\
+\n\
+void applyFog(inout vec4 color, vec4 groundAtmosphereColor, vec3 lightDirection, float distanceToCamera) {\n\
+\n\
+    vec3 fogColor = groundAtmosphereColor.rgb;\n\
+\n\
+    // If there is dynamic lighting, apply that to the fog.\n\
+    const float NONE = 0.0;\n\
+    if (czm_atmosphereDynamicLighting != NONE) {\n\
+        float darken = clamp(dot(normalize(czm_viewerPositionWC), lightDirection), czm_fogMinimumBrightness, 1.0);\n\
+        fogColor *= darken;\n\
+    }\n\
+\n\
+    // Tonemap if HDR rendering is disabled\n\
+    #ifndef HDR\n\
+        fogColor.rgb = czm_acesTonemapping(fogColor.rgb);\n\
+        fogColor.rgb = czm_inverseGamma(fogColor.rgb);\n\
+    #endif\n\
+\n\
+    // Matches the constant in GlobeFS.glsl. This makes the fog falloff\n\
+    // more gradual.\n\
+    const float fogModifier = 0.15;\n\
+    vec3 withFog = czm_fog(distanceToCamera, color.rgb, fogColor, fogModifier);\n\
+    color = vec4(withFog, color.a);\n\
+}\n\
+\n\
+void atmosphereStage(inout vec4 color, in ProcessedAttributes attributes) {\n\
+    vec3 rayleighColor;\n\
+    vec3 mieColor;\n\
+    float opacity;\n\
+\n\
+    vec3 positionWC;\n\
+    vec3 lightDirection;\n\
+\n\
+    // When the camera is in space, compute the position per-fragment for\n\
+    // more accurate ground atmosphere. All other cases will use\n\
+    //\n\
+    // The if condition will be added in https://github.com/CesiumGS/cesium/issues/11717\n\
+    if (false) {\n\
+        positionWC = computeEllipsoidPositionWC(attributes.positionMC);\n\
+        lightDirection = czm_getDynamicAtmosphereLightDirection(positionWC, czm_atmosphereDynamicLighting);\n\
+\n\
+        // The fog color is derived from the ground atmosphere color\n\
+        czm_computeGroundAtmosphereScattering(\n\
+            positionWC,\n\
+            lightDirection,\n\
+            rayleighColor,\n\
+            mieColor,\n\
+            opacity\n\
+        );\n\
+    } else {\n\
+        positionWC = attributes.positionWC;\n\
+        lightDirection = czm_getDynamicAtmosphereLightDirection(positionWC, czm_atmosphereDynamicLighting);\n\
+        rayleighColor = v_atmosphereRayleighColor;\n\
+        mieColor = v_atmosphereMieColor;\n\
+        opacity = v_atmosphereOpacity;\n\
+    }\n\
+\n\
+    //color correct rayleigh and mie colors\n\
+    const bool ignoreBlackPixels = true;\n\
+    rayleighColor = czm_applyHSBShift(rayleighColor, czm_atmosphereHsbShift, ignoreBlackPixels);\n\
+    mieColor = czm_applyHSBShift(mieColor, czm_atmosphereHsbShift, ignoreBlackPixels);\n\
+\n\
+    vec4 groundAtmosphereColor = czm_computeAtmosphereColor(positionWC, lightDirection, rayleighColor, mieColor, opacity);\n\
+\n\
+    if (u_isInFog) {\n\
+        float distanceToCamera = length(attributes.positionEC);\n\
+        applyFog(color, groundAtmosphereColor, lightDirection, distanceToCamera);\n\
+    } else {\n\
+        // Ground atmosphere\n\
+    }\n\
+}\n\
+";

package/Source/Shaders/Model/AtmosphereStageVS.glsl

@@ -0,0 +1,12 @@
+void atmosphereStage(ProcessedAttributes attributes) {
+    vec3 lightDirection = czm_getDynamicAtmosphereLightDirection(v_positionWC, czm_atmosphereDynamicLighting);
+
+    czm_computeGroundAtmosphereScattering(
+        // This assumes the geometry stage came before this.
+        v_positionWC,
+        lightDirection,
+        v_atmosphereRayleighColor,
+        v_atmosphereMieColor,
+        v_atmosphereOpacity
+    );
+}

package/Source/Shaders/Model/AtmosphereStageVS.js

@@ -0,0 +1,14 @@
+//This file is automatically rebuilt by the Cesium build process.
+export default "void atmosphereStage(ProcessedAttributes attributes) {\n\
+    vec3 lightDirection = czm_getDynamicAtmosphereLightDirection(v_positionWC, czm_atmosphereDynamicLighting);\n\
+\n\
+    czm_computeGroundAtmosphereScattering(\n\
+        // This assumes the geometry stage came before this.\n\
+        v_positionWC,\n\
+        lightDirection,\n\
+        v_atmosphereRayleighColor,\n\
+        v_atmosphereMieColor,\n\
+        v_atmosphereOpacity\n\
+    );\n\
+}\n\
+";

package/Source/Shaders/Model/GeometryStageFS.glsl

@@ -3,7 +3,7 @@ void geometryStage(out ProcessedAttributes attributes)
   attributes.positionMC = v_positionMC;
   attributes.positionEC = v_positionEC;
 
-  #ifdef COMPUTE_POSITION_WC_CUSTOM_SHADER
+  #if defined(COMPUTE_POSITION_WC_CUSTOM_SHADER) || defined(COMPUTE_POSITION_WC_STYLE) || defined(COMPUTE_POSITION_WC_ATMOSPHERE)
   attributes.positionWC = v_positionWC;
   #endif
 

package/Source/Shaders/Model/GeometryStageFS.js

@@ -4,7 +4,7 @@ export default "void geometryStage(out ProcessedAttributes attributes)\n\
   attributes.positionMC = v_positionMC;\n\
   attributes.positionEC = v_positionEC;\n\
 \n\
-  #ifdef COMPUTE_POSITION_WC_CUSTOM_SHADER\n\
+  #if defined(COMPUTE_POSITION_WC_CUSTOM_SHADER) || defined(COMPUTE_POSITION_WC_STYLE) || defined(COMPUTE_POSITION_WC_ATMOSPHERE)\n\
   attributes.positionWC = v_positionWC;\n\
   #endif\n\
 \n\

package/Source/Shaders/Model/GeometryStageVS.glsl

@@ -1,4 +1,4 @@
-vec4 geometryStage(inout ProcessedAttributes attributes, mat4 modelView, mat3 normal) 
+vec4 geometryStage(inout ProcessedAttributes attributes, mat4 modelView, mat3 normal)
 {
     vec4 computedPosition;
 
@@ -16,7 +16,7 @@ vec4 geometryStage(inout ProcessedAttributes attributes, mat4 modelView, mat3 no
     #endif
 
     // Sometimes the custom shader and/or style needs this
-    #if defined(COMPUTE_POSITION_WC_CUSTOM_SHADER) || defined(COMPUTE_POSITION_WC_STYLE)
+    #if defined(COMPUTE_POSITION_WC_CUSTOM_SHADER) || defined(COMPUTE_POSITION_WC_STYLE) || defined(COMPUTE_POSITION_WC_ATMOSPHERE)
     // Note that this is a 32-bit position which may result in jitter on small
     // scales.
     v_positionWC = (czm_model * vec4(positionMC, 1.0)).xyz;
@@ -27,7 +27,7 @@ vec4 geometryStage(inout ProcessedAttributes attributes, mat4 modelView, mat3 no
     #endif
 
     #ifdef HAS_TANGENTS
-    v_tangentEC = normalize(normal * attributes.tangentMC);    
+    v_tangentEC = normalize(normal * attributes.tangentMC);
     #endif
 
     #ifdef HAS_BITANGENTS
@@ -37,6 +37,6 @@ vec4 geometryStage(inout ProcessedAttributes attributes, mat4 modelView, mat3 no
     // All other varyings need to be dynamically generated in
     // GeometryPipelineStage
     setDynamicVaryings(attributes);
-    
+
     return computedPosition;
 }

package/Source/Shaders/Model/GeometryStageVS.js

@@ -1,5 +1,5 @@
 //This file is automatically rebuilt by the Cesium build process.
-export default "vec4 geometryStage(inout ProcessedAttributes attributes, mat4 modelView, mat3 normal) \n\
+export default "vec4 geometryStage(inout ProcessedAttributes attributes, mat4 modelView, mat3 normal)\n\
 {\n\
     vec4 computedPosition;\n\
 \n\
@@ -17,7 +17,7 @@ export default "vec4 geometryStage(inout ProcessedAttributes attributes, mat4 mo
     #endif\n\
 \n\
     // Sometimes the custom shader and/or style needs this\n\
-    #if defined(COMPUTE_POSITION_WC_CUSTOM_SHADER) || defined(COMPUTE_POSITION_WC_STYLE)\n\
+    #if defined(COMPUTE_POSITION_WC_CUSTOM_SHADER) || defined(COMPUTE_POSITION_WC_STYLE) || defined(COMPUTE_POSITION_WC_ATMOSPHERE)\n\
     // Note that this is a 32-bit position which may result in jitter on small\n\
     // scales.\n\
     v_positionWC = (czm_model * vec4(positionMC, 1.0)).xyz;\n\
@@ -28,7 +28,7 @@ export default "vec4 geometryStage(inout ProcessedAttributes attributes, mat4 mo
     #endif\n\
 \n\
     #ifdef HAS_TANGENTS\n\
-    v_tangentEC = normalize(normal * attributes.tangentMC);    \n\
+    v_tangentEC = normalize(normal * attributes.tangentMC);\n\
     #endif\n\
 \n\
     #ifdef HAS_BITANGENTS\n\
@@ -38,7 +38,7 @@ export default "vec4 geometryStage(inout ProcessedAttributes attributes, mat4 mo
     // All other varyings need to be dynamically generated in\n\
     // GeometryPipelineStage\n\
     setDynamicVaryings(attributes);\n\
-    \n\
+\n\
     return computedPosition;\n\
 }\n\
 ";

package/Source/Shaders/Model/ModelColorStageFS.glsl

@@ -4,4 +4,4 @@ void modelColorStage(inout czm_modelMaterial material)
     float highlight = ceil(model_colorBlend);
     material.diffuse *= mix(model_color.rgb, vec3(1.0), highlight);
     material.alpha *= model_color.a;
-}
+}

package/Source/Shaders/Model/ModelColorStageFS.js

@@ -5,4 +5,5 @@ export default "void modelColorStage(inout czm_modelMaterial material)\n\
     float highlight = ceil(model_colorBlend);\n\
     material.diffuse *= mix(model_color.rgb, vec3(1.0), highlight);\n\
     material.alpha *= model_color.a;\n\
-}";
+}\n\
+";

package/Source/Shaders/Model/ModelFS.glsl

@@ -79,5 +79,9 @@ void main()
     silhouetteStage(color);
     #endif
 
+    #ifdef HAS_ATMOSPHERE
+    atmosphereStage(color, attributes);
+    #endif
+
     out_FragColor = color;
 }

package/Source/Shaders/Model/ModelFS.js

@@ -79,6 +79,10 @@ void main()\n\
     #if defined(HAS_SILHOUETTE) && defined(HAS_NORMALS)\n\
     silhouetteStage(color);\n\
     #endif\n\
+\n\
+    #ifdef HAS_ATMOSPHERE\n\
+    atmosphereStage(color, attributes);\n\
+    #endif\n\
 \n\
     out_FragColor = color;\n\
 }\n\

package/Source/Shaders/Model/ModelVS.glsl

@@ -7,7 +7,7 @@ czm_modelVertexOutput defaultVertexOutput(vec3 positionMC) {
     return vsOutput;
 }
 
-void main() 
+void main()
 {
     // Initialize the attributes struct with all
     // attributes except quantized ones.
@@ -66,14 +66,14 @@ void main()
     // Update the position for this instance in place
     #ifdef HAS_INSTANCING
 
-        // The legacy instance stage is used when rendering i3dm models that 
+        // The legacy instance stage is used when rendering i3dm models that
         // encode instances transforms in world space, as opposed to glTF models
         // that use EXT_mesh_gpu_instancing, where instance transforms are encoded
         // in object space.
         #ifdef USE_LEGACY_INSTANCING
         mat4 instanceModelView;
         mat3 instanceModelViewInverseTranspose;
-        
+
         legacyInstancingStage(attributes, instanceModelView, instanceModelViewInverseTranspose);
 
         modelView = instanceModelView;
@@ -93,6 +93,10 @@ void main()
     MetadataStatistics metadataStatistics;
     metadataStage(metadata, metadataClass, metadataStatistics, attributes);
 
+    #ifdef HAS_VERTICAL_EXAGGERATION
+    verticalExaggerationStage(attributes);
+    #endif
+
     #ifdef HAS_CUSTOM_VERTEX_SHADER
     czm_modelVertexOutput vsOutput = defaultVertexOutput(attributes.positionMC);
     customShaderStage(vsOutput, attributes, featureIds, metadata, metadataClass, metadataStatistics);
@@ -100,7 +104,12 @@ void main()
 
     // Compute the final position in each coordinate system needed.
     // This returns the value that will be assigned to gl_Position.
-    vec4 positionClip = geometryStage(attributes, modelView, normal);    
+    vec4 positionClip = geometryStage(attributes, modelView, normal);
+
+    // This must go after the geometry stage as it needs v_positionWC
+    #ifdef HAS_ATMOSPHERE
+    atmosphereStage(attributes);
+    #endif
 
     #ifdef HAS_SILHOUETTE
     silhouetteStage(attributes, positionClip);

package/Source/Shaders/Model/ModelVS.js

@@ -8,7 +8,7 @@ czm_modelVertexOutput defaultVertexOutput(vec3 positionMC) {\n\
     return vsOutput;\n\
 }\n\
 \n\
-void main() \n\
+void main()\n\
 {\n\
     // Initialize the attributes struct with all\n\
     // attributes except quantized ones.\n\
@@ -67,14 +67,14 @@ void main() \n\
     // Update the position for this instance in place\n\
     #ifdef HAS_INSTANCING\n\
 \n\
-        // The legacy instance stage is used when rendering i3dm models that \n\
+        // The legacy instance stage is used when rendering i3dm models that\n\
         // encode instances transforms in world space, as opposed to glTF models\n\
         // that use EXT_mesh_gpu_instancing, where instance transforms are encoded\n\
         // in object space.\n\
         #ifdef USE_LEGACY_INSTANCING\n\
         mat4 instanceModelView;\n\
         mat3 instanceModelViewInverseTranspose;\n\
-        \n\
+\n\
         legacyInstancingStage(attributes, instanceModelView, instanceModelViewInverseTranspose);\n\
 \n\
         modelView = instanceModelView;\n\
@@ -93,6 +93,10 @@ void main() \n\
     MetadataClass metadataClass;\n\
     MetadataStatistics metadataStatistics;\n\
     metadataStage(metadata, metadataClass, metadataStatistics, attributes);\n\
+\n\
+    #ifdef HAS_VERTICAL_EXAGGERATION\n\
+    verticalExaggerationStage(attributes);\n\
+    #endif\n\
 \n\
     #ifdef HAS_CUSTOM_VERTEX_SHADER\n\
     czm_modelVertexOutput vsOutput = defaultVertexOutput(attributes.positionMC);\n\
@@ -101,7 +105,12 @@ void main() \n\
 \n\
     // Compute the final position in each coordinate system needed.\n\
     // This returns the value that will be assigned to gl_Position.\n\
-    vec4 positionClip = geometryStage(attributes, modelView, normal);    \n\
+    vec4 positionClip = geometryStage(attributes, modelView, normal);\n\
+\n\
+    // This must go after the geometry stage as it needs v_positionWC\n\
+    #ifdef HAS_ATMOSPHERE\n\
+    atmosphereStage(attributes);\n\
+    #endif\n\
 \n\
     #ifdef HAS_SILHOUETTE\n\
     silhouetteStage(attributes, positionClip);\n\

package/Source/Shaders/Model/VerticalExaggerationStageVS.glsl

@@ -0,0 +1,39 @@
+void verticalExaggerationStage(
+  inout ProcessedAttributes attributes
+) {
+  // Compute the distance from the camera to the local center of curvature.
+  vec4 vertexPositionENU = czm_modelToEnu * vec4(attributes.positionMC, 1.0);
+  vec2 vertexAzimuth = normalize(vertexPositionENU.xy);
+  // Curvature = 1 / radius of curvature.
+  float azimuthalCurvature = dot(vertexAzimuth * vertexAzimuth, czm_eyeEllipsoidCurvature);
+  float eyeToCenter = 1.0 / azimuthalCurvature + czm_eyeHeight;
+
+  // Compute the approximate ellipsoid normal at the vertex position.
+  // Uses a circular approximation for the Earth curvature along the geodesic.
+  vec3 vertexPositionEC = (czm_modelView * vec4(attributes.positionMC, 1.0)).xyz;
+  vec3 centerToVertex = eyeToCenter * czm_eyeEllipsoidNormalEC + vertexPositionEC;
+  vec3 vertexNormal = normalize(centerToVertex);
+
+  // Estimate the (sine of the) angle between the camera direction and the vertex normal
+  float verticalDistance = dot(vertexPositionEC, czm_eyeEllipsoidNormalEC);
+  float horizontalDistance = length(vertexPositionEC - verticalDistance * czm_eyeEllipsoidNormalEC);
+  float sinTheta = horizontalDistance / (eyeToCenter + verticalDistance);
+  bool isSmallAngle = clamp(sinTheta, 0.0, 0.05) == sinTheta;
+
+  // Approximate the change in height above the ellipsoid, from camera to vertex position.
+  float exactVersine = 1.0 - dot(czm_eyeEllipsoidNormalEC, vertexNormal);
+  float smallAngleVersine = 0.5 * sinTheta * sinTheta;
+  float versine = isSmallAngle ? smallAngleVersine : exactVersine;
+  float dHeight = dot(vertexPositionEC, vertexNormal) - eyeToCenter * versine;
+  float vertexHeight = czm_eyeHeight + dHeight;
+
+  // Transform the approximate vertex normal to model coordinates.
+  vec3 vertexNormalMC = (czm_inverseModelView * vec4(vertexNormal, 0.0)).xyz;
+  vertexNormalMC = normalize(vertexNormalMC);
+
+  // Compute the exaggeration and apply it along the approximate vertex normal.
+  float stretch = u_verticalExaggerationAndRelativeHeight.x;
+  float shift = u_verticalExaggerationAndRelativeHeight.y;
+  float exaggeration = (vertexHeight - shift) * (stretch - 1.0);
+  attributes.positionMC += exaggeration * vertexNormalMC;
+}

package/Source/Shaders/Model/VerticalExaggerationStageVS.js

@@ -0,0 +1,41 @@
+//This file is automatically rebuilt by the Cesium build process.
+export default "void verticalExaggerationStage(\n\
+  inout ProcessedAttributes attributes\n\
+) {\n\
+  // Compute the distance from the camera to the local center of curvature.\n\
+  vec4 vertexPositionENU = czm_modelToEnu * vec4(attributes.positionMC, 1.0);\n\
+  vec2 vertexAzimuth = normalize(vertexPositionENU.xy);\n\
+  // Curvature = 1 / radius of curvature.\n\
+  float azimuthalCurvature = dot(vertexAzimuth * vertexAzimuth, czm_eyeEllipsoidCurvature);\n\
+  float eyeToCenter = 1.0 / azimuthalCurvature + czm_eyeHeight;\n\
+\n\
+  // Compute the approximate ellipsoid normal at the vertex position.\n\
+  // Uses a circular approximation for the Earth curvature along the geodesic.\n\
+  vec3 vertexPositionEC = (czm_modelView * vec4(attributes.positionMC, 1.0)).xyz;\n\
+  vec3 centerToVertex = eyeToCenter * czm_eyeEllipsoidNormalEC + vertexPositionEC;\n\
+  vec3 vertexNormal = normalize(centerToVertex);\n\
+\n\
+  // Estimate the (sine of the) angle between the camera direction and the vertex normal\n\
+  float verticalDistance = dot(vertexPositionEC, czm_eyeEllipsoidNormalEC);\n\
+  float horizontalDistance = length(vertexPositionEC - verticalDistance * czm_eyeEllipsoidNormalEC);\n\
+  float sinTheta = horizontalDistance / (eyeToCenter + verticalDistance);\n\
+  bool isSmallAngle = clamp(sinTheta, 0.0, 0.05) == sinTheta;\n\
+\n\
+  // Approximate the change in height above the ellipsoid, from camera to vertex position.\n\
+  float exactVersine = 1.0 - dot(czm_eyeEllipsoidNormalEC, vertexNormal);\n\
+  float smallAngleVersine = 0.5 * sinTheta * sinTheta;\n\
+  float versine = isSmallAngle ? smallAngleVersine : exactVersine;\n\
+  float dHeight = dot(vertexPositionEC, vertexNormal) - eyeToCenter * versine;\n\
+  float vertexHeight = czm_eyeHeight + dHeight;\n\
+\n\
+  // Transform the approximate vertex normal to model coordinates.\n\
+  vec3 vertexNormalMC = (czm_inverseModelView * vec4(vertexNormal, 0.0)).xyz;\n\
+  vertexNormalMC = normalize(vertexNormalMC);\n\
+\n\
+  // Compute the exaggeration and apply it along the approximate vertex normal.\n\
+  float stretch = u_verticalExaggerationAndRelativeHeight.x;\n\
+  float shift = u_verticalExaggerationAndRelativeHeight.y;\n\
+  float exaggeration = (vertexHeight - shift) * (stretch - 1.0);\n\
+  attributes.positionMC += exaggeration * vertexNormalMC;\n\
+}\n\
+";

package/Source/Shaders/SkyAtmosphereCommon.glsl

@@ -8,16 +8,6 @@ float interpolateByDistance(vec4 nearFarScalar, float distance)
     return mix(startValue, endValue, t);
 }
 
-vec3 getLightDirection(vec3 positionWC)
-{
-    float lightEnum = u_radiiAndDynamicAtmosphereColor.z;
-    vec3 lightDirection =
-        positionWC * float(lightEnum == 0.0) +
-        czm_lightDirectionWC * float(lightEnum == 1.0) +
-        czm_sunDirectionWC * float(lightEnum == 2.0);
-    return normalize(lightDirection);
-}
-
 void computeAtmosphereScattering(vec3 positionWC, vec3 lightDirection, out vec3 rayleighColor, out vec3 mieColor, out float opacity, out float underTranslucentGlobe)
 {
     float ellipsoidRadiiDifference = czm_ellipsoidRadii.x - czm_ellipsoidRadii.z;
@@ -28,7 +18,7 @@ void computeAtmosphereScattering(vec3 positionWC, vec3 lightDirection, out vec3
     float distanceAdjustModifier = ellipsoidRadiiDifference / 2.0;
     float distanceAdjust = distanceAdjustModifier * clamp((czm_eyeHeight - distanceAdjustMin) / (distanceAdjustMax - distanceAdjustMin), 0.0, 1.0);
 
-    // Since atmosphere scattering assumes the atmosphere is a spherical shell, we compute an inner radius of the atmosphere best fit 
+    // Since atmosphere scattering assumes the atmosphere is a spherical shell, we compute an inner radius of the atmosphere best fit
     // for the position on the ellipsoid.
     float radiusAdjust = (ellipsoidRadiiDifference / 4.0) + distanceAdjust;
     float atmosphereInnerRadius = (length(czm_viewerPositionWC) - czm_eyeHeight) - radiusAdjust;
@@ -46,7 +36,7 @@ void computeAtmosphereScattering(vec3 positionWC, vec3 lightDirection, out vec3
         // Check for intersection with the inner radius of the atmopshere.
         czm_raySegment primaryRayEarthIntersect = czm_raySphereIntersectionInterval(primaryRay, vec3(0.0), atmosphereInnerRadius + radiusAdjust);
         if (primaryRayEarthIntersect.start > 0.0 && primaryRayEarthIntersect.stop > 0.0) {
-            
+
             // Compute position on globe.
             vec3 direction = normalize(positionWC);
             czm_ray ellipsoidRay = czm_ray(positionWC, -direction);
@@ -62,7 +52,7 @@ void computeAtmosphereScattering(vec3 positionWC, vec3 lightDirection, out vec3
             vec3 nearColor = vec3(0.0);
 
             rayleighColor = mix(nearColor, horizonColor, exp(-angle) * opacity);
-            
+
             // Set the traslucent flag to avoid alpha adjustment in computeFinalColor funciton.
             underTranslucentGlobe = 1.0;
             return;

package/Source/Shaders/SkyAtmosphereCommon.js

@@ -9,16 +9,6 @@ export default "float interpolateByDistance(vec4 nearFarScalar, float distance)\
     return mix(startValue, endValue, t);\n\
 }\n\
 \n\
-vec3 getLightDirection(vec3 positionWC)\n\
-{\n\
-    float lightEnum = u_radiiAndDynamicAtmosphereColor.z;\n\
-    vec3 lightDirection =\n\
-        positionWC * float(lightEnum == 0.0) +\n\
-        czm_lightDirectionWC * float(lightEnum == 1.0) +\n\
-        czm_sunDirectionWC * float(lightEnum == 2.0);\n\
-    return normalize(lightDirection);\n\
-}\n\
-\n\
 void computeAtmosphereScattering(vec3 positionWC, vec3 lightDirection, out vec3 rayleighColor, out vec3 mieColor, out float opacity, out float underTranslucentGlobe)\n\
 {\n\
     float ellipsoidRadiiDifference = czm_ellipsoidRadii.x - czm_ellipsoidRadii.z;\n\
@@ -29,7 +19,7 @@ void computeAtmosphereScattering(vec3 positionWC, vec3 lightDirection, out vec3
     float distanceAdjustModifier = ellipsoidRadiiDifference / 2.0;\n\
     float distanceAdjust = distanceAdjustModifier * clamp((czm_eyeHeight - distanceAdjustMin) / (distanceAdjustMax - distanceAdjustMin), 0.0, 1.0);\n\
 \n\
-    // Since atmosphere scattering assumes the atmosphere is a spherical shell, we compute an inner radius of the atmosphere best fit \n\
+    // Since atmosphere scattering assumes the atmosphere is a spherical shell, we compute an inner radius of the atmosphere best fit\n\
     // for the position on the ellipsoid.\n\
     float radiusAdjust = (ellipsoidRadiiDifference / 4.0) + distanceAdjust;\n\
     float atmosphereInnerRadius = (length(czm_viewerPositionWC) - czm_eyeHeight) - radiusAdjust;\n\
@@ -47,7 +37,7 @@ void computeAtmosphereScattering(vec3 positionWC, vec3 lightDirection, out vec3
         // Check for intersection with the inner radius of the atmopshere.\n\
         czm_raySegment primaryRayEarthIntersect = czm_raySphereIntersectionInterval(primaryRay, vec3(0.0), atmosphereInnerRadius + radiusAdjust);\n\
         if (primaryRayEarthIntersect.start > 0.0 && primaryRayEarthIntersect.stop > 0.0) {\n\
-            \n\
+\n\
             // Compute position on globe.\n\
             vec3 direction = normalize(positionWC);\n\
             czm_ray ellipsoidRay = czm_ray(positionWC, -direction);\n\
@@ -63,7 +53,7 @@ void computeAtmosphereScattering(vec3 positionWC, vec3 lightDirection, out vec3
             vec3 nearColor = vec3(0.0);\n\
 \n\
             rayleighColor = mix(nearColor, horizonColor, exp(-angle) * opacity);\n\
-            \n\
+\n\
             // Set the traslucent flag to avoid alpha adjustment in computeFinalColor funciton.\n\
             underTranslucentGlobe = 1.0;\n\
             return;\n\

package/Source/Shaders/SkyAtmosphereFS.glsl

@@ -11,8 +11,9 @@ in float v_translucent;
 
 void main (void)
 {
-    vec3 lightDirection = getLightDirection(v_outerPositionWC);
-   
+    float lightEnum = u_radiiAndDynamicAtmosphereColor.z;
+    vec3 lightDirection = czm_getDynamicAtmosphereLightDirection(v_outerPositionWC, lightEnum);
+
     vec3 mieColor;
     vec3 rayleighColor;
     float opacity;
@@ -42,14 +43,8 @@ void main (void)
     #endif
 
     #ifdef COLOR_CORRECT
-        // Convert rgb color to hsb
-        vec3 hsb = czm_RGBToHSB(color.rgb);
-        // Perform hsb shift
-        hsb.x += u_hsbShift.x; // hue
-        hsb.y = clamp(hsb.y + u_hsbShift.y, 0.0, 1.0); // saturation
-        hsb.z = hsb.z > czm_epsilon7 ? hsb.z + u_hsbShift.z : 0.0; // brightness
-        // Convert shifted hsb back to rgb
-        color.rgb = czm_HSBToRGB(hsb);
+        const bool ignoreBlackPixels = true;
+        color.rgb = czm_applyHSBShift(color.rgb, u_hsbShift, ignoreBlackPixels);
     #endif
 
     // For the parts of the sky atmosphere that are not behind a translucent globe,

package/Source/Shaders/SkyAtmosphereFS.js

@@ -12,8 +12,9 @@ in float v_translucent;\n\
 \n\
 void main (void)\n\
 {\n\
-    vec3 lightDirection = getLightDirection(v_outerPositionWC);\n\
-   \n\
+    float lightEnum = u_radiiAndDynamicAtmosphereColor.z;\n\
+    vec3 lightDirection = czm_getDynamicAtmosphereLightDirection(v_outerPositionWC, lightEnum);\n\
+\n\
     vec3 mieColor;\n\
     vec3 rayleighColor;\n\
     float opacity;\n\
@@ -43,14 +44,8 @@ void main (void)\n\
     #endif\n\
 \n\
     #ifdef COLOR_CORRECT\n\
-        // Convert rgb color to hsb\n\
-        vec3 hsb = czm_RGBToHSB(color.rgb);\n\
-        // Perform hsb shift\n\
-        hsb.x += u_hsbShift.x; // hue\n\
-        hsb.y = clamp(hsb.y + u_hsbShift.y, 0.0, 1.0); // saturation\n\
-        hsb.z = hsb.z > czm_epsilon7 ? hsb.z + u_hsbShift.z : 0.0; // brightness\n\
-        // Convert shifted hsb back to rgb\n\
-        color.rgb = czm_HSBToRGB(hsb);\n\
+        const bool ignoreBlackPixels = true;\n\
+        color.rgb = czm_applyHSBShift(color.rgb, u_hsbShift, ignoreBlackPixels);\n\
     #endif\n\
 \n\
     // For the parts of the sky atmosphere that are not behind a translucent globe,\n\

package/Source/Shaders/SkyAtmosphereVS.glsl

@@ -12,7 +12,8 @@ out float v_translucent;
 void main(void)
 {
     vec4 positionWC = czm_model * position;
-    vec3 lightDirection = getLightDirection(positionWC.xyz);
+    float lightEnum = u_radiiAndDynamicAtmosphereColor.z;
+    vec3 lightDirection = czm_getDynamicAtmosphereLightDirection(positionWC.xyz, lightEnum);
 
     #ifndef PER_FRAGMENT_ATMOSPHERE
         computeAtmosphereScattering(
@@ -24,7 +25,7 @@ void main(void)
             v_translucent
         );
     #endif
-    
+
     v_outerPositionWC = positionWC.xyz;
     gl_Position = czm_modelViewProjection * position;
 }

package/Source/Shaders/SkyAtmosphereVS.js

@@ -13,7 +13,8 @@ out float v_translucent;\n\
 void main(void)\n\
 {\n\
     vec4 positionWC = czm_model * position;\n\
-    vec3 lightDirection = getLightDirection(positionWC.xyz);\n\
+    float lightEnum = u_radiiAndDynamicAtmosphereColor.z;\n\
+    vec3 lightDirection = czm_getDynamicAtmosphereLightDirection(positionWC.xyz, lightEnum);\n\
 \n\
     #ifndef PER_FRAGMENT_ATMOSPHERE\n\
         computeAtmosphereScattering(\n\
@@ -25,7 +26,7 @@ void main(void)\n\
             v_translucent\n\
         );\n\
     #endif\n\
-    \n\
+\n\
     v_outerPositionWC = positionWC.xyz;\n\
     gl_Position = czm_modelViewProjection * position;\n\
 }\n\