deeptwins-engine-3d 0.1.49 → 0.1.50

package/dist/esm/drawCommand/viewShed/config.d.ts

@@ -0,0 +1,20 @@
+/**
+ * 公共着色器
+ * @returns {ShaderSource}
+ */
+export declare const sensorComm = "\nuniform vec4 u_intersectionColor;\nuniform float u_intersectionWidth;\nuniform vec4 u_lineColor;\n\nbool inSensorShadow(vec3 coneVertexWC, vec3 pointWC)\n{    \n    // Diagonal matrix from the unscaled ellipsoid space to the scaled space.\n    vec3 D = czm_ellipsoidInverseRadii;    \n    // Sensor vertex in the scaled ellipsoid space    \n    vec3 q = D * coneVertexWC;    \n    float qMagnitudeSquared = dot(q, q);    \n    float test = qMagnitudeSquared - 1.0;        \n    // Sensor vertex to fragment vector in the ellipsoid's scaled space \n    vec3 temp = D * pointWC - q;    \n    float d = dot(temp, q);        \n    // Behind silhouette plane and inside silhouette cone  \n    return (d < -test) && (d / length(temp) < -sqrt(test));\n}\n\nvec4 getLineColor()\n{ \n    return u_lineColor;\n}\n\nvec4 getIntersectionColor()\n{  \n    return u_intersectionColor;\n}\n\nfloat getIntersectionWidth()\n{ \n    return u_intersectionWidth;\n}\n\nvec2 sensor2dTextureCoordinates(float sensorRadius, vec3 pointMC)\n{   \n    // (s, t) both in the range [0, 1] \n    float t = pointMC.z / sensorRadius;  \n    float s = 1.0 + (atan(pointMC.y, pointMC.x) / czm_twoPi); \n    s = s - floor(s);      \n    return vec2(s, t);\n}\n\n";
+/**
+ * VS着色器
+ * @returns {ShaderSource}
+ */
+export declare const sensorVS = "\nin vec4 position;\nin vec3 normal;\nout vec3 v_position;\nout vec3 v_positionWC;\nout vec3 v_positionEC;\nout vec3 v_normalEC;\nvoid main()\n{   \n    gl_Position = czm_modelViewProjection * position;   \n    v_position = vec3(position);    \n    v_positionWC = (czm_model * position).xyz;   \n    v_positionEC = (czm_modelView * position).xyz;  \n    v_normalEC = czm_normal * normal;\n}";
+/**
+ * FS着色器
+ * @returns {ShaderSource}
+ */
+export declare const sensorFS = "\nuniform bool u_showIntersection;\nuniform bool u_showThroughEllipsoid;\nuniform float u_radius;\nuniform float u_xHalfAngle;\nuniform float u_yHalfAngle;\nuniform float u_normalDirection;\nuniform float u_type;\nin vec3 v_position;\nin vec3 v_positionWC;\nin vec3 v_positionEC;\nin vec3 v_normalEC;\nvec4 getColor(float sensorRadius, vec3 pointEC)\n{   \n    czm_materialInput materialInput;   \n    vec3 pointMC = (czm_inverseModelView * vec4(pointEC, 1.0)).xyz; \n    materialInput.st = sensor2dTextureCoordinates(sensorRadius, pointMC);   \n    materialInput.str = pointMC / sensorRadius;   \n\n    vec3 positionToEyeEC = -v_positionEC;   \n    materialInput.positionToEyeEC = positionToEyeEC;    \n\n    vec3 normalEC = normalize(v_normalEC);  \n    materialInput.normalEC = u_normalDirection * normalEC; \n\n    czm_material material = czm_getMaterial(materialInput);   \n\n     return mix(czm_phong(normalize(positionToEyeEC), material,czm_lightDirectionEC),vec4(material.diffuse, material.alpha), 0.4);\n}\n\nbool isOnBoundary(float value, float epsilon)\n{    \n    float width = getIntersectionWidth();  \n    float tolerance = width * epsilon;\n    #ifdef GL_OES_standard_derivatives    \n    float delta = max(abs(dFdx(value)), abs(dFdy(value))); \n    float pixels = width * delta;   \n    float temp = abs(value);    \n    // There are a couple things going on here.   \n    // First we test the value at the current fragment to see if it is within the tolerance.   \n    // We also want to check if the value of an adjacent pixel is within the tolerance,    \n    // but we don't want to admit points that are obviously not on the surface.    \n    // For example, if we are looking for \"value\" to be close to 0, but value is 1 and the adjacent value is 2,  \n    // then the delta would be 1 and \"temp - delta\" would be \"1 - 1\" which is zero even though neither of   \n    // the points is close to zero.   \n    return temp < tolerance && temp < pixels || (delta < 10.0 * tolerance && temp - delta < tolerance && temp < pixels);\n    #else    \n    return abs(value) < tolerance;\n    #endif\n}\n\n\nvec4 shade(bool isOnBoundary)\n{   \n\n    if (u_showIntersection && isOnBoundary) \n        {      \n            return getIntersectionColor();  \n        }   \n    if(u_type == 1.0)\n        {       \n            return getLineColor(); \n        }   \n\n    return getColor(u_radius, v_positionEC);\n\n}\n                   \nfloat ellipsoidSurfaceFunction(vec3 point)\n{   \n    vec3 scaled = czm_ellipsoidInverseRadii * point;  \n    return dot(scaled, scaled) - 1.0;\n}\n\nvoid main()\n{    \n    vec3 sensorVertexWC = czm_model[3].xyz;      // (0.0, 0.0, 0.0) in model coordinates   \n    vec3 sensorVertexEC = czm_modelView[3].xyz;  // (0.0, 0.0, 0.0) in model coordinates   \n    //vec3 pixDir = normalize(v_position);    \n    float positionX = v_position.x;  \n    float positionY = v_position.y;  \n    float positionZ = v_position.z;  \n    vec3 zDir = vec3(0.0, 0.0, 1.0); \n    vec3 lineX = vec3(positionX, 0 ,positionZ);  \n    vec3 lineY = vec3(0, positionY, positionZ);  \n    float resX = dot(normalize(lineX), zDir);  \n    if(resX < cos(u_xHalfAngle)-0.00001)\n    {       \n        discard; \n    }    \n    float resY = dot(normalize(lineY), zDir);  \n    if(resY < cos(u_yHalfAngle)-0.00001)\n    {    \n        discard;  \n    } \n    float ellipsoidValue = ellipsoidSurfaceFunction(v_positionWC);\n    // Occluded by the ellipsoid?\n    if (!u_showThroughEllipsoid)\n  {   \n    // Discard if in the ellipsoid  \n    // PERFORMANCE_IDEA: A coarse check for ellipsoid intersection could be done on the CPU first. \n    if (ellipsoidValue < 0.0) \n    {       \n        discard;  \n    }   \n    // Discard if in the sensor's shadow   \n    if (inSensorShadow(sensorVertexWC, v_positionWC)) \n    {      \n        discard;  \n    }   \n  }   \n    // Notes: Each surface functions should have an associated tolerance based on the floating point error. \n    bool isOnEllipsoid = isOnBoundary(ellipsoidValue, czm_epsilon3);   \n    //isOnEllipsoid = false;   \n    // if((resX >= 0.8 && resX <= 0.81)||(resY >= 0.8 && resY <= 0.81))\n    // {     \n    //     out_FragColor = shade(false);\n    // }\n    out_FragColor = shade(isOnEllipsoid);\n}\n";
+/**
+ * 扫描面FS着色器
+ * @returns {ShaderSource}
+ */
+export declare const scanPlaneFS = "\nuniform bool u_showIntersection;\nuniform bool u_showThroughEllipsoid;\nuniform float u_radius;\nuniform float u_xHalfAngle;\nuniform float u_yHalfAngle;\nuniform float u_normalDirection;\nuniform vec4 u_color;\nin vec3 v_position;\nin vec3 v_positionWC;\nin vec3 v_positionEC;\nin vec3 v_normalEC;\nvec4 getColor(float sensorRadius, vec3 pointEC)\n{    \n    czm_materialInput materialInput;    \n    vec3 pointMC = (czm_inverseModelView * vec4(pointEC, 1.0)).xyz;  \n    materialInput.st = sensor2dTextureCoordinates(sensorRadius, pointMC);\n    materialInput.str = pointMC / sensorRadius;    \n    vec3 positionToEyeEC = -v_positionEC;    \n    materialInput.positionToEyeEC = positionToEyeEC;    \n    vec3 normalEC = normalize(v_normalEC);   \n    materialInput.normalEC = u_normalDirection * normalEC;    \n    czm_material material = czm_getMaterial(materialInput);   \n    material.diffuse = u_color.rgb;    \n    material.alpha = u_color.a;    \n    return mix(czm_phong(normalize(positionToEyeEC), material,czm_lightDirectionEC), vec4(material.diffuse, material.alpha), 0.4);\n}\n\nbool isOnBoundary(float value, float epsilon)\n{    \n    float width = getIntersectionWidth();   \n    float tolerance = width * epsilon;\n    #ifdef GL_OES_standard_derivatives    \n    float delta = max(abs(dFdx(value)), abs(dFdy(value)));    \n    float pixels = width * delta;    \n    float temp = abs(value);    \n    // There are a couple things going on here.    \n    //  First we test the value at the current fragment to see if it is within the tolerance.    \n    //  We also want to check if the value of an adjacent pixel is within the tolerance,   \n    //  but we don't want to admit points that are obviously not on the surface.    \n    //  For example, if we are looking for value to be close to 0, but value is 1 and the adjacent value is 2,    \n    //  then the delta would be 1 and temp - delta would be 1 - 1 which is zero even though neither of    \n    //  the points is close to zero.    \n    return temp < tolerance && temp < pixels || (delta < 10.0 * tolerance && temp - delta < tolerance && temp < pixels);\n    #else    \n    return abs(value) < tolerance;\n    #endif\n}\n\n\nvec4 shade(bool isOnBoundary)\n    {   \n    if (u_showIntersection && isOnBoundary)  \n    {    \n        return getIntersectionColor(); \n    }   \n    return getColor(u_radius, v_positionEC);\n}\n\n\nfloat ellipsoidSurfaceFunction(vec3 point)\n{    \n    vec3 scaled = czm_ellipsoidInverseRadii * point;   \n    return dot(scaled, scaled) - 1.0;\n}\n\n\nvoid main()\n{  \n    vec3 sensorVertexWC = czm_model[3].xyz;      // (0.0, 0.0, 0.0) in model coordinates   \n    vec3 sensorVertexEC = czm_modelView[3].xyz;  // (0.0, 0.0, 0.0) in model coordinates   \n    //vec3 pixDir = normalize(v_position);    \n    float positionX = v_position.x;    \n    float positionY = v_position.y;    \n    float positionZ = v_position.z;    \n    vec3 zDir = vec3(0.0, 0.0, 1.0);    \n    vec3 lineX = vec3(positionX, 0 ,positionZ);    \n    vec3 lineY = vec3(0, positionY, positionZ);    \n    float resX = dot(normalize(lineX), zDir);    \n    if(resX < cos(u_xHalfAngle) - 0.0001)\n    {       \n        discard;  \n    }    \n    float resY = dot(normalize(lineY), zDir);  \n    if(resY < cos(u_yHalfAngle)- 0.0001)\n    {   \n        discard;  \n    }   \n    float ellipsoidValue = ellipsoidSurfaceFunction(v_positionWC);  \n    // Occluded by the ellipsoid\n    if (!u_showThroughEllipsoid)\n  {   \n    // Discard if in the ellipsoid   \n    // PERFORMANCE_IDEA: A coarse check for ellipsoid intersection could be done on the CPU first.  \n    if (ellipsoidValue < 0.0)  \n    {      \n        discard;  \n    }   \n    // Discard if in the sensor's shadow \n    if (inSensorShadow(sensorVertexWC, v_positionWC)) \n    {    \n        discard;  \n    }   \n  }   \n    // Notes: Each surface functions should have an associated tolerance based on the floating point error.  \n    bool isOnEllipsoid = isOnBoundary(ellipsoidValue, czm_epsilon3);  \n    out_FragColor = shade(isOnEllipsoid);\n}\n    ";

