knowledge(요구): 요구 557 implicit 시험 결과 — 20만 점이면 매끈 + 세부 지킴 · 누락 = 외곽 마스크(DJI 지면점 없는 띠 36 %) · 빈 띠 메우기는 정답 없음 · 외부 값과 비교 됨(폭)

Co-Authored-By: Claude Opus 5.5 <noreply@anthropic.com>
2026-11-01 19:35:20 +09:00
parent f27587b823
commit a1c33a2d9d
1 changed files with 101 additions and 0 deletions
@@ -0,0 +1,101 @@
+# 557 결과 — implicit 서피스 품질 · 누락 · 외부 무료 라이브러리 값과 비교
+
+조사만 · 앱 코드 · storage · git 손대지 않음. 앱 implicit 코드(`build_implicit` · `TerrainContext` · 설정 · CSF)는 랩탑 메인(origin/main_laptop_1)과 같음 확인 → 이 폴더 앱 함수를 그대로 불러 씀. 자료 = 용화.las(DJI L1 · 4,900.7 만 점). 모두 같은 1 m 격자(dtm_classification x · y). 비교군 A0 · A1 · B1 · B2 · C2 는 553 결과 그대로(다시 안 돌림).
+
+- I1 = implicit(업체 분류 지면점) · I2 = implicit(앱 내장 CSF 지면점) · 뒤 숫자 = RBF 에 쓰는 점 수(앱 기본 10k · 50k · 200k · all = 전부).
+- 검사점(업체 지면 20 %)은 RBF 점에서 뺌(553 B 방법과 같게) · 외곽(footprint)은 앱처럼 지면점 전체로.
+- 범위 둘 — **업체 외곽** = 553 기준 범위(업체 지면점 외곽 · 23.2 만 칸) · **LAS 점 범위** = 점이 있는 곳 전체(같은 외곽 짓기 · 36.1 만 칸).
+
+## ① 결론 다섯 줄
+
+1. **implicit 를 기본으로 써도 됨 — 다만 점 수를 올려서.** 앱 기본 1 만 점은 너무 적어 깎기 비탈 끝 · 골 바닥을 1.4 m 까지 뭉갬. 20 만 점이면 등고선은 여전히 매끄럽고(꺾임 19/km — 지금 DTM 188 · TIN 18) 업체 지면점 · 외부 TIN 과 거의 같음(B1 과 RMS 0.18 m). 시간은 +25 초뿐(짓기 42 초 · 1 m 셈 47 초).
+2. **누락의 주인은 implicit 가 아니라 「외곽 마스크」.** 업체 외곽 안에서는 implicit 빈칸 0 % · 계산 실패 0. 그러나 LAS 점이 있는 범위의 **36 %** 가 비어 있음 — 업체(DJI Terra)가 지면점을 하나도 안 준 띠(주로 북 · 동 가장자리 숲 · 칸당 점 37 개로 적지 않음). 지금 DTM 도 똑같이 36 % 빔. 앱 CSF 로 가르면 11 % 로 줆.
+3. **메우는 길 — 「자료 없음」 표시가 가장 정직.** 빈 띠를 RBF 로 끝까지 세거나(a) 앱 DTM 으로 메워도(b) 그 자리에선 외부 도구끼리도 값이 크게 갈림(SMRF − WBT 평균 1.9 m · RMS 4.3 m). 어느 방법이든 그 띠 칸의 6~13 % 가 그 칸 가장 낮은 레이저점보다 0.5 m 넘게 떠 있음(= 확실히 땅 위). 메우면 「추정」으로 따로 칠해야 함.
+4. **외부 라이브러리와 비교 — 됨.** 아래 표 ③ 칸이 각 방법 − B1(업체 분류 + PDAL TIN) · − B2(PDAL SMRF) · − C2(WhiteboxTools)의 평균 · RMS · 95 %. 단 외부 셋끼리도 0.55~1.4 m RMS 로 다름 → 「외부 값」은 정답 하나가 아니라 폭. I1 20 만 점은 그 폭 안(B1 과 0.18 · B2 와 0.57 · C2 와 1.39 m).
+5. **권함** — 분류 있는 LAS: implicit + 점 수 20 만(설정 키 `implicit_max_points_per_tile` 하나) · 빈 띠는 「자료 없음」. 분류 없는 LAS(CSF + implicit): implicit 는 매끈하게 덮지만 **CSF 가 잘못 고른 점을 그대로 따라감**(점을 늘릴수록 혹이 늘어 꺾임 6 → 67/km) → 「등고선이 매끄러우니 품질만 보면 됨」은 1 만 점일 때의 착시. CSF 바꾸기(553 의 SMRF) 가 먼저.
+
+## ② 잣대 표 — 업체 외곽 안(553 기준 범위)
+
+「정답」 순위 아님. 잣대마다 기준이 다름 — ②는 업체 분류를 정답으로 둔 참고값 · ③은 외부 도구 값과의 거리 · ④는 꼴.
+
+| 방법 | ① 빈칸 % 업체 외곽 / LAS 점 범위 | ② 검사점 RMS · 평균 차(참고) m | ③ − B1 평균 · RMS · 95 % | ③ − B2 평균 · RMS · 95 % | ③ − C2 평균 · RMS · 95 % | ④ 등고선 km · 꺾임/km | ④ 기울기 p50 / p90 % · 빨강 % | ⑤ 시간 s · 메모리 GB |
+| --- | --- | --- | --- | --- | --- | --- | --- | --- |
+| A0 지금 DTM | 0 / 35.6 | 0.767 · +0.601 | −0.49 · 0.83 · 1.59 | −0.66 · 1.08 · 2.08 | −0.82 · 1.71 · 3.32 | 189 · 188 | 69 / 120 · 32.0 | 앱 WF1 안 0.5 |
+| A1 지금 평활본 | 2.2※ / 37.0 | 0.747 · +0.602 | −0.50 · 0.68 · 1.26 | −0.67 · 0.97 · 1.82 | −0.83 · 1.65 · 3.32 | 163 · 70 | 68 / 104 · 25.7 | — |
+| B1 업체 분류 + PDAL TIN | 0.02 / 0.04 | 0.117 · −0.008 | 0 | −0.17 · 0.55 · 1.12 | −0.33 · 1.37 · 2.73 | 159 · 18 | 70 / 89 · 16.3 | 20 · 3.6 |
+| B2 PDAL SMRF + TIN | 0.03 / 0.02 | 0.265 · −0.108 | +0.17 · 0.55 · 1.12 | 0 | −0.16 · 1.36 · 2.66 | 169 · 70 | 70 / 98 · 22.5 | 58 · 6.7 |
+| C2 WBT 경사 필터 + TIN | 0.02 / 0.02 | 0.257 · +0.056 | +0.33 · 1.37 · 2.73 | +0.16 · 1.36 · 2.66 | 0 | 200 · 100 | 72 / 130 · 27.5 | 453 · 5.4 |
+| **I1 10k(앱 기본)** | 0 / 35.6 | 0.273 · +0.003 | −0.03 · 0.40 · 0.83 | −0.20 · 0.66 · 1.35 | −0.36 · 1.45 · 2.96 | 157 · 4 | 69 / 86 · 13.3 | 29 + 31 · 1.45 |
+| I1 50k | 0 / 35.6 | 0.160 · −0.002 | −0.02 · 0.26 · 0.49 | −0.19 · 0.59 · 1.19 | −0.35 · 1.41 · 2.82 | 158 · 10 | 70 / 88 · 15.5 | 36 + 42 · 1.49 |
+| **I1 200k** | 0 / 35.6 | 0.118 · −0.004 | −0.01 · 0.18 · 0.35 | −0.17 · 0.57 · 1.15 | −0.33 · 1.39 · 2.75 | 160 · 19 | 70 / 90 · 17.2 | 42 + 47 · 1.52 |
+| I1 all(91 만) | 0 / 35.6 | 0.101 · −0.007 | +0.01 · 0.13 · 0.22 | −0.16 · 0.57 · 1.14 | −0.32 · 1.38 · 2.71 | 162 · 29 | 70 / 91 · 18.8 | 44 + 48 · 1.58 |
+| **I2 10k(앱 CSF · 기본)** | 0.17 / 10.9 | 0.330 · +0.014 | +0.01 · 0.50 · 1.01 | −0.16 · 0.68 · 1.42 | −0.32 · 1.42 · 2.83 | 157 · 6 | 69 / 87 · 14.9 | 33 + 32 · 1.46 |
+| I2 50k | 0.17 / 10.9 | 0.223 · +0.011 | +0.04 · 0.47 · 0.85 | −0.13 · 0.64 · 1.32 | −0.29 · 1.34 · 2.63 | 161 · 22 | 70 / 91 · 18.7 | 40 + 43 · 1.49 |
+| I2 200k | 0.17 / 10.9 | 0.178 · +0.009 | +0.07 · 0.53 · 0.90 | −0.10 · 0.65 · 1.28 | −0.26 · 1.27 · 2.45 | 166 · 47 | 71 / 95 · 22.2 | 43 + 51 · 1.51 |
+| I2 all(64 만) | 0.17 / 10.9 | 0.165 · +0.009 | +0.10 · 0.58 · 1.02 | −0.07 · 0.66 · 1.29 | −0.22 · 1.22 · 2.31 | 172 · 67 | 71 / 99 · 23.9 | 44 + 48 · 1.57 |
+
+- ① 빈칸 까닭(LAS 점 범위 기준) — I1: 외곽 마스크 · 칸 둘레에 업체 지면점 없음 20.3 % + 지면점 드문 곳을 외곽 다듬기(틈 메우기 1 m · 안쪽 1 m 깎기 · 가장 큰 덩이만)가 버림 15.3 % · 계산 실패(NaN) 0 · 등고선 화면만의 추가 빈칸(뽑힌 점 범위로 격자 축소) 0~0.01 %. I2: 4.9 % + 6.1 % · 나머지 0.
+- ⑤ 시간 = implicit 짓기(앱 `build_implicit` · 미리보기 격자 셈 포함) + 1 m 격자 셈(종횡단 sampler 꼴). 메모리 = 프로세스 최대(점 4,900 만 개를 미리 올린 1.4 GB 포함 · implicit 몫은 +0.1~0.2 GB). B · C 는 553 기록(LAS 읽기부터 따로 실행 · 그대로 비교 못함).
+- 등고선 = 1 m 간격 총 길이 · 꺾임 = 2 m 간격으로 다시 찍어 30° 넘게 꺾이는 곳 수 / km(작을수록 매끄러움 · 너무 작으면 세부를 지운 것).
+- 빨강 = 앱 급경사 셋째 범위 경계 83.9 % 이상 칸 비율.
+- ※ A1 2.2 % 는 평활 격자(1.078 m)를 1 m 로 옮길 때 가장자리 빠짐 — 실제 빈칸 아님.
+
+### 단면 셋 — 업체 지면점 원값 대비(±0.5 m 띠) · 최대 차 / RMS m
+
+| 단면 | A0 | A1 | B1 | B2 | C2 | I1 10k | I2 10k | I1 200k |
+| --- | --- | --- | --- | --- | --- | --- | --- | --- |
+| S1 임도 노면(40 m) | 4.20 / 1.63 | 4.05 / 1.60 | 0.50 / 0.11 | 0.46 / 0.10 | 1.49 / 0.36 | 0.77 / 0.15 | 2.11 / 0.39 | 0.33 / 0.08 |
+| S2 깎기 비탈(34 m) | 3.30 / 1.36 | 3.40 / 1.33 | 0.54 / 0.11 | 0.46 / 0.10 | 1.07 / 0.26 | 1.40 / 0.28 | 0.74 / 0.18 | 0.41 / 0.11 |
+| S3 숲 골(65 m · 깊이 14 m) | 1.81 / 0.57 | 1.45 / 0.53 | 0.33 / 0.13 | 1.75 / 0.29 | 1.66 / 0.27 | 1.41 / 0.43 | 1.97 / 0.54 | 0.58 / 0.14 |
+
+### 메우기 셋 — 업체 외곽 밖 · LAS 점 범위 안(빈 띠 12.8 만 칸)만
+
+외곽 안 값은 셋 다 그대로(② · ③ 같음). 빈 띠엔 정답이 없음 → 외부 도구와의 차 · 「뜬 칸」(그 칸 가장 낮은 레이저점보다 0.5 m 넘게 높음 = 땅 위가 확실) 으로만 봄.
+
+| 방법 | 덮임 % | − B2 평균 · RMS m | − C2 평균 · RMS m | 뜬 칸 % |
+| --- | --- | --- | --- | --- |
+| (c) 메우지 않음 · I1 | 0 | — | — | — |
+| (a) I1 10k RBF 끝까지 | 100 | −0.29 · 3.22 | −2.17 · 5.06 | 12.8 |
+| (b) I1 10k + 앱 DTM(가까운 값) | 100 | −0.34 · 3.06 | −2.22 · 4.83 | 6.2 |
+| (a) I1 200k RBF 끝까지 | 100 | 0.00 · 2.59 | −1.88 · 4.43 | 11.1 |
+| (c) I2 10k 메우지 않음(CSF 가 이미 덮은 곳) | 69.7 | +0.26 · 3.23 | −1.28 · 3.84 | 15.0 |
+| (a) I2 10k RBF 끝까지 | 100 | +0.35 · 3.25 | −1.53 · 4.17 | 13.4 |
+| (b) I2 10k + 앱 CSF DTM | 100 | +0.30 · 3.53 | −1.58 · 4.04 | 10.9 |
+| 참고: B1(업체 분류 TIN 이 50 m 까지 건너 이음) | 99.9 | +0.03 · 2.53 | −1.85 · 4.35 | 10.8 |
+| 참고: B2 · C2 자기 지면 | 100 | B2 − C2 = −1.88 · 4.32 | | 11.0 · 11.3 |
+
+→ 빈 띠는 어느 방법이든 외부 값과 RMS 2.5~5 m. 메움 b(가까운 값)가 뜬 칸은 가장 적으나 등고선이 계단(꺾임 94/km). 메움 a 는 매끈하지만 그럴듯한 틀린 값(뜬 칸 11~13 %).
+
+## ③ 그림(`figs557/` · `storage/_shared/_lab_557/`)
+
+| 파일 | 보이는 것 |
+| --- | --- |
+| `557_누락_자리.png` | LAS 점 범위 전체. 빨강 = 업체 지면점 없음으로 외곽에서 잘린 칸 · 분홍 = 지면점은 드물게 있으나 외곽 다듬기로 버린 칸 · 파란 선 = 업체 외곽. I1 은 북 · 동 가장자리 넓은 띠가 빔(36 %). I2(앱 CSF)는 띠 대부분을 덮으나 노랑(업체 외곽 안에서 B1 과 1 m 넘게 다른 곳 · CSF 가 고른 점 차이)이 점점이 남. 검은 선 = 단면 셋. |
+| `557_외부와_차이.png` | 6 방법 × 외부 셋(B1 · B2 · C2) 차이 지도(±1.5 m). A0 · A1 은 전면 파랑(0.5 m 낮음). I1 10k 는 흰 바탕에 혹 · 웅덩이 점(점이 성긴 자리) · I1 200k 는 B1 과 거의 흰색. C2 열은 모든 방법에서 파랑 덩이 — C2 자체가 그 자리에 나무 위 점을 땅으로 넣음(C2 쪽 흠). |
+| `557_단면_셋.png` | 검은 점 = 업체 지면점 원값 · 주황 점 = 앱 CSF 지면점. S1 · S2 에서 A0 · A1 은 비탈 전체가 1~4 m 낮게 그려짐 · I1 10k(빨강)는 깎기 비탈 끝 모서리를 둥글게 뭉갬(S2 최대 1.4 m) · I1 200k(보라)는 B1 과 겹침. S3 골은 I1 10k 가 바닥 · 둔덕을 곧은 선으로 지나감 · SMRF · WBT 는 비탈 중간에 덤불 혹(1.7 m). |
+| `557_등고선_가장자리.png` | 임도가 지나는 100 m 창의 1 m 등고선. A0 계단 · 줄무늬 · A1 흐려진 줄무늬 · B1 깔끔 · I1 10k 가장 매끈(임도 모서리도 같이 무뎌짐) · I1 200k 는 B1 과 거의 같은 꼴 · C2 는 둥근 혹 등고선. |
+| `557_메우기.png` | 빈 띠를 메운 값 − B2. a · b 모두 빈 띠에서 ±3 m 넘는 빨강 · 파랑 얼룩 · 오른쪽 빨강 칸 = 레이저점보다 떠 있는 칸(띠 전체에 흩어짐). |
+
+## ④ 점 수 감도(설정 키 `implicit_max_points_per_tile` · 앱엔 타일 나눔 없음 — 이름만 「per tile」)
+
+| 점 수 | 뽑힌 점 사이 가운데 · 95 % m | 칸 → 가장 가까운 점 95 % · 최대 m | I1 검사점 RMS | I1 − B1 RMS | I1 꺾임/km | I2 꺾임/km | 짓기 + 셈 s |
+| --- | --- | --- | --- | --- | --- | --- | --- |
+| 10k(기본) | 1.28 · 6.1 | 8.0 · 16.6 | 0.273 | 0.40 | 4 | 6 | 60 |
+| 50k | 0.45 · 2.6 | 4.4 · 10.1 | 0.160 | 0.26 | 10 | 22 | 78 |
+| 200k | 0.20 · 1.1 | 2.9 · 7.2 | 0.118 | 0.18 | 19 | 47 | 89 |
+| all | 0.09 · 0.4 | 2.1 · 7.0 | 0.101 | 0.13 | 29 | 67 | 92 |
+
+- 점 뽑기가 **파일 순서대로 일정 간격**(자리 고르게가 아님) → 1 만 점이면 점이 비행 줄을 따라 몰려 이웃 사이 가운데 1.3 m 인데 칸에서 가장 가까운 점이 최대 16.6 m · 그 사이를 곡면이 곧게 건넘(골 바닥 · 비탈 끝이 뭉개지는 까닭).
+- I1(업체 분류): 점을 늘리면 꾸준히 좋아지고 200k 부터 B1 과 비슷 · 시간은 거의 안 늚(곡면 셈이 이웃 64 개만 씀).
+- I2(앱 CSF): 점을 늘리면 검사점은 좋아지나 꺾임 · 빨강이 급히 늚(CSF 가 넣은 덤불 · 그루터기 점까지 충실히 따라감).
+- 매끄러움 계수(`implicit_smoothing` 0.1)는 이번에 바꾸지 않음.
+
+## ⑤ 위험
+
+- 사업 하나(용화 · 숲 98 % · 급경사)만 잼. 평지 · 다른 DJI 기종 · 분류 없는 LAS 는 따로.
+- ② 검사점 · ④ 단면 원값은 업체 분류를 정답으로 둔 값 — 참고 잣대. 업체가 놓친 땅 · 잘못 넣은 점은 가르지 못함.
+- 사용자가 화면에서 본 「implicit 누락」이 업체 외곽 **안쪽**이라면 이 측정과 다른 까닭(다른 사업 · 화면 쪽) — 이 자료 · 이 코드로는 외곽 안 빈칸 0. 랩탑 쪽 코드 조사와 맞춰 볼 것.
+- 빈 띠(36 %)에는 정답이 없음 — 「메운 값」의 질은 외부 도구끼리의 차로만 짐작. 설계선이 그 띠를 지나면 어떤 서피스든 m 단위 오차.
+- 1 m 격자에서 잰 꼴 · 등고선 꺾임은 이 시험의 정의 · 앱 화면의 등고선(미리보기 격자 1.075 m)과 조금 다름.
+- I 판은 검사점 22.8 만 개를 RBF 점에서 뺌 — 앱(전부 씀)보다 아주 조금 불리.
+- 점 수를 올리면 저장 파일이 커짐(압축 npz · 10k 0.08 MB → 200k 1.3 MB → all 5 MB) · 종횡단 · 등고선마다 곡면을 다시 세우는 시간(초 단위)은 이번에 따로 안 잼.