@cesium/engine 21.0.1 → 22.1.0

package/Source/Shaders/BillboardCollectionVS.js

@@ -1,5 +1,6 @@
 //This file is automatically rebuilt by the Cesium build process.
-export default "#ifdef INSTANCED\n\
+export default "uniform float u_threePointDepthTestDistance;\n\
+#ifdef INSTANCED\n\
 in vec2 direction;\n\
 #endif\n\
 in vec4 positionHighAndScale;\n\
@@ -10,10 +11,10 @@ in vec4 compressedAttribute2;                       // label horizontal origin,
 in vec4 eyeOffset;                                  // eye offset in meters, 4 bytes free (texture range)\n\
 in vec4 scaleByDistance;                            // near, nearScale, far, farScale\n\
 in vec4 pixelOffsetScaleByDistance;                 // near, nearScale, far, farScale\n\
-in vec4 compressedAttribute3;                       // distance display condition near, far, disableDepthTestDistance, dimensions\n\
+in vec4 compressedAttribute3;                       // distance display condition near, far, disableDepthTestDistanceSq, dimensions\n\
 in vec2 sdf;                                        // sdf outline color (rgb) and width (w)\n\
 in float splitDirection;                            // splitDirection\n\
-#if defined(VERTEX_DEPTH_CHECK) || defined(FRAGMENT_DEPTH_CHECK)\n\
+#if defined(VS_THREE_POINT_DEPTH_CHECK) || defined(FS_THREE_POINT_DEPTH_CHECK)\n\
 in vec4 textureCoordinateBoundsOrLabelTranslate;    // the min and max x and y values for the texture coordinates\n\
 #endif\n\
 #ifdef VECTOR_TILE\n\
@@ -21,12 +22,12 @@ in float a_batchId;\n\
 #endif\n\
 \n\
 out vec2 v_textureCoordinates;\n\
-#ifdef FRAGMENT_DEPTH_CHECK\n\
+#ifdef FS_THREE_POINT_DEPTH_CHECK\n\
 out vec4 v_textureCoordinateBounds;\n\
 out vec4 v_originTextureCoordinateAndTranslate;\n\
-out vec4 v_compressed;                                 // x: eyeDepth, y: applyTranslate & enableDepthCheck, z: dimensions, w: imageSize\n\
 out mat2 v_rotationMatrix;\n\
 #endif\n\
+out vec4 v_compressed;                                 // x: eyeDepth, y: applyTranslate & enableDepthCheck, z: dimensions, w: imageSize\n\
 \n\
 out vec4 v_pickColor;\n\
 out vec4 v_color;\n\
@@ -93,7 +94,7 @@ vec4 addScreenSpaceOffset(vec4 positionEC, vec2 imageSize, float scale, vec2 dir
     return positionEC;\n\
 }\n\
 \n\