package/dist/esm/drawCommand/viewShed/config.js

@@ -0,0 +1,21 @@
+/**
+ * 公共着色器
+ * @returns {ShaderSource}
+ */
+export var sensorComm = "\nuniform vec4 u_intersectionColor;\nuniform float u_intersectionWidth;\nuniform vec4 u_lineColor;\n\nbool inSensorShadow(vec3 coneVertexWC, vec3 pointWC)\n{    \n    // Diagonal matrix from the unscaled ellipsoid space to the scaled space.\n    vec3 D = czm_ellipsoidInverseRadii;    \n    // Sensor vertex in the scaled ellipsoid space    \n    vec3 q = D * coneVertexWC;    \n    float qMagnitudeSquared = dot(q, q);    \n    float test = qMagnitudeSquared - 1.0;        \n    // Sensor vertex to fragment vector in the ellipsoid's scaled space \n    vec3 temp = D * pointWC - q;    \n    float d = dot(temp, q);        \n    // Behind silhouette plane and inside silhouette cone  \n    return (d < -test) && (d / length(temp) < -sqrt(test));\n}\n\nvec4 getLineColor()\n{ \n    return u_lineColor;\n}\n\nvec4 getIntersectionColor()\n{  \n    return u_intersectionColor;\n}\n\nfloat getIntersectionWidth()\n{ \n    return u_intersectionWidth;\n}\n\nvec2 sensor2dTextureCoordinates(float sensorRadius, vec3 pointMC)\n{   \n    // (s, t) both in the range [0, 1] \n    float t = pointMC.z / sensorRadius;  \n    float s = 1.0 + (atan(pointMC.y, pointMC.x) / czm_twoPi); \n    s = s - floor(s);      \n    return vec2(s, t);\n}\n\n";
+/**
+ * VS着色器
+ * @returns {ShaderSource}
+ */
+export var sensorVS = "\nin vec4 position;\nin vec3 normal;\nout vec3 v_position;\nout vec3 v_positionWC;\nout vec3 v_positionEC;\nout vec3 v_normalEC;\nvoid main()\n{   \n    gl_Position = czm_modelViewProjection * position;   \n    v_position = vec3(position);    \n    v_positionWC = (czm_model * position).xyz;   \n    v_positionEC = (czm_modelView * position).xyz;  \n    v_normalEC = czm_normal * normal;\n}";
+
+/**
+ * FS着色器
+ * @returns {ShaderSource}
+ */
+export var sensorFS = "\nuniform bool u_showIntersection;\nuniform bool u_showThroughEllipsoid;\nuniform float u_radius;\nuniform float u_xHalfAngle;\nuniform float u_yHalfAngle;\nuniform float u_normalDirection;\nuniform float u_type;\nin vec3 v_position;\nin vec3 v_positionWC;\nin vec3 v_positionEC;\nin vec3 v_normalEC;\nvec4 getColor(float sensorRadius, vec3 pointEC)\n{   \n    czm_materialInput materialInput;   \n    vec3 pointMC = (czm_inverseModelView * vec4(pointEC, 1.0)).xyz; \n    materialInput.st = sensor2dTextureCoordinates(sensorRadius, pointMC);   \n    materialInput.str = pointMC / sensorRadius;   \n\n    vec3 positionToEyeEC = -v_positionEC;   \n    materialInput.positionToEyeEC = positionToEyeEC;    \n\n    vec3 normalEC = normalize(v_normalEC);  \n    materialInput.normalEC = u_normalDirection * normalEC; \n\n    czm_material material = czm_getMaterial(materialInput);   \n\n     return mix(czm_phong(normalize(positionToEyeEC), material,czm_lightDirectionEC),vec4(material.diffuse, material.alpha), 0.4);\n}\n\nbool isOnBoundary(float value, float epsilon)\n{    \n    float width = getIntersectionWidth();  \n    float tolerance = width * epsilon;\n    #ifdef GL_OES_standard_derivatives    \n    float delta = max(abs(dFdx(value)), abs(dFdy(value))); \n    float pixels = width * delta;   \n    float temp = abs(value);    \n    // There are a couple things going on here.   \n    // First we test the value at the current fragment to see if it is within the tolerance.   \n    // We also want to check if the value of an adjacent pixel is within the tolerance,    \n    // but we don't want to admit points that are obviously not on the surface.    \n    // For example, if we are looking for \"value\" to be close to 0, but value is 1 and the adjacent value is 2,  \n    // then the delta would be 1 and \"temp - delta\" would be \"1 - 1\" which is zero even though neither of   \n    // the points is close to zero.   \n    return temp < tolerance && temp < pixels || (delta < 10.0 * tolerance && temp - delta < tolerance && temp < pixels);\n    #else    \n    return abs(value) < tolerance;\n    #endif\n}\n\n\nvec4 shade(bool isOnBoundary)\n{   \n\n    if (u_showIntersection && isOnBoundary) \n        {      \n            return getIntersectionColor();  \n        }   \n    if(u_type == 1.0)\n        {       \n            return getLineColor(); \n        }   \n\n    return getColor(u_radius, v_positionEC);\n\n}\n                   \nfloat ellipsoidSurfaceFunction(vec3 point)\n{   \n    vec3 scaled = czm_ellipsoidInverseRadii * point;  \n    return dot(scaled, scaled) - 1.0;\n}\n\nvoid main()\n{    \n    vec3 sensorVertexWC = czm_model[3].xyz;      // (0.0, 0.0, 0.0) in model coordinates   \n    vec3 sensorVertexEC = czm_modelView[3].xyz;  // (0.0, 0.0, 0.0) in model coordinates   \n    //vec3 pixDir = normalize(v_position);    \n    float positionX = v_position.x;  \n    float positionY = v_position.y;  \n    float positionZ = v_position.z;  \n    vec3 zDir = vec3(0.0, 0.0, 1.0); \n    vec3 lineX = vec3(positionX, 0 ,positionZ);  \n    vec3 lineY = vec3(0, positionY, positionZ);  \n    float resX = dot(normalize(lineX), zDir);  \n    if(resX < cos(u_xHalfAngle)-0.00001)\n    {       \n        discard; \n    }    \n    float resY = dot(normalize(lineY), zDir);  \n    if(resY < cos(u_yHalfAngle)-0.00001)\n    {    \n        discard;  \n    } \n    float ellipsoidValue = ellipsoidSurfaceFunction(v_positionWC);\n    // Occluded by the ellipsoid?\n    if (!u_showThroughEllipsoid)\n  {   \n    // Discard if in the ellipsoid  \n    // PERFORMANCE_IDEA: A coarse check for ellipsoid intersection could be done on the CPU first. \n    if (ellipsoidValue < 0.0) \n    {       \n        discard;  \n    }   \n    // Discard if in the sensor's shadow   \n    if (inSensorShadow(sensorVertexWC, v_positionWC)) \n    {      \n        discard;  \n    }   \n  }   \n    // Notes: Each surface functions should have an associated tolerance based on the floating point error. \n    bool isOnEllipsoid = isOnBoundary(ellipsoidValue, czm_epsilon3);   \n    //isOnEllipsoid = false;   \n    // if((resX >= 0.8 && resX <= 0.81)||(resY >= 0.8 && resY <= 0.81))\n    // {     \n    //     out_FragColor = shade(false);\n    // }\n    out_FragColor = shade(isOnEllipsoid);\n}\n";
+/**
+ * 扫描面FS着色器
+ * @returns {ShaderSource}
+ */
+export var scanPlaneFS = "\nuniform bool u_showIntersection;\nuniform bool u_showThroughEllipsoid;\nuniform float u_radius;\nuniform float u_xHalfAngle;\nuniform float u_yHalfAngle;\nuniform float u_normalDirection;\nuniform vec4 u_color;\nin vec3 v_position;\nin vec3 v_positionWC;\nin vec3 v_positionEC;\nin vec3 v_normalEC;\nvec4 getColor(float sensorRadius, vec3 pointEC)\n{    \n    czm_materialInput materialInput;    \n    vec3 pointMC = (czm_inverseModelView * vec4(pointEC, 1.0)).xyz;  \n    materialInput.st = sensor2dTextureCoordinates(sensorRadius, pointMC);\n    materialInput.str = pointMC / sensorRadius;    \n    vec3 positionToEyeEC = -v_positionEC;    \n    materialInput.positionToEyeEC = positionToEyeEC;    \n    vec3 normalEC = normalize(v_normalEC);   \n    materialInput.normalEC = u_normalDirection * normalEC;    \n    czm_material material = czm_getMaterial(materialInput);   \n    material.diffuse = u_color.rgb;    \n    material.alpha = u_color.a;    \n    return mix(czm_phong(normalize(positionToEyeEC), material,czm_lightDirectionEC), vec4(material.diffuse, material.alpha), 0.4);\n}\n\nbool isOnBoundary(float value, float epsilon)\n{    \n    float width = getIntersectionWidth();   \n    float tolerance = width * epsilon;\n    #ifdef GL_OES_standard_derivatives    \n    float delta = max(abs(dFdx(value)), abs(dFdy(value)));    \n    float pixels = width * delta;    \n    float temp = abs(value);    \n    // There are a couple things going on here.    \n    //  First we test the value at the current fragment to see if it is within the tolerance.    \n    //  We also want to check if the value of an adjacent pixel is within the tolerance,   \n    //  but we don't want to admit points that are obviously not on the surface.    \n    //  For example, if we are looking for value to be close to 0, but value is 1 and the adjacent value is 2,    \n    //  then the delta would be 1 and temp - delta would be 1 - 1 which is zero even though neither of    \n    //  the points is close to zero.    \n    return temp < tolerance && temp < pixels || (delta < 10.0 * tolerance && temp - delta < tolerance && temp < pixels);\n    #else    \n    return abs(value) < tolerance;\n    #endif\n}\n\n\nvec4 shade(bool isOnBoundary)\n    {   \n    if (u_showIntersection && isOnBoundary)  \n    {    \n        return getIntersectionColor(); \n    }   \n    return getColor(u_radius, v_positionEC);\n}\n\n\nfloat ellipsoidSurfaceFunction(vec3 point)\n{    \n    vec3 scaled = czm_ellipsoidInverseRadii * point;   \n    return dot(scaled, scaled) - 1.0;\n}\n\n\nvoid main()\n{  \n    vec3 sensorVertexWC = czm_model[3].xyz;      // (0.0, 0.0, 0.0) in model coordinates   \n    vec3 sensorVertexEC = czm_modelView[3].xyz;  // (0.0, 0.0, 0.0) in model coordinates   \n    //vec3 pixDir = normalize(v_position);    \n    float positionX = v_position.x;    \n    float positionY = v_position.y;    \n    float positionZ = v_position.z;    \n    vec3 zDir = vec3(0.0, 0.0, 1.0);    \n    vec3 lineX = vec3(positionX, 0 ,positionZ);    \n    vec3 lineY = vec3(0, positionY, positionZ);    \n    float resX = dot(normalize(lineX), zDir);    \n    if(resX < cos(u_xHalfAngle) - 0.0001)\n    {       \n        discard;  \n    }    \n    float resY = dot(normalize(lineY), zDir);  \n    if(resY < cos(u_yHalfAngle)- 0.0001)\n    {   \n        discard;  \n    }   \n    float ellipsoidValue = ellipsoidSurfaceFunction(v_positionWC);  \n    // Occluded by the ellipsoid\n    if (!u_showThroughEllipsoid)\n  {   \n    // Discard if in the ellipsoid   \n    // PERFORMANCE_IDEA: A coarse check for ellipsoid intersection could be done on the CPU first.  \n    if (ellipsoidValue < 0.0)  \n    {      \n        discard;  \n    }   \n    // Discard if in the sensor's shadow \n    if (inSensorShadow(sensorVertexWC, v_positionWC)) \n    {    \n        discard;  \n    }   \n  }   \n    // Notes: Each surface functions should have an associated tolerance based on the floating point error.  \n    bool isOnEllipsoid = isOnBoundary(ellipsoidValue, czm_epsilon3);  \n    out_FragColor = shade(isOnEllipsoid);\n}\n    ";

