ĀjgSdZddlmZddlZddlZddlZddlZddlmZddl Z e jdddl m Z ddlmZddlZddlZddlmZddlmZdd lmZd d ge j4d <d e j4d<dddZej8gdej:ZddZd dZ d!dZ!d"dZ"d#dZ# d$dZ$d%dZ% d&dZ& d'dZ'd(dZ(d)dZ)e*dk(re)yy)*a+Match STEP/CAD, process XLS data, measured PLY, and defect CSV by run. The exported arrays are diagnostic supervision scaffolds. The supplied CSVs contain a point but no defect width/height, and scan-to-CAD registration is not independently validated, so ``training_valid`` is deliberately false. ) annotationsN)PathAgg)LineCollection) coo_matrix)dijkstra)cKDTreezApple SD Gothic Neoz DejaVu Sansz font.familyFzaxes.unicode_minus 16h_mirror18h01)16hr )g@g$@g4@cVttjjdS)N)r__file__resolveparents전처리/match_targets.py root_defaultr!s > ! ! # + +A ..rcF||ddddfjz|dddfzS)N)T)pointsmatrixs r transformr%s4 F2A2rr6N$$ $vbqb!e} 44rct||d|dff||fj}tjt |d|tj S)Nrr)shapeF)directedindices)rtocsrnpasarrayrfloat32) node_count edge_index edge_lengthseedgraphs rgeodesic_distancer))sP  jmZ]%CD(*5 77(>(@U3s8_dW   5 4s C((C1c8tjjt|}|j |j tj jj}|jtj jjj||jtjj|jt j"}|dj%}|dj%}||z }|t j&j)|djt j"t!j*d||jt j"fS)Nrprimitive_normalsr)axisij,ij->i)o3dioread_triangle_meshstrremove_duplicated_verticescompute_triangle_normalstgeometryRaycastingScene add_triangles TriangleMesh from_legacycompute_closest_pointscoreTensorastyper!r#numpylinalgnormeinsum) mesh_pathrmeshscenehitclosestnormalsdeltas r closest_cadrd>s 66 $ $S^ 4D##%t'D'D'F EENN * * ,E 33??EF  & &sxxv}}RZZ7P'Q RC(m!!#G%&,,.G W E BIINN5qN188DbiiE7G$$*F2::$6 77rct|Dcic]$\}}t|jddd|&}}}tjdtj tj }tjdtj tj }tddD]5}d|fd|ffD](\} } | |d} | |vs||| s||| | |dz <*7||fScc}}w) Nz [rru 상부온도u 하부온도02d) enumeraterLsplitr!fullnanr#range) namesvaluesvalidinamelookuptopbottomzoneprefixtargetkeys r heater_valuesrzKsrf) r=r!ptpfloorminmaxrXint32clip)cadrtruscalervs rprovisional_heater_projectionris "&&QT# $E 88SAYQT0Ct4DDrI J Q QRTRZRZ [D 774B D TF4L3<*> >>rc |jjddtjdddd\}}|j }t j d|jdd}t j||d|fddddgf||d|fddddgfgd}t j|| g}|t j|}|dj|d d t|jt|j }|djddgd d g|djt!dt!dd|dj#d|dj%d|j'||ddtt j(|}t!dD]}t!dD]t}|||f}t j|r;|dk(r| n| |}|r|ddn|d}|dj+|||ddd||krdndZ|dj+||dddddv|dddft j,dk\}t/dd t|jt|j }|d| d!f|d| d"ffD]\}} }!|j0|dddf|dddffd#| i|}"|j1||df||dfdd$d%d&'|j1|d|dd(d)d%*|j%|!d+z|j'|"|d|t j|}#tt3|#r1tt jt j4|#nd,d,}$|dj1|dddf|dddf|dd-|$ |$.}"|dj1||df||dfdd$d%d&'|dj1|d|dd(d)d/*|dj%d0|j'|"|dd1t j6dt3|dz td2t3|t83}%|d4j1||%df||%df||%d5d6dtt j:|d7.}|j'||d4d8|d4j%d9|d:j=t?|d;d|d:j1|dddf|dddf|dddfdd?dd.|d:j1|d|dd@d)dA*|d:jA|d:j%dB|ddD]5}|jCdC|j#dD|jEdE7|jG|dF|z|jI|dGHtjJ|y)INTrexist_okr|r)rf )figsizeconstrained_layoutrrautoinferno)aspectcmapvminvmaxuupper / 상부ulower / 하부rfrgzheater zone numberu52×17 sandwich heater map (* = linearly interpolated)utemperature (°C))axlabelz.1f*z.0fcenterwhiteblack)havafontsizecolorNAcyang@)srrrzupper heater overlayzlower heater overlaycnonez#00e5ffgffffff?)r facecolors edgecolors linewidthsFx)rmarkerrz cyan = weak 10 mm AOIg?coolwarm)rrrrrz#111u0upper − lower heater overlay AOI shown in cyanutemperature difference (°C)i:)dtypeg?viridis_z!nearest CAD surface distance (mm)zCregistered measured PLY deformation proxy, not tracked displacementgQ?z#55768bg?)rcolorsalphamagma<z#00bcd4zBweak 10 mm geodesic area label assumed radius, not observed extentequalz CAD X (mm)z CAD Z (mm)u: STEP + XLS heaters matched to PLY + defect CSV CAD overlay assumes heater 01→17 uniformly ordered along CAD X; physical footprints are required )dpi)&parentmkdirpltsubplotsravelr!arangerstackrimshowr=rr set_yticks set_xticksrm set_xlabel set_titlecolorbarnanmeantextexpdictscatterr;abslinspaceint percentileadd_collectionr autoscale set_aspect set_ylabelsuptitlesavefigclose)&r@runredgesscanscan_absweakanchorrtru top_imputedbottom_imputednode_top node_bottom node_deltanotefigaxestakesegtempfiniteimagemidpointrCcolvalue was_imputedraoicommonrrotitledraw delta_finite delta_limshows& rpreviewrrs  KKdT2 Q8MIC ::RZZ'(HQx9CcNE{{5!.1Qh{NCP +65+QuSkQ S#uhQR.3h.>7GMQ S#tXPQY_ ` q!t*s #C CieFJJL.AfjjlH[ \F#Aw2HI#Aw 5KLNFErzz#ad)SAYC&CFC 3sAv;CF rf'B  8 6!9fQi2cY G U667 Tb(; <Nbkk*56L3|;LE"&& !567RTVXYI 7??3q!t9c!Q$i::!*   ##iWZ[\GOOC1Is1a4yDAJ!'PQXYOZGOOF1IvayBsiOHGa**+pq12h g l ;R]]<=XLLC5xxKK#K #rc|dz }|dz }|dz }tj|dz j}tj|dz j}tj|d5}|j D cic]} | ||  } } dddtj|d5}|j D cic]} | ||  } } ddd|j dd |d z } | j d |j dd  d | d | d}}} | ddddf}g}td|dz jdD}tjD]\}}ttj d|k(dd}|d|}tj|dtj}|d|dz }tj|d5}|djtj }dddt#|jtj }t%||\}}}t|j|dk(rdnd}t'|}tj||dtj } t#| d|djtj }!tt)| j+|!d}"t-t/| |||"}#|#dddft0dk}$tj2d|#dddft0dz d zzjtj }%t5| d!| d"|| d#|\}&}'t7|&\}(})t7|'\}*}+t9| |(|*\},}-}.}/t)|}0|0j+| d $\}1}2||2| z }3tj:d%|3|jtj }4t=|j?d&d}5|d'duxr|d(du}6tj@t/| t<}7d)|5d*|6d+}8||d,z }9tjB|9fidtj|d-| d.|d/|d"| d"|d#| d#|d!| d!d0|d1|d2|d3|d4|d5|1jtj d6|4d7tj|d8d9tj|d:tj d;tj|d'xstjD|d(xstjDgtj d<| d=|!d>tj|"tjFd?|#d@t0dA|$dB|%dC|(dD|*dE|&dF|'dGtjH|&dHtjH|'dI|)dJ|+dK|,dL|-dM|.dN|/dOtjddP|dQtjtjJjM|ddRddRfdStj|5dTtj|6dU|7dVdWdXgdYdZgd[d\d]d^d_d`dadbdcdddedf dg}:|9jOdhjQtjR|:dd idjk| |dlz };tU|;|| ||||%|!|(|*|)|+|-|.|/|8|jWidm|dn|do|9jXdpdq|;jXdrt[|j]|dsdtdut[|j]|dvt[|j]|dwt/|dxt/| dy|"dzt_tjJjM|ddRddRfdS|5dT|6d{t|7jadOdd|d}d~dd|8ddddddddddddddd|gdd}<|dz jQtjR|zprepare..sTE1Esdatas_2z 3. CSV File*run_idvariantsrzanalysis/results/z_all_merge_pointcloud_full.npzrr z*16*.csvz*18*.csv raw_xyz_mmrrr|process_feature_names process_xprocess_feature_valid)workersrHindependently_validatedr4r5z%registration independently validated=z; CSV defect extent observed=z6; all targets diagnostic/weak and training_valid=falsez_cad_ply_targets.npz cad_pos_mm cad_normalcad_edge_indexscan_pos_raw_mmscan_pos_cad_mmscan_closest_cad_mmscan_surface_distance_mmscan_normal_residual_mmcad_to_scan_distance_mmcad_to_scan_normal_residual_mm defect_classr6defect_pixel_xyr3defect_width_height_pxdefect_anchor_raw_mmdefect_anchor_cad_mmdefect_anchor_cad_node#geodesic_distance_to_defect_seed_mmweak_area_radii_mmweak_defect_maskweak_defect_probability heater_top_cheater_bottom_cheater_top_c_rawheater_bottom_c_rawheater_top_observedheater_bottom_observedheater_top_imputedheater_bottom_imputedheater_zone_index_zero_based%provisional_heater_top_c_per_cad_node(provisional_heater_bottom_c_per_cad_node'provisional_heater_delta_c_per_cad_nodeheater_projection_validtransform_raw_ply_to_cadtransform_determinantr$registration_independently_validateddefect_extent_observedtraining_validzone physical production runz(CAD node/edge geometry derived from STEPz%XLS process values and validity masksz#fused all_merge_pointcloud_full PLYz defect CSVz/T1/T2/TOP are registration and coverage QA onlyz3PLY point to nearest CAD surface; deformation proxyz5oriented normal component of nearest-surface residualz$CAD node to nearest fused scan pointz0CAD-node geodesic dilation at weak_area_radii_mmz.Gaussian weak label around point-only CSV seedz?upper/lower heater channels after internal linear interpolationzsource values with NaN retainedzKtrue only where an internal NaN was interpolated between valid zone numbersz/uniform CAD-X projection for visualization onlyz@false until registration and spatial target extent are validated) rrrr(r)zheater_top_c/heater_bottom_czheater_*_c_rawzheater_*_imputed!provisional_heater_*_per_cad_noder;) sample_unit model_inputpost_process_target_sourcescomponent_scansarraysz .schema.json) ensure_asciiindentzutf-8)r,z_cad_ply_targets.pngrposenpzrzvisualizations/process_source cad_sourcezdatas_2/AC28_DT RR UPR STD.stp ply_source defect_source scan_points cad_nodesdefect_anchor_noderegistration_determinanttraining_valid_nodesheater_projectionu5provisional uniform zone 01→17 ordering along CAD Xznlinear interpolation only for missing zones bracketed by valid neighboring zone numbers; no edge extrapolationzBall_merge_pointcloud_full only: fused post-process surface outcomez]T1_DOWN/T1_UP/T2_LEFT/T2_RIGHT/TOP are registration, coverage, and sensor-consistency QA only)heater_missing_value_methodprimary_scan_rolecomponent_scan_rolestatuszgSTEP design + XLS process conditions -> CAD-domain defect risk/shape segmentation and deformation fieldzmpost-process supervision source during training/validation; not an inference input for pre-process predictionzTall_merge_pointcloud_full.ply (represented by the existing deterministic 2 mm cache)zVT1/T2/TOP are not independent samples or inference inputs; retain for pose/coverage QA)primaryr@zVmatch CAD and fused PLY geometry, then verify residuals by region and component sensorz-a physical rigid transform has determinant +1zcurrent best-looking mapping has determinant -1; treat only as coordinate-handedness/reflection diagnostic until scanner axes are confirmed)geometry_featuresproper_rigid_requirement 16h_statuszWgeodesic dilation of a CSV seed at assumed 5/10/20 mm radii; not observed defect extentzPregistered scan-to-nearest-CAD surface residual; not material-point displacementzWnormal component of nearest-surface residual; sign depends on oriented CAD tessellationzKdiagnostic uniform X-bin overlay; not a validated physical heater footprint)r(rrr<)z5confirm CSV image origin and pixel-to-PLY calibrationzUobtain physically valid independently verified PLY-to-STEP registration for every runzFcollect defect polygons/masks or W/H instead of point-only coordinateszNcollect many OK and NG runs/design revisions and split by run and CAD revisionz[provide heater-zone coordinates/footprints and their relation to the STEP coordinate system)problem role_of_plyscan_selectionregistration_policytarget_semanticsentriesblocking_validationz index.jsonr+c3K|]^}dtj|ddtj|dd|dd|ddtj|dd `y w) z

