Files
Aislo/B04_PreProcess/B04_PreProcess_Schema.py
T
eomsangdonandClaude Opus 5 ea254a741f fix(B04): 지면 필터의 고정 상수를 걷어내 LAS 기복·경사와 무관하게 지면을 잡는다
용화.las(기복 251.9m) 전처리에서 지면점이 0.06%까지 사라졌다. 원인은 현장 하나에
맞춰 박아둔 상수였다.

- CSF 하강 예산이 0.3185m x 150회 = 47.8m로 고정돼, 기복이 그보다 큰 산악지에서
  천이 지면에 닿지 못했다. 반복 수를 collision_grid 기준 필요 하강량에서 뽑고
  SURFACE_CSF_MAX_ITERATIONS로 상한만 둔다.
- CSF 6단계 수목 필터는 비교 피연산자가 뒤집혀 조건이 항상 참이었다. 셀 지면
  후보 대비 높이를 보도록 순서를 바로잡는다.
- grid_min_z의 3x3 minimum_filter가 급경사에서 기준면을 경사만큼 파고들어
  (중앙 1.88m > 임계 1.5m) 지면점을 떨궜다. 제거한다.
- 기준면 배열이 float32라 셀 최저점이 반올림으로 자기 기준면보다 낮아져
  탈락했다. 원본 좌표와 같은 float64로 둔다.
- LAS가 지면분류(class 2)를 싣고 오면 filter_classification을 자동으로 붙인다.
  미분류 LAS에서는 조용히 빠진다.
- 필터별 ground_ratio를 로그에 남기고, 하한 미만이면 WARNING을 띄운다.
  지금까지는 0.06%가 나와도 "계산 완료"로만 보였다.

실측 (build_ground_masks 재계산):
  용화       grid_min_z 0.52% -> 4.91%  (class2 회수율 17.9% -> 93.3%)
             csf        0.06% -> 1.57%  (class2 회수율  1.3% -> 55.1%)
  cloud_merged grid_min_z 11.60% -> 17.51%
             csf         8.84% -> 10.68%

Co-Authored-By: Claude Opus 5 (1M context) <noreply@anthropic.com>
2026-09-01 11:10:45 +09:00

160 lines
4.9 KiB
Python

"""B04 지표면 분석 요청·응답 검증 모델."""
from typing import Any
from pydantic import BaseModel, ConfigDict, Field
from config.config_system import (
SURFACE_MODEL_PRECOMPUTE,
SURFACE_MODEL_SOURCE_FILTERS,
)
# classification은 LAS가 지면분류를 싣고 있을 때 엔진이 자동으로 붙이지만,
# 재계산 요청에서 이름으로 지정하는 것도 허용한다.
_ALLOWED_FILTERS = set(SURFACE_MODEL_SOURCE_FILTERS) | {"ransac", "classification"}
_ALLOWED_METHODS = set(SURFACE_MODEL_PRECOMPUTE)
class SurfaceAnalyzeRequest(BaseModel):
"""지표면 분석 실행 요청."""
model_config = ConfigDict(extra="forbid")
input_file_id: int = Field(gt=0, description="구조화 대상 원본 LAS input_files.id")
source_filters: list[str] | None = Field(
default=None, description="실행할 지면 필터 (미지정 시 config 기본값)"
)
methods: list[str] | None = Field(
default=None, description="생성할 지표면 표현 (미지정 시 config 기본값)"
)
force: bool = Field(default=False, description="캐시 무시 후 강제 재계산")
def resolved_filters(self) -> list[str]:
filters = self.source_filters or list(SURFACE_MODEL_SOURCE_FILTERS)
invalid = [f for f in filters if f not in _ALLOWED_FILTERS]
if invalid:
raise ValueError(f"허용되지 않은 지면 필터입니다: {invalid}")
return filters
def resolved_methods(self) -> list[str]:
methods = self.methods or list(SURFACE_MODEL_PRECOMPUTE)
invalid = [m for m in methods if m not in _ALLOWED_METHODS]
if invalid:
raise ValueError(f"허용되지 않은 지표면 표현입니다: {invalid}")
return methods
class SurfaceConfirmRequest(BaseModel):
"""사용자가 선택한 지표면 모델 확정 요청."""
model_config = ConfigDict(extra="forbid")
model_id: int = Field(gt=0, description="확정할 surface_models.id")
smooth: bool | None = Field(default=None, description="확정 시 적용한 스무딩 여부")
contour_interval_m: float | None = Field(
default=None,
gt=0,
description="확정 시 적용한 등고선 간격(m)",
)
class SurfaceConfirmResponse(BaseModel):
"""지표면 모델 확정 결과."""
status: str = "success"
project_id: str
model_id: int
confirmed: bool = True
class SurfaceModelSummary(BaseModel):
"""저장된 지표면 모델 요약."""
id: int
model_type: str
status: str
resolution_m: float | None = None
model_file_path: str | None = None
generation_params: dict[str, Any] | None = None
created_at: str | None = None
class SurfaceInputFileSummary(BaseModel):
"""WF1 분석 대상으로 선택 가능한 LAS/LAZ 입력 파일 요약."""
id: int
file_type: str
original_filename: str
raw_file_path: str
file_size_mb: float | None = None
crs_epsg: int | None = None
status: str | None = None
created_at: str | None = None
class SurfaceInputFileListResponse(BaseModel):
"""프로젝트 지표면 분석 입력 파일 목록 응답."""
status: str = "success"
project_id: str
files: list[SurfaceInputFileSummary]
class SurfacePointCloudSampleResponse(BaseModel):
"""B04 3D 미리보기용 포인트클라우드 샘플."""
status: str = "success"
project_id: str
point_count: int
sampled_count: int
bounds: dict[str, float]
points: list[list[float]]
rgb: list[list[int]] | None = None
class SurfaceConfirmedResponse(BaseModel):
"""확정 지표면 요약 — 화면 진입 판정·준비화면·B05가 공통으로 쓰는 단일 출처.
포인트 배열 없이 확정값과 지형 가장자리만 담아 수 KB로 유지한다.
signature는 확정 구성이 바뀌었는지 프론트가 한 줄로 비교하기 위한 값이다.
"""
status: str = "success"
project_id: str
model_id: int | None = None
source_filter: str | None = None
method: str | None = None
smooth: bool | None = None
contour_interval_m: float | None = None
signature: str
point_count: int | None = None
bounds: dict[str, float] | None = None
# 계획노선(B03 CSV)의 평면 범위. 지도 초기 화면을 도로 기준으로 맞출 때 쓴다.
route_bounds: dict[str, float] | None = None
class SurfaceGroundStatsResponse(BaseModel):
"""필터별 지면 포인트 통계 응답."""
status: str = "success"
project_id: str
filters: dict[str, dict[str, Any]]
class SurfaceAnalyzeResponse(BaseModel):
"""지표면 분석 실행 결과."""
status: str = "success"
project_id: str
ground_summary: dict[str, Any]
manifest_status: str
surface_model_ids: list[int]
class SurfaceModelListResponse(BaseModel):
"""프로젝트 지표면 모델 목록 응답."""
status: str = "success"
project_id: str
models: list[SurfaceModelSummary]