-#ifdef VERTEX_DEPTH_CHECK\n\
+#ifdef VS_THREE_POINT_DEPTH_CHECK\n\
 float getGlobeDepth(vec4 positionEC)\n\
 {\n\
     vec4 posWC = czm_eyeToWindowCoordinates(positionEC);\n\
@@ -138,7 +139,7 @@ void main()\n\
     origin.y = floor(compressed * SHIFT_RIGHT3);\n\
     compressed -= origin.y * SHIFT_LEFT3;\n\
 \n\
-#ifdef FRAGMENT_DEPTH_CHECK\n\
+#ifdef FS_THREE_POINT_DEPTH_CHECK\n\
     vec2 depthOrigin = origin.xy;\n\
 #endif\n\
     origin -= vec2(1.0);\n\
@@ -177,10 +178,10 @@ void main()\n\
 \n\
     vec2 imageSize = vec2(floor(temp), temp2);\n\
 \n\
-#ifdef FRAGMENT_DEPTH_CHECK\n\
+#ifdef FS_THREE_POINT_DEPTH_CHECK\n\
     float labelHorizontalOrigin = floor(compressedAttribute2.w - (temp2 * SHIFT_LEFT2));\n\
     float applyTranslate = 0.0;\n\
-    if (labelHorizontalOrigin != 0.0) // is a billboard, so set apply translate to false\n\
+    if (labelHorizontalOrigin != 0.0) // is a label, so set apply translate to true\n\
     {\n\
         applyTranslate = 1.0;\n\
         labelHorizontalOrigin -= 2.0;\n\
@@ -201,7 +202,7 @@ void main()\n\
     translucencyByDistance.w = ((temp - floor(temp)) * SHIFT_LEFT8) / 255.0;\n\
 #endif\n\
 \n\
-#if defined(VERTEX_DEPTH_CHECK) || defined(FRAGMENT_DEPTH_CHECK)\n\
+#if defined(VS_THREE_POINT_DEPTH_CHECK) || defined(FS_THREE_POINT_DEPTH_CHECK)\n\
     temp = compressedAttribute3.w;\n\
     temp = temp * SHIFT_RIGHT12;\n\
 \n\
@@ -250,10 +251,6 @@ void main()\n\
 \n\
     vec4 p = czm_translateRelativeToEye(positionHigh, positionLow);\n\
     vec4 positionEC = czm_modelViewRelativeToEye * p;\n\
-\n\
-#if defined(FRAGMENT_DEPTH_CHECK) || defined(VERTEX_DEPTH_CHECK)\n\
-    float eyeDepth = positionEC.z;\n\
-#endif\n\
 \n\
     positionEC = czm_eyeOffset(positionEC, eyeOffset.xyz);\n\
     positionEC.xyz *= show;\n\
@@ -312,35 +309,56 @@ void main()\n\
     mat2 rotationMatrix;\n\
     float mpp;\n\
 \n\
+    float enableDepthCheck = 1.0;\n\
 #ifdef DISABLE_DEPTH_DISTANCE\n\
-    float disableDepthTestDistance = compressedAttribute3.z;\n\
+    float disableDepthTestDistanceSq = compressedAttribute3.z;\n\
+    if (disableDepthTestDistanceSq == 0.0 && czm_minimumDisableDepthTestDistance != 0.0)\n\
+    {\n\
+        disableDepthTestDistanceSq = czm_minimumDisableDepthTestDistance;\n\
+    }\n\
+\n\
+    if (lengthSq < disableDepthTestDistanceSq || disableDepthTestDistanceSq < 0.0)\n\
+    {\n\
+        enableDepthCheck = 0.0;\n\
+    }\n\
 #endif\n\
 \n\
-#ifdef VERTEX_DEPTH_CHECK\n\
-if (lengthSq < disableDepthTestDistance) {\n\
-    float depthsilon = 10.0;\n\
+    v_compressed.y = enableDepthCheck;\n\
 \n\
-    vec2 labelTranslate = textureCoordinateBoundsOrLabelTranslate.xy;\n\
-    vec4 pEC1 = addScreenSpaceOffset(positionEC, dimensions, scale, vec2(0.0), origin, labelTranslate, pixelOffset, alignedAxis, validAlignedAxis, rotation, sizeInMeters, rotationMatrix, mpp);\n\
+#ifdef VS_THREE_POINT_DEPTH_CHECK\n\
+if (lengthSq < (u_threePointDepthTestDistance * u_threePointDepthTestDistance) && (enableDepthCheck == 1.0)) {\n\
+    float depthsilon = 10.0;\n\
+    vec2 depthOrigin;\n\
+    // Horizontal origin for labels comes from a special attribute. If that value is 0, this is a billboard - use the regular origin. \n\
+    // Otherwise, transform the label origin to -1, 0, 1 (right, center, left).\n\
+    depthOrigin.x = floor(compressedAttribute2.w - (temp2 * SHIFT_LEFT2));\n\
+    depthOrigin.x = czm_branchFreeTernary(depthOrigin.x == 0.0, origin.x, depthOrigin.x - 2.0);\n\
+    depthOrigin.y = origin.y;\n\
+\n\
+    vec4 pEC1 = addScreenSpaceOffset(positionEC, dimensions, scale, vec2(0.0), depthOrigin, vec2(0.0), pixelOffset, alignedAxis, validAlignedAxis, rotation, sizeInMeters, rotationMatrix, mpp);\n\
     float globeDepth1 = getGlobeDepth(pEC1);\n\
 \n\
     if (globeDepth1 != 0.0 && pEC1.z + depthsilon < globeDepth1)\n\
     {\n\
-        vec4 pEC2 = addScreenSpaceOffset(positionEC, dimensions, scale, vec2(0.0, 1.0), origin, labelTranslate, pixelOffset, alignedAxis, validAlignedAxis, rotation, sizeInMeters, rotationMatrix, mpp);\n\
+        vec4 pEC2 = addScreenSpaceOffset(positionEC, dimensions, scale, vec2(0.0, 1.0), depthOrigin, vec2(0.0), pixelOffset, alignedAxis, validAlignedAxis, rotation, sizeInMeters, rotationMatrix, mpp);\n\
         float globeDepth2 = getGlobeDepth(pEC2);\n\
 \n\
         if (globeDepth2 != 0.0 && pEC2.z + depthsilon < globeDepth2)\n\
         {\n\
-            vec4 pEC3 = addScreenSpaceOffset(positionEC, dimensions, scale, vec2(1.0), origin, labelTranslate, pixelOffset, alignedAxis, validAlignedAxis, rotation, sizeInMeters, rotationMatrix, mpp);\n\
+            vec4 pEC3 = addScreenSpaceOffset(positionEC, dimensions, scale, vec2(1.0), depthOrigin, vec2(0.0), pixelOffset, alignedAxis, validAlignedAxis, rotation, sizeInMeters, rotationMatrix, mpp);\n\
             float globeDepth3 = getGlobeDepth(pEC3);\n\
             if (globeDepth3 != 0.0 && pEC3.z + depthsilon < globeDepth3)\n\
             {\n\
+                // \"Discard\" this vertex, as three key points fail depth test.\n\
                 positionEC.xyz = vec3(0.0);\n\
             }\n\
         }\n\
     }\n\
 }\n\
 #endif\n\
+    // Write out the eyespace depth before applying the screen space offset, but after potentially \"discarding\" the vertex\n\
+    // by setting its eyespace position to zero, via the three-point depth test above.\n\
+    v_compressed.x = positionEC.z;\n\
 \n\
     positionEC = addScreenSpaceOffset(positionEC, imageSize, scale, direction, origin, translate, pixelOffset, alignedAxis, validAlignedAxis, rotation, sizeInMeters, rotationMatrix, mpp);\n\
     gl_Position = czm_projection * positionEC;\n\
@@ -351,17 +369,14 @@ if (lengthSq < disableDepthTestDistance) {\n\
 #endif\n\
 \n\
 #ifdef DISABLE_DEPTH_DISTANCE\n\
-    if (disableDepthTestDistance == 0.0 && czm_minimumDisableDepthTestDistance != 0.0)\n\
-    {\n\
-        disableDepthTestDistance = czm_minimumDisableDepthTestDistance;\n\
-    }\n\
 \n\
-    if (disableDepthTestDistance != 0.0)\n\
+    if (disableDepthTestDistanceSq != 0.0)\n\
     {\n\
         // Don't try to \"multiply both sides\" by w.  Greater/less-than comparisons won't work for negative values of w.\n\
         float zclip = gl_Position.z / gl_Position.w;\n\
         bool clipped = (zclip < -1.0 || zclip > 1.0);\n\
-        if (!clipped && (disableDepthTestDistance < 0.0 || (lengthSq > 0.0 && lengthSq < disableDepthTestDistance)))\n\
+        // disableDepthTestDistanceSq can be less than zero if it's explicitly set to -1 in JS (as a sentinel value equivalent to infinity)\n\
+        if (!clipped && (disableDepthTestDistanceSq < 0.0 || (lengthSq > 0.0 && lengthSq < disableDepthTestDistanceSq)))\n\
         {\n\
             // Position z on the near plane.\n\
             gl_Position.z = -gl_Position.w;\n\
@@ -372,7 +387,7 @@ if (lengthSq < disableDepthTestDistance) {\n\
     }\n\
 #endif\n\
 \n\
-#ifdef FRAGMENT_DEPTH_CHECK\n\
+#ifdef FS_THREE_POINT_DEPTH_CHECK\n\
     if (sizeInMeters) {\n\
         translate /= mpp;\n\
         dimensions /= mpp;\n\
@@ -384,21 +399,13 @@ if (lengthSq < disableDepthTestDistance) {\n\
 #else\n\
     v_rotationMatrix = mat2(1.0, 0.0, 0.0, 1.0);\n\
 #endif\n\
-\n\
-    float enableDepthCheck = 0.0;\n\
-    if (lengthSq < disableDepthTestDistance)\n\
-    {\n\
-        enableDepthCheck = 1.0;\n\
-    }\n\
-\n\
     float dw = floor(clamp(dimensions.x, 0.0, SHIFT_LEFT12));\n\
     float dh = floor(clamp(dimensions.y, 0.0, SHIFT_LEFT12));\n\
 \n\
     float iw = floor(clamp(imageSize.x, 0.0, SHIFT_LEFT12));\n\
     float ih = floor(clamp(imageSize.y, 0.0, SHIFT_LEFT12));\n\
 \n\
-    v_compressed.x = eyeDepth;\n\
-    v_compressed.y = applyTranslate * SHIFT_LEFT1 + enableDepthCheck;\n\
+    v_compressed.y += applyTranslate * SHIFT_LEFT1;\n\
     v_compressed.z = dw * SHIFT_LEFT12 + dh;\n\
     v_compressed.w = iw * SHIFT_LEFT12 + ih;\n\
     v_originTextureCoordinateAndTranslate.xy = depthOrigin;\n\

package/Source/Shaders/Builtin/Constants/passCesium3DTile.glsl

@@ -6,4 +6,4 @@
  *
  * @see czm_pass
  */
-const float czm_passCesium3DTile = 4.0;
+const float czm_passCesium3DTile = 5.0;

package/Source/Shaders/Builtin/Constants/passCesium3DTile.js

@@ -7,5 +7,5 @@ export default "/**\n\
  *\n\
  * @see czm_pass\n\
  */\n\
-const float czm_passCesium3DTile = 4.0;\n\
+const float czm_passCesium3DTile = 5.0;\n\
 ";

package/Source/Shaders/Builtin/Constants/passCesium3DTileClassification.glsl

@@ -6,4 +6,4 @@
  *
  * @see czm_pass
  */
-const float czm_passCesium3DTileClassification = 5.0;
+const float czm_passCesium3DTileClassification = 6.0;

package/Source/Shaders/Builtin/Constants/passCesium3DTileClassification.js

@@ -7,5 +7,5 @@ export default "/**\n\
  *\n\
  * @see czm_pass\n\
  */\n\
-const float czm_passCesium3DTileClassification = 5.0;\n\
+const float czm_passCesium3DTileClassification = 6.0;\n\
 ";

package/Source/Shaders/Builtin/Constants/passCesium3DTileClassificationIgnoreShow.glsl

@@ -6,4 +6,4 @@
  *
  * @see czm_pass
  */
-const float czm_passCesium3DTileClassificationIgnoreShow = 6.0;
+const float czm_passCesium3DTileClassificationIgnoreShow = 7.0;

package/Source/Shaders/Builtin/Constants/passCesium3DTileClassificationIgnoreShow.js

@@ -7,5 +7,5 @@ export default "/**\n\
  *\n\
  * @see czm_pass\n\
  */\n\
-const float czm_passCesium3DTileClassificationIgnoreShow = 6.0;\n\
+const float czm_passCesium3DTileClassificationIgnoreShow = 7.0;\n\
 ";

package/Source/Shaders/Builtin/Constants/passCesium3DTileEdges.glsl

@@ -0,0 +1,10 @@
+/**
+ * The automatic GLSL constant for {@link Pass#CESIUM_3D_TILE_EDGES}
+ *
+ * @name czm_passCesium3DTileEdges
+ * @glslConstant
+ *
+ * @see czm_pass
+ */
+const float czm_passCesium3DTileEdges = 4.0;
+

package/Source/Shaders/Builtin/Constants/passCesium3DTileEdges.js

@@ -0,0 +1,12 @@
+//This file is automatically rebuilt by the Cesium build process.
+export default "/**\n\
+ * The automatic GLSL constant for {@link Pass#CESIUM_3D_TILE_EDGES}\n\
+ *\n\
+ * @name czm_passCesium3DTileEdges\n\
+ * @glslConstant\n\
+ *\n\
+ * @see czm_pass\n\
+ */\n\
+const float czm_passCesium3DTileEdges = 4.0;\n\
+\n\
+";

package/Source/Shaders/Builtin/Constants/passGaussianSplats.glsl

@@ -6,4 +6,4 @@
  *
  * @see czm_pass
  */
-const float czm_passGaussianSplats = 10.0;
+const float czm_passGaussianSplats = 11.0;

package/Source/Shaders/Builtin/Constants/passGaussianSplats.js

@@ -7,5 +7,5 @@ export default "/**\n\
  *\n\
  * @see czm_pass\n\
  */\n\
-const float czm_passGaussianSplats = 10.0;\n\
+const float czm_passGaussianSplats = 11.0;\n\
 ";

package/Source/Shaders/Builtin/Constants/passOpaque.glsl

@@ -6,4 +6,4 @@
  *
  * @see czm_pass
  */
-const float czm_passOpaque = 7.0;
+const float czm_passOpaque = 8.0;

package/Source/Shaders/Builtin/Constants/passOpaque.js

@@ -7,5 +7,5 @@ export default "/**\n\
  *\n\
  * @see czm_pass\n\
  */\n\
-const float czm_passOpaque = 7.0;\n\
+const float czm_passOpaque = 8.0;\n\
 ";

package/Source/Shaders/Builtin/Constants/passOverlay.glsl

@@ -6,4 +6,4 @@
  *
  * @see czm_pass
  */
-const float czm_passOverlay = 11.0;
+const float czm_passOverlay = 12.0;

package/Source/Shaders/Builtin/Constants/passOverlay.js

@@ -7,5 +7,5 @@ export default "/**\n\
  *\n\
  * @see czm_pass\n\
  */\n\
-const float czm_passOverlay = 11.0;\n\
+const float czm_passOverlay = 12.0;\n\
 ";

package/Source/Shaders/Builtin/Constants/passTranslucent.glsl

@@ -6,4 +6,4 @@
  *
  * @see czm_pass
  */
-const float czm_passTranslucent = 8.0;
+const float czm_passTranslucent = 9.0;

package/Source/Shaders/Builtin/Constants/passTranslucent.js

@@ -7,5 +7,5 @@ export default "/**\n\
  *\n\
  * @see czm_pass\n\
  */\n\
-const float czm_passTranslucent = 8.0;\n\
+const float czm_passTranslucent = 9.0;\n\
 ";

package/Source/Shaders/Builtin/Constants/passVoxels.glsl

@@ -6,4 +6,4 @@
  *
  * @see czm_pass
  */
-const float czm_passVoxels = 9.0;
+const float czm_passVoxels = 10.0;

package/Source/Shaders/Builtin/Constants/passVoxels.js

@@ -7,5 +7,5 @@ export default "/**\n\
  *\n\
  * @see czm_pass\n\
  */\n\
-const float czm_passVoxels = 9.0;\n\
+const float czm_passVoxels = 10.0;\n\
 ";

package/Source/Shaders/Builtin/CzmBuiltins.js

@@ -14,6 +14,7 @@ import czm_oneOverTwoPi from './Constants/oneOverTwoPi.js'
 import czm_passCesium3DTile from './Constants/passCesium3DTile.js'
 import czm_passCesium3DTileClassification from './Constants/passCesium3DTileClassification.js'
 import czm_passCesium3DTileClassificationIgnoreShow from './Constants/passCesium3DTileClassificationIgnoreShow.js'
+import czm_passCesium3DTileEdges from './Constants/passCesium3DTileEdges.js'
 import czm_passClassification from './Constants/passClassification.js'
 import czm_passCompute from './Constants/passCompute.js'
 import czm_passEnvironment from './Constants/passEnvironment.js'
@@ -153,6 +154,7 @@ export default {
     czm_passCesium3DTile : czm_passCesium3DTile,
     czm_passCesium3DTileClassification : czm_passCesium3DTileClassification,
     czm_passCesium3DTileClassificationIgnoreShow : czm_passCesium3DTileClassificationIgnoreShow,
+    czm_passCesium3DTileEdges : czm_passCesium3DTileEdges,
     czm_passClassification : czm_passClassification,
     czm_passCompute : czm_passCompute,
     czm_passEnvironment : czm_passEnvironment,

package/Source/Shaders/GlobeVS.glsl

@@ -141,7 +141,7 @@ void main()
 #elif defined(INCLUDE_WEB_MERCATOR_Y)
     float webMercatorT = czm_decompressTextureCoordinates(compressed0.w).x;
     float encodedNormal = 0.0;
-#elif defined(ENABLE_VERTEX_LIGHTING) || defined(GENERATE_POSITION_AND_NORMAL)
+#elif defined(ENABLE_VERTEX_LIGHTING) || defined(GENERATE_POSITION_AND_NORMAL) || defined(APPLY_MATERIAL)
     float webMercatorT = textureCoordinates.y;
     float encodedNormal = compressed0.w;
 #else

package/Source/Shaders/GlobeVS.js

@@ -142,7 +142,7 @@ void main()\n\
 #elif defined(INCLUDE_WEB_MERCATOR_Y)\n\
     float webMercatorT = czm_decompressTextureCoordinates(compressed0.w).x;\n\
     float encodedNormal = 0.0;\n\
-#elif defined(ENABLE_VERTEX_LIGHTING) || defined(GENERATE_POSITION_AND_NORMAL)\n\
+#elif defined(ENABLE_VERTEX_LIGHTING) || defined(GENERATE_POSITION_AND_NORMAL) || defined(APPLY_MATERIAL)\n\
     float webMercatorT = textureCoordinates.y;\n\
     float encodedNormal = compressed0.w;\n\
 #else\n\

package/Source/Shaders/Model/EdgeDetectionStageFS.glsl

@@ -0,0 +1,58 @@
+void edgeDetectionStage(inout vec4 color, inout FeatureIds featureIds) {
+    if (u_isEdgePass) {
+        return;
+    }
+
+    vec2 screenCoord = gl_FragCoord.xy / czm_viewport.zw;
+
+    vec4 edgeColor = texture(czm_edgeColorTexture, screenCoord);
+    vec4 edgeId = texture(czm_edgeIdTexture, screenCoord);
+
+    // Packed window-space depth from edge pass (0..1)
+    float edgeDepthWin = czm_unpackDepth(texture(czm_edgeDepthTexture, screenCoord));
+
+    // Near / far for current frustum
+    float n = czm_currentFrustum.x;
+    float f = czm_currentFrustum.y;
+
+    // geometry depth in eye coordinate
+    vec4 geomEC = czm_windowToEyeCoordinates(gl_FragCoord);
+    float geomDepthLinear = -geomEC.z;
+
+    // Convert edge depth to linear depth
+    float z_ndc_edge = edgeDepthWin * 2.0 - 1.0;
+    float edgeDepthLinear = (2.0 * n * f) / (f + n - z_ndc_edge * (f - n));
+
+    float d = abs(edgeDepthLinear - geomDepthLinear);
+
+    // Adaptive epsilon using linear depth fwidth for robustness
+    float pixelStepLinear = fwidth(geomDepthLinear);
+    float rel = geomDepthLinear * 0.0005;
+    float eps = max(n * 1e-4, max(pixelStepLinear * 1.5, rel));
+
+    // If Edge isn't behind any geometry and the pixel has edge data
+    if (d < eps && edgeId.r > 0.0) {
+#ifdef HAS_EDGE_FEATURE_ID
+        float edgeFeatureId    = edgeId.g;
+        float currentFeatureId = float(featureIds.featureId_0);
+#endif
+        float globeDepth = czm_unpackDepth(texture(czm_globeDepthTexture, screenCoord));
+        // Background / sky / globe: always show edge
+        bool isBackground = geomDepthLinear > globeDepth;
+        bool drawEdge = isBackground;
+
+#ifdef HAS_EDGE_FEATURE_ID
+        bool hasEdgeFeature = edgeFeatureId > 0.0;
+        bool hasCurrentFeature = currentFeatureId > 0.0;
+        bool featuresMatch = edgeFeatureId == currentFeatureId;
+
+        drawEdge = drawEdge || !hasEdgeFeature || !hasCurrentFeature || featuresMatch;
+#else
+        drawEdge = true;
+#endif
+
+        if (drawEdge) {
+            color = edgeColor;
+        }
+    }
+}

package/Source/Shaders/Model/EdgeDetectionStageFS.js

@@ -0,0 +1,60 @@
+//This file is automatically rebuilt by the Cesium build process.
+export default "void edgeDetectionStage(inout vec4 color, inout FeatureIds featureIds) {\n\
+    if (u_isEdgePass) {\n\
+        return;\n\
+    }\n\
+\n\
+    vec2 screenCoord = gl_FragCoord.xy / czm_viewport.zw;\n\
+\n\
+    vec4 edgeColor = texture(czm_edgeColorTexture, screenCoord);\n\
+    vec4 edgeId = texture(czm_edgeIdTexture, screenCoord);\n\
+\n\
+    // Packed window-space depth from edge pass (0..1)\n\
+    float edgeDepthWin = czm_unpackDepth(texture(czm_edgeDepthTexture, screenCoord));\n\
+\n\
+    // Near / far for current frustum\n\
+    float n = czm_currentFrustum.x;\n\
+    float f = czm_currentFrustum.y;\n\
+\n\
+    // geometry depth in eye coordinate\n\
+    vec4 geomEC = czm_windowToEyeCoordinates(gl_FragCoord);\n\
+    float geomDepthLinear = -geomEC.z;\n\
+\n\
+    // Convert edge depth to linear depth\n\
+    float z_ndc_edge = edgeDepthWin * 2.0 - 1.0;\n\
+    float edgeDepthLinear = (2.0 * n * f) / (f + n - z_ndc_edge * (f - n));\n\
+\n\
+    float d = abs(edgeDepthLinear - geomDepthLinear);\n\
+\n\
+    // Adaptive epsilon using linear depth fwidth for robustness\n\
+    float pixelStepLinear = fwidth(geomDepthLinear);\n\
+    float rel = geomDepthLinear * 0.0005;\n\
+    float eps = max(n * 1e-4, max(pixelStepLinear * 1.5, rel));\n\
+\n\
+    // If Edge isn't behind any geometry and the pixel has edge data\n\
+    if (d < eps && edgeId.r > 0.0) {\n\
+#ifdef HAS_EDGE_FEATURE_ID\n\
+        float edgeFeatureId    = edgeId.g;\n\
+        float currentFeatureId = float(featureIds.featureId_0);\n\
+#endif\n\
+        float globeDepth = czm_unpackDepth(texture(czm_globeDepthTexture, screenCoord));\n\
+        // Background / sky / globe: always show edge\n\
+        bool isBackground = geomDepthLinear > globeDepth;\n\
+        bool drawEdge = isBackground;\n\
+\n\
+#ifdef HAS_EDGE_FEATURE_ID\n\
+        bool hasEdgeFeature = edgeFeatureId > 0.0;\n\
+        bool hasCurrentFeature = currentFeatureId > 0.0;\n\
+        bool featuresMatch = edgeFeatureId == currentFeatureId;\n\
+\n\
+        drawEdge = drawEdge || !hasEdgeFeature || !hasCurrentFeature || featuresMatch;\n\
+#else\n\
+        drawEdge = true;\n\
+#endif\n\
+\n\
+        if (drawEdge) {\n\
+            color = edgeColor;\n\
+        }\n\
+    }\n\
+}\n\
+";

package/Source/Shaders/Model/EdgeVisibilityStageFS.glsl

@@ -0,0 +1,83 @@
+// CESIUM_REDIRECTED_COLOR_OUTPUT flag is used to avoid color attachment conflicts
+// when shaders are processed by different rendering pipelines (e.g., OIT).
+// Only declare MRT outputs when not in a derived shader context.
+#if defined(HAS_EDGE_VISIBILITY_MRT) && !defined(CESIUM_REDIRECTED_COLOR_OUTPUT)
+layout(location = 1) out vec4 out_id;        // edge id / metadata
+layout(location = 2) out vec4 out_edgeDepth; // packed depth
+#endif
+
+void edgeVisibilityStage(inout vec4 color, inout FeatureIds featureIds)
+{
+#ifdef HAS_EDGE_VISIBILITY
+
+    if (!u_isEdgePass) {
+        return;
+    }
+    
+    float edgeTypeInt = v_edgeType * 255.0;
+    
+    // Color code different edge types
+    vec4 edgeColor = vec4(0.0);
+    
+    if (edgeTypeInt < 0.5) { // HIDDEN (0)
+        edgeColor = vec4(0.0, 0.0, 0.0, 0.0); // Transparent for hidden edges
+    }
+    else if (edgeTypeInt > 0.5 && edgeTypeInt < 1.5) { // SILHOUETTE (1) - Conditional visibility
+        // Proper silhouette detection using face normals
+        vec3 normalA = normalize(v_faceNormalAView);
+        vec3 normalB = normalize(v_faceNormalBView);
+        
+        // Calculate view direction using existing eye-space position varying (v_positionEC)
+        vec3 viewDir = -normalize(v_positionEC);
+        
+        // Calculate dot products to determine triangle facing
+        float dotA = dot(normalA, viewDir);
+        float dotB = dot(normalB, viewDir);
+        
+        const float eps = 1e-3;
+        bool frontA = dotA > eps;
+        bool backA  = dotA < -eps;
+        bool frontB = dotB > eps;
+        bool backB  = dotB < -eps;
+        
+        // True silhouette: one triangle front-facing, other back-facing
+        bool oppositeFacing = (frontA && backB) || (backA && frontB);
+        
+        // Exclude edges where both triangles are nearly grazing (perpendicular to view)
+        // This handles the top-view cylinder case where both normals are ~horizontal
+        bool bothNearGrazing = (abs(dotA) <= eps && abs(dotB) <= eps);
+        
+        if (!(oppositeFacing && !bothNearGrazing)) {
+            discard; // Not a true silhouette edge
+        } else {
+            // True silhouette
+            edgeColor = vec4(1.0, 0.0, 0.0, 1.0);
+        }
+    }
+    else if (edgeTypeInt > 1.5 && edgeTypeInt < 2.5) { // HARD (2) - BRIGHT GREEN
+        edgeColor = vec4(0.0, 1.0, 0.0, 1.0); // Extra bright green
+    }
+    else if (edgeTypeInt > 2.5 && edgeTypeInt < 3.5) { // REPEATED (3)
+        edgeColor = vec4(0.0, 0.0, 1.0, 1.0);
+    } else {
+        edgeColor = vec4(0.0, 0.0, 0.0, 0.0);
+    }
+
+    // Temporary color: white
+    edgeColor = vec4(1.0, 1.0, 1.0, 1.0);
+    color = edgeColor;
+    
+    #if defined(HAS_EDGE_VISIBILITY_MRT) && !defined(CESIUM_REDIRECTED_COLOR_OUTPUT)
+        // Write edge metadata
+        out_id = vec4(0.0);
+        out_id.r = edgeTypeInt;                    // Edge type (0-3)
+#ifdef HAS_EDGE_FEATURE_ID
+        out_id.g = float(featureIds.featureId_0); // Feature ID if available
+#else
+        out_id.g = 0.0;
+#endif
+        // Pack depth into separate MRT attachment
+        out_edgeDepth = czm_packDepth(gl_FragCoord.z);
+    #endif
+#endif
+}