rz

rSz

raw NPZ · metadata

matched CAD PLY targets for z
N)htmlescape)r es rr zprepare..TsA  AeH-.ht{{1X;7O6PQxj!y\N"EdkkRSTYRZF[E\\h jsA$A&z index.htmluCAD–PLY matched targets

CAD–PLY matched supervision

PLY is an outcome used to build labels. At prediction time the intended inputs are STEP geometry and XLS process conditions.

zcad_ply_target_summary.json)2jsonloads read_textr!loadfilesrnextglobRUN_POSEitemsrwherer"float64rXr#rrdrDr queryr)r; WEAK_RADII_MMrrzrrr\rr>rsavez_compressedrlrrrZdet with_suffix write_textdumpsrappendrrrL relative_tor=sumjoin)>rootout result_dir graph_path table_pathr]metarrryr(tablevisualrrrlengthr]csv_dirrrDrC pose_metar scan_pathscan_rawscan_cad scan_closestr scan_signed defect_pathdefect anchor_raw anchor_cad anchor_nodegeo weak_masksweak_probabilityheater_top_rawheater_bottom_raw heater_topr0 heater_bottomr1 heater_zonerrr scan_tree cad_scan_abscad_scan_index cad_deltacad_scan_signedregistration_validextent_observedr;rSoutputschemapngmanifestcardss> rpreparersl ??JIIJ99I ::tJJUUW XD**d%YYddfgK % 0C*-))4)3c#h)4 1 % 0C*-))4)3c#h)4 1IIdTI*S;K5KFV\\cg\MhTD1x%/5;NC ; 1 %FGTti/55nETTG^^% T"((5?c1215a89$T* Ik2BJJ?.se3QRR WWYU 3s8}++BJJ7H4Xv.55bjjA.9)X.N+ h 7< K z$/8;BB2::N '#,,,Z8;< C%ED\]4%88 66#QW d8K)KPQ(Q"QRYYZ\ZdZde,9 ) *E+,>s,CUKbEcdgEh-j)))D^)T& &/JK\/], ,9V ]:,6 X{JH% '0sA'F$ n^,s2 ))J 6BII"**U!)--0I5"QR ,D8\VK=PX\=\#c(D134F3GH'((^ ` #233   ::c? /2 ?E ! -2+,>s,C #((?"@"E  #((?"@  %  7?  !- HP %0 %1$7$7 $C ,; F7O4 GIjjQWXbQcegeoeoFp $&::vj/A/KRVVVT_M`Mjdfdjdj.kmomwmw#x ", CM $&::k288#D 14 -! @J! "%5# $$% $6C% &,' &BS' (!# N ;) (UWT_T_`qTr) * 2+ *J_+ ,*5- .3;/ 06A1 25?3 4%'JJu$55 6&,7 6DF::biimm\bcedecegihigici\jNkCl7 824; :P^; @9FHop,QS_+`P,a+b+Q$V+[0q"C$q5f"d   $ >*55 JJvE! ? AMOo  #i33D9: =LSQ\QhQhimQnMo  3x=  +6s3x  BVWb  'biimmF2A2rr6N.K(L  34F   %o 8NsSaSeSeSgOh  &u  !X ,\!e$C  }&^}Gmw "z(Wh !z(z(A1~   'H6 <##DJJxeTU$Vah#i GG  E  <## cfk  k  $//;;DJJx^clmD$))"3"3"54 (( 4d44Cll5d%55GtS'*H $**c(9"56s3xHW\ ]^r__main__)returnr)r np.ndarrayrrrr) r$rr%rr&rr'rrr)r@rrdict[str, object])r]rrrrz)tuple[np.ndarray, np.ndarray, np.ndarray])rnrrorrprrtuple[np.ndarray, np.ndarray])rorrr)rrrtrrurrz5tuple[np.ndarray, np.ndarray, np.ndarray, np.ndarray])"r@rrrLrrrrrrrrrrrrrtrrurrrrrrrrrrrrrLrNone)ryrrzrr{rrr)rr)+__doc__ __future__rrr9r`rcpathlibr matplotlibusematplotlib.pyplotpyplotrmatplotlib.collectionsrrYr!open3drI scipy.sparserscipy.sparse.csgraphr scipy.spatialr rcParamsrjr"r#rorrr)rDrdrzrrrrr__name__rrrrs #  u1#)!4mD ]%* !"' 2 ,bjj9 /5Q#-Q58Q=GQ  7 # (E "?+5?:o?D/D/(2D/:DD/D/%/D/9CD/$D/6@D/! D/0: D/HR D/  D/  D/NhV _ zFr