package/dist/esm/graphicLayer/PolygonGroundPrimitiveImageMaterialInstance.js

@@ -46,7 +46,7 @@ var PolygonGroundPrimitiveImageMaterialInstance = /*#__PURE__*/function (_Polygo
           position: new Cesium.GeometryAttribute({
             componentDatatype: Cesium.ComponentDatatype.DOUBLE,
             componentsPerAttribute: 3,
-            values: Cesium.Cartesian3.packArray(positions)
+            values: Cesium.Cartesian3.packArray(utils.buildPolygonHierarchy(positions).positions)
           })
         },
         indices: new Uint32Array(indices),

package/dist/esm/index.d.ts

@@ -70,6 +70,7 @@ export declare class DeepTwinsEngine3D {
     static PoiLayer: any;
     static FlightPlanning: any;
     static FlightRiskEvaluation: any;
-    static BuildClampedPolygonGrid: any;
+    static BuildClampedPolygonGridUnified: any;
+    static ViewShedAnalyserLayer: any;
 }
 export default DeepTwinsEngine3D;

package/dist/esm/index.js

@@ -42,7 +42,7 @@ import ArcGisLayer from "./tileLayer/ArcGisLayer";
 import RasterLayer from "./tileLayer/RasterLayer";
 import WmsLayer from "./tileLayer/WmsLayer";
 import WmtsLayer from "./tileLayer/WmtsLayer";
