keyence analysis

%23%20%2F%2F%2F%20script%0A%23%20requires-python%20%3D%20%22%3E%3D3.11%22%0A%23%20dependencies%20%3D%20%5B%0A%23%20%20%20%22marimo%3E%3D0.13%22%2C%0A%23%20%20%20%22pandas%3E%3D2.0%22%2C%0A%23%20%20%20%22matplotlib%3E%3D3.7%22%2C%0A%23%20%20%20%22anywidget%3D%3D0.11.0%22%2C%0A%23%20%20%20%22traitlets%3D%3D5.16.1%22%2C%0A%23%20%5D%0A%23%20%2F%2F%2F%0A%0Aimport%20marimo%0A%0A__generated_with%20%3D%20%220.24.0%22%0Aapp%20%3D%20marimo.App(width%3D%22medium%22)%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20marimo%20as%20mo%0A%20%20%20%20import%20numpy%20as%20np%0A%20%20%20%20import%20pandas%20as%20pd%0A%20%20%20%20import%20matplotlib.pyplot%20as%20plt%0A%20%20%20%20from%20pathlib%20import%20Path%0A%20%20%20%20import%20anywidget%2C%20traitlets%0A%0A%0A%20%20%20%20return%20anywidget%2C%20mo%2C%20np%2C%20pd%2C%20plt%2C%20traitlets%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%20Keyence%202D%20Scan%20%E2%80%94%20Channel%20Height%20Analysis%0A%0A%20%20%20%20Analyze%20height%20measurements%20from%20a%20Keyence%20confocal%20sensor%202D%20scan.%0A%20%20%20%20The%20workflow%20has%20five%20steps%3A%20load%20data%20%E2%86%92%20review%20the%20scan%20grid%20%E2%86%92%20filter%20and%20crop%0A%20%20%20%20to%20a%20region%20of%20interest%20%E2%86%92%20assign%20regions%20to%20sections%20%E2%86%92%20export%20results.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Step%201%20%E2%80%94%20Load%20Data%0A%0A%20%20%20%20**New%20session%3A**%20Upload%20a%20raw%20Keyence%20CSV%20export%20file%20below.%0A%20%20%20%20The%20file%20should%20have%20columns%3A%20%60time%60%2C%20%60measurement%60%20(mm)%2C%20%60x%60%20(mm)%2C%20%60y%60%20(mm).%0A%20%20%20%20If%20the%20file%20contains%20two%20datasets%20separated%20by%20blank%20rows%2C%20only%20the%20first%20is%20loaded.%0A%0A%20%20%20%20**Restore%20session%3A**%20Upload%20the%20three%20session%20files%20saved%20from%20a%20previous%20session%0A%20%20%20%20(%60avg_grid.csv%60%2C%20%60roi_points.csv%60%2C%20%60session_params.csv%60).%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20file_upload%20%3D%20mo.ui.file(%0A%20%20%20%20%20%20%20%20filetypes%3D%5B%22.csv%22%5D%2C%0A%20%20%20%20%20%20%20%20kind%3D%22area%22%2C%0A%20%20%20%20%20%20%20%20label%3D%22**Upload%20Keyence%20CSV%20file**%20(drag-and-drop%20or%20click%20to%20browse)%22%2C%0A%20%20%20%20)%0A%20%20%20%20file_upload%0A%20%20%20%20return%20(file_upload%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(anywidget%2C%20traitlets)%3A%0A%20%20%20%20class%20SessionLoader(anywidget.AnyWidget)%3A%0A%20%20%20%20%20%20%20%20loaded%20%3D%20traitlets.Dict(%7B%7D).tag(sync%3DTrue)%0A%20%20%20%20%20%20%20%20_esm%20%3D%20'%5Cnasync%20function%20render(%7B%20model%2C%20el%20%7D)%20%7B%5Cn%20%20el.style.cssText%20%3D%20%22font-family%3Asans-serif%3Bdisplay%3Aflex%3Bflex-direction%3Acolumn%3Bgap%3A6px%3B%22%3B%5Cn%5Cn%20%20const%20input%20%3D%20document.createElement(%22input%22)%3B%5Cn%20%20input.type%20%3D%20%22file%22%3B%5Cn%20%20input.webkitdirectory%20%3D%20true%3B%5Cn%20%20input.style.display%20%3D%20%22none%22%3B%5Cn%20%20el.appendChild(input)%3B%5Cn%5Cn%20%20const%20btn%20%3D%20document.createElement(%22button%22)%3B%5Cn%20%20btn.textContent%20%3D%20%22%F0%9F%93%82%20Open%20session%20folder%22%3B%5Cn%20%20btn.style.cssText%20%3D%20%22padding%3A6px%2014px%3Bfont-size%3A14px%3Bcursor%3Apointer%3Bborder-radius%3A4px%3Bborder%3A1px%20solid%20%23888%3Bwidth%3Afit-content%3B%22%3B%5Cn%20%20const%20statusDiv%20%3D%20document.createElement(%22div%22)%3B%5Cn%20%20statusDiv.style.cssText%20%3D%20%22font-size%3A13px%3Bmin-height%3A20px%3B%22%3B%5Cn%20%20el.appendChild(btn)%3B%5Cn%20%20el.appendChild(statusDiv)%3B%5Cn%5Cn%20%20btn.onclick%20%3D%20()%20%3D%3E%20%7B%20input.value%20%3D%20%22%22%3B%20input.click()%3B%20%7D%3B%5Cn%5Cn%20%20input.onchange%20%3D%20async%20()%20%3D%3E%20%7B%5Cn%20%20%20%20const%20files%20%3D%20Array.from(input.files)%3B%5Cn%20%20%20%20if%20(!files.length)%20%7B%20statusDiv.textContent%20%3D%20%22%22%3B%20return%3B%20%7D%5Cn%20%20%20%20statusDiv.style.color%20%3D%20%22%22%3B%5Cn%20%20%20%20statusDiv.textContent%20%3D%20%22%E2%8F%B3%20Reading%20files%E2%80%A6%22%3B%5Cn%5Cn%20%20%20%20const%20required%20%3D%20%5B%22avg_grid.csv%22%2C%20%22roi_points.csv%22%2C%20%22session_params.csv%22%5D%3B%5Cn%20%20%20%20const%20loaded%20%3D%20%7B%7D%3B%5Cn%20%20%20%20for%20(const%20name%20of%20required)%20%7B%5Cn%20%20%20%20%20%20const%20f%20%3D%20files.find(f%20%3D%3E%20f.name%20%3D%3D%3D%20name)%3B%5Cn%20%20%20%20%20%20loaded%5Bname%5D%20%3D%20f%20%3F%20await%20f.text()%20%3A%20null%3B%5Cn%20%20%20%20%7D%5Cn%5Cn%20%20%20%20const%20missing%20%3D%20required.filter(n%20%3D%3E%20!loaded%5Bn%5D)%3B%5Cn%20%20%20%20if%20(missing.length%20%3D%3D%3D%200)%20%7B%5Cn%20%20%20%20%20%20statusDiv.textContent%20%3D%20%22%E2%9C%85%20Session%20files%20loaded%20%E2%80%94%20restore%20mode%20active.%22%3B%5Cn%20%20%20%20%20%20statusDiv.style.color%20%3D%20%22green%22%3B%5Cn%20%20%20%20%7D%20else%20%7B%5Cn%20%20%20%20%20%20statusDiv.textContent%20%3D%20%22%E2%9A%A0%EF%B8%8F%20Missing%20in%20selected%20folder%3A%20%22%20%2B%20missing.join(%22%2C%20%22)%3B%5Cn%20%20%20%20%20%20statusDiv.style.color%20%3D%20%22darkorange%22%3B%5Cn%20%20%20%20%7D%5Cn%20%20%20%20model.set(%22loaded%22%2C%20loaded)%3B%5Cn%20%20%20%20model.save_changes()%3B%5Cn%20%20%7D%3B%5Cn%7D%5Cnexport%20default%20%7B%20render%20%7D%3B%5Cn'%0A%0A%0A%20%20%20%20return%20(SessionLoader%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(SessionLoader%2C%20mo)%3A%0A%20%20%20%20session_loader%20%3D%20mo.ui.anywidget(SessionLoader())%0A%20%20%20%20session_loader%0A%20%20%20%20return%20(session_loader%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo%2C%20pd%2C%20session_loader)%3A%0A%20%20%20%20import%20io%20as%20_restore_io%0A%0A%20%20%20%20_loaded%20%3D%20session_loader.value.get(%22loaded%22%2C%20%7B%7D)%0A%0A%20%20%20%20restore_mode%20%3D%20bool(%0A%20%20%20%20%20%20%20%20_loaded.get(%22avg_grid.csv%22)%0A%20%20%20%20%20%20%20%20and%20_loaded.get(%22roi_points.csv%22)%0A%20%20%20%20%20%20%20%20and%20_loaded.get(%22session_params.csv%22)%0A%20%20%20%20)%0A%0A%20%20%20%20if%20restore_mode%3A%0A%20%20%20%20%20%20%20%20_sp%20%3D%20pd.read_csv(_restore_io.StringIO(_loaded%5B%22session_params.csv%22%5D))%0A%20%20%20%20%20%20%20%20session_params%20%20%20%20%3D%20%7Brow%5B%22param%22%5D%3A%20float(row%5B%22value%22%5D)%20for%20_%2C%20row%20in%20_sp.iterrows()%7D%0A%20%20%20%20%20%20%20%20avg_grid_content%20%20%3D%20_loaded%5B%22avg_grid.csv%22%5D%0A%20%20%20%20%20%20%20%20roi_points_content%20%3D%20_loaded%5B%22roi_points.csv%22%5D%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20session_params%20%20%20%20%20%3D%20%7B%7D%0A%20%20%20%20%20%20%20%20avg_grid_content%20%20%20%3D%20%22%22%0A%20%20%20%20%20%20%20%20roi_points_content%20%3D%20%22%22%0A%0A%20%20%20%20_any_loaded%20%3D%20any(v%20is%20not%20None%20for%20v%20in%20_loaded.values())%20if%20_loaded%20else%20False%0A%20%20%20%20_missing%20%3D%20%5Bf%20for%20f%20in%20%5B%22avg_grid.csv%22%2C%20%22roi_points.csv%22%2C%20%22session_params.csv%22%5D%20if%20not%20_loaded.get(f)%5D%0A%20%20%20%20(%0A%20%20%20%20%20%20%20%20mo.callout(mo.md(f%22**Missing%20session%20files%3A**%20%7B'%2C%20'.join(_missing)%7D%22)%2C%20kind%3D%22danger%22)%0A%20%20%20%20%20%20%20%20if%20(_any_loaded%20and%20_missing)%0A%20%20%20%20%20%20%20%20else%20mo.md(%22%22)%0A%20%20%20%20)%0A%20%20%20%20return%20avg_grid_content%2C%20restore_mode%2C%20roi_points_content%2C%20session_params%0A%0A%0A%40app.cell%0Adef%20_(file_upload%2C%20pd)%3A%0A%20%20%20%20def%20load_first_dataset(content%3A%20str)%20-%3E%20pd.DataFrame%3A%0A%20%20%20%20%20%20%20%20%22%22%22Read%20the%20first%20dataset%20from%20a%20Keyence%20CSV%2C%20stopping%20at%20the%20first%20blank%20row.%22%22%22%0A%20%20%20%20%20%20%20%20rows%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20header%20%3D%20None%0A%20%20%20%20%20%20%20%20for%20line%20in%20content.splitlines()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20stripped%20%3D%20line.strip()%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20not%20stripped%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20rows%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20break%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20continue%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20header%20is%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20header%20%3D%20stripped.split(%22%2C%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20continue%0A%20%20%20%20%20%20%20%20%20%20%20%20rows.append(stripped.split(%22%2C%22))%0A%20%20%20%20%20%20%20%20df%20%3D%20pd.DataFrame(rows%2C%20columns%3Dheader)%0A%20%20%20%20%20%20%20%20df%5B%22measurement%22%5D%20%3D%20pd.to_numeric(df%5B%22measurement%22%5D%2C%20errors%3D%22coerce%22)%0A%20%20%20%20%20%20%20%20df%5B%22x%22%5D%20%3D%20pd.to_numeric(df%5B%22x%22%5D%2C%20errors%3D%22coerce%22)%0A%20%20%20%20%20%20%20%20df%5B%22y%22%5D%20%3D%20pd.to_numeric(df%5B%22y%22%5D%2C%20errors%3D%22coerce%22)%0A%20%20%20%20%20%20%20%20return%20df%0A%0A%20%20%20%20if%20file_upload.value%3A%0A%20%20%20%20%20%20%20%20df_raw%20%3D%20load_first_dataset(file_upload.contents().decode(%22utf-8%22))%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20df_raw%20%3D%20pd.DataFrame(columns%3D%5B%22time%22%2C%20%22measurement%22%2C%20%22x%22%2C%20%22y%22%5D)%0A%0A%20%20%20%20df_raw%0A%20%20%20%20return%20(df_raw%2C)%0A%0A%0A%40app.cell%0Adef%20_(df_raw%2C%20mo)%3A%0A%20%20%20%20mo.md(f%22%22%22%0A%20%20%20%20%23%23%20Dataset%20summary%0A%0A%20%20%20%20-%20**Rows%20loaded%3A**%20%7Blen(df_raw)%3A%2C%7D%0A%20%20%20%20-%20**x%20range%3A**%20%7Bdf_raw%5B'x'%5D.min()%3A.3f%7D%20%E2%80%93%20%7Bdf_raw%5B'x'%5D.max()%3A.3f%7D%20mm%0A%20%20%20%20-%20**y%20range%3A**%20%7Bdf_raw%5B'y'%5D.min()%3A.3f%7D%20%E2%80%93%20%7Bdf_raw%5B'y'%5D.max()%3A.3f%7D%20mm%0A%20%20%20%20-%20**Measurement%20range%3A**%20%7Bdf_raw%5B'measurement'%5D.min()%3A.4f%7D%20%E2%80%93%20%7Bdf_raw%5B'measurement'%5D.max()%3A.4f%7D%20mm%0A%20%20%20%20-%20**Invalid%20readings%20(%E2%88%9299.9999)%3A**%20%7B(df_raw%5B'measurement'%5D%20%3D%3D%20-99.9999).sum()%3A%2C%7D%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(df_raw%2C%20mo)%3A%0A%20%20%20%20%23%20Filter%20out%20invalid%20readings%20for%20valid-data%20summary%0A%20%20%20%20INVALID_VALUE%20%3D%20-99.9999%0A%20%20%20%20df_valid%20%3D%20df_raw%5Bdf_raw%5B%22measurement%22%5D%20!%3D%20INVALID_VALUE%5D.copy()%0A%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20f%22%22%22%0A%20%20%20%20%20%20%20%20%23%23%23%20Valid%20measurements%20only%0A%0A%20%20%20%20%20%20%20%20-%20**Valid%20rows%3A**%20%7Blen(df_valid)%3A%2C%7D%20of%20%7Blen(df_raw)%3A%2C%7D%20(%7B100*len(df_valid)%2Flen(df_raw)%3A.1f%7D%25)%0A%20%20%20%20%20%20%20%20-%20**Measurement%20range%3A**%20%7Bdf_valid%5B'measurement'%5D.min()%3A.4f%7D%20%E2%80%93%20%7Bdf_valid%5B'measurement'%5D.max()%3A.4f%7D%20mm%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20)%20if%20not%20df_valid.empty%20else%20mo.md(%22No%20valid%20measurements%20found.%22)%0A%20%20%20%20return%20(df_valid%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(avg_grid_content%2C%20df_valid%2C%20mo%2C%20pd%2C%20restore_mode)%3A%0A%20%20%20%20mo.stop(%0A%20%20%20%20%20%20%20%20df_valid.empty%20and%20not%20restore_mode%2C%0A%20%20%20%20%20%20%20%20mo.callout(mo.md(%22**Upload%20a%20Keyence%20CSV%20file%20above%2C%20or%20restore%20a%20previous%20session.**%22)%2C%20kind%3D%22info%22)%0A%20%20%20%20)%0A%0A%20%20%20%20if%20restore_mode%3A%0A%20%20%20%20%20%20%20%20import%20io%20as%20_avg_io%0A%20%20%20%20%20%20%20%20df_avg%20%3D%20pd.read_csv(_avg_io.StringIO(avg_grid_content))%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20_MERGE_TOL%20%3D%200.002%0A%20%20%20%20%20%20%20%20_y_sorted%20%3D%20sorted(df_valid%5B%22y%22%5D.unique())%0A%20%20%20%20%20%20%20%20_y_map%2C%20_cluster%20%3D%20%7B%7D%2C%20%5B_y_sorted%5B0%5D%5D%0A%20%20%20%20%20%20%20%20for%20_v%20in%20_y_sorted%5B1%3A%5D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_v%20-%20_cluster%5B-1%5D%20%3C%3D%20_MERGE_TOL%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_cluster.append(_v)%0A%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_rep%20%3D%20round(sum(_cluster)%20%2F%20len(_cluster)%2C%206)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20for%20_y%20in%20_cluster%3A%20_y_map%5B_y%5D%20%3D%20_rep%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_cluster%20%3D%20%5B_v%5D%0A%20%20%20%20%20%20%20%20_rep%20%3D%20round(sum(_cluster)%20%2F%20len(_cluster)%2C%206)%0A%20%20%20%20%20%20%20%20for%20_y%20in%20_cluster%3A%20_y_map%5B_y%5D%20%3D%20_rep%0A%0A%20%20%20%20%20%20%20%20df_avg%20%3D%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20df_valid%0A%20%20%20%20%20%20%20%20%20%20%20%20.assign(y_bin%3Ddf_valid%5B%22y%22%5D.map(_y_map))%0A%20%20%20%20%20%20%20%20%20%20%20%20.groupby(%5B%22x%22%2C%20%22y_bin%22%5D%2C%20as_index%3DFalse)%5B%22measurement%22%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20.mean()%0A%20%20%20%20%20%20%20%20%20%20%20%20.rename(columns%3D%7B%22y_bin%22%3A%20%22y%22%2C%20%22measurement%22%3A%20%22avg_height_mm%22%7D)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20df_avg%0A%20%20%20%20return%20(df_avg%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Step%202%20%E2%80%94%20Data%20Overview%0A%0A%20%20%20%20The%20table%20below%20shows%20any%20gaps%20in%20the%20scan%20grid%20%E2%80%94%20x%20or%20y%20positions%20where%20the%0A%20%20%20%20sensor%20skipped%20one%20or%20more%20steps.%20These%20appear%20as%20stretched%20pixels%20in%20the%20full%20map.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(df_avg%2C%20mo)%3A%0A%20%20%20%20def%20dedupe_positions(values%2C%20merge_tol%3D0.002)%3A%0A%20%20%20%20%20%20%20%20%22%22%22Merge%20positions%20within%20merge_tol%20into%20a%20single%20representative%20value.%22%22%22%0A%20%20%20%20%20%20%20%20vals%20%3D%20sorted(values)%0A%20%20%20%20%20%20%20%20clusters%2C%20group%20%3D%20%5B%5D%2C%20%5Bvals%5B0%5D%5D%0A%20%20%20%20%20%20%20%20for%20v%20in%20vals%5B1%3A%5D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20v%20-%20group%5B-1%5D%20%3C%3D%20merge_tol%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20group.append(v)%0A%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20clusters.append(sum(group)%20%2F%20len(group))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20group%20%3D%20%5Bv%5D%0A%20%20%20%20%20%20%20%20clusters.append(sum(group)%20%2F%20len(group))%0A%20%20%20%20%20%20%20%20return%20clusters%0A%0A%20%20%20%20def%20find_gaps(values%2C%20label%2C%20merge_tol%3D0.002)%3A%0A%20%20%20%20%20%20%20%20positions%20%3D%20dedupe_positions(values%2C%20merge_tol)%0A%20%20%20%20%20%20%20%20steps%20%3D%20%5Bround(b%20-%20a%2C%206)%20for%20a%2C%20b%20in%20zip(positions%2C%20positions%5B1%3A%5D)%5D%0A%20%20%20%20%20%20%20%20nominal_step%20%3D%20min(steps)%0A%20%20%20%20%20%20%20%20rows%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20a%2C%20b%2C%20step%20in%20zip(positions%2C%20positions%5B1%3A%5D%2C%20steps)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20step%20%3E%3D%20nominal_step%20*%202%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20missing%20%3D%20round((step%20-%20nominal_step)%20%2F%20nominal_step)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20rows.append(%7B%22axis%22%3A%20label%2C%20%22from%20(mm)%22%3A%20round(a%2C%204)%2C%20%22to%20(mm)%22%3A%20round(b%2C%204)%2C%20%22gap%20(mm)%22%3A%20round(step%2C%204)%2C%20%22missing%20steps%22%3A%20int(missing)%7D)%0A%20%20%20%20%20%20%20%20return%20rows%2C%20nominal_step%0A%0A%20%20%20%20x_gaps%2C%20x_step%20%3D%20find_gaps(df_avg%5B%22x%22%5D.unique()%2C%20%22x%22)%0A%20%20%20%20y_gaps%2C%20y_step%20%3D%20find_gaps(df_avg%5B%22y%22%5D.unique()%2C%20%22y%22)%0A%20%20%20%20all_gaps%20%3D%20x_gaps%20%2B%20y_gaps%0A%0A%20%20%20%20mo.vstack(%5B%0A%20%20%20%20%20%20%20%20mo.md(f%22**Nominal%20grid%20spacing**%20%E2%80%94%20x%3A%20%7Bx_step%3A.4f%7D%20mm%2C%20y%3A%20%7By_step%3A.4f%7D%20mm%22)%2C%0A%20%20%20%20%20%20%20%20mo.md(%22**Gaps%20in%20scan%20grid**%20(regions%20with%20no%20measurements)%3A%22)%2C%0A%20%20%20%20%20%20%20%20mo.ui.table(all_gaps)%20if%20all_gaps%20else%20mo.md(%22_No%20gaps%20found._%22)%2C%0A%20%20%20%20%5D)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Step%203%20%E2%80%94%20Filter%20and%20Define%20Region%20of%20Interest%0A%0A%20%20%20%20**Height%20threshold%3A**%20measurements%20above%20this%20value%20are%20excluded%20(e.g.%2C%20set%20to%20a%20little%20larger%20than%20the%20maximum%0A%20%20%20%20channel%20depth%20to%20hide%20the%20surrounding%20substrate).%0A%0A%20%20%20%20**ROI%20bounds%3A**%20the%20four%20inputs%20below%20define%20the%20x%2Fy%20window%20shown%20in%20the%20right-hand%0A%20%20%20%20plot.%20Adjust%20the%20bounds%20so%20the%20right%20plot%20frames%20the%20channel%20region%20of%20interest.%0A%20%20%20%20The%20left%20plot%20shows%20the%20full%20map%20with%20the%20ROI%20indicated%20by%20black%20lines.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo%2C%20session_params)%3A%0A%20%20%20%20height_threshold%20%3D%20mo.ui.number(%0A%20%20%20%20%20%20%20%20start%3D0.0%2C%20stop%3D2.0%2C%20step%3D0.01%2C%0A%20%20%20%20%20%20%20%20value%3Dsession_params.get(%22height_threshold%22%2C%200.5)%2C%0A%20%20%20%20%20%20%20%20label%3D%22Max%20height%20threshold%20(mm)%22%2C%0A%20%20%20%20)%0A%20%20%20%20height_threshold%0A%20%20%20%20return%20(height_threshold%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(df_avg%2C%20mo%2C%20session_params)%3A%0A%20%20%20%20if%20session_params%3A%0A%20%20%20%20%20%20%20%20_xmin_def%20%3D%20session_params%5B%22x_min%22%5D%0A%20%20%20%20%20%20%20%20_xmax_def%20%3D%20session_params%5B%22x_max%22%5D%0A%20%20%20%20%20%20%20%20_ymin_def%20%3D%20session_params%5B%22y_min%22%5D%0A%20%20%20%20%20%20%20%20_ymax_def%20%3D%20session_params%5B%22y_max%22%5D%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20_x0%2C%20_x1%20%3D%20df_avg%5B%22x%22%5D.min()%2C%20df_avg%5B%22x%22%5D.max()%0A%20%20%20%20%20%20%20%20_y0%2C%20_y1%20%3D%20df_avg%5B%22y%22%5D.min()%2C%20df_avg%5B%22y%22%5D.max()%0A%20%20%20%20%20%20%20%20_xr%2C%20_yr%20%3D%20_x1%20-%20_x0%2C%20_y1%20-%20_y0%0A%20%20%20%20%20%20%20%20_xmin_def%20%3D%20round(_x0%20%2B%20_xr%20%2F%203%2C%204)%0A%20%20%20%20%20%20%20%20_xmax_def%20%3D%20round(_x0%20%2B%202%20*%20_xr%20%2F%203%2C%204)%0A%20%20%20%20%20%20%20%20_ymin_def%20%3D%20round(_y0%2C%204)%0A%20%20%20%20%20%20%20%20_ymax_def%20%3D%20round(_y0%20%2B%202%20*%20_yr%20%2F%203%2C%204)%0A%0A%20%20%20%20x_min%20%3D%20mo.ui.number(start%3D-10.0%2C%20stop%3D10.0%2C%20step%3D0.001%2C%20value%3D_xmin_def%2C%20label%3D%22x%20min%20(mm)%22)%0A%20%20%20%20x_max%20%3D%20mo.ui.number(start%3D-10.0%2C%20stop%3D10.0%2C%20step%3D0.001%2C%20value%3D_xmax_def%2C%20label%3D%22x%20max%20(mm)%22)%0A%20%20%20%20y_min%20%3D%20mo.ui.number(start%3D-10.0%2C%20stop%3D10.0%2C%20step%3D0.001%2C%20value%3D_ymin_def%2C%20label%3D%22y%20min%20(mm)%22)%0A%20%20%20%20y_max%20%3D%20mo.ui.number(start%3D-10.0%2C%20stop%3D10.0%2C%20step%3D0.001%2C%20value%3D_ymax_def%2C%20label%3D%22y%20max%20(mm)%22)%0A%20%20%20%20mo.vstack(%5B%0A%20%20%20%20%20%20%20%20mo.hstack(%5Bx_min%2C%20x_max%5D)%2C%0A%20%20%20%20%20%20%20%20mo.hstack(%5By_min%2C%20y_max%5D)%2C%0A%20%20%20%20%5D)%0A%20%20%20%20return%20x_max%2C%20x_min%2C%20y_max%2C%20y_min%0A%0A%0A%40app.cell%0Adef%20_(df_avg%2C%20height_threshold%2C%20plt%2C%20x_max%2C%20x_min%2C%20y_max%2C%20y_min)%3A%0A%20%20%20%20df_filtered%20%3D%20df_avg%5Bdf_avg%5B%22avg_height_mm%22%5D%20%3C%20height_threshold.value%5D%0A%20%20%20%20pivot_full%20%3D%20df_filtered.pivot_table(index%3D%22y%22%2C%20columns%3D%22x%22%2C%20values%3D%22avg_height_mm%22%2C%20aggfunc%3D%22mean%22)%0A%0A%20%20%20%20%23%20Expand%20ROI%20selection%20by%20one%20grid%20point%20beyond%20each%20boundary%20so%20pcolormesh%0A%20%20%20%20%23%20computes%20correct%20cell%20edges%20at%20the%20ROI%20limits%3B%20xlim%2Fylim%20clips%20the%20display.%0A%20%20%20%20all_x%20%3D%20sorted(df_filtered%5B%22x%22%5D.unique())%0A%20%20%20%20all_y%20%3D%20sorted(df_filtered%5B%22y%22%5D.unique())%0A%20%20%20%20x_lo%20%3D%20max((x%20for%20x%20in%20all_x%20if%20x%20%3C%20x_min.value)%2C%20default%3Dx_min.value)%0A%20%20%20%20x_hi%20%3D%20min((x%20for%20x%20in%20all_x%20if%20x%20%3E%20x_max.value)%2C%20default%3Dx_max.value)%0A%20%20%20%20y_lo%20%3D%20max((y%20for%20y%20in%20all_y%20if%20y%20%3C%20y_min.value)%2C%20default%3Dy_min.value)%0A%20%20%20%20y_hi%20%3D%20min((y%20for%20y%20in%20all_y%20if%20y%20%3E%20y_max.value)%2C%20default%3Dy_max.value)%0A%0A%20%20%20%20df_roi_expanded%20%3D%20df_filtered%5B%0A%20%20%20%20%20%20%20%20(df_filtered%5B%22x%22%5D%20%3E%3D%20x_lo)%20%26%20(df_filtered%5B%22x%22%5D%20%3C%3D%20x_hi)%20%26%0A%20%20%20%20%20%20%20%20(df_filtered%5B%22y%22%5D%20%3E%3D%20y_lo)%20%26%20(df_filtered%5B%22y%22%5D%20%3C%3D%20y_hi)%0A%20%20%20%20%5D%0A%20%20%20%20pivot_roi%20%3D%20df_roi_expanded.pivot_table(index%3D%22y%22%2C%20columns%3D%22x%22%2C%20values%3D%22avg_height_mm%22%2C%20aggfunc%3D%22mean%22)%0A%20%20%20%20n_roi%20%3D%20df_filtered%5B%0A%20%20%20%20%20%20%20%20(df_filtered%5B%22x%22%5D%20%3E%3D%20x_min.value)%20%26%20(df_filtered%5B%22x%22%5D%20%3C%3D%20x_max.value)%20%26%0A%20%20%20%20%20%20%20%20(df_filtered%5B%22y%22%5D%20%3E%3D%20y_min.value)%20%26%20(df_filtered%5B%22y%22%5D%20%3C%3D%20y_max.value)%0A%20%20%20%20%5D.shape%5B0%5D%0A%0A%20%20%20%20vmax%20%3D%20height_threshold.value%0A%20%20%20%20fig%2C%20axes%20%3D%20plt.subplots(1%2C%202%2C%20figsize%3D(6%2C%205)%2C%20constrained_layout%3DTrue)%0A%0A%20%20%20%20mesh%20%3D%20axes%5B0%5D.pcolormesh(%0A%20%20%20%20%20%20%20%20pivot_full.columns%2C%20pivot_full.index%2C%20pivot_full.values%2C%0A%20%20%20%20%20%20%20%20cmap%3D%22viridis%22%2C%20shading%3D%22auto%22%2C%20vmin%3D0%2C%20vmax%3Dvmax%2C%0A%20%20%20%20)%0A%20%20%20%20axes%5B0%5D.axvline(x_min.value%2C%20color%3D%22black%22%2C%20linewidth%3D1.2)%0A%20%20%20%20axes%5B0%5D.axvline(x_max.value%2C%20color%3D%22black%22%2C%20linewidth%3D1.2)%0A%20%20%20%20axes%5B0%5D.axhline(y_min.value%2C%20color%3D%22black%22%2C%20linewidth%3D1.2)%0A%20%20%20%20axes%5B0%5D.axhline(y_max.value%2C%20color%3D%22black%22%2C%20linewidth%3D1.2)%0A%20%20%20%20axes%5B0%5D.set_xlabel(%22x%20(mm)%22)%0A%20%20%20%20axes%5B0%5D.set_ylabel(%22y%20(mm)%22)%0A%20%20%20%20axes%5B0%5D.set_title(f%22Full%20map%20%20(n%3D%7Blen(df_filtered)%3A%2C%7D)%22)%0A%20%20%20%20axes%5B0%5D.set_aspect(%22equal%22)%0A%0A%20%20%20%20if%20not%20pivot_roi.empty%3A%0A%20%20%20%20%20%20%20%20axes%5B1%5D.pcolormesh(%0A%20%20%20%20%20%20%20%20%20%20%20%20pivot_roi.columns%2C%20pivot_roi.index%2C%20pivot_roi.values%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20cmap%3D%22viridis%22%2C%20shading%3D%22auto%22%2C%20vmin%3D0%2C%20vmax%3Dvmax%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20axes%5B1%5D.set_xlim(x_min.value%2C%20x_max.value)%0A%20%20%20%20axes%5B1%5D.set_ylim(y_min.value%2C%20y_max.value)%0A%20%20%20%20axes%5B1%5D.set_xlabel(%22x%20(mm)%22)%0A%20%20%20%20axes%5B1%5D.set_ylabel(%22y%20(mm)%22)%0A%20%20%20%20axes%5B1%5D.set_title(f%22ROI%20%20(n%3D%7Bn_roi%3A%2C%7D)%22)%0A%20%20%20%20axes%5B1%5D.set_aspect(%22equal%22)%0A%0A%20%20%20%20fig.colorbar(mesh%2C%20ax%3Daxes%2C%20label%3D%22Avg%20height%20(mm)%22%2C%20shrink%3D0.8)%0A%20%20%20%20fig%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Step%204%20%E2%80%94%20Section%20Analysis%0A%0A%20%20%20%20Paint%20cells%20in%20the%20grid%20below%20to%20assign%20them%20to%20one%20of%20six%20sections%20(S1%E2%80%93S6).%0A%20%20%20%20Click%20a%20button%20at%20the%20top%20to%20choose%20the%20active%20section%20colour%2C%20then%20click%20or%20drag%0A%20%20%20%20over%20grid%20cells%20to%20assign%20them.%20Clicking%20an%20already-assigned%20cell%20removes%20it.%0A%20%20%20%20The%203D%20scatter%20updates%20live%2C%20and%20the%20statistics%20table%20below%20shows%20mean%20and%20SD%0A%20%20%20%20height%20per%20section%20in%20%C2%B5m.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(df_avg%2C%20height_threshold%2C%20np%2C%20x_max%2C%20x_min%2C%20y_max%2C%20y_min)%3A%0A%20%20%20%20%23%20ROI%20grid%20data%20for%20the%20section-selector%20widget%0A%20%20%20%20df_roi_sel%20%3D%20df_avg%5B%0A%20%20%20%20%20%20%20%20(df_avg%5B%22avg_height_mm%22%5D%20%3C%20height_threshold.value)%20%26%0A%20%20%20%20%20%20%20%20(df_avg%5B%22x%22%5D%20%3E%3D%20x_min.value)%20%26%20(df_avg%5B%22x%22%5D%20%3C%3D%20x_max.value)%20%26%0A%20%20%20%20%20%20%20%20(df_avg%5B%22y%22%5D%20%3E%3D%20y_min.value)%20%26%20(df_avg%5B%22y%22%5D%20%3C%3D%20y_max.value)%0A%20%20%20%20%5D%0A%20%20%20%20pivot_sel%20%3D%20df_roi_sel.pivot_table(index%3D%22y%22%2C%20columns%3D%22x%22%2C%20values%3D%22avg_height_mm%22%2C%20aggfunc%3D%22mean%22)%0A%20%20%20%20xs_sel%20%3D%20pivot_sel.columns.values%0A%20%20%20%20ys_sel%20%3D%20pivot_sel.index.values%0A%20%20%20%20Z_sel%20%20%3D%20pivot_sel.values%0A%0A%20%20%20%20XX%2C%20YY%20%3D%20np.meshgrid(xs_sel%2C%20ys_sel)%0A%20%20%20%20_mask%20%3D%20~np.isnan(Z_sel)%0A%20%20%20%20flat_x_sel%20%20%3D%20XX%5B_mask%5D%0A%20%20%20%20flat_y_sel%20%20%3D%20YY%5B_mask%5D%0A%20%20%20%20flat_z_sel%20%20%3D%20Z_sel%5B_mask%5D%20*%201000%20%20%23%20%C2%B5m%0A%0A%20%20%20%20_ny%2C%20_nx%20%3D%20Z_sel.shape%0A%20%20%20%20_c2i%20%3D%20np.full(_ny%20*%20_nx%2C%20-1%2C%20dtype%3Dint)%0A%20%20%20%20_c2i%5B_mask.flatten()%5D%20%3D%20np.arange(int(_mask.sum()))%0A%20%20%20%20x_vals_sel%20%20%20%3D%20xs_sel.tolist()%0A%20%20%20%20y_vals_sel%20%20%20%3D%20ys_sel.tolist()%0A%20%20%20%20cell_to_idx_sel%20%3D%20_c2i.tolist()%0A%0A%20%20%20%20print(f%22ROI%20grid%3A%20%7BZ_sel.shape%7D%20%20valid%20pts%3A%20%7Bint(_mask.sum())%7D%20%20%22%0A%20%20%20%20%20%20%20%20%20%20f%22NaN%3A%20%7Bnp.isnan(Z_sel).mean()%3A.1%25%7D%22)%0A%20%20%20%20print(f%22Height%20range%3A%20%5B%7Bnp.nanmin(Z_sel)*1000%3A.1f%7D%2C%20%7Bnp.nanmax(Z_sel)*1000%3A.1f%7D%5D%20%C2%B5m%22)%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20cell_to_idx_sel%2C%0A%20%20%20%20%20%20%20%20flat_x_sel%2C%0A%20%20%20%20%20%20%20%20flat_y_sel%2C%0A%20%20%20%20%20%20%20%20flat_z_sel%2C%0A%20%20%20%20%20%20%20%20x_vals_sel%2C%0A%20%20%20%20%20%20%20%20y_vals_sel%2C%0A%20%20%20%20)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(anywidget%2C%20traitlets)%3A%0A%20%20%20%20class%20ChannelSectionSelector(anywidget.AnyWidget)%3A%0A%20%20%20%20%20%20%20%20x_vals%20%20%20%20%20%20%3D%20traitlets.List().tag(sync%3DTrue)%0A%20%20%20%20%20%20%20%20y_vals%20%20%20%20%20%20%3D%20traitlets.List().tag(sync%3DTrue)%0A%20%20%20%20%20%20%20%20cell_to_idx%20%3D%20traitlets.List().tag(sync%3DTrue)%0A%20%20%20%20%20%20%20%20x_pts%20%20%20%20%20%20%20%3D%20traitlets.List().tag(sync%3DTrue)%0A%20%20%20%20%20%20%20%20y_pts%20%20%20%20%20%20%20%3D%20traitlets.List().tag(sync%3DTrue)%0A%20%20%20%20%20%20%20%20z_pts%20%20%20%20%20%20%20%3D%20traitlets.List().tag(sync%3DTrue)%0A%20%20%20%20%20%20%20%20value%20%20%20%20%20%20%20%3D%20traitlets.List().tag(sync%3DTrue)%0A%20%20%20%20%20%20%20%20sec_colors%20%20%3D%20traitlets.List().tag(sync%3DTrue)%0A%0A%20%20%20%20%20%20%20%20_esm%20%3D%20'%5Cnfunction%20ensurePlotly()%20%7B%5Cn%20%20return%20new%20Promise(resolve%20%3D%3E%20%7B%5Cn%20%20%20%20if%20(window.Plotly)%20%7B%20resolve()%3B%20return%3B%20%7D%5Cn%20%20%20%20const%20s%20%3D%20document.createElement(%22script%22)%3B%5Cn%20%20%20%20s.src%20%3D%20%22https%3A%2F%2Fcdn.plot.ly%2Fplotly-2.35.2.min.js%22%3B%5Cn%20%20%20%20s.onload%20%3D%20resolve%3B%5Cn%20%20%20%20document.head.appendChild(s)%3B%5Cn%20%20%7D)%3B%5Cn%7D%5Cn%5Cnasync%20function%20render(%7B%20model%2C%20el%20%7D)%20%7B%5Cn%20%20await%20ensurePlotly()%3B%5Cn%20%20const%20Plotly%20%3D%20window.Plotly%3B%5Cn%5Cn%20%20el.style.cssText%20%3D%20%22display%3Aflex%3Bflex-direction%3Acolumn%3Bgap%3A6px%3B%22%3B%5Cn%5Cn%20%20const%20bar%20%3D%20document.createElement(%22div%22)%3B%5Cn%20%20bar.style.cssText%20%3D%20%22display%3Aflex%3Bgap%3A6px%3Balign-items%3Acenter%3Bpadding%3A4px%200%3Bflex-wrap%3Awrap%3B%22%3B%5Cn%20%20const%20lbl%20%3D%20document.createElement(%22b%22)%3B%20lbl.textContent%20%3D%20%22Active%20section%3A%22%3B%5Cn%20%20lbl.style.marginRight%20%3D%20%224px%22%3B%20bar.appendChild(lbl)%3B%5Cn%20%20const%20sc%20%3D%20model.get(%22sec_colors%22)%3B%5Cn%20%20let%20activeSec%20%3D%200%3B%5Cn%20%20const%20btns%20%3D%20sc.map((c%2C%20i)%20%3D%3E%20%7B%5Cn%20%20%20%20const%20b%20%3D%20document.createElement(%22button%22)%3B%5Cn%20%20%20%20b.textContent%20%3D%20%22S%22%20%2B%20(i%20%2B%201)%3B%5Cn%20%20%20%20b.style.cssText%20%3D%20%22background%3A%22%20%2B%20c%20%2B%20%22%3Bcolor%3Awhite%3Bborder%3A3px%20solid%20transparent%3B%22%20%2B%5Cn%20%20%20%20%20%20%22border-radius%3A4px%3Bpadding%3A3px%2010px%3Bcursor%3Apointer%3Bfont-weight%3Abold%3Bfont-size%3A14px%3B%22%3B%5Cn%20%20%20%20b.onclick%20%3D%20()%20%3D%3E%20%7B%5Cn%20%20%20%20%20%20activeSec%20%3D%20i%3B%5Cn%20%20%20%20%20%20btns.forEach((bb%2C%20ii)%20%3D%3E%20%7B%20bb.style.borderColor%20%3D%20ii%20%3D%3D%3D%20i%20%3F%20%22%23000%22%20%3A%20%22transparent%22%3B%20%7D)%3B%5Cn%20%20%20%20%7D%3B%5Cn%20%20%20%20bar.appendChild(b)%3B%5Cn%20%20%20%20return%20b%3B%5Cn%20%20%7D)%3B%5Cn%20%20btns%5B0%5D.style.borderColor%20%3D%20%22%23000%22%3B%5Cn%20%20el.appendChild(bar)%3B%5Cn%5Cn%20%20const%20hint%20%3D%20document.createElement(%22div%22)%3B%5Cn%20%20hint.style.cssText%20%3D%20%22font-size%3A12px%3Bcolor%3A%23666%3Bpadding%3A2px%204px%3B%22%3B%5Cn%20%20hint.textContent%20%3D%20%22Click%20or%20drag%20to%20assign%20cells%20to%20the%20active%20section.%20Clicking%20an%20assigned%20cell%20removes%20it.%22%3B%5Cn%20%20el.appendChild(hint)%3B%5Cn%5Cn%20%20const%20row%20%3D%20document.createElement(%22div%22)%3B%5Cn%20%20row.style.cssText%20%3D%20%22display%3Aflex%3Bgap%3A8px%3Balign-items%3Aflex-start%3B%22%3B%5Cn%20%20el.appendChild(row)%3B%5Cn%5Cn%20%20const%20wrapper%20%3D%20document.createElement(%22div%22)%3B%5Cn%20%20wrapper.style.cssText%20%3D%20%22overflow%3Aauto%3Bmax-height%3A580px%3Bborder%3A1px%20solid%20%23555%3B%22%20%2B%5Cn%20%20%20%20%22border-radius%3A4px%3Bflex-shrink%3A0%3Bcursor%3Acrosshair%3B%22%3B%5Cn%20%20row.appendChild(wrapper)%3B%5Cn%20%20const%20canvas%20%3D%20document.createElement(%22canvas%22)%3B%5Cn%20%20canvas.style.display%20%3D%20%22block%22%3B%5Cn%20%20wrapper.appendChild(canvas)%3B%5Cn%5Cn%20%20const%20div3d%20%3D%20document.createElement(%22div%22)%3B%5Cn%20%20div3d.style.cssText%20%3D%20%22flex%3A1%3Bheight%3A580px%3Bmin-width%3A0%3B%22%3B%5Cn%20%20row.appendChild(div3d)%3B%5Cn%5Cn%20%20const%20xv%20%3D%20model.get(%22x_vals%22)%2C%20yv%20%3D%20model.get(%22y_vals%22)%3B%5Cn%20%20const%20cidx%20%3D%20model.get(%22cell_to_idx%22)%3B%5Cn%20%20const%20nx%20%3D%20xv.length%2C%20ny%20%3D%20yv.length%3B%5Cn%20%20let%20colors%20%3D%20model.get(%22value%22).slice()%3B%5Cn%5Cn%20%20const%20UNSET%20%3D%20%22%23aaaaaa%22%3B%5Cn%20%20const%20CW%20%3D%2044%2C%20CH%20%3D%209%3B%5Cn%20%20const%20PL%20%3D%2052%2C%20PR%20%3D%208%2C%20PT%20%3D%208%2C%20PB%20%3D%2038%3B%5Cn%20%20canvas.width%20%20%3D%20PL%20%2B%20nx%20*%20CW%20%2B%20PR%3B%5Cn%20%20canvas.height%20%3D%20PT%20%2B%20ny%20*%20CH%20%2B%20PB%3B%5Cn%20%20const%20ctx2d%20%3D%20canvas.getContext(%222d%22)%3B%5Cn%5Cn%20%20function%20cellTop(yi)%20%7B%20return%20PT%20%2B%20(ny%20-%201%20-%20yi)%20*%20CH%3B%20%7D%5Cn%5Cn%20%20function%20drawFullGrid()%20%7B%5Cn%20%20%20%20ctx2d.fillStyle%20%3D%20%22%23111%22%3B%5Cn%20%20%20%20ctx2d.fillRect(0%2C%200%2C%20canvas.width%2C%20canvas.height)%3B%5Cn%20%20%20%20for%20(let%20yi%20%3D%200%3B%20yi%20%3C%20ny%3B%20yi%2B%2B)%20%7B%5Cn%20%20%20%20%20%20for%20(let%20xi%20%3D%200%3B%20xi%20%3C%20nx%3B%20xi%2B%2B)%20%7B%5Cn%20%20%20%20%20%20%20%20const%20k%20%3D%20cidx%5Byi%20*%20nx%20%2B%20xi%5D%3B%5Cn%20%20%20%20%20%20%20%20ctx2d.fillStyle%20%3D%20k%20%3C%200%20%3F%20%22%232a2a2a%22%20%3A%20colors%5Bk%5D%3B%5Cn%20%20%20%20%20%20%20%20ctx2d.fillRect(PL%20%2B%20xi%20*%20CW%20%2B%201%2C%20cellTop(yi)%20%2B%201%2C%20CW%20-%202%2C%20CH%20-%202)%3B%5Cn%20%20%20%20%20%20%7D%5Cn%20%20%20%20%7D%5Cn%20%20%20%20ctx2d.fillStyle%20%3D%20%22%23999%22%3B%20ctx2d.font%20%3D%20%229px%20monospace%22%3B%20ctx2d.textAlign%20%3D%20%22right%22%3B%5Cn%20%20%20%20const%20ystep%20%3D%20Math.max(1%2C%20Math.round(ny%20%2F%2028))%3B%5Cn%20%20%20%20for%20(let%20yi%20%3D%200%3B%20yi%20%3C%20ny%3B%20yi%20%2B%3D%20ystep)%5Cn%20%20%20%20%20%20ctx2d.fillText(yv%5Byi%5D.toFixed(3)%2C%20PL%20-%203%2C%20cellTop(yi)%20%2B%20CH%20%2F%202%20%2B%203)%3B%5Cn%20%20%20%20ctx2d.textAlign%20%3D%20%22center%22%3B%5Cn%20%20%20%20for%20(let%20xi%20%3D%200%3B%20xi%20%3C%20nx%3B%20xi%2B%2B)%5Cn%20%20%20%20%20%20ctx2d.fillText(xv%5Bxi%5D.toFixed(3)%2C%20PL%20%2B%20xi%20*%20CW%20%2B%20CW%20%2F%202%2C%20PT%20%2B%20ny%20*%20CH%20%2B%2016)%3B%5Cn%20%20%20%20ctx2d.fillStyle%20%3D%20%22%23777%22%3B%20ctx2d.font%20%3D%20%2210px%20sans-serif%22%3B%20ctx2d.save()%3B%5Cn%20%20%20%20ctx2d.translate(11%2C%20PT%20%2B%20(ny%20*%20CH)%20%2F%202)%3B%20ctx2d.rotate(-Math.PI%20%2F%202)%3B%5Cn%20%20%20%20ctx2d.textAlign%20%3D%20%22center%22%3B%20ctx2d.fillText(%22y%20(mm)%22%2C%200%2C%200)%3B%20ctx2d.restore()%3B%5Cn%20%20%20%20ctx2d.textAlign%20%3D%20%22center%22%3B%5Cn%20%20%20%20ctx2d.fillText(%22x%20(mm)%22%2C%20PL%20%2B%20(nx%20*%20CW)%20%2F%202%2C%20PT%20%2B%20ny%20*%20CH%20%2B%2030)%3B%5Cn%20%20%7D%5Cn%5Cn%20%20let%20painting%20%3D%20false%2C%20dirty%20%3D%20false%2C%20paintMode%20%3D%20%22paint%22%3B%5Cn%5Cn%20%20function%20applyCell(xi%2C%20yi)%20%7B%5Cn%20%20%20%20const%20k%20%3D%20cidx%5Byi%20*%20nx%20%2B%20xi%5D%3B%5Cn%20%20%20%20if%20(k%20%3C%200)%20return%20false%3B%5Cn%20%20%20%20const%20target%20%3D%20paintMode%20%3D%3D%3D%20%22erase%22%20%3F%20UNSET%20%3A%20sc%5BactiveSec%5D%3B%5Cn%20%20%20%20if%20(colors%5Bk%5D%20%3D%3D%3D%20target)%20return%20false%3B%5Cn%20%20%20%20colors%5Bk%5D%20%3D%20target%3B%5Cn%20%20%20%20ctx2d.fillStyle%20%3D%20target%3B%5Cn%20%20%20%20ctx2d.fillRect(PL%20%2B%20xi%20*%20CW%20%2B%201%2C%20cellTop(yi)%20%2B%201%2C%20CW%20-%202%2C%20CH%20-%202)%3B%5Cn%20%20%20%20return%20true%3B%5Cn%20%20%7D%5Cn%5Cn%20%20function%20hitTest(e)%20%7B%5Cn%20%20%20%20const%20r%20%3D%20canvas.getBoundingClientRect()%3B%5Cn%20%20%20%20const%20sx%20%3D%20canvas.width%20%2F%20r.width%2C%20sy%20%3D%20canvas.height%20%2F%20r.height%3B%5Cn%20%20%20%20const%20mx%20%3D%20(e.clientX%20-%20r.left)%20*%20sx%2C%20my%20%3D%20(e.clientY%20-%20r.top)%20*%20sy%3B%5Cn%20%20%20%20const%20xi%20%3D%20Math.floor((mx%20-%20PL)%20%2F%20CW)%3B%5Cn%20%20%20%20const%20yi%20%3D%20ny%20-%201%20-%20Math.floor((my%20-%20PT)%20%2F%20CH)%3B%5Cn%20%20%20%20if%20(xi%20%3C%200%20%7C%7C%20xi%20%3E%3D%20nx%20%7C%7C%20yi%20%3C%200%20%7C%7C%20yi%20%3E%3D%20ny)%20return%20null%3B%5Cn%20%20%20%20return%20%5Bxi%2C%20yi%5D%3B%5Cn%20%20%7D%5Cn%5Cn%20%20canvas.addEventListener(%22mousedown%22%2C%20e%20%3D%3E%20%7B%5Cn%20%20%20%20painting%20%3D%20true%3B%5Cn%20%20%20%20const%20h%20%3D%20hitTest(e)%3B%5Cn%20%20%20%20if%20(!h)%20return%3B%5Cn%20%20%20%20const%20k%20%3D%20cidx%5Bh%5B1%5D%20*%20nx%20%2B%20h%5B0%5D%5D%3B%5Cn%20%20%20%20paintMode%20%3D%20(k%20%3E%3D%200%20%26%26%20colors%5Bk%5D%20%3D%3D%3D%20sc%5BactiveSec%5D)%20%3F%20%22erase%22%20%3A%20%22paint%22%3B%5Cn%20%20%20%20if%20(applyCell(h%5B0%5D%2C%20h%5B1%5D))%20dirty%20%3D%20true%3B%5Cn%20%20%7D)%3B%5Cn%20%20canvas.addEventListener(%22mousemove%22%2C%20e%20%3D%3E%20%7B%5Cn%20%20%20%20if%20(!painting)%20return%3B%5Cn%20%20%20%20const%20h%20%3D%20hitTest(e)%3B%5Cn%20%20%20%20if%20(h%20%26%26%20applyCell(h%5B0%5D%2C%20h%5B1%5D))%20dirty%20%3D%20true%3B%5Cn%20%20%7D)%3B%5Cn%20%20function%20commit()%20%7B%5Cn%20%20%20%20painting%20%3D%20false%3B%5Cn%20%20%20%20if%20(!dirty)%20return%3B%5Cn%20%20%20%20dirty%20%3D%20false%3B%5Cn%20%20%20%20Plotly.restyle(div3d%2C%20%7B%22marker.color%22%3A%20%5Bcolors.slice()%5D%7D)%3B%5Cn%20%20%20%20model.set(%22value%22%2C%20colors.slice())%3B%5Cn%20%20%20%20model.save_changes()%3B%5Cn%20%20%7D%5Cn%20%20canvas.addEventListener(%22mouseup%22%2C%20%20%20%20commit)%3B%5Cn%20%20canvas.addEventListener(%22mouseleave%22%2C%20commit)%3B%5Cn%5Cn%20%20await%20Plotly.newPlot(div3d%2C%20%5B%7B%5Cn%20%20%20%20x%3A%20model.get(%22x_pts%22)%2C%20y%3A%20model.get(%22y_pts%22)%2C%20z%3A%20model.get(%22z_pts%22)%2C%5Cn%20%20%20%20type%3A%20%22scatter3d%22%2C%20mode%3A%20%22markers%22%2C%5Cn%20%20%20%20marker%3A%20%7B%20color%3A%20colors.slice()%2C%20size%3A%205%2C%20opacity%3A%200.85%20%7D%2C%5Cn%20%20%7D%5D%2C%20%7B%5Cn%20%20%20%20uirevision%3A%20%22camera-fixed%22%2C%5Cn%20%20%20%20scene%3A%20%7B%5Cn%20%20%20%20%20%20xaxis%3A%20%7B%20title%3A%20%22x%20(mm)%22%20%7D%2C%20yaxis%3A%20%7B%20title%3A%20%22y%20(mm)%22%20%7D%2C%5Cn%20%20%20%20%20%20zaxis%3A%20%7B%20title%3A%20%22Height%20(%C2%B5m)%22%20%7D%2C%5Cn%20%20%20%20%20%20aspectmode%3A%20%22manual%22%2C%20aspectratio%3A%20%7B%20x%3A%200.5%2C%20y%3A%202.5%2C%20z%3A%201%20%7D%2C%5Cn%20%20%20%20%20%20camera%3A%20%7B%5Cn%20%20%20%20%20%20%20%20eye%3A%20%20%20%20%7B%20x%3A%201.5835825282309526%2C%20%20y%3A%20-2.567896581879239%2C%20%20%20z%3A%201.0965461841144974%20%20%7D%2C%5Cn%20%20%20%20%20%20%20%20up%3A%20%20%20%20%20%7B%20x%3A%20-0.007445266559157256%2C%20y%3A%20-0.006403405783215493%2C%20z%3A%200.9999517810375851%20%7D%2C%5Cn%20%20%20%20%20%20%20%20center%3A%20%7B%20x%3A%20-0.4315610370627384%2C%20%20y%3A%20-0.3781042286667127%2C%20%20z%3A%20-0.33805537739416563%20%7D%2C%5Cn%20%20%20%20%20%20%7D%2C%5Cn%20%20%20%20%7D%2C%5Cn%20%20%20%20title%3A%20%7B%20text%3A%20%223D%20view%22%2C%20font%3A%20%7B%20size%3A%2013%20%7D%20%7D%2C%5Cn%20%20%20%20margin%3A%20%7B%20l%3A%200%2C%20r%3A%200%2C%20t%3A%2045%2C%20b%3A%200%20%7D%2C%20height%3A%20580%2C%5Cn%20%20%7D%2C%20%7B%20responsive%3A%20true%20%7D)%3B%5Cn%5Cn%20%20drawFullGrid()%3B%5Cn%5Cn%20%20model.on(%22change%3Avalue%22%2C%20()%20%3D%3E%20%7B%5Cn%20%20%20%20if%20(painting)%20return%3B%5Cn%20%20%20%20colors%20%3D%20model.get(%22value%22).slice()%3B%5Cn%20%20%20%20drawFullGrid()%3B%5Cn%20%20%20%20Plotly.restyle(div3d%2C%20%7B%22marker.color%22%3A%20%5Bcolors.slice()%5D%7D)%3B%5Cn%20%20%7D)%3B%5Cn%7D%5Cn%5Cnexport%20default%20%7B%20render%20%7D%3B%5Cn'%0A%0A%0A%20%20%20%20return%20(ChannelSectionSelector%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20ChannelSectionSelector%2C%0A%20%20%20%20cell_to_idx_sel%2C%0A%20%20%20%20flat_x_sel%2C%0A%20%20%20%20flat_y_sel%2C%0A%20%20%20%20flat_z_sel%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20pd%2C%0A%20%20%20%20restore_mode%2C%0A%20%20%20%20roi_points_content%2C%0A%20%20%20%20x_vals_sel%2C%0A%20%20%20%20y_vals_sel%2C%0A)%3A%0A%20%20%20%20_SEC_COLORS%20%3D%20%5B%22%23e41a1c%22%2C%22%23377eb8%22%2C%22%234daf4a%22%2C%22%23984ea3%22%2C%22%23ff7f00%22%2C%22%23a65628%22%5D%0A%0A%20%20%20%20if%20restore_mode%20and%20roi_points_content%3A%0A%20%20%20%20%20%20%20%20import%20io%20as%20_roi_io2%0A%20%20%20%20%20%20%20%20_sec_to_color%20%3D%20%7Bf%22S%7Bi%2B1%7D%22%3A%20c%20for%20i%2C%20c%20in%20enumerate(_SEC_COLORS)%7D%0A%20%20%20%20%20%20%20%20_df_roi%20%3D%20pd.read_csv(_roi_io2.StringIO(roi_points_content))%0A%20%20%20%20%20%20%20%20_point_color%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20(round(float(r.x_mm)%2C%206)%2C%20round(float(r.y_mm)%2C%206))%3A%20_sec_to_color.get(r.section%2C%20%22%23aaaaaa%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20_%2C%20r%20in%20_df_roi.iterrows()%0A%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20%20%20_init_colors%20%3D%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20_point_color.get((round(float(x)%2C%206)%2C%20round(float(y)%2C%206))%2C%20%22%23aaaaaa%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20x%2C%20y%20in%20zip(flat_x_sel%2C%20flat_y_sel)%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20_n_restored%20%3D%20sum(1%20for%20c%20in%20_init_colors%20if%20c%20!%3D%20%22%23aaaaaa%22)%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20_init_colors%20%3D%20%5B%22%23aaaaaa%22%5D%20*%20len(flat_x_sel)%0A%20%20%20%20%20%20%20%20_n_restored%20%3D%200%0A%0A%20%20%20%20_raw%20%3D%20ChannelSectionSelector(%0A%20%20%20%20%20%20%20%20x_vals%3Dx_vals_sel%2C%0A%20%20%20%20%20%20%20%20y_vals%3Dy_vals_sel%2C%0A%20%20%20%20%20%20%20%20cell_to_idx%3Dcell_to_idx_sel%2C%0A%20%20%20%20%20%20%20%20x_pts%3Dflat_x_sel.tolist()%2C%0A%20%20%20%20%20%20%20%20y_pts%3Dflat_y_sel.tolist()%2C%0A%20%20%20%20%20%20%20%20z_pts%3Dflat_z_sel.tolist()%2C%0A%20%20%20%20%20%20%20%20value%3D_init_colors%2C%0A%20%20%20%20%20%20%20%20sec_colors%3D_SEC_COLORS%2C%0A%20%20%20%20)%0A%20%20%20%20selector%20%3D%20mo.ui.anywidget(_raw)%0A%0A%20%20%20%20mo.vstack(%5B%0A%20%20%20%20%20%20%20%20selector%2C%0A%20%20%20%20%20%20%20%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.callout(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mo.md(f%22**Section%20assignments%20restored**%20from%20%60roi_points.csv%60%20(%7B_n_restored%7D%20assigned%20points).%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20kind%3D%22success%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_n_restored%20%3E%200%20else%20mo.md(%22%22)%0A%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%5D)%0A%20%20%20%20return%20(selector%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(flat_z_sel%2C%20mo%2C%20np%2C%20selector)%3A%0A%20%20%20%20_SEC_COLORS%20%3D%20%5B%22%23e41a1c%22%2C%22%23377eb8%22%2C%22%234daf4a%22%2C%22%23984ea3%22%2C%22%23ff7f00%22%2C%22%23a65628%22%5D%0A%20%20%20%20_pt_colors%20%3D%20selector.value%5B%22value%22%5D%0A%20%20%20%20_rows%20%3D%20%5B%22%7C%20Section%20%7C%20N%20pts%20%7C%20Mean%20(%C2%B5m)%20%7C%20SD%20(%C2%B5m)%20%7C%22%2C%20%22%7C%3A---%7C---%3A%7C---%3A%7C---%3A%7C%22%5D%0A%20%20%20%20for%20_i%2C%20_c%20in%20enumerate(_SEC_COLORS)%3A%0A%20%20%20%20%20%20%20%20_mask%20%3D%20np.array(%5Bv%20%3D%3D%20_c%20for%20v%20in%20_pt_colors%5D)%0A%20%20%20%20%20%20%20%20_z%20%3D%20flat_z_sel%5B_mask%5D%0A%20%20%20%20%20%20%20%20if%20len(_z)%20%3E%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_rows.append(f%22%7C%20S%7B_i%2B1%7D%20%7C%20%7Blen(_z)%7D%20%7C%20%7Bnp.mean(_z)%3A.1f%7D%20%7C%20%7Bnp.std(_z)%3A.1f%7D%20%7C%22)%0A%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_rows.append(f%22%7C%20S%7B_i%2B1%7D%20%7C%20%E2%80%94%20%7C%20%E2%80%94%20%7C%20%E2%80%94%20%7C%22)%0A%20%20%20%20mo.md(%22%23%23%23%20Section%20height%20estimates%20(%C2%B5m)%5Cn%5Cn%22%20%2B%20%22%5Cn%22.join(_rows))%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Step%205%20%E2%80%94%20Save%20Session%0A%0A%20%20%20%20Enter%20a%20session%20name%20and%20click%20**%F0%9F%92%BE%20Save%20session%20(.zip)**%20to%20download%20a%20ZIP%20file%0A%20%20%20%20containing%20all%20four%20session%20files.%20Unzip%20it%20to%20get%20the%20session%20folder%2C%20then%20use%0A%20%20%20%20**%F0%9F%93%82%20Open%20session%20folder**%20in%20Step%201%20to%20restore%20the%20session.%0A%0A%20%20%20%20%7C%20File%20%7C%20Contents%20%7C%0A%20%20%20%20%7C---%7C---%7C%0A%20%20%20%20%7C%20**avg_grid.csv**%20%7C%20Full%20averaged%20height%20grid%20(x%2C%20y%2C%20avg_height_mm)%20%7C%0A%20%20%20%20%7C%20**roi_points.csv**%20%7C%20ROI%20points%20with%20section%20labels%20and%20height%20in%20mm%20and%20%C2%B5m%20%7C%0A%20%20%20%20%7C%20**section_stats.csv**%20%7C%20Per-section%20statistics%20(N%2C%20mean%2C%20SD%2C%20min%2C%20max%20in%20%C2%B5m)%20%7C%0A%20%20%20%20%7C%20**session_params.csv**%20%7C%20Height%20threshold%20and%20ROI%20bounds%20for%20session%20restore%20%7C%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(anywidget%2C%20traitlets)%3A%0A%20%20%20%20class%20SessionSaver(anywidget.AnyWidget)%3A%0A%20%20%20%20%20%20%20%20files%20%3D%20traitlets.Dict(%7B%7D).tag(sync%3DTrue)%0A%20%20%20%20%20%20%20%20label%20%3D%20traitlets.Unicode('').tag(sync%3DTrue)%0A%20%20%20%20%20%20%20%20_esm%20%3D%20'%5Cnfunction%20crc32(data)%20%7B%5Cn%20%20const%20t%20%3D%20new%20Uint32Array(256)%3B%5Cn%20%20for%20(let%20i%20%3D%200%3B%20i%20%3C%20256%3B%20i%2B%2B)%20%7B%5Cn%20%20%20%20let%20c%20%3D%20i%3B%5Cn%20%20%20%20for%20(let%20j%20%3D%200%3B%20j%20%3C%208%3B%20j%2B%2B)%20c%20%3D%20(c%20%26%201)%20%3F%20(0xEDB88320%20%5E%20(c%20%3E%3E%3E%201))%20%3A%20(c%20%3E%3E%3E%201)%3B%5Cn%20%20%20%20t%5Bi%5D%20%3D%20c%3B%5Cn%20%20%7D%5Cn%20%20let%20crc%20%3D%200xFFFFFFFF%3B%5Cn%20%20for%20(const%20b%20of%20data)%20crc%20%3D%20t%5B(crc%20%5E%20b)%20%26%200xFF%5D%20%5E%20(crc%20%3E%3E%3E%208)%3B%5Cn%20%20return%20(crc%20%5E%200xFFFFFFFF)%20%3E%3E%3E%200%3B%5Cn%7D%5Cn%5Cnfunction%20makeZip(fileList)%20%7B%5Cn%20%20const%20enc%20%3D%20new%20TextEncoder()%3B%5Cn%20%20const%20locals%20%3D%20%5B%5D%2C%20cd%20%3D%20%5B%5D%3B%5Cn%20%20let%20offset%20%3D%200%3B%5Cn%20%20for%20(const%20%7B%20name%2C%20data%20%7D%20of%20fileList)%20%7B%5Cn%20%20%20%20const%20nb%20%3D%20enc.encode(name)%3B%5Cn%20%20%20%20const%20crc%20%3D%20crc32(data)%3B%5Cn%20%20%20%20const%20lh%20%3D%20new%20DataView(new%20ArrayBuffer(30%20%2B%20nb.length))%3B%5Cn%20%20%20%20lh.setUint32(0%2C%200x04034b50%2C%20true)%3B%20lh.setUint16(4%2C%2020%2C%20true)%3B%5Cn%20%20%20%20lh.setUint16(6%2C%200%2C%20true)%3B%20lh.setUint16(8%2C%200%2C%20true)%3B%5Cn%20%20%20%20lh.setUint16(10%2C%200%2C%20true)%3B%20lh.setUint16(12%2C%200%2C%20true)%3B%5Cn%20%20%20%20lh.setUint32(14%2C%20crc%2C%20true)%3B%20lh.setUint32(18%2C%20data.length%2C%20true)%3B%5Cn%20%20%20%20lh.setUint32(22%2C%20data.length%2C%20true)%3B%20lh.setUint16(26%2C%20nb.length%2C%20true)%3B%5Cn%20%20%20%20lh.setUint16(28%2C%200%2C%20true)%3B%5Cn%20%20%20%20new%20Uint8Array(lh.buffer%2C%2030).set(nb)%3B%5Cn%20%20%20%20const%20ce%20%3D%20new%20DataView(new%20ArrayBuffer(46%20%2B%20nb.length))%3B%5Cn%20%20%20%20ce.setUint32(0%2C%200x02014b50%2C%20true)%3B%20ce.setUint16(4%2C%2020%2C%20true)%3B%5Cn%20%20%20%20ce.setUint16(6%2C%2020%2C%20true)%3B%20ce.setUint16(8%2C%200%2C%20true)%3B%20ce.setUint16(10%2C%200%2C%20true)%3B%5Cn%20%20%20%20ce.setUint16(12%2C%200%2C%20true)%3B%20ce.setUint16(14%2C%200%2C%20true)%3B%5Cn%20%20%20%20ce.setUint32(16%2C%20crc%2C%20true)%3B%20ce.setUint32(20%2C%20data.length%2C%20true)%3B%5Cn%20%20%20%20ce.setUint32(24%2C%20data.length%2C%20true)%3B%20ce.setUint16(28%2C%20nb.length%2C%20true)%3B%5Cn%20%20%20%20ce.setUint16(30%2C%200%2C%20true)%3B%20ce.setUint16(32%2C%200%2C%20true)%3B%5Cn%20%20%20%20ce.setUint16(34%2C%200%2C%20true)%3B%20ce.setUint16(36%2C%200%2C%20true)%3B%5Cn%20%20%20%20ce.setUint32(38%2C%200%2C%20true)%3B%20ce.setUint32(42%2C%20offset%2C%20true)%3B%5Cn%20%20%20%20new%20Uint8Array(ce.buffer%2C%2046).set(nb)%3B%5Cn%20%20%20%20locals.push(new%20Uint8Array(lh.buffer)%2C%20data)%3B%5Cn%20%20%20%20cd.push(new%20Uint8Array(ce.buffer))%3B%5Cn%20%20%20%20offset%20%2B%3D%2030%20%2B%20nb.length%20%2B%20data.length%3B%5Cn%20%20%7D%5Cn%20%20const%20cdSize%20%3D%20cd.reduce((s%2C%20b)%20%3D%3E%20s%20%2B%20b.length%2C%200)%3B%5Cn%20%20const%20eocd%20%3D%20new%20DataView(new%20ArrayBuffer(22))%3B%5Cn%20%20eocd.setUint32(0%2C%200x06054b50%2C%20true)%3B%5Cn%20%20eocd.setUint16(8%2C%20fileList.length%2C%20true)%3B%20eocd.setUint16(10%2C%20fileList.length%2C%20true)%3B%5Cn%20%20eocd.setUint32(12%2C%20cdSize%2C%20true)%3B%20eocd.setUint32(16%2C%20offset%2C%20true)%3B%5Cn%20%20const%20parts%20%3D%20%5B...locals%2C%20...cd%2C%20new%20Uint8Array(eocd.buffer)%5D%3B%5Cn%20%20const%20total%20%3D%20parts.reduce((s%2C%20b)%20%3D%3E%20s%20%2B%20b.length%2C%200)%3B%5Cn%20%20const%20out%20%3D%20new%20Uint8Array(total)%3B%5Cn%20%20let%20pos%20%3D%200%3B%5Cn%20%20for%20(const%20p%20of%20parts)%20%7B%20out.set(p%2C%20pos)%3B%20pos%20%2B%3D%20p.length%3B%20%7D%5Cn%20%20return%20out%3B%5Cn%7D%5Cn%5Cnasync%20function%20render(%7B%20model%2C%20el%20%7D)%20%7B%5Cn%20%20el.style.cssText%20%3D%20%22font-family%3Asans-serif%3Bdisplay%3Aflex%3Bflex-direction%3Acolumn%3Bgap%3A8px%3B%22%3B%5Cn%20%20const%20labelRow%20%3D%20document.createElement(%22div%22)%3B%5Cn%20%20labelRow.style.cssText%20%3D%20%22display%3Aflex%3Bgap%3A8px%3Balign-items%3Acenter%3B%22%3B%5Cn%20%20const%20labelLbl%20%3D%20document.createElement(%22label%22)%3B%5Cn%20%20labelLbl.textContent%20%3D%20%22Session%20name%3A%22%3B%5Cn%20%20labelLbl.style.cssText%20%3D%20%22font-size%3A13px%3Bwhite-space%3Anowrap%3B%22%3B%5Cn%20%20const%20labelInput%20%3D%20document.createElement(%22input%22)%3B%5Cn%20%20labelInput.type%20%3D%20%22text%22%3B%5Cn%20%20labelInput.value%20%3D%20model.get(%22label%22)%3B%5Cn%20%20labelInput.style.cssText%20%3D%20%22font-size%3A13px%3Bpadding%3A3px%208px%3Bborder%3A1px%20solid%20%23aaa%3Bborder-radius%3A4px%3Bwidth%3A220px%3B%22%3B%5Cn%20%20labelInput.oninput%20%3D%20()%20%3D%3E%20%7B%20model.set(%22label%22%2C%20labelInput.value)%3B%20model.save_changes()%3B%20%7D%3B%5Cn%20%20labelRow.appendChild(labelLbl)%3B%5Cn%20%20labelRow.appendChild(labelInput)%3B%5Cn%20%20el.appendChild(labelRow)%3B%5Cn%20%20const%20btnRow%20%3D%20document.createElement(%22div%22)%3B%5Cn%20%20btnRow.style.cssText%20%3D%20%22display%3Aflex%3Bgap%3A10px%3Balign-items%3Acenter%3B%22%3B%5Cn%20%20const%20btn%20%3D%20document.createElement(%22button%22)%3B%5Cn%20%20btn.textContent%20%3D%20%22%F0%9F%92%BE%20Save%20session%20(.zip)%22%3B%5Cn%20%20btn.style.cssText%20%3D%20%22padding%3A6px%2016px%3Bfont-size%3A14px%3Bcursor%3Apointer%3Bborder-radius%3A4px%3Bborder%3A1px%20solid%20%23888%3Bwidth%3Afit-content%3B%22%3B%5Cn%20%20const%20statusSpan%20%3D%20document.createElement(%22span%22)%3B%5Cn%20%20statusSpan.style.fontSize%20%3D%20%2213px%22%3B%5Cn%20%20btnRow.appendChild(btn)%3B%5Cn%20%20btnRow.appendChild(statusSpan)%3B%5Cn%20%20el.appendChild(btnRow)%3B%5Cn%20%20const%20hint%20%3D%20document.createElement(%22div%22)%3B%5Cn%20%20hint.style.cssText%20%3D%20%22font-size%3A12px%3Bcolor%3A%23888%3B%22%3B%5Cn%20%20hint.textContent%20%3D%20%22Downloads%20a%20.zip%20file.%20Unzip%20it%20to%20get%20the%20session%20folder%2C%20then%20use%20%5C'Open%20session%20folder%5C'%20above%20to%20restore.%22%3B%5Cn%20%20el.appendChild(hint)%3B%5Cn%20%20btn.onclick%20%3D%20()%20%3D%3E%20%7B%5Cn%20%20%20%20try%20%7B%5Cn%20%20%20%20%20%20const%20label%20%3D%20(model.get(%22label%22)%20%7C%7C%20%22session%22).replace(%2F%5B%5C%5C%2F%5C%5C%5C%5C%3A*%3F%22%3C%3E%7C%5D%2Fg%2C%20%22_%22)%3B%5Cn%20%20%20%20%20%20statusSpan.style.color%20%3D%20%22%22%3B%5Cn%20%20%20%20%20%20statusSpan.textContent%20%3D%20%22%E2%8F%B3%20Building%20ZIP%E2%80%A6%22%3B%5Cn%20%20%20%20%20%20const%20enc%20%3D%20new%20TextEncoder()%3B%5Cn%20%20%20%20%20%20const%20fileList%20%3D%20Object.entries(model.get(%22files%22)).map((%5Bname%2C%20content%5D)%20%3D%3E%20(%7B%5Cn%20%20%20%20%20%20%20%20name%2C%20data%3A%20enc.encode(content)%5Cn%20%20%20%20%20%20%7D))%3B%5Cn%20%20%20%20%20%20const%20zipData%20%3D%20makeZip(fileList)%3B%5Cn%20%20%20%20%20%20const%20blob%20%3D%20new%20Blob(%5BzipData%5D%2C%20%7B%20type%3A%20%22application%2Fzip%22%20%7D)%3B%5Cn%20%20%20%20%20%20const%20url%20%3D%20URL.createObjectURL(blob)%3B%5Cn%20%20%20%20%20%20const%20a%20%3D%20document.createElement(%22a%22)%3B%5Cn%20%20%20%20%20%20a.href%20%3D%20url%3B%20a.download%20%3D%20label%20%2B%20%22.zip%22%3B%5Cn%20%20%20%20%20%20document.body.appendChild(a)%3B%20a.click()%3B%20document.body.removeChild(a)%3B%5Cn%20%20%20%20%20%20setTimeout(()%20%3D%3E%20URL.revokeObjectURL(url)%2C%201000)%3B%5Cn%20%20%20%20%20%20statusSpan.textContent%20%3D%20%22%E2%9C%85%20Downloaded%20%22%20%2B%20label%20%2B%20%22.zip%22%3B%5Cn%20%20%20%20%20%20statusSpan.style.color%20%3D%20%22green%22%3B%5Cn%20%20%20%20%7D%20catch(e)%20%7B%5Cn%20%20%20%20%20%20statusSpan.textContent%20%3D%20%22%E2%9D%8C%20%22%20%2B%20e.message%3B%5Cn%20%20%20%20%20%20statusSpan.style.color%20%3D%20%22red%22%3B%5Cn%20%20%20%20%7D%5Cn%20%20%7D%3B%5Cn%7D%5Cnexport%20default%20%7B%20render%20%7D%3B%5Cn'%0A%0A%0A%20%20%20%20return%20(SessionSaver%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20SessionSaver%2C%0A%20%20%20%20df_avg%2C%0A%20%20%20%20flat_x_sel%2C%0A%20%20%20%20flat_y_sel%2C%0A%20%20%20%20flat_z_sel%2C%0A%20%20%20%20height_threshold%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20np%2C%0A%20%20%20%20pd%2C%0A%20%20%20%20selector%2C%0A%20%20%20%20x_max%2C%0A%20%20%20%20x_min%2C%0A%20%20%20%20y_max%2C%0A%20%20%20%20y_min%2C%0A)%3A%0A%20%20%20%20import%20datetime%0A%0A%20%20%20%20_SEC_COLORS%20%3D%20%5B%22%23e41a1c%22%2C%22%23377eb8%22%2C%22%234daf4a%22%2C%22%23984ea3%22%2C%22%23ff7f00%22%2C%22%23a65628%22%5D%0A%20%20%20%20_color_to_sec%20%3D%20%7Bc%3A%20f%22S%7Bi%2B1%7D%22%20for%20i%2C%20c%20in%20enumerate(_SEC_COLORS)%7D%0A%20%20%20%20_pt_colors%20%3D%20selector.value%5B%22value%22%5D%0A%0A%20%20%20%20_roi_df%20%3D%20pd.DataFrame(%7B%0A%20%20%20%20%20%20%20%20%22x_mm%22%3A%20flat_x_sel%2C%20%22y_mm%22%3A%20flat_y_sel%2C%0A%20%20%20%20%20%20%20%20%22height_mm%22%3A%20flat_z_sel%20%2F%201000%2C%20%22height_um%22%3A%20flat_z_sel%2C%0A%20%20%20%20%20%20%20%20%22section%22%3A%20%5B_color_to_sec.get(c%2C%20%22unassigned%22)%20for%20c%20in%20_pt_colors%5D%2C%0A%20%20%20%20%7D)%0A%20%20%20%20_stats_rows%20%3D%20%5B%5D%0A%20%20%20%20for%20_i%2C%20_c%20in%20enumerate(_SEC_COLORS)%3A%0A%20%20%20%20%20%20%20%20_mask%20%3D%20np.array(%5Bv%20%3D%3D%20_c%20for%20v%20in%20_pt_colors%5D)%0A%20%20%20%20%20%20%20%20_z%20%3D%20flat_z_sel%5B_mask%5D%0A%20%20%20%20%20%20%20%20_stats_rows.append(%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22section%22%3A%20f%22S%7B_i%2B1%7D%22%2C%20%22n_pts%22%3A%20int(len(_z))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22mean_um%22%3A%20round(float(np.mean(_z))%2C%202)%20if%20len(_z)%20else%20None%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22sd_um%22%3A%20%20%20round(float(np.std(_z))%2C%20%202)%20if%20len(_z)%20else%20None%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22min_um%22%3A%20%20round(float(np.min(_z))%2C%20%202)%20if%20len(_z)%20else%20None%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22max_um%22%3A%20%20round(float(np.max(_z))%2C%20%202)%20if%20len(_z)%20else%20None%2C%0A%20%20%20%20%20%20%20%20%7D)%0A%0A%20%20%20%20session_saver%20%3D%20mo.ui.anywidget(SessionSaver(%0A%20%20%20%20%20%20%20%20files%3D%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22avg_grid.csv%22%3A%20%20%20%20%20%20%20df_avg.to_csv(index%3DFalse)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22roi_points.csv%22%3A%20%20%20%20%20_roi_df.to_csv(index%3DFalse)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22section_stats.csv%22%3A%20%20pd.DataFrame(_stats_rows).to_csv(index%3DFalse)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22session_params.csv%22%3A%20pd.DataFrame(%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%22param%22%3A%20%22height_threshold%22%2C%20%22value%22%3A%20height_threshold.value%7D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%22param%22%3A%20%22x_min%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%22value%22%3A%20x_min.value%7D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%22param%22%3A%20%22x_max%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%22value%22%3A%20x_max.value%7D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%22param%22%3A%20%22y_min%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%22value%22%3A%20y_min.value%7D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%22param%22%3A%20%22y_max%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%22value%22%3A%20y_max.value%7D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D).to_csv(index%3DFalse)%2C%0A%20%20%20%20%20%20%20%20%7D%2C%0A%20%20%20%20%20%20%20%20label%3Ddatetime.datetime.now().strftime(%22%25Y-%25m-%25d_%25H%25M%25S%22)%2C%0A%20%20%20%20))%0A%20%20%20%20session_saver%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20return%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A