dimine-python-sdk 0.2.1__py3-none-any.whl → 0.2.2__py3-none-any.whl

dimine_python_sdk/__init__.py

@@ -3,7 +3,7 @@ from pathlib import Path
 from .conn import ConnClient
 
 DOCS_PATH = Path(__file__).parent / "docs"
-__VERSION__ = "0.2.1"  # 由 build.bat 中的 generate_version.py 自动更新
+__VERSION__ = "0.2.2"  # 由 build.bat 中的 generate_version.py 自动更新
 
 
 __all__ = [

dimine_python_sdk/cli/commands/file.py

@@ -195,13 +195,13 @@ def register_commands(subparsers, format_parent) -> None:
     p.add_argument("--file-id", default=None, help="文件ID，多个文件打开时必填")
 
     p = subparsers.add_parser("file-open", parents=[format_parent], help="打开文件到 Dimine 场景")
-    p.add_argument("--path", required=True, help="要打开的 .dmf/.dwg/.dxf 等文件路径")
+    p.add_argument("--path", "--file", required=True, help="要打开的 .dmf/.dwg/.dxf 等文件路径")
 
     p = subparsers.add_parser("file-close", parents=[format_parent], help="关闭 Dimine 中的文件")
     p.add_argument("--file-id", required=True, help="要关闭的文件ID")
 
     p = subparsers.add_parser("file-create", parents=[format_parent], help="在 Dimine 中创建新的空文件")
-    p.add_argument("--path", required=True, help="新文件的路径")
+    p.add_argument("--path", "--file", required=True, help="新文件的路径")
 
     p = subparsers.add_parser("file-visibility", parents=[format_parent], help="设置文件显示/隐藏状态")
     p.add_argument("--file-id", required=True, help="文件ID")

dimine_python_sdk/docs/guides/lib/exploitation/contour_modeling.md

@@ -331,7 +331,105 @@ for i, merged in enumerate(merged_results):
     print(f"  实体[{i}]: 顶点={merged.points.shape[0]}, 面片={merged.faces.shape[0]}")
 ```
 
-## 6. 算法类型
+## 6. 地表建模
+
+从多条等高线生成三角网地表模型（DTM）。输入为一组等高线，输出为三角网地表。
+
+### 基础用法
+
+```python
+import numpy as np
+from dimine_python_sdk.lib.exploitation import SurfaceModeling
+
+# 多条等高线（不同高程的封闭轮廓线）
+contour1 = [[0, 0, 100], [10, 0, 100], [10, 10, 100], [0, 10, 100], [0, 0, 100]]
+contour2 = [[2, 2, 200], [8, 2, 200], [8, 8, 200], [2, 8, 200], [2, 2, 200]]
+contour3 = [[4, 4, 300], [6, 4, 300], [6, 6, 300], [4, 6, 300], [4, 4, 300]]
+
+# 地表建模
+surfaces = SurfaceModeling.create([contour1, contour2, contour3], mode=0)
+
+print(f"生成地表数量: {len(surfaces)}")
+for i, tin in enumerate(surfaces):
+    print(f"  地表[{i}] 顶点数: {tin.points.shape[0]}, 三角面片数: {tin.faces.shape[0]}")
+```
+
+### 支持的输入格式
+
+`SurfaceModeling.create(contour_lines, mode=0)` 的每条等高线支持三种输入格式：
+
+```python
+import numpy as np
+from dimine_python_sdk.models.types import Line
+from dimine_python_sdk.lib.exploitation import SurfaceModeling
+
+# 方式一：List[List[List[float]]] 纯列表
+result = SurfaceModeling.create([
+    [[0, 0, 100], [10, 0, 100], [10, 10, 100], [0, 10, 100], [0, 0, 100]],
+    [[2, 2, 200], [8, 2, 200], [8, 8, 200], [2, 8, 200], [2, 2, 200]],
+])
+
+# 方式二：List[np.ndarray]
+arr1 = np.array([[0, 0, 100], [10, 0, 100], [10, 10, 100], [0, 10, 100], [0, 0, 100]])
+arr2 = np.array([[2, 2, 200], [8, 2, 200], [8, 8, 200], [2, 8, 200], [2, 2, 200]])
+result = SurfaceModeling.create([arr1, arr2])
+
+# 方式三：List[Line Pydantic 模型]
+line1 = Line(geometry=[[0, 0, 100], [10, 0, 100], [10, 10, 100], [0, 10, 100], [0, 0, 100]])
+line2 = Line(geometry=[[2, 2, 200], [8, 2, 200], [8, 8, 200], [2, 8, 200], [2, 2, 200]])
+result = SurfaceModeling.create([line1, line2])
+```
+
+### 从 DMF 文件读取等高线进行地表建模
+
+```python
+from dimine_python_sdk.lib.io import DmfFile
+from dimine_python_sdk.lib.exploitation import SurfaceModeling
+
+# 从 DMF 文件加载等高线
+dmf = DmfFile("path/to/contours.dmf")
+
+# 提取所有多段线实体作为等高线
+contour_lines = []
+for layer in dmf.layers:
+    for entity in layer.entities:
+        if hasattr(entity, 'geometry'):
+            contour_lines.append(entity.geometry)
+
+# 执行地表建模
+surfaces = SurfaceModeling.create(contour_lines, mode=0)
+
+print(f"从 {len(contour_lines)} 条等高线生成 {len(surfaces)} 个地表模型")
+```
+
+### 参数说明
+
+| 参数 | 类型 | 默认值 | 说明 |
+|------|------|--------|------|
+| `contour_lines` | `List[Line]` / `List[np.ndarray]` / `List[List[List[float]]]` | — | 等高线列表，每条为三维点序列 |
+| `mode` | `int` | `0` | 建模模式，默认 0 |
+| `vertical_angle_threshold` | `float` | `25.0` | 伪影判定阈值（度）：三角面法向量与水平面夹角小于此值则判定为伪影 |
+| `vertical_iterations` | `int` | `10` | 伪影剥离最大迭代次数，某轮无新伪影则提前退出 |
+
+### 伪影移除原理
+
+Delaunay 三角网在凸包边界处会生成近似垂直的伪影三角形。`SurfaceModeling.create`
+自动通过迭代剥离移除这些伪影：
+
+1. 找出严格边界三角形（包含仅被一个三角形引用的边）
+2. 扩展为边界区域（边界三角形 + 与之共享边的邻接三角形）
+3. 对边界区域内三角面做法向量检查：\|法向量Z\| < sin(threshold) → 移除
+4. 移除后形成新边界，重复步骤 1-3，直到某轮无新伪影或达到最大迭代次数
+
+```python
+# 默认参数（推荐）
+surfaces = SurfaceModeling.create(contour_lines, mode=0)
+
+# 调整阈值（更激进 / 更保守）
+surfaces = SurfaceModeling.create(contour_lines, mode=0, vertical_angle_threshold=30.0)
+```
+
+## 7. 算法类型
 
 通过 `ContourAlgorithm` 枚举指定两轮廓线之间的连接算法：
 
@@ -354,19 +452,21 @@ models = ContourModeling.two_contour_modeling(
 )
 ```
 
-## 7. 异常处理
+## 8. 异常处理
 
 轮廓线建模定义了层次化异常，建议调用时捕获具体异常：
 
 ```python
 from dimine_python_sdk.lib.exploitation import (
     ContourModeling,
+    SurfaceModeling,
     ContourAlgorithm,
     ModelingError,
     TwoContourModelingError,
     SingleContourModelingError,
     ContourExtrapolateModelingError,
     ParallelContourModelingError,
+    SurfaceModelingError,
 )
 
 try:
@@ -377,6 +477,12 @@ except TwoContourModelingError as e:
     print(f"两轮廓线建模失败: {e}")
 except ModelingError as e:
     print(f"轮廓线建模错误: {e}")
+
+# 地表建模异常处理
+try:
+    surfaces = SurfaceModeling.create(contour_lines, mode=0)
+except SurfaceModelingError as e:
+    print(f"地表建模失败: {e}")
 ```
 
 异常继承关系：
@@ -387,10 +493,11 @@ RuntimeError
       ├── TwoContourModelingError         (两轮廓线建模失败)
       ├── SingleContourModelingError      (单轮廓线封闭为面失败)
       ├── ContourExtrapolateModelingError (轮廓线外推建模失败)
-      └── ParallelContourModelingError    (多平行轮廓线建模失败)
+      ├── ParallelContourModelingError    (多平行轮廓线建模失败)
+      └── SurfaceModelingError            (地表建模失败)
 ```
 
-## 8. 提交到 Dimine 场景
+## 9. 提交到 Dimine 场景
 
 生成的三角网模型可以通过 `conn` 模块提交到 Dimine 场景中显示：
 
@@ -435,7 +542,7 @@ if __name__ == "__main__":
     asyncio.run(main())
 ```
 
-## 9. 完整示例
+## 10. 完整示例
 
 ```python
 import numpy as np

dimine_python_sdk/lib/exploitation/__init__.py

@@ -8,6 +8,8 @@ from dimine_python_sdk.lib.exploitation.contour_modeling import (
     ModelingError,
     ParallelContourModelingError,
     SingleContourModelingError,
+    SurfaceModeling,
+    SurfaceModelingError,
     TwoContourModelingError,
 )
 from dimine_python_sdk.lib.exploitation.implicit_modeling import (
@@ -31,6 +33,7 @@ from dimine_python_sdk.lib.exploitation.models import (
     ImplicitStratumModelingParam,
     LanewayParam,
     SectionParam,
+    SurfaceModelingParam,
 )
 
 __all__ = [
@@ -45,10 +48,12 @@ __all__ = [
     "TwoContourModelingError",
     "SingleContourModelingError",
     "ContourExtrapolateModelingError",
+    "SurfaceModelingError",
     "ParallelContourModelingError",
     "ContourAlgorithm",
     "ExtrapolateMode",
     "ExtrapolateDirection",
+    "SurfaceModeling",
     "ImplicitModeling",
     "ImplicitModelingError",
     "ImplicitStratumModelingError",
@@ -56,4 +61,5 @@ __all__ = [
     "ImplicitGridParam",
     "ImplicitEstimationParam",
     "ImplicitStratumModelingParam",
+    "SurfaceModelingParam",
 ]

dimine_python_sdk/lib/exploitation/contour_modeling.py

@@ -23,6 +23,7 @@ from dimine_python_sdk.lib.native.exploitation import (
     contour_extrapolate_modeling as _native_contour_extrapolate_modeling,
     multi_parallel_contour_modeling as _native_multi_parallel_contour_modeling,
     single_contour_modeling as _native_single_contour_modeling,
+    surface_modeling as _native_surface_modeling,
     two_contour_modeling as _native_two_contour_modeling,
     unite_polydata as _native_unite_polydata,
 )
@@ -35,12 +36,15 @@ __all__ = [
     "SingleContourModelingError",
     "ContourExtrapolateModelingError",
     "ParallelContourModelingError",
+    "SurfaceModelingError",
     "ContourModeling",
+    "SurfaceModeling",
     "two_contour_modeling_impl",
     "single_contour_modeling_impl",
     "contour_extrapolate_modeling_impl",
     "unite_polydata_impl",
     "multi_parallel_contour_modeling_impl",
+    "surface_modeling_impl",
 ]
 
 
@@ -64,11 +68,124 @@ class ParallelContourModelingError(ModelingError):
     """多平行轮廓线建模失败"""
 
 
+class SurfaceModelingError(ModelingError):
+    """地表建模失败"""
+
+
 def _native_tin_to_tingeometry(tin: NativeTin) -> TINGeometry:
     """将 NativeTin 转换为 TINGeometry"""
     return TINGeometry(points=tin.points, faces=tin.faces)
 
 
+def _remove_vertical_triangles(
+    tin: NativeTin,
+    angle_threshold: float = 25.0,
+    iterations: int = 10,
+) -> NativeTin:
+    """移除地表模型中凸包边界附近的伪影三角形。
+
+    这些三角形出现在 Delaunay 三角网凸包边界附近，由边界边连接远距离
+    顶点形成近似垂直的三角面，并非真实地形表面的一部分。
+
+    每次迭代：
+    1. 找出严格边界三角形（包含只被一个三角形引用的边）
+    2. 边界三角形 + 与之共享边的邻接三角形共同构成「边界区域」
+    3. 对边界区域内的三角形做法向量角度检查（默认 15°）
+    4. 移除不合规三角形，形成新的边界
+    5. 重复 1-4，逐步向内剥离伪影层
+
+    Parameters
+    ----------
+    tin : NativeTin
+        待处理的三角网模型。
+    angle_threshold : float
+        法向量与水平面的夹角阈值（度），默认 15.0。
+    iterations : int
+        最大迭代次数，默认 3。若某轮未移除任何三角形则提前终止。
+
+    Returns
+    -------
+    NativeTin
+        剔除伪影后的三角网模型，顶点已重新编号。
+    """
+    if len(tin.faces) == 0:
+        return tin
+
+    points = np.asarray(tin.points, dtype=float)
+    faces = np.asarray(tin.faces, dtype=int)
+    min_nz = np.sin(np.radians(angle_threshold))
+
+    for _ in range(iterations):
+        n_faces = len(faces)
+        if n_faces == 0:
+            break
+
+        # --- 构建边→三角形映射，找出边界边和严格边界三角形 ---
+        edge_to_tris: dict[tuple, list] = {}
+        for fi, f in enumerate(faces):
+            i0, i1, i2 = int(f[0]), int(f[1]), int(f[2])
+            for e in [(min(i0, i1), max(i0, i1)),
+                      (min(i1, i2), max(i1, i2)),
+                      (min(i2, i0), max(i2, i0))]:
+                edge_to_tris.setdefault(e, []).append(fi)
+
+        boundary_tri_set: set = set()
+        for tris in edge_to_tris.values():
+            if len(tris) == 1:
+                boundary_tri_set.update(tris)
+
+        if not boundary_tri_set:
+            break
+
+        # --- 扩展边界区域：纳入与边界三角形共享边的邻接三角形 ---
+        boundary_region: set = set(boundary_tri_set)
+        for tri_idx in list(boundary_tri_set):
+            f = faces[tri_idx]
+            i0, i1, i2 = int(f[0]), int(f[1]), int(f[2])
+            for e in [(min(i0, i1), max(i0, i1)),
+                      (min(i1, i2), max(i1, i2)),
+                      (min(i2, i0), max(i2, i0))]:
+                for neighbor in edge_to_tris.get(e, []):
+                    boundary_region.add(neighbor)
+
+        region_indices = np.array(sorted(boundary_region), dtype=int)
+
+        # --- 对边界区域做法向量检查 ---
+        v0 = points[faces[region_indices, 0]]
+        v1 = points[faces[region_indices, 1]]
+        v2 = points[faces[region_indices, 2]]
+
+        normals = np.cross(v1 - v0, v2 - v0)
+        mag = np.linalg.norm(normals, axis=1, keepdims=True)
+        mag[mag == 0] = 1e-30
+        normals = normals / mag
+
+        region_keep = np.abs(normals[:, 2]) >= min_nz  # |nz| < sin(θ) → 法向量接近水平 → 三角面接近垂直
+        keep_mask = np.ones(n_faces, dtype=bool)
+        keep_mask[region_indices] = region_keep
+
+        kept_faces = faces[keep_mask]
+
+        if len(kept_faces) == n_faces:
+            break  # 本轮没有移除任何三角形，提前终止
+
+        faces = kept_faces
+
+    if len(faces) == 0:
+        return tin
+
+    # --- 顶点重新编号 ---
+    used_indices = np.unique(faces)
+    reindex_map = {old: new for new, old in enumerate(used_indices)}
+    reindexed_faces = np.array(
+        [[reindex_map[int(i)] for i in face] for face in faces],
+        dtype=int,
+    )
+    kept_points = points[used_indices]
+
+    return NativeTin(points=kept_points, faces=reindexed_faces)
+
+
 def _tingeometry_to_native_tin(tin: TINGeometry) -> NativeTin:
     """将 TINGeometry 转换为 NativeTin"""
     return NativeTin(
@@ -325,3 +442,76 @@ def multi_parallel_contour_modeling_impl(
     except Exception as exc:
         raise ParallelContourModelingError(f"多平行轮廓线建模失败: {exc}") from exc
     return [_native_tin_to_tingeometry(item) for item in result]
+
+
+class SurfaceModeling:
+    """
+    地表建模 — 从等高线生成三角网地表（DTM）
+
+    推荐用法（无需实例化）
+    --------------------
+    >>> lines = [line1, line2, line3]  # 等高线列表
+    >>> SurfaceModeling.create(lines, mode=0)
+    """
+
+    @staticmethod
+    def create(contour_lines, mode: int = 0, *,
+               vertical_angle_threshold: float = 25.0,
+               vertical_iterations: int = 10) -> List[TINGeometry]:
+        """
+        地表建模：从等高线生成三角网地表
+
+        Parameters
+        ----------
+        contour_lines : List[Line] or List[np.ndarray] or List[List[List[float]]]
+            等高线列表，每条为三维点序列。
+        mode : int
+            建模模式，默认 0。
+        vertical_angle_threshold : float
+            判定伪影三角形的法向量夹角阈值（度），默认 25.0。
+        vertical_iterations : int
+            伪影剥离最大迭代次数，默认 10。某轮无新伪影则提前退出。
+
+        Returns
+        -------
+        List[TINGeometry]
+            生成的三角网地表列表。
+        """
+        return surface_modeling_impl(
+            contour_lines, mode,
+            vertical_angle_threshold=vertical_angle_threshold,
+            vertical_iterations=vertical_iterations,
+        )
+
+
+def surface_modeling_impl(contour_lines, mode: int = 0, *,
+                         vertical_angle_threshold: float = 25.0,
+                         vertical_iterations: int = 10) -> List[TINGeometry]:
+    """地表建模：输入等高线，输出三角网地表
+
+    Parameters
+    ----------
+    contour_lines : list
+        等高线列表。
+    mode : int
+        建模模式，默认 0。
+    vertical_angle_threshold : float
+        判定伪影三角形的法向量夹角阈值（度），默认 25.0。
+    vertical_iterations : int
+        伪影剥离最大迭代次数，默认 10。某轮无新伪影则提前退出。
+    """
+    normalized = [_normalize_line(line) for line in contour_lines]
+    try:
+        result = _native_surface_modeling(normalized, mode=mode)
+    except Exception as exc:
+        raise SurfaceModelingError(f"地表建模失败: {exc}") from exc
+
+    tins: List[TINGeometry] = []
+    for item in result:
+        item = _remove_vertical_triangles(
+            item,
+            angle_threshold=vertical_angle_threshold,
+            iterations=vertical_iterations,
+        )
+        tins.append(_native_tin_to_tingeometry(item))
+    return tins

dimine_python_sdk/lib/exploitation/models.py

@@ -261,3 +261,19 @@ def _normalize_polylines_set(data) -> List[List[List[float]]]:
     raise TypeError(
         f"不支持的 polylines_set 类型: {type(data)}，支持 List[Line]、List[np.ndarray]、List[List[List[float]]]"
     )
+
+
+class SurfaceModelingParam(ParamModel):
+    """
+    地表建模参数
+
+    示例
+    ------
+    >>> SurfaceModelingParam(mode=0)
+    """
+
+    mode: int = Field(
+        default=0,
+        ge=0,
+        description="建模模式，0=默认模式",
+    )

dimine_python_sdk/lib/io/dmf_file.py

@@ -221,17 +221,14 @@ def _model_entity_to_native(entity: Any) -> NativeEntity:
             properties=props,
         )
     if isinstance(entity, TextInfo):
-        text_props = [p for p in props if p.name != "text"]
-        text_props.append(
-            NativeProperty(name="text", type="string", value=entity.text)
-        )
-        return NativeEntity(
-            entity_type="text",
-            name=name,
-            feature_name=feature_name,
-            geometry=np.asarray(entity.position, dtype=float),
+        normal = getattr(entity, "normal", [0, 0, 1]) or [0, 0, 1]
+        return NativeText(
+            position=np.asarray(entity.position, dtype=float),
+            text=entity.text,
             color=color,
-            properties=text_props,
+            height=entity.size,
+            normal=np.asarray(normal, dtype=float),
+            rotation=getattr(entity, "rotate", 0.0) or 0.0,
         )
     raise DmfInsertError(f"不支持的实体类型: {type(entity).__name__}")
 
@@ -271,6 +268,10 @@ def _model_layer_to_native(layer: DmfLayer) -> NativeLayer:
     entities = []
     for idx, entity in enumerate(layer.entities):
         ne = _model_entity_to_native(entity)
+        # NativeText 没有 name/feature_name 属性，跳过命名逻辑
+        if isinstance(ne, NativeText):
+            entities.append(ne)
+            continue
         # 兼容旧行为：未命名实体使用 "feature-index" 作为实体名称
         if not ne.name:
             ne.name = f"{ne.feature_name or '0'}-{idx}"
@@ -557,11 +558,12 @@ class DmfFile:
         *,
         name: str = "",
         color: Optional[List[int]] = None,
+        size: float = 20,
         properties: Optional[List[PropertyValue]] = None,
     ) -> None:
         """向指定图层插入文本标注"""
         _, lyr = self._resolve_layer(layer)
-        lyr.insert_text(position, text, name=name, color=color)
+        lyr.insert_text(position, text, name=name, color=color, size=size)
         if properties:
             lyr.entities[-1].properties = list(properties)
 

dimine_python_sdk/lib/io/models.py

@@ -300,7 +300,7 @@ class DmfLayer(LayerInfo):
         tin = TINGeometry(points=pts[:, :3], faces=fcs)
         self.entities.append(Shell(feature=name, geometry=tin, color=color))
 
-    def insert_text(self, position: np.ndarray, text: str, *, name: str = "", color: Optional[List[int]] = None) -> None:
+    def insert_text(self, position: np.ndarray, text: str, *, name: str = "", color: Optional[List[int]] = None, size: float = 20) -> None:
         """向本图层插入文本标注"""
         from dimine_python_sdk.lib.io.dmf_file import DmfInsertError
 
@@ -308,4 +308,4 @@ class DmfLayer(LayerInfo):
         if pos.shape[0] < 3:
             raise DmfInsertError(f"文本位置至少需要3个分量，当前: {pos.shape[0]}")
 
-        self.entities.append(TextInfo(feature=name, text=text, position=pos[:3].tolist(), color=color or [255, 255, 255]))
+        self.entities.append(TextInfo(feature=name, text=text, position=pos[:3].tolist(), color=color or [255, 255, 255], size=size))

dimine_python_sdk/lib/native/__init__.py

@@ -82,6 +82,7 @@ from dimine_python_sdk.lib.native.exploitation import (
     laneway_modeling,
     multi_parallel_contour_modeling,
     single_contour_modeling,
+    surface_modeling,
     two_contour_modeling,
     unite_polydata,
 )
@@ -161,4 +162,5 @@ __all__ = [
     "contour_extrapolate_modeling",
     "multi_parallel_contour_modeling",
     "unite_polydata",
+    "surface_modeling",
 ]

dimine_python_sdk/lib/native/exploitation.py

@@ -11,6 +11,7 @@ from typing import Any
 from dimine_python_sdk.lib.native._core import Dm
 from dimine_python_sdk.lib.native._base import (
     NativeExploitationError,
+    ensure_points,
     parse_json_response,
 )
 from dimine_python_sdk.lib.native.models import NativeLithologyModel, NativeTin
@@ -159,6 +160,66 @@ def implicit_stratum_modeling(param: dict | str) -> list[NativeLithologyModel]:
     ]
 
 
+def surface_modeling(contour_lines: list, mode: int = 0) -> list[NativeTin]:
+    """
+    地表建模：从等高线生成三角网地表（DTM）。
+
+    将输入的等高线写入临时 DMF 文件，调用 C++ SurfaceModeling 生成地表
+    三角网，读取结果后清理临时文件。
+
+    Parameters
+    ----------
+    contour_lines : list of np.ndarray
+        等高线列表，每条为 shape (N, 3) 的 numpy 数组。
+    mode : int
+        建模模式，默认 0。
+
+    Returns
+    -------
+    list[NativeTin]
+        生成的三角网地表列表。
+    """
+    import json
+    import os
+    import shutil
+    import tempfile
+
+    import numpy as np
+
+    from dimine_python_sdk.lib.native.io import write_dmf
+    from dimine_python_sdk.lib.native.models import NativeEntity, NativeLayer
+
+    # 构建 NativeLayer，每条等高线为一个 polyline 实体
+    entities: list[NativeEntity] = []
+    for i, line in enumerate(contour_lines):
+        pts = np.asarray(line, dtype=float)
+        if pts.ndim != 2 or pts.shape[1] < 2:
+            raise NativeExploitationError(
+                f"等高线 {i} 格式无效: shape={pts.shape}，期望 (N, 3)"
+            )
+        entities.append(
+            NativeEntity(
+                entity_type="polyline",
+                name=f"contour_{i}",
+                geometry=ensure_points(pts),
+            )
+        )
+    layer = NativeLayer(name="contours", entities=entities)
+
+    tmp_dir = tempfile.mkdtemp(prefix="dm_surface_")
+    tmp_file = os.path.join(tmp_dir, "surface_input.dmf")
+    try:
+        write_dmf(tmp_file, [layer])
+        param = json.dumps({"file": tmp_file, "mode": int(mode)}, ensure_ascii=False)
+        try:
+            result = Dm.SurfaceModeling(param)
+        except Exception as exc:
+            raise NativeExploitationError(f"地表建模失败: {exc}") from exc
+        return _convert_polydata_list(result, "地表建模")
+    finally:
+        shutil.rmtree(tmp_dir, ignore_errors=True)
+
+
 __all__ = [
     "get_section_contour",
     "laneway_modeling",
@@ -168,4 +229,5 @@ __all__ = [
     "multi_parallel_contour_modeling",
     "unite_polydata",
     "implicit_stratum_modeling",
+    "surface_modeling",
 ]