-import BuildClampedPolygonGrid from "./tool/BuildClampedPolygonGrid";
+import BuildClampedPolygonGridUnified from "./tool/BuildClampedPolygonGrid";
 import Compass from "./tool/Compass";
 import TimerInterval from "./tool/TimerInterval";
 import { ClampToGround, ShapeSection } from "./tool/common";
@@ -53,17 +53,19 @@ import Frustum from "./visualization/Frustum";
 import Heatmap2d from "./visualization/Heatmap2d";
 import Heatmap3d from "./visualization/Heatmap3d";
 import PointCluster from "./visualization/PointCluster";
+// drawCommand
+import { ViewShedAnalyserLayer } from "./drawCommand";
 window.Cesium = Cesium;
 // 将DEEP_TWINS_BASE_URL赋值给CESIUM_BASE_URL
 DEEP_TWINS_BASE_URL && (window.CESIUM_BASE_URL = DEEP_TWINS_BASE_URL);
 // 全局加载css
 loadCss(Cesium.buildModuleUrl('Widgets/widgets.css'));
 // 打印版本信息
-console.log('DeepTwinsEngine3D Version:', "0.1.49");
+console.log('DeepTwinsEngine3D Version:', "0.1.50");
 export var DeepTwinsEngine3D = /*#__PURE__*/_createClass(function DeepTwinsEngine3D() {
   _classCallCheck(this, DeepTwinsEngine3D);
 });
-_defineProperty(DeepTwinsEngine3D, "Version", "0.1.49");
+_defineProperty(DeepTwinsEngine3D, "Version", "0.1.50");
 _defineProperty(DeepTwinsEngine3D, "Map", Map);
 _defineProperty(DeepTwinsEngine3D, "GroundSkyBox", GroundSkyBox);
 _defineProperty(DeepTwinsEngine3D, "RasterLayer", RasterLayer);
@@ -133,5 +135,6 @@ _defineProperty(DeepTwinsEngine3D, "VisitLayer", VisitLayer);
 _defineProperty(DeepTwinsEngine3D, "PoiLayer", PoiLayer);
 _defineProperty(DeepTwinsEngine3D, "FlightPlanning", FlightPlanning);
 _defineProperty(DeepTwinsEngine3D, "FlightRiskEvaluation", FlightRiskEvaluation);
-_defineProperty(DeepTwinsEngine3D, "BuildClampedPolygonGrid", BuildClampedPolygonGrid);
+_defineProperty(DeepTwinsEngine3D, "BuildClampedPolygonGridUnified", BuildClampedPolygonGridUnified);
+_defineProperty(DeepTwinsEngine3D, "ViewShedAnalyserLayer", ViewShedAnalyserLayer);
 export default DeepTwinsEngine3D;

package/dist/esm/map/Map.js

