j)dZddlZddlZddlZddlmZddlZejdddlm Z ddl Z ddl Z ddlmZddlmZddlmZejdeeejddlmZmZmZmZmZd Zd Z d Z!d Z"ed z Z#dZ$dZ%ddZ&dZ'dZ(dZ)dZ*dZ+e,dkr e*dSdS)aNoisy versions of the report's meshes, topology kept: ../meshes/_noise.ply for sigma = 0.2, 0.6, 1 % of the bounding-box diagonal, and ../report/noisy.html (a validation table and renders). Each mesh is first subdivided to edges of at most MAX_EDGE of the diagonal (same shape: a coarse mesh would otherwise only move its few big faces). Every vertex x then moves by u(x), a smooth random 3D vector field: per component, white noise on a grid Gaussian-filtered to a correlation length of CORR sigma, read at the vertices (trilinear), scaled to a standard deviation of sigma (so the displacement along any direction, the normal included, has std sigma). The correlation length is raised until the field's largest gradient (Frobenius norm over the grid, an upper bound of its Lipschitz constant) is below LIP: then x -> x + u(x) is one-to-one in all of space, so the surface cannot cross itself or turn a face over, and a thin wall moves as a whole. Faces are unchanged (same connectivity, watertight). Checked: no face turned over, and the exact sections at the 128 layer centres along up are simple rings. Tried first and dropped: per-vertex noise along the normal folds faces at this size, and noise along the normal (even smooth) makes surfaces cross at sharp edges: a wall moving in and the floor next to it moving up meet. env -u PYTHONPATH PYTHONWARNINGS=ignore .venv/bin/python ~/code/data/meshes/fcp_2710_tests/scripts/make_noisy.py N)PathAgg)Poly3DCollection)ndimage) LinearRing)CSSOBJREPORTROOTUP)gMb`?g~jtx?g{Gz?g{Gzt?g@g?meshesc8d|zdddS)Ndg.p)replace)ss @/home/joao/code/data/meshes/fcp_2710_tests/scripts/make_noisy.pylabelr*s Ag>> ! !#s + ++c(tj||dddf||dddfz ||dddf||dddfz }|tjtj|dddz S)NrT)axiskeepdimsgYn)npcrossmaximumlinalgnorm)VFns r face_normalsr%.s 1QQQT7a!!!Q$j(!AaaadG*q111a4z*ABBA rz")..T.BBFKK KKrctj|t}|j|z}||}}|tj|dz||z z|z z}|||z|||z }dx}} |D]@} | | jD]3} t| dkr|dz }| t| j z } 4A|| fS)zU(rings, rings that self-cross) over the exact sections at the layer centres along up.g?) plane_origin plane_normalheightsrNr) rasarrayfloatverticesminmaxarangesection_multiplanediscretelenr is_simple) meshuplayershlohizssecr$badrLs r sections_okr@3s B  B A UUWWaeeggB ry  3&273f< z vector_field..Lsa      !4!4U!;!;!B!B2:!N!NV\_`V`gm n n n   rcDg|]}||z zS)std)rKrrSs rrPz vector_field..Ns*(((cAEEGGm (((rc3RK|]!}tj|D] }|dzV "dS)rN)rgradient)rKrgrr9s r zvector_field..Os@KKQ9J9JKK2B!GKKKKKKKrcTg|]$}tj|z z jd%S)r)order)rmap_coordinatesT)rKrr"r9r:s rrPz vector_field..Ps6SSSA')!q2vl-=QGGGSSSr) r/tuplerceilr0rIintranger-sqrtsumstack) r"rSrMrNGlipur9r:rOs ```` @@@r vector_fieldrfFs@  A qAJ B "'15588a&j025:;;BB3GG!K L LE       ((   A((((a(((A KKKKQKKKKKLLPPRR S SC SSSSSSQRSSSUVWWA c6Mrc jt|z|z} t|j||z||\}}|tkrn|dz}/ddlm}ddlm}|j|j |j } } } | tj t|t} } t|j|z| | zddkx}rbd| | |<| | | dddf| | dddfz}||tjt||dddf|dddffft|fdz d \}}| }tj|| D]A}tj| ||kzd}| |d||<Bt|j|z| | zddkx}b|j|z}t/t|| | zddk}t1j|| d |||z ||fS) zo(noisy mesh, displacements u, correlation length used / diag, the field's largest gradient, faces turned over).Tg?r) coo_matrix)connected_componentsrNr)rOF)directedprocess)CORRrfr.LIP scipy.sparserhscipy.sparse.csgraphrifacesr% edges_uniquecopyrzerosr4boolraanyravelonesuniquenonzeromeanr^trimeshTrimesh)r6rEdiagrNrMrerdrhrir#n0Eu0SterLlabcvsr"turneds rnoisyrTs E\D FdmUT\63GG3 99 $  ('''''999999z4,d.?1rA FFHHbhs1vvt,,Bdma/33b8==a@@AE E1 J J L L#!A$**,, a!!!Q$j1Qqqq!tW:% &%%jj"'#a&&//AaaadGQqqqRStWCU1V_bcd_e_e^gjk^k&l&l&lw|}}}3 GGII3q6"" # #AAN++A.BrFKKNNAbEEdma/33b8==a@@AE E1 J J L L# A <1%%*//22a7<<>> ? ?F ?1a / / /FTM3 NNrc tjdt|zdfd}tjgd}|tj|z}|dj}|d| d}}||zdz ||z dz } }tt||D]h\} \} } | dt|| dzd d } | j }d d tj||zzz}| t!| jtj|||fgdzdd| |d| z |d| z| |d| z |d| z| |d| z |d| z| d| | dd| | dj|||tj|dS)Nr+g@r)figsizedpi)g?g333333g?rrr3dortho) projection proj_typeg??)g333333?g= ףp=?rnone) facecolors edgecolors linewidths)rrri )fontsize)pltfigurer4rarrayr r!r.r/r0 enumeratezip add_subplotr%absadd_collection3dr trianglesc_set_xlimset_ylimset_zlimset_box_aspect set_axis_off view_init set_title tight_layoutsavefigclose)r titlespathfiglightallvr:r;mrkr6raxfnshades rrenderrqs_ *a#f++os3 = = =C H''' ( (E%29>>%3H3H*H% !9 D XXa[[$((1++B Gq=27--//A-qA!#ff"5"566$$ 9D! __QF QUtw_ W W  tbfR%Z0000 ,T^eUZ\aNaHbewewewHw8>1NNN O O O AaD1HadQh'''QqTAXqtax)H)H)H"++VWXYVZ]^V^`abc`dgh`hJiJiJi )$$$boo&7&7&7b#9N9N9N Q ####KKIcNNNNNrctdg}tjD](\}}t j}t j|d}ttj |j }t j |j|jt"|zd\}}t j||d}|jr|js J||j} t-j|d} t1|| \} } |gd t3|jd d g}} t5t6D]\}}t9|||tjd |z\}}}}}|t|d tA|dz t1|| \}}|| z!d}tE||t3|jtG|jt|j$|j$z t|%|z ||||| ||  }|&|| &||&dd|zddd|dzddtO|ddd|zddd|dzddd|zdd|dd |d!|d"|d#| d"| d$|d%d&d'tQ| |tRd(|d)z tOd*|d+t j|z dd,d'*td-z }d.|D|*r2fd/tWj,|-Dng}||z}|.tWj/|d0ta|dS)1NT)exist_okr6)force)max_iterrk)rrrzclean (subdivided, ,z faces)i_noisez.plyr) namerErq watertight volume_ratiorScorrrdrrings rings_clean bad_rings bad_cleanuσ = rrz % diag (along the normal rS.2fz %)16sz sigma z4.1fz% normal std .3fz% corr .1fz% lip z turned z self-crossing rings /z (clean z) vol r.4f)flushz img/noisy_z.pngz z: z sznoisy_meshes.jsonch|] }|d S)rrRrKrs r zmain..s $ $ $!AfI $ $ $rc@g|]}|dv |ddk|S)rloaf_panrR)rKrdones rrPzmain..s7 f f fAfIT4I4IaPVi[eNeNe1NeNeNer)indent)1OUTmkdirr itemstime perf_counterr|loadr-rr r!extentsremeshsubdivide_to_sizer.rqMAX_EDGEr} is_watertightis_winding_consistentvertex_normalsr getr@r4rSIGMASrrandom default_rngexportrradictruvolumerSappendprintrr existsjsonloads read_text write_textdumps write_html)rowsrrt0rr~vfbaser$r7r0b0r rirrrerrdrrr>drowjsoldrs @rmainrsuIItI D9;;IIa    L& ) ) )RY^^AI..//~// AGXPT__a/bb1q!T222!Fd&@FF$FF@   VD) $ $T2&&B"RDJ"R"R"R"R!Sf%% @ @DAq&+D!T29;P;PQUXYQY;Z;Z&[&[ #Aq$V HHSd88%((8888 9 9 9$Q++JE3Q AADQW$qJ_J_$)!(T[*@$A$AuQUUWWW[^G\G\cgmp$ErS\^```C KK    MM!    MM`#'```sSQVZGW```` a a a Ts33cAgs33#E BRs33`cfj`js33j33+133IL33OT33^`33ce33^,233:> @ @ @ @ @ vvv(?T(?(?(??@@@ ;4;;4,..3;;;;4HHHHH " "B $ $t $ $ $Djljsjsjuju } f f f fdj00 f f f f{}C :DMM$*T!,,,---trcdd|D}ddtd|DD}dtddtzdd t dd t dd tjd d |d|d}tdz |dS)Nc3ZK|]}d|ddd|dzddd|dzddd|d zd d|d d d|d dd|drdndd|dd|d|dkrdndd|dd|dd|dd|dd|dddVdS) zrz rrErz %rSrrrrdrrqrryesz norzz / rrrrz NrRrs rrWzwrite_html..snnde _1V9 _ _sQwZ'7 B _ _cAeHn B _ _ehklmsktet B _ _x D _ _&'j D _ _FG o=fUU[f _ _qrs{q| _ _ ~+>>eeB _ _JKK. _ _]^_f]g _ _~ _ _#$M"2 _ _>?~=N^ _ _ _nnnnnnrc3.K|]}d|d|d|dVdS)z

NrR)rKr$s rrWzwrite_html..sWII\]WaWW1WWqWWWIIIIIIrc3&K|] }|dV dS)rNrRrs rrWzwrite_html..sJoHoH}~opqwoxoHoHoHoHoHoHrz Noisy meshes

← report

Noisy meshes

Files: meshes/<name>_noise0p2.ply, _noise0p6, _noise1.ply; script scripts/make_noisy.py. Each mesh is subdivided to edges ≤ rru % of the diagonal (same shape), then every vertex x moves by u(x), a smooth random 3D vector field (seeded): each component Gaussian with std σ = 0.2 / 0.6 / 1 % of the diagonal and a Gaussian correlation length of at least u8 σ, raised until the field's largest gradient is below u. Then x → x + u(x) is one-to-one, so the surface cannot cross itself or turn a face over, whatever the wall thickness or the edges. Faces are unchanged (same connectivity, watertight). "along the normal": the achieved displacement along the vertex normals, over the diagonal. Self-crossing rings: the exact sections at the 128 layer centres along up that are not simple (red: more than the clean mesh). Generated z%Y-%m-%d %H:%Mu.

z
objectσalong the normalcorrelation lengthmax gradientfaces watertightfaces turned overself-crossing ringscleanvolume / clean
z z noisy.html) joinrfromkeysrrrmrnrstrftimer r)rtrimgspages rrrsw nnim nnn n nB 77IIaeananoHoHCGoHoHoHbHbHIII I ID fi   '*Hn        HK   cgboqAcBcB  uw     D l&&t,,,,,r__main__)r&)-__doc__rsysrpathlibr matplotlibusematplotlib.pyplotpyplotrnumpyrr|mpl_toolkits.mplot3d.art3drscipyrshapely.geometryrrinsertstr__file__parent parts_reportrr r r r rrrmrnrrr%r@rfrrrr__name__rRrrrs"  u777777''''''33ttH~~,--...33333333333333     Xo,,,LLL &   OOO:(###L---2 zDFFFFFr