@@ -151,7 +151,7 @@ var Map = /*#__PURE__*/function (_Cesium$Viewer) {
     value: function _createVideoContainer() {
       var el = document.createElement('div');
       el.id = constant.VIDEO_CONTAINER_ID;
-      el.style.cssText = 'position: absolute; left: 0; top: 0; right: 0; bottom: 0; pointer-events: none; overflow: hidden; display: none';
+      el.style.cssText = 'position: absolute; left: 0; top: 0; pointer-events: none; opacity: 0';
       this.container.appendChild(el);
     }
 

package/dist/esm/material/primitive/ImageMaterialAppearance.d.ts

@@ -2,7 +2,7 @@ import BaseMaterialAppearance from './BaseMaterialAppearance';
 declare class ImageMaterialAppearance extends BaseMaterialAppearance {
     private readonly currMaterial;
     private readonly currAppearance;
-    constructor(material: any, appearance: any);
+    constructor(material: any, appearance?: any);
     static handleMaterial(style: any): any;
     private _createMaterial;
     private _createAppearance;

package/dist/esm/tool/BuildClampedPolygonGrid.d.ts

@@ -1,11 +1,4 @@
-import * as Cesium from 'deeptwins-cesium';
-/**
- * 构建贴地的多边形网格
- */
-export default function BuildClampedPolygonGrid(map: any, geoJson: any, gridSize?: number, heightBuffer?: number): Promise<{
+export default function BuildClampedPolygonGridUnified(map: any, geoJson: any, gridSize?: number, heightBuffer?: number): Promise<{
     type: string;
-    source: {
-        positions: Cesium.Cartesian3[];
-        indices: number[];
-    }[];
+    source: any[];
 }>;

package/dist/esm/tool/BuildClampedPolygonGrid.js

@@ -1,11 +1,40 @@
 function _typeof(o) { "@babel/helpers - typeof"; return _typeof = "function" == typeof Symbol && "symbol" == typeof Symbol.iterator ? function (o) { return typeof o; } : function (o) { return o && "function" == typeof Symbol && o.constructor === Symbol && o !== Symbol.prototype ? "symbol" : typeof o; }, _typeof(o); }
 function _regeneratorRuntime() { "use strict"; /*! regenerator-runtime -- Copyright (c) 2014-present, Facebook, Inc. -- license (MIT): https://github.com/facebook/regenerator/blob/main/LICENSE */ _regeneratorRuntime = function _regeneratorRuntime() { return e; }; var t, e = {}, r = Object.prototype, n = r.hasOwnProperty, o = Object.defineProperty || function (t, e, r) { t[e] = r.value; }, i = "function" == typeof Symbol ? Symbol : {}, a = i.iterator || "@@iterator", c = i.asyncIterator || "@@asyncIterator", u = i.toStringTag || "@@toStringTag"; function define(t, e, r) { return Object.defineProperty(t, e, { value: r, enumerable: !0, configurable: !0, writable: !0 }), t[e]; } try { define({}, ""); } catch (t) { define = function define(t, e, r) { return t[e] = r; }; } function wrap(t, e, r, n) { var i = e && e.prototype instanceof Generator ? e : Generator, a = Object.create(i.prototype), c = new Context(n || []); return o(a, "_invoke", { value: makeInvokeMethod(t, r, c) }), a; } function tryCatch(t, e, r) { try { return { type: "normal", arg: t.call(e, r) }; } catch (t) { return { type: "throw", arg: t }; } } e.wrap = wrap; var h = "suspendedStart", l = "suspendedYield", f = "executing", s = "completed", y = {}; function Generator() {} function GeneratorFunction() {} function GeneratorFunctionPrototype() {} var p = {}; define(p, a, function () { return this; }); var d = Object.getPrototypeOf, v = d && d(d(values([]))); v && v !== r && n.call(v, a) && (p = v); var g = GeneratorFunctionPrototype.prototype = Generator.prototype = Object.create(p); function defineIteratorMethods(t) { ["next", "throw", "return"].forEach(function (e) { define(t, e, function (t) { return this._invoke(e, t); }); }); } function AsyncIterator(t, e) { function invoke(r, o, i, a) { var c = tryCatch(t[r], t, o); if ("throw" !== c.type) { var u = c.arg, h = u.value; return h && "object" == _typeof(h) && n.call(h, "__await") ? e.resolve(h.__await).then(function (t) { invoke("next", t, i, a); }, function (t) { invoke("throw", t, i, a); }) : e.resolve(h).then(function (t) { u.value = t, i(u); }, function (t) { return invoke("throw", t, i, a); }); } a(c.arg); } var r; o(this, "_invoke", { value: function value(t, n) { function callInvokeWithMethodAndArg() { return new e(function (e, r) { invoke(t, n, e, r); }); } return r = r ? r.then(callInvokeWithMethodAndArg, callInvokeWithMethodAndArg) : callInvokeWithMethodAndArg(); } }); } function makeInvokeMethod(e, r, n) { var o = h; return function (i, a) { if (o === f) throw new Error("Generator is already running"); if (o === s) { if ("throw" === i) throw a; return { value: t, done: !0 }; } for (n.method = i, n.arg = a;;) { var c = n.delegate; if (c) { var u = maybeInvokeDelegate(c, n); if (u) { if (u === y) continue; return u; } } if ("next" === n.method) n.sent = n._sent = n.arg;else if ("throw" === n.method) { if (o === h) throw o = s, n.arg; n.dispatchException(n.arg); } else "return" === n.method && n.abrupt("return", n.arg); o = f; var p = tryCatch(e, r, n); if ("normal" === p.type) { if (o = n.done ? s : l, p.arg === y) continue; return { value: p.arg, done: n.done }; } "throw" === p.type && (o = s, n.method = "throw", n.arg = p.arg); } }; } function maybeInvokeDelegate(e, r) { var n = r.method, o = e.iterator[n]; if (o === t) return r.delegate = null, "throw" === n && e.iterator.return && (r.method = "return", r.arg = t, maybeInvokeDelegate(e, r), "throw" === r.method) || "return" !== n && (r.method = "throw", r.arg = new TypeError("The iterator does not provide a '" + n + "' method")), y; var i = tryCatch(o, e.iterator, r.arg); if ("throw" === i.type) return r.method = "throw", r.arg = i.arg, r.delegate = null, y; var a = i.arg; return a ? a.done ? (r[e.resultName] = a.value, r.next = e.nextLoc, "return" !== r.method && (r.method = "next", r.arg = t), r.delegate = null, y) : a : (r.method = "throw", r.arg = new TypeError("iterator result is not an object"), r.delegate = null, y); } function pushTryEntry(t) { var e = { tryLoc: t[0] }; 1 in t && (e.catchLoc = t[1]), 2 in t && (e.finallyLoc = t[2], e.afterLoc = t[3]), this.tryEntries.push(e); } function resetTryEntry(t) { var e = t.completion || {}; e.type = "normal", delete e.arg, t.completion = e; } function Context(t) { this.tryEntries = [{ tryLoc: "root" }], t.forEach(pushTryEntry, this), this.reset(!0); } function values(e) { if (e || "" === e) { var r = e[a]; if (r) return r.call(e); if ("function" == typeof e.next) return e; if (!isNaN(e.length)) { var o = -1, i = function next() { for (; ++o < e.length;) if (n.call(e, o)) return next.value = e[o], next.done = !1, next; return next.value = t, next.done = !0, next; }; return i.next = i; } } throw new TypeError(_typeof(e) + " is not iterable"); } return GeneratorFunction.prototype = GeneratorFunctionPrototype, o(g, "constructor", { value: GeneratorFunctionPrototype, configurable: !0 }), o(GeneratorFunctionPrototype, "constructor", { value: GeneratorFunction, configurable: !0 }), GeneratorFunction.displayName = define(GeneratorFunctionPrototype, u, "GeneratorFunction"), e.isGeneratorFunction = function (t) { var e = "function" == typeof t && t.constructor; return !!e && (e === GeneratorFunction || "GeneratorFunction" === (e.displayName || e.name)); }, e.mark = function (t) { return Object.setPrototypeOf ? Object.setPrototypeOf(t, GeneratorFunctionPrototype) : (t.__proto__ = GeneratorFunctionPrototype, define(t, u, "GeneratorFunction")), t.prototype = Object.create(g), t; }, e.awrap = function (t) { return { __await: t }; }, defineIteratorMethods(AsyncIterator.prototype), define(AsyncIterator.prototype, c, function () { return this; }), e.AsyncIterator = AsyncIterator, e.async = function (t, r, n, o, i) { void 0 === i && (i = Promise); var a = new AsyncIterator(wrap(t, r, n, o), i); return e.isGeneratorFunction(r) ? a : a.next().then(function (t) { return t.done ? t.value : a.next(); }); }, defineIteratorMethods(g), define(g, u, "Generator"), define(g, a, function () { return this; }), define(g, "toString", function () { return "[object Generator]"; }), e.keys = function (t) { var e = Object(t), r = []; for (var n in e) r.push(n); return r.reverse(), function next() { for (; r.length;) { var t = r.pop(); if (t in e) return next.value = t, next.done = !1, next; } return next.done = !0, next; }; }, e.values = values, Context.prototype = { constructor: Context, reset: function reset(e) { if (this.prev = 0, this.next = 0, this.sent = this._sent = t, this.done = !1, this.delegate = null, this.method = "next", this.arg = t, this.tryEntries.forEach(resetTryEntry), !e) for (var r in this) "t" === r.charAt(0) && n.call(this, r) && !isNaN(+r.slice(1)) && (this[r] = t); }, stop: function stop() { this.done = !0; var t = this.tryEntries[0].completion; if ("throw" === t.type) throw t.arg; return this.rval; }, dispatchException: function dispatchException(e) { if (this.done) throw e; var r = this; function handle(n, o) { return a.type = "throw", a.arg = e, r.next = n, o && (r.method = "next", r.arg = t), !!o; } for (var o = this.tryEntries.length - 1; o >= 0; --o) { var i = this.tryEntries[o], a = i.completion; if ("root" === i.tryLoc) return handle("end"); if (i.tryLoc <= this.prev) { var c = n.call(i, "catchLoc"), u = n.call(i, "finallyLoc"); if (c && u) { if (this.prev < i.catchLoc) return handle(i.catchLoc, !0); if (this.prev < i.finallyLoc) return handle(i.finallyLoc); } else if (c) { if (this.prev < i.catchLoc) return handle(i.catchLoc, !0); } else { if (!u) throw new Error("try statement without catch or finally"); if (this.prev < i.finallyLoc) return handle(i.finallyLoc); } } } }, abrupt: function abrupt(t, e) { for (var r = this.tryEntries.length - 1; r >= 0; --r) { var o = this.tryEntries[r]; if (o.tryLoc <= this.prev && n.call(o, "finallyLoc") && this.prev < o.finallyLoc) { var i = o; break; } } i && ("break" === t || "continue" === t) && i.tryLoc <= e && e <= i.finallyLoc && (i = null); var a = i ? i.completion : {}; return a.type = t, a.arg = e, i ? (this.method = "next", this.next = i.finallyLoc, y) : this.complete(a); }, complete: function complete(t, e) { if ("throw" === t.type) throw t.arg; return "break" === t.type || "continue" === t.type ? this.next = t.arg : "return" === t.type ? (this.rval = this.arg = t.arg, this.method = "return", this.next = "end") : "normal" === t.type && e && (this.next = e), y; }, finish: function finish(t) { for (var e = this.tryEntries.length - 1; e >= 0; --e) { var r = this.tryEntries[e]; if (r.finallyLoc === t) return this.complete(r.completion, r.afterLoc), resetTryEntry(r), y; } }, catch: function _catch(t) { for (var e = this.tryEntries.length - 1; e >= 0; --e) { var r = this.tryEntries[e]; if (r.tryLoc === t) { var n = r.completion; if ("throw" === n.type) { var o = n.arg; resetTryEntry(r); } return o; } } throw new Error("illegal catch attempt"); }, delegateYield: function delegateYield(e, r, n) { return this.delegate = { iterator: values(e), resultName: r, nextLoc: n }, "next" === this.method && (this.arg = t), y; } }, e; }
+function ownKeys(e, r) { var t = Object.keys(e); if (Object.getOwnPropertySymbols) { var o = Object.getOwnPropertySymbols(e); r && (o = o.filter(function (r) { return Object.getOwnPropertyDescriptor(e, r).enumerable; })), t.push.apply(t, o); } return t; }
+function _objectSpread(e) { for (var r = 1; r < arguments.length; r++) { var t = null != arguments[r] ? arguments[r] : {}; r % 2 ? ownKeys(Object(t), !0).forEach(function (r) { _defineProperty(e, r, t[r]); }) : Object.getOwnPropertyDescriptors ? Object.defineProperties(e, Object.getOwnPropertyDescriptors(t)) : ownKeys(Object(t)).forEach(function (r) { Object.defineProperty(e, r, Object.getOwnPropertyDescriptor(t, r)); }); } return e; }
+function _defineProperty(obj, key, value) { key = _toPropertyKey(key); if (key in obj) { Object.defineProperty(obj, key, { value: value, enumerable: true, configurable: true, writable: true }); } else { obj[key] = value; } return obj; }
+function _toPropertyKey(t) { var i = _toPrimitive(t, "string"); return "symbol" == _typeof(i) ? i : String(i); }
+function _toPrimitive(t, r) { if ("object" != _typeof(t) || !t) return t; var e = t[Symbol.toPrimitive]; if (void 0 !== e) { var i = e.call(t, r || "default"); if ("object" != _typeof(i)) return i; throw new TypeError("@@toPrimitive must return a primitive value."); } return ("string" === r ? String : Number)(t); }
+function _createForOfIteratorHelper(o, allowArrayLike) { var it = typeof Symbol !== "undefined" && o[Symbol.iterator] || o["@@iterator"]; if (!it) { if (Array.isArray(o) || (it = _unsupportedIterableToArray(o)) || allowArrayLike && o && typeof o.length === "number") { if (it) o = it; var i = 0; var F = function F() {}; return { s: F, n: function n() { if (i >= o.length) return { done: true }; return { done: false, value: o[i++] }; }, e: function e(_e) { throw _e; }, f: F }; } throw new TypeError("Invalid attempt to iterate non-iterable instance.\nIn order to be iterable, non-array objects must have a [Symbol.iterator]() method."); } var normalCompletion = true, didErr = false, err; return { s: function s() { it = it.call(o); }, n: function n() { var step = it.next(); normalCompletion = step.done; return step; }, e: function e(_e2) { didErr = true; err = _e2; }, f: function f() { try { if (!normalCompletion && it.return != null) it.return(); } finally { if (didErr) throw err; } } }; }
+function _unsupportedIterableToArray(o, minLen) { if (!o) return; if (typeof o === "string") return _arrayLikeToArray(o, minLen); var n = Object.prototype.toString.call(o).slice(8, -1); if (n === "Object" && o.constructor) n = o.constructor.name; if (n === "Map" || n === "Set") return Array.from(o); if (n === "Arguments" || /^(?:Ui|I)nt(?:8|16|32)(?:Clamped)?Array$/.test(n)) return _arrayLikeToArray(o, minLen); }
+function _arrayLikeToArray(arr, len) { if (len == null || len > arr.length) len = arr.length; for (var i = 0, arr2 = new Array(len); i < len; i++) arr2[i] = arr[i]; return arr2; }
 function asyncGeneratorStep(gen, resolve, reject, _next, _throw, key, arg) { try { var info = gen[key](arg); var value = info.value; } catch (error) { reject(error); return; } if (info.done) { resolve(value); } else { Promise.resolve(value).then(_next, _throw); } }
 function _asyncToGenerator(fn) { return function () { var self = this, args = arguments; return new Promise(function (resolve, reject) { var gen = fn.apply(self, args); function _next(value) { asyncGeneratorStep(gen, resolve, reject, _next, _throw, "next", value); } function _throw(err) { asyncGeneratorStep(gen, resolve, reject, _next, _throw, "throw", err); } _next(undefined); }); }; }
+import { Delaunay } from 'd3-delaunay';
 import * as Cesium from 'deeptwins-cesium';
 import BaseSource from "../graphicLayer/BaseSource";
 import { translateSource } from "../tool/common";
 import * as utils from "./utils";
+function distanceToSegment(p, a, b) {
+  var abx = b.x - a.x;
+  var aby = b.y - a.y;
+  var apx = p.x - a.x;
+  var apy = p.y - a.y;
+  var abLen2 = abx * abx + aby * aby;
+  var t = abLen2 === 0 ? 0 : Math.max(0, Math.min(1, (apx * abx + apy * aby) / abLen2));
+  var dx = a.x + t * abx - p.x;
+  var dy = a.y + t * aby - p.y;
+  return Math.sqrt(dx * dx + dy * dy);
+}
+function distanceToPolygonEdge(p, polygon) {
+  var min = Infinity;
+  for (var i = 0; i < polygon.length; i++) {
+    var a = polygon[i];
+    var b = polygon[(i + 1) % polygon.length];
+    min = Math.min(min, distanceToSegment(p, a, b));
+  }
+  return min;
+}
 function pointInPolygon(pt, polygon) {
   var inside = false;
   for (var i = 0, j = polygon.length - 1; i < polygon.length; j = i++) {
@@ -19,38 +48,15 @@ function pointInPolygon(pt, polygon) {
   return inside;
 }
 
-// 修改 buildGridIndices：不删顶点，只生成 polygon 内部三角形索引
-function buildGridIndices(cols, rows, polygon, gridPoints2D, step) {
-  var indices = [];
-  for (var y = 0; y < rows - 1; y++) {
-    for (var x = 0; x < cols - 1; x++) {
-      var i0 = y * cols + x;
-      var i1 = i0 + 1;
-      var i2 = i0 + cols;
-      var i3 = i2 + 1;
-
-      // 三角形中心点，判断是否在 polygon 内
-      var cx = (x + 0.5) * step + gridPoints2D[0].x; // 这里偏移加上最小X
-      var cy = (y + 0.5) * step + gridPoints2D[0].y; // 偏移加上最小Y
-      if (pointInPolygon(new Cesium.Cartesian2(cx, cy), polygon)) {
-        indices.push(i0, i1, i2);
-        indices.push(i1, i3, i2);
-      }
-    }
-  }
-  return indices;
-}
-
-/**
- * 构建贴地的多边形网格
- */
-export default function BuildClampedPolygonGrid(_x, _x2) {
-  return _BuildClampedPolygonGrid.apply(this, arguments);
+// 构建贴模型/贴地网格
+export default function BuildClampedPolygonGridUnified(_x, _x2) {
+  return _BuildClampedPolygonGridUnified.apply(this, arguments);
 }
-function _BuildClampedPolygonGrid() {
-  _BuildClampedPolygonGrid = _asyncToGenerator( /*#__PURE__*/_regeneratorRuntime().mark(function _callee(map, geoJson) {
+function _BuildClampedPolygonGridUnified() {
+  _BuildClampedPolygonGridUnified = _asyncToGenerator( /*#__PURE__*/_regeneratorRuntime().mark(function _callee(map, geoJson) {
     var gridSize,
       heightBuffer,
+      scene,
       translatedSource,
       source,
       finalSource,
@@ -62,90 +68,206 @@ function _BuildClampedPolygonGrid() {
         case 0:
           gridSize = _args2.length > 2 && _args2[2] !== undefined ? _args2[2] : 10;
           heightBuffer = _args2.length > 3 && _args2[3] !== undefined ? _args2[3] : 0.5;
-          if (!(!map.terrainProvider || map.terrainProvider instanceof Cesium.EllipsoidTerrainProvider)) {
-            _context2.next = 4;
-            break;
-          }
-          throw new Error('请设置地形');
-        case 4:
+          scene = map.scene; // 地形至少加载一次（给 fallback 用）
+          _context2.next = 5;
+          return utils.isTerrainReady(map);
+        case 5:
           translatedSource = translateSource(geoJson);
           source = BaseSource.geoJsonToGeoCartesian3Array(translatedSource);
           finalSource = [];
           _loop = /*#__PURE__*/_regeneratorRuntime().mark(function _loop() {
-            var positions, ring, center, plane, poly2D, minX, minY, maxX, maxY, cols, rows, gridPoints2D, y, x, indices, positions3D, cartographics, clamped, finalPositions;
+            var positions, ring, center, plane, poly2D, minX, minY, maxX, maxY, _iterator, _step, p, gridPoints2D, edgeFactor, edgeThreshold, y, x, center2D, dist, step, j, a, b, edgeLen, n, k, t, delaunay, indices, positions3D, modelClamped, needTerrain, terrainClamped, cartos, sampled, finalPositions;
             return _regeneratorRuntime().wrap(function _loop$(_context) {
               while (1) switch (_context.prev = _context.next) {
                 case 0:
-                  positions = source[i].positions; // 1️⃣ 去闭合
-                  ring = utils.buildPolygonHierarchy(positions).positions; // 2️⃣ 局部切平面
+                  positions = source[i].positions;
+                  ring = utils.buildPolygonHierarchy(positions).positions;
+                  /* ---------------------------------- */
+                  /* 1. 局部切平面 */
+                  /* ---------------------------------- */
                   center = Cesium.BoundingSphere.fromPoints(ring).center;
-                  plane = new Cesium.EllipsoidTangentPlane(center, Cesium.Ellipsoid.WGS84); // 3️⃣ 投影到 2D
+                  plane = new Cesium.EllipsoidTangentPlane(center, Cesium.Ellipsoid.WGS84);
                   poly2D = ring.map(function (p) {
                     return plane.projectPointOntoPlane(p);
-                  }); // 4️⃣ 包围盒
-                  minX = Infinity, minY = Infinity;
-                  maxX = -Infinity, maxY = -Infinity;
-                  poly2D.forEach(function (p) {
-                    minX = Math.min(minX, p.x);
-                    minY = Math.min(minY, p.y);
-                    maxX = Math.max(maxX, p.x);
-                    maxY = Math.max(maxY, p.y);
                   });
-                  cols = Math.floor((maxX - minX) / gridSize) + 1;
-                  rows = Math.floor((maxY - minY) / gridSize) + 1; // 5️⃣ 生成完整规则网格点（不删点）
+                  /* ---------------------------------- */
+                  /* 2. 包围盒 */
+                  /* ---------------------------------- */
+                  minX = Infinity, minY = Infinity, maxX = -Infinity, maxY = -Infinity;
+                  _iterator = _createForOfIteratorHelper(poly2D);
+                  try {
+                    for (_iterator.s(); !(_step = _iterator.n()).done;) {
+                      p = _step.value;
+                      minX = Math.min(minX, p.x);
+                      minY = Math.min(minY, p.y);
+                      maxX = Math.max(maxX, p.x);
+                      maxY = Math.max(maxY, p.y);
+                    }
+
+                    /* ---------------------------------- */
+                    /* 3. 自适应采样 */
+                    /* ---------------------------------- */
+                  } catch (err) {
+                    _iterator.e(err);
+                  } finally {
+                    _iterator.f();
+                  }
                   gridPoints2D = [];
-                  for (y = 0; y < rows; y++) {
-                    for (x = 0; x < cols; x++) {
-                      gridPoints2D.push(new Cesium.Cartesian2(minX + x * gridSize, minY + y * gridSize));
+                  edgeFactor = 0.3;
+                  edgeThreshold = gridSize;
+                  y = minY;
+                case 12:
+                  if (!(y < maxY)) {
+                    _context.next = 28;
+                    break;
+                  }
+                  x = minX;
+                case 14:
+                  if (!(x < maxX)) {
+                    _context.next = 25;
+                    break;
+                  }
+                  center2D = new Cesium.Cartesian2(x + gridSize * 0.5, y + gridSize * 0.5);
+                  if (pointInPolygon(center2D, poly2D)) {
+                    _context.next = 19;
+                    break;
+                  }
+                  x += gridSize;
+                  return _context.abrupt("continue", 14);
+                case 19:
+                  dist = distanceToPolygonEdge(center2D, poly2D);
+                  step = dist < edgeThreshold ? gridSize * edgeFactor : gridSize;
+                  gridPoints2D.push(center2D);
+                  x += step;
+                  _context.next = 14;
+                  break;
+                case 25:
+                  y += gridSize;
+                  _context.next = 12;
+                  break;
+                case 28:
+                  /* ---------------------------------- */
+                  /* 4. 边界插值 */
+                  /* ---------------------------------- */
+
+                  for (j = 0; j < poly2D.length; j++) {
+                    a = poly2D[j];
+                    b = poly2D[(j + 1) % poly2D.length];
+                    edgeLen = Cesium.Cartesian2.distance(a, b);
+                    n = Math.ceil(edgeLen / (gridSize * edgeFactor));
+                    for (k = 0; k <= n; k++) {
+                      t = k / n;
+                      gridPoints2D.push(new Cesium.Cartesian2(Cesium.Math.lerp(a.x, b.x, t), Cesium.Math.lerp(a.y, b.y, t)));
                     }
                   }
 
-                  // 6️⃣ 生成三角形索引（只保留 polygon 内部三角形）
-                  indices = buildGridIndices(cols, rows, poly2D, gridPoints2D, gridSize); // 7️⃣ 回投 3D
+                  /* ---------------------------------- */
+                  /* 5. Delaunay */
+                  /* ---------------------------------- */
+                  delaunay = Delaunay.from(gridPoints2D.map(function (p) {
+                    return [p.x, p.y];
+                  }));
+                  indices = Array.from(delaunay.triangles);
+                  /* ---------------------------------- */
+                  /* 6. 回投 3D（无高度） */
+                  /* ---------------------------------- */
                   positions3D = gridPoints2D.map(function (p) {
                     return plane.projectPointOntoEllipsoid(p);
-                  }); // 8️⃣ 贴地
-                  cartographics = positions3D.map(function (p) {
+                  });
+                  /* ---------------------------------- */
+                  /* 7. 贴模型（优先） */
+                  /* ---------------------------------- */
+                  // await scene.render();
+                  modelClamped = null;
+                  _context.prev = 33;
+                  _context.next = 36;
+                  return scene.clampToHeightMostDetailed(JSON.parse(JSON.stringify(positions3D)));
+                case 36:
+                  modelClamped = _context.sent;
+                  _context.next = 42;
+                  break;
+                case 39:
+                  _context.prev = 39;
+                  _context.t0 = _context["catch"](33);
+                  modelClamped = null;
+                case 42:
+                  /* ---------------------------------- */
+                  /* 8. terrain fallback（只给失败点） */
+                  /* ---------------------------------- */
+                  needTerrain = !modelClamped || modelClamped.some(function (p) {
+                    return !Cesium.defined(p);
+                  });
+                  terrainClamped = null;
+                  if (!needTerrain) {
+                    _context.next = 50;
+                    break;
+                  }
+                  cartos = JSON.parse(JSON.stringify(positions3D)).map(function (p) {
                     return Cesium.Cartographic.fromCartesian(p);
                   });
-                  _context.next = 17;
-                  return Cesium.sampleTerrainMostDetailed(map.terrainProvider, cartographics);
-                case 17:
-                  clamped = _context.sent;
-                  finalPositions = clamped.map(function (c) {
-                    return Cesium.Cartesian3.fromRadians(c.longitude, c.latitude, c.height + heightBuffer);
+                  _context.next = 48;
+                  return Cesium.sampleTerrainMostDetailed(map.terrainProvider, cartos);
+                case 48:
+                  sampled = _context.sent;
+                  terrainClamped = sampled.map(function (c) {
+                    return Cesium.Cartesian3.fromRadians(c.longitude, c.latitude, c.height);
+                  });
+                case 50:
+                  /* ---------------------------------- */
+                  /* 9. 合并结果（模型优先，terrain 兜底） */
+                  /* ---------------------------------- */
+                  finalPositions = positions3D.map(function (_, idx) {
+                    var p = modelClamped && Cesium.defined(modelClamped[idx]) ? modelClamped[idx] : terrainClamped ? terrainClamped[idx] : null;
+                    if (!p) return null;
+                    var carto = Cesium.Cartographic.fromCartesian(p);
+                    return Cesium.Cartesian3.fromRadians(carto.longitude, carto.latitude, carto.height + heightBuffer);
                   });
-                  finalSource.push({
+                  if (!finalPositions.some(function (p) {
+                    return !p;
+                  })) {
+                    _context.next = 54;
+                    break;
+                  }
+                  console.warn('存在无法贴地/贴模型的点，已跳过该 polygon');
+                  return _context.abrupt("return", 1);
+                case 54:
+                  finalSource.push(_objectSpread(_objectSpread({}, source[i]), {}, {
                     positions: finalPositions,
                     indices: indices
-                  });
-                case 20:
+                  }));
+                case 55:
                 case "end":
                   return _context.stop();
               }
-            }, _loop);
+            }, _loop, null, [[33, 39]]);
           });
           i = 0;
-        case 9:
+        case 10:
           if (!(i < source.length)) {
+            _context2.next = 17;
+            break;
+          }
+          return _context2.delegateYield(_loop(), "t0", 12);
+        case 12:
+          if (!_context2.t0) {
             _context2.next = 14;
             break;
           }
-          return _context2.delegateYield(_loop(), "t0", 11);
-        case 11:
+          return _context2.abrupt("continue", 14);
+        case 14:
           i++;
-          _context2.next = 9;
+          _context2.next = 10;
           break;
-        case 14:
+        case 17:
           return _context2.abrupt("return", {
             type: 'clampedPolygonGrid',
             source: finalSource
           });
-        case 15:
+        case 18:
         case "end":
           return _context2.stop();
       }
     }, _callee);
   }));
-  return _BuildClampedPolygonGrid.apply(this, arguments);
+  return _BuildClampedPolygonGridUnified.apply(this, arguments);
 }

package/dist/esm/tool/utils.d.ts

@@ -71,3 +71,4 @@ export declare function getMidpoint(position1: Cartesian3, position2: Cartesian3
     alt: number;
 };
 export declare function densifyLine(p1: Cartesian3, p2: Cartesian3, count?: number): Cesium.Cartesian3[];
+export declare function isTerrainReady(map: any): Promise<void>;

package/dist/esm/tool/utils.js

@@ -477,4 +477,25 @@ export function densifyLine(p1, p2) {
     result.push(Cesium.Cartesian3.lerp(p1, p2, i / count, new Cesium.Cartesian3()));
   }
   return result;
+}
+
+// 判断地形是否加载完成
+export function isTerrainReady(map) {
+  return new Promise(function (resolve) {
+    if (!map.scene.globe || !map.scene.terrainProvider) {
+      resolve();
+      return;
+    }
+    if (map.scene.globe.tilesLoaded) {
+      resolve();
+      return;
+    }
+    var remove = map.scene.postRender.addEventListener(function () {
+      if (map.scene.globe.tilesLoaded) {
+        remove();
+        resolve();
+      }
+    });
+    map.scene.requestRender();
+  });
 }

package/dist/esm/videoFusion/BaseVideo.js

@@ -157,8 +157,10 @@ var BaseVideo = /*#__PURE__*/function () {
   }, {
     key: "_drawCanvas",
     value: function _drawCanvas(options) {
-      var width = options.width,
-        height = options.height;
+      var _options$videoOptions = options.videoOptions,
+        width = _options$videoOptions.width,
+        height = _options$videoOptions.height,
+        edgeFeather = options.style.edgeFeather;
       var canvas = document.createElement('canvas');
       canvas.width = width;
       canvas.height = height;
@@ -169,18 +171,20 @@ var BaseVideo = /*#__PURE__*/function () {
       ctx.drawImage(this.videoElement, 0, 0, width, height);
 
       // 2. 应用羽化遮罩
-      ctx.globalCompositeOperation = 'destination-in'; // 使用遮罩影响透明度
+      if (edgeFeather) {
+        ctx.globalCompositeOperation = 'destination-in'; // 使用遮罩影响透明度
 
-      var maskCanvas = document.createElement('canvas');
-      maskCanvas.width = width;
-      maskCanvas.height = height;
-      var maskCtx = maskCanvas.getContext('2d');
-      if (maskCtx && this._maskImage) {
-        // 绘制并拉伸遮罩图片到整个画布
-        maskCtx.drawImage(this._maskImage, 0, 0, width, height);
-        ctx.drawImage(maskCanvas, 0, 0);
+        var maskCanvas = document.createElement('canvas');
+        maskCanvas.width = width;
+        maskCanvas.height = height;
+        var maskCtx = maskCanvas.getContext('2d');
+        if (maskCtx && this._maskImage) {
+          // 绘制并拉伸遮罩图片到整个画布
+          maskCtx.drawImage(this._maskImage, 0, 0, width, height);
+          ctx.drawImage(maskCanvas, 0, 0);
+        }
+        ctx.globalCompositeOperation = 'source-over'; // 还原默认混合模式
       }
-      ctx.globalCompositeOperation = 'source-over'; // 还原默认混合模式
       return canvas;
     }
 
@@ -213,11 +217,6 @@ var BaseVideo = /*#__PURE__*/function () {
             _this2.onErrorFun && _this2.onErrorFun(e);
           });
         }
-        if (this.videoPolygonLayer) {
-          this.videoPolygonLayer.setStyle({
-            material: this.videoElement
-          });
-        }
       }
       // 暂停后继续播放
       if (this.synchronizer && this.videoType === 'hls') {
@@ -250,6 +249,7 @@ var BaseVideo = /*#__PURE__*/function () {
         });
         this.getMap().clock.shouldAnimate = true;
       }
+      this.videoElement.style.display = 'none';
     }
 
     // 监听错误