rm-graphical-computing 1.0.6 → 1.0.8

package/dist/index.cjs

@@ -1 +1 @@
-Object.defineProperty(exports,Symbol.toStringTag,{value:`Module`});var e=Object.create,t=Object.defineProperty,n=Object.getOwnPropertyDescriptor,r=Object.getOwnPropertyNames,i=Object.getPrototypeOf,a=Object.prototype.hasOwnProperty,o=(e,i,o,s)=>{if(i&&typeof i==`object`||typeof i==`function`)for(var c=r(i),l=0,u=c.length,d;l<u;l++)d=c[l],!a.call(e,d)&&d!==o&&t(e,d,{get:(e=>i[e]).bind(null,d),enumerable:!(s=n(i,d))||s.enumerable});return e},s=(n,r,a)=>(a=n==null?{}:e(i(n)),o(r||!n||!n.__esModule?t(a,`default`,{value:n,enumerable:!0}):a,n));let c=require(`three`);c=s(c);function l(e,t=!1){let n=e[0].index!==null,r=new Set(Object.keys(e[0].attributes)),i=new Set(Object.keys(e[0].morphAttributes)),a={},o={},s=e[0].morphTargetsRelative,l=new c.BufferGeometry,d=0;for(let c=0;c<e.length;++c){let u=e[c],f=0;if(n!==(u.index!==null))return console.error(`THREE.BufferGeometryUtils: .mergeGeometries() failed with geometry at index `+c+`. All geometries must have compatible attributes; make sure index attribute exists among all geometries, or in none of them.`),null;for(let e in u.attributes){if(!r.has(e))return console.error(`THREE.BufferGeometryUtils: .mergeGeometries() failed with geometry at index `+c+`. All geometries must have compatible attributes; make sure "`+e+`" attribute exists among all geometries, or in none of them.`),null;a[e]===void 0&&(a[e]=[]),a[e].push(u.attributes[e]),f++}if(f!==r.size)return console.error(`THREE.BufferGeometryUtils: .mergeGeometries() failed with geometry at index `+c+`. Make sure all geometries have the same number of attributes.`),null;if(s!==u.morphTargetsRelative)return console.error(`THREE.BufferGeometryUtils: .mergeGeometries() failed with geometry at index `+c+`. .morphTargetsRelative must be consistent throughout all geometries.`),null;for(let e in u.morphAttributes){if(!i.has(e))return console.error(`THREE.BufferGeometryUtils: .mergeGeometries() failed with geometry at index `+c+`.  .morphAttributes must be consistent throughout all geometries.`),null;o[e]===void 0&&(o[e]=[]),o[e].push(u.morphAttributes[e])}if(t){let e;if(n)e=u.index.count;else if(u.attributes.position!==void 0)e=u.attributes.position.count;else return console.error(`THREE.BufferGeometryUtils: .mergeGeometries() failed with geometry at index `+c+`. The geometry must have either an index or a position attribute`),null;l.addGroup(d,e,c),d+=e}}if(n){let t=0,n=[];for(let r=0;r<e.length;++r){let i=e[r].index;for(let e=0;e<i.count;++e)n.push(i.getX(e)+t);t+=e[r].attributes.position.count}l.setIndex(n)}for(let e in a){let t=u(a[e]);if(!t)return console.error(`THREE.BufferGeometryUtils: .mergeGeometries() failed while trying to merge the `+e+` attribute.`),null;l.setAttribute(e,t)}for(let e in o){let t=o[e][0].length;if(t===0)break;l.morphAttributes=l.morphAttributes||{},l.morphAttributes[e]=[];for(let n=0;n<t;++n){let t=[];for(let r=0;r<o[e].length;++r)t.push(o[e][r][n]);let r=u(t);if(!r)return console.error(`THREE.BufferGeometryUtils: .mergeGeometries() failed while trying to merge the `+e+` morphAttribute.`),null;l.morphAttributes[e].push(r)}}return l}function u(e){let t,n,r,i=-1,a=0;for(let o=0;o<e.length;++o){let s=e[o];if(t===void 0&&(t=s.array.constructor),t!==s.array.constructor)return console.error(`THREE.BufferGeometryUtils: .mergeAttributes() failed. BufferAttribute.array must be of consistent array types across matching attributes.`),null;if(n===void 0&&(n=s.itemSize),n!==s.itemSize)return console.error(`THREE.BufferGeometryUtils: .mergeAttributes() failed. BufferAttribute.itemSize must be consistent across matching attributes.`),null;if(r===void 0&&(r=s.normalized),r!==s.normalized)return console.error(`THREE.BufferGeometryUtils: .mergeAttributes() failed. BufferAttribute.normalized must be consistent across matching attributes.`),null;if(i===-1&&(i=s.gpuType),i!==s.gpuType)return console.error(`THREE.BufferGeometryUtils: .mergeAttributes() failed. BufferAttribute.gpuType must be consistent across matching attributes.`),null;a+=s.count*n}let o=new t(a),s=new c.BufferAttribute(o,n,r),l=0;for(let t=0;t<e.length;++t){let r=e[t];if(r.isInterleavedBufferAttribute){let e=l/n;for(let t=0,i=r.count;t<i;t++)for(let i=0;i<n;i++){let n=r.getComponent(t,i);s.setComponent(t+e,i,n)}}else o.set(r.array,l);l+=r.count*n}return i!==void 0&&(s.gpuType=i),s}var d=class extends c.Object3D{constructor(e=document.createElement(`div`)){super(),this.isCSS2DObject=!0,this.element=e,this.element.style.position=`absolute`,this.element.style.userSelect=`none`,this.element.setAttribute(`draggable`,!1),this.center=new c.Vector2(.5,.5),this.addEventListener(`removed`,function(){this.traverse(function(e){e.element&&e.element instanceof e.element.ownerDocument.defaultView.Element&&e.element.parentNode!==null&&e.element.remove()})})}copy(e,t){return super.copy(e,t),this.element=e.element.cloneNode(!0),this.center=e.center,this}};new c.Vector3,new c.Matrix4,new c.Matrix4,new c.Vector3,new c.Vector3;var f=(e,t,n,r,i=1e-8)=>{let a=t.x-e.x,o=t.y-e.y,s=r.x-n.x,l=r.y-n.y,u=a*l-o*s;if(Math.abs(u)<i)return null;let d=((n.x-e.x)*l-(n.y-e.y)*s)/u,f=((n.x-e.x)*o-(n.y-e.y)*a)/u;return d<-i||d>1+i||f<-i||f>1+i?null:{point:new c.Vector3(e.x+d*a,e.y+d*o,e.z)}},p=(e,t,n={})=>{let{fontSize:r=36,fontFamily:i=`sans-serif`,backgroundColor:a=`rgba(0,0,0,0)`,textColor:o=`#ffffff`,padding:s=10,maxWidth:l=256}=n,u=document.createElement(`canvas`).getContext(`2d`);u.font=`bold ${r}px ${i}`;let d=u.measureText(e).width,f=Math.min(d+s*2,l),p=r+s*2,h=document.createElement(`canvas`);h.width=f,h.height=p;let g=h.getContext(`2d`);g.fillStyle=a,g.fillRect(0,0,f,p),g.fillStyle=o,g.font=`bold ${r}px ${i}`,g.textAlign=`center`,g.textBaseline=`middle`,d>l-s*2?m(g,e,f/2,p/2,l,r):g.fillText(e,f/2,p/2);let _=new c.CanvasTexture(h),v=new c.SpriteMaterial({map:_,depthTest:!1,transparent:!0}),y=new c.Sprite(v);y.position.set(t.x,t.y,t.z);let b=.005;return y.scale.set(f*b,p*b,1),y};function m(e,t,n,r,i,a){let o=t.split(``),s=``,c=a*1.2,l=r;for(let t=0;t<o.length;t++){let r=s+o[t];e.measureText(r).width>i&&t>0?(e.fillText(s,n,l),s=o[t],l+=c):s=r}e.fillText(s,n,l)}var h=(e,t,n=.01)=>{let r=Math.sqrt(e.x*e.x+e.y*e.y+e.z*e.z),i=Math.sqrt(t.x*t.x+t.y*t.y+t.z*t.z);if(r<1e-10||i<1e-10)return{parallel:!1,sameDirection:null,angle:NaN};let a={x:e.x/r,y:e.y/r,z:e.z/r},o={x:t.x/i,y:t.y/i,z:t.z/i},s=a.x*o.x+a.y*o.y+a.z*o.z;s=Math.max(-1,Math.min(1,s));let c=Math.acos(Math.abs(s)),l=c<=n;return{parallel:l,sameDirection:l?s>0:null,angle:c}};function g(e,t){let n=new c.BufferGeometry,r=[],i=[];e.forEach(e=>r.push(e.x,e.y,e.z)),t.forEach(e=>r.push(e.x,e.y,e.z)),i.push(0,2,1,0,3,2),i.push(4,5,6,4,6,7),i.push(0,1,5,0,5,4),i.push(1,2,6,1,6,5),i.push(2,3,7,2,7,6),i.push(3,0,4,3,4,7),n.setAttribute(`position`,new c.Float32BufferAttribute(r,3)),n.setIndex(i),n.computeVertexNormals();let a=new c.MeshBasicMaterial({color:`green`,transparent:!0,opacity:.5,side:c.DoubleSide});return new c.Mesh(n,a)}var _=class{constructor(e={}){this.config={value:e.value||`0.0000m`,position:e.position||new c.Vector3(0,0,0),style:{fontFamily:e.style?.fontFamily||`Arial, sans-serif`,fontSize:e.style?.fontSize||`16px`,fontWeight:e.style?.fontWeight||`normal`,fontStyle:e.style?.fontStyle||`normal`},backgroundOpacity:e.backgroundOpacity||.7},this.label=null,this.labelElement=null,this.createLabel()}createLabel(){return this.labelElement=document.createElement(`div`),this.labelElement.textContent=this.config.value,this.applyBlackWhiteStyle(),this.label=new d(this.labelElement),this.label.position.copy(this.config.position),this.label}applyBlackWhiteStyle(){let e=this.labelElement.style,t=this.config;e.fontFamily=t.style.fontFamily,e.fontSize=t.style.fontSize,e.fontWeight=t.style.fontWeight,e.fontStyle=t.style.fontStyle,e.color=`#FFFFFF`,e.textShadow=`0 1px 2px rgba(0, 0, 0, 0.5)`,e.backgroundColor=`rgba(0, 0, 0, ${t.backgroundOpacity})`,e.padding=`6px 12px`,e.borderRadius=`4px`,e.position=`absolute`,e.whiteSpace=`nowrap`,e.userSelect=`none`,e.pointerEvents=`none`,e.zIndex=`1000`,e.textAlign=`center`,e.lineHeight=`1.2`}updateValue(e){this.labelElement&&(this.labelElement.textContent=e)}updatePosition(e){this.label&&this.label.position.copy(e)}updateFontStyle(e){if(this.labelElement){let t=this.labelElement.style;e.fontFamily&&(t.fontFamily=e.fontFamily),e.fontSize&&(t.fontSize=e.fontSize),e.fontWeight&&(t.fontWeight=e.fontWeight),e.fontStyle&&(t.fontStyle=e.fontStyle)}}getLabel(){return this.label}getLabelElement(){return this.labelElement}dispose(){this.labelElement&&this.labelElement.parentNode&&this.labelElement.parentNode.removeChild(this.labelElement),this.label=null,this.labelElement=null}},v=e=>{let t=new c.Group;if(e&&e.length){let n=`green.red.pink.white.azure.yellow.bisque.grey.blanchedalmond.blue.blueviolet.brown.burlywood.cadetblue.chartreuse.chocolate.coral.cornflowerblue.cornsilk.crimson.cyan.darkblue.darkcyan.darkgoldenrod.darkgray.darkgreen.darkgrey.darkkhaki.darkmagenta.darkolivegreen.darkorange.darkorchid.darkred.darksalmon.darkseagreen.darkslateblue.darkslategray.darkslategrey.darkturquoise.darkviolet.deeppink.deepskyblue.dimgray.dimgrey.dodgerblue.firebrick.floralwhite.forestgreen.fuchsia.gainsboro.ghostwhite.gold.goldenrod.gray.aliceblue.greenyellow.black.honeydew.hotpink.indianred.indigo.ivory.khaki.lavender.lavenderblush.lawngreen.lemonchiffon.lightblue.lightcoral.lightcyan.lightgoldenrodyellow.lightgray.lightgreen.lightgrey.lightpink.lightsalmon.lightseagreen.lightskyblue.lightslategray.lightslategrey.lightsteelblue.lightyellow.lime.limegreen.linen.magenta.maroon.mediumaquamarine.mediumblue.mediumorchid.mediumpurple.mediumseagreen.mediumslateblue.mediumspringgreen.mediumturquoise.mediumvioletred.midnightblue.mintcream.mistyrose.moccasin.navajowhite.navy.oldlace.olive.olivedrab.orange.orangered.orchid.palegoldenrod.palegreen.paleturquoise.palevioletred.papayawhip.peachpuff.peru.pink.plum.powderblue.purple.rebeccapurple.yellow.rosybrown.royalblue.saddlebrown.salmon.sandybrown.seagreen.seashell.sienna.silver.skyblue.slateblue.slategray.slategrey.snow.springgreen.steelblue.tan.teal.thistle.tomato.turquoise.violet.wheat.white.whitesmoke.yellow.yellowgreen`.split(`.`);e.forEach((e,r)=>{let i=n[r],a=new c.Group,{pcdPoints:o,category:s,name:l,center:u,box:d,obbBox:f,contour:p,nline:m,coordinatesByWidth:h,coordinatesByArea:v}=e,y=o.length;if(s!=`door`)return;s===`balcony railing`&&console.log(`pcdObj`,e);let b=new c.Box3(new c.Vector3(d.min.x,d.min.y,d.min.z),new c.Vector3(d.max.x,d.max.y,d.max.z)),x=new c.Box3Helper(b,i);a.add(x);let S=new _({value:l,position:new c.Vector3(u.x,u.y,u.z)}).getLabel();a.add(S);let C=new c.InstancedMesh(new c.SphereGeometry(.007,8,8),new c.MeshBasicMaterial({color:i}),y);a.add(C);let w=new c.Matrix4,T=0;if(o.forEach(e=>{w.setPosition(e.x,e.y,e.z),C.setMatrixAt(T,w),T++}),C.instanceMatrix.needsUpdate=!0,m){let e=new c.BufferGeometry;e.setFromPoints([new c.Vector3(m.start.x,m.start.y,m.start.z),new c.Vector3(m.end.x,m.end.y,m.end.z)]);let t=new c.LineBasicMaterial({color:`red`,depthTest:!1,depthWrite:!1});new c.Line(e,t);let n=new c.Vector3().copy(new c.Vector3(m.start.x,m.start.y,m.start.z)).add(new c.Vector3(m.end.x,m.end.y,m.end.z)).multiplyScalar(.5),r=new c.Vector3(m.centerDir.x,m.centerDir.y,m.centerDir.z);new c.ArrowHelper(r,n,1,`red`)}if(p){let e=[];p.forEach((t,n)=>{e.push(new c.Vector3(t.x,t.y,t.z))});let n=new c.BufferGeometry;n.setFromPoints(e);let r=new c.LineBasicMaterial({color:`yellow`,depthTest:!1,depthWrite:!1}),i=new c.Line(n,r);t.add(i)}if(v){let e=[],n=[];v.coordinates.forEach((t,r)=>{e.push(new c.Vector3(t.x,t.y,t.z)),n.length<4&&n.push(new c.Vector3(t.x,t.y,t.z))});let r=new c.BufferGeometry;r.setFromPoints(e);let i=new c.LineBasicMaterial({color:`green`,depthTest:!1,depthWrite:!1}),a=new c.Line(r,i);t.add(a);let o=g(n.map(e=>new c.Vector3(e.x,e.y,v.heightData.minZ)),n.map(e=>new c.Vector3(e.x,e.y,e.z)));t.add(o);let[s,l,u]=e,d=s.distanceTo(l),f=new c.Vector3().copy(s).add(l).multiplyScalar(.5);new _({value:String(d.toFixed(3)),position:f}).getLabel();let p=l.distanceTo(u),m=new c.Vector3().copy(l).add(u).multiplyScalar(.5);new _({value:String(p.toFixed(3)),position:m}).getLabel()}if(f){let e=new c.BoxGeometry(f.size.x,f.size.y,f.size.z),n=new c.LineBasicMaterial({color:i,linewidth:2}),r=new c.EdgesGeometry(e),a=new c.LineSegments(r,n),o=new Vector3(f.center.x,f.center.y,f.center.z),s=new Quaternion(f.rotation.x,f.rotation.y,f.rotation.z,f.rotation.w);a.position.copy(o),a.quaternion.copy(s),t.add(a)}})}return t},y=(e,t,n,r)=>{let i=e,a=[];if(r&&r.length&&(a=r.map(e=>e.clone())),!i||i.length<4)return{};let o=new c.Vector3().subVectors(i[1],i[0]).normalize(),s=new c.Vector3().subVectors(i[3],i[0]).normalize(),l=o.clone().cross(s).normalize(),u=1/0,d=-1/0,f=1/0,p=-1/0;for(let e of i){let t=e.dot(o),n=e.dot(s);t<u&&(u=t),t>d&&(d=t),n<f&&(f=n),n>p&&(p=n)}let m=i[0].dot(l),g=o.clone().multiplyScalar(u).add(s.clone().multiplyScalar(f)).add(l.clone().multiplyScalar(m)),_=o.clone().multiplyScalar(d).add(s.clone().multiplyScalar(f)).add(l.clone().multiplyScalar(m)),v=o.clone().multiplyScalar(d).add(s.clone().multiplyScalar(p)).add(l.clone().multiplyScalar(m)),y=o.clone().multiplyScalar(u).add(s.clone().multiplyScalar(p)).add(l.clone().multiplyScalar(m)),b=[g,_,v,y],x=h(new c.Vector3().subVectors(g,_),new c.Vector3(0,0,1));x&&x.parallel&&(b=[g,y,v,_]);let S=0;if(r.length>0){let e=new c.Vector3().subVectors(b[1],b[0]).normalize(),r=e.clone().negate(),i=new c.Vector3().subVectors(b[3],b[0]).normalize(),o=i.clone().negate(),s=b[0].distanceTo(b[1]),l=b[0].distanceTo(b[3]),u=s/t,d=l/n,f=.001;for(let r=0;r<u;r++){let r=b[0],o=0;for(let s=0;s<d;s++){let s=r.clone().addScaledVector(e,t),l=s.clone().addScaledVector(i,n),u=r.clone().addScaledVector(i,n),d=new c.Vector3().add(r).add(s).add(l).add(u).multiplyScalar(.25);for(let e=0;e<a.length;e++)if(a[e].distanceTo(d)<f){o++;break}r=u}if(o>=d/2)break;S+=o,b[0]=b[0].addScaledVector(e,t),b[3]=b[3].addScaledVector(e,t)}s=b[0].distanceTo(b[1]),u=s/t;for(let e=0;e<u;e++){let e=b[1],o=0;for(let s=0;s<d;s++){let s=e.clone().addScaledVector(r,t),l=s.clone().addScaledVector(i,n),u=e.clone().addScaledVector(i,n),d=new c.Vector3().add(e).add(s).add(l).add(u).multiplyScalar(.25);for(let e=0;e<a.length;e++)if(a[e].distanceTo(d)<f){o++;break}e=u}if(o>=d/2)break;S+=o,b[1]=b[1].addScaledVector(r,t),b[2]=b[2].addScaledVector(r,t)}s=b[0].distanceTo(b[1]),u=s/t;for(let r=0;r<d;r++){let r=b[3],i=0;for(let s=0;s<u;s++){let s=r.clone().addScaledVector(e,t),l=s.clone().addScaledVector(o,n),u=r.clone().addScaledVector(o,n),d=new c.Vector3().add(r).add(s).add(l).add(u).multiplyScalar(.25);for(let e=0;e<a.length;e++)if(a[e].distanceTo(d)<f){i++;break}r=s}if(i>=u/2)break;S+=i,b[2]=b[2].addScaledVector(o,n),b[3]=b[3].addScaledVector(o,n)}}let C=b[0].distanceTo(b[1])*b[0].distanceTo(b[3]);return{facePoints:b,boxArea:C,totalInPlaneNum:S}},b=(e,t,n,r,i=c.MathUtils.degToRad(10))=>{let a=new c.Vector3,o=new c.Vector3;a.subVectors(new c.Vector3(t.x,t.y,0),new c.Vector3(e.x,e.y,0)).normalize(),o.subVectors(new c.Vector3(r.x,r.y,0),new c.Vector3(n.x,n.y,0)).normalize();let s=a.dot(o),l=Math.acos(Math.min(Math.abs(s),1)),u=1-Math.min(l/i,1),d=s>0?`same`:`opposite`;return{parallelism:u,angle:l,angleDeg:c.MathUtils.radToDeg(l),isParallel:l<i,direction:d,vectors:{v1:a,v2:o}}},x=(e,t,n,r)=>{let i=new c.Vector3().subVectors(t,e),a=new c.Vector3().subVectors(r,n),o=i.length(),s=a.length();if(o===0||s===0)return{rate:NaN,angle:NaN,isPerpendicular:!1,isParallel:!1};i.normalize(),a.normalize();let l=Math.abs(i.dot(a)),u=Math.min(1,Math.max(0,l)),d=Math.acos(u),f=c.MathUtils.radToDeg(d),p=1-u;return{rate:p,percent:`${(p*100).toFixed(2)}%`,angle:f,isPerpendicular:Math.abs(f-90)<.01,isParallel:f<.01}},S=(e,t,n=!1)=>new c.Vector3((e.x+t.x)/2,(e.y+t.y)/2,n?0:(e.z+t.z)/2),C=(e,t,n,r,i={})=>{let{parallelAngleMax:a=c.MathUtils.degToRad(10),collinearAngleMax:o=c.MathUtils.degToRad(7),distanceThreshold:s=2,overlapThreshold:l=.01}=i,u=new c.Vector3().subVectors(t,e),d=new c.Vector3().subVectors(r,n),f=u.clone().normalize(),p=d.clone().normalize(),m=c.MathUtils.clamp(f.dot(p),-1,1),h=Math.acos(Math.abs(m)),g=c.MathUtils.radToDeg(h);if(h>a)return{type:`not_parallel`,angleDeg:g,avgPerpendicularDistance:null,maxPerpendicularDistance:null,gap:null,overlap:null,closestDistance:null};m<0&&p.negate();let _=new c.Vector3().addVectors(f,p).normalize(),v=new c.Vector3().addVectors(e,t).multiplyScalar(.5),y=w(n,v,f),b=w(r,v,f),x=(y+b)/2,S=Math.max(y,b),C=e.dot(_),E=t.dot(_),D=n.dot(_),O=r.dot(_),k=Math.min(C,E),A=Math.max(C,E),j=Math.min(D,O),M=Math.max(D,O),N=A-k,P=M-j,F=Math.min(N,P),I=Math.min(A,M)-Math.max(k,j),L=I<0?-I:0,R=I>0?I:0,z=F>0?R/F:0,B=T(e,t,n,r),V=h<=o&&x<s,H;return V?H=z>.5?`collinear_overlap`:`collinear_gap`:h<=a&&(H=`parallel_offset`),{type:H,angleDeg:g,avgPerpendicularDistance:x,maxPerpendicularDistance:S,gap:L,overlap:R,closestDistance:B}};function w(e,t,n){let r=new c.Vector3().subVectors(e,t),i=r.dot(n),a=n.clone().multiplyScalar(i);return new c.Vector3().subVectors(r,a).length()}function T(e,t,n,r){let i=new c.Vector3().subVectors(t,e),a=new c.Vector3().subVectors(r,n),o=new c.Vector3().subVectors(n,e),s=i.dot(i),l=i.dot(a),u=a.dot(a),d=i.dot(o),f=a.dot(o),p=s*u-l*l,m,h;p<1e-10?(m=0,h=u===0?0:f/u):(m=(l*f-u*d)/p,h=(s*f-l*d)/p),m=c.MathUtils.clamp(m,0,1),h=u===0?0:c.MathUtils.clamp((l*m+f)/u,0,1),m=s===0?0:c.MathUtils.clamp((l*h-d)/s,0,1);let g=e.clone().add(i.clone().multiplyScalar(m)),_=n.clone().add(a.clone().multiplyScalar(h));return g.distanceTo(_)}var E=(e,t,n,r=.05,i=.02)=>{let[a,o,s,l]=e,u=new c.Vector3().subVectors(o,a),d=new c.Vector3().subVectors(l,a),f=new c.Vector3().crossVectors(u,d).normalize(),p=u.length(),m=d.length(),h=u.clone().normalize(),g=d.clone().normalize(),_=n-r,v=n+r,y=new c.Vector3,b=[],x=Math.ceil(p/i),S=Math.ceil(m/i),C=p/x,w=m/S,T=new Uint8Array(x*S);for(let e=0;e<t.length;e++){y.subVectors(t[e],a);let n=Math.abs(y.dot(f));if(n<_||n>v)continue;let r=y.dot(h),i=y.dot(g);if(r<0||r>p||i<0||i>m)continue;b.push(new c.Vector3(t[e].x,t[e].y,t[e].z));let o=Math.min(Math.floor(r/C),x-1),s=Math.min(Math.floor(i/w),S-1);T[s*x+o]=1}let E=0;for(let e=0;e<T.length;e++)T[e]&&E++;let D=x*S,O=p*m,k=E==0||D==0?0:E/D*O;return{filteredPoints:b,coveragePercent:E==0||D==0?0:E/D*100,quadArea:O,coveredArea:k}},D=(e,t)=>Math.round(Math.max(5,Math.min(350,t*Math.sqrt(e)))),O=e=>{let t=new c.Vector3(e.start.x,e.start.y,e.start.z),n=new c.Vector3(e.end.x,e.end.y,e.end.z),r=t.distanceTo(n),i=e.rooftopPz-e.start.z,a=D(r,100),o=D(i,70),s=i/o,l=new Map;for(let e=0;e<o;e++)l.set((e+1)*s,{count:0,minZ:1/0,maxZ:-1/0});let u=new c.Line3(t,n),d=e.originalPoints,f=new c.Vector3;for(let e=0;e<d.length;e++){let t=d[e],n=u.closestPointToPoint(t,!0,f).distanceTo(t),r=Math.floor(n/s);if(r<o){let e=l.get((r+1)*s);e.count++,t.z<e.minZ&&(e.minZ=t.z),t.z>e.maxZ&&(e.maxZ=t.z)}}let p,m,h=a/2;for(let[e,t]of l)if(t.count>=h){p=t.minZ;break}for(let[e,t]of[...l.entries()].reverse())if(t.count>=h){m=t.maxZ;break}return{minZ:p,maxZ:m}},k=(e,t,n)=>{let r=new c.Plane,i=new c.Vector3(e.x,e.y,e.z),a=new c.Vector3(t.x,t.y,t.z),o=new c.Vector3(t.x,t.y,t.z+5);r.setFromCoplanarPoints(i,a,o);let s=[];n.forEach(e=>{let t=new c.Vector3;r.projectPoint(e,t),t.z=i.z,s.push(t)});let l=new c.Vector3,u=new c.Vector3;r.projectPoint(i,l),r.projectPoint(a,u);let d=new c.Line3(l,u),f=[],p=[],m=new c.Vector3;return s.forEach((e,t)=>{d.closestPointToPoint(e,!0,m).distanceTo(e)<.001?f.push(n[t]):p.push(e)}),{newOriginalPoints:f,newProjectPoints:p}},A=e=>{function t(e,t,n,r,i,a,o,s,l){let u=new c.Line3(e,t),d=new Map,f=new c.Vector3,p=e.clone(),m=t.clone(),h=0;for(;;){let g=new Map,_=0;for(let e=0;e<l.checkResults.length;e++)if(!(l.checkResults[e].allCenterPoints.length/s*100<1))for(let t=0;t<l.checkResults[e].allCenterPoints.length;t++){if(d.has(e)&&d.get(e).removePtsIndex.includes(t))continue;let n=l.checkResults[e].allCenterPoints[t];n=new c.Vector3(n.x,n.y,n.z);let r=u.closestPointToPoint(n,!0,f),a=n.distanceTo(r);Math.abs(a-i)<.01&&(_++,g.has(e)?g.get(e).removePtsIndex.push(t):g.set(e,{index:e,removePtsIndex:[t]}))}if(h==o){e=p,t=m;break}if(_<a/20*19||_==0){g.clear(),o-h<5&&(e=p,t=m,d.clear());break}else{h++;for(let[e,t]of g)d.has(e)?d.get(e).removePtsIndex.push(...t.removePtsIndex):d.set(e,t);g.clear(),e=e.addScaledVector(n,r),t=t.addScaledVector(n,r),u.set(e,t)}}if(d.size>0&&h!=o)for(let[e,t]of d){let n=[],r=[];for(let i=0;i<l.checkResults[e].allCenterPoints.length;i++)t.removePtsIndex.includes(i)||(n.push(l.checkResults[e].allCenterPoints[i]),r.push(l.checkResults[e].originalVertices[i]));l.checkResults[e].allCenterPoints=n,l.checkResults[e].originalVertices=r}}for(let n=0;n<e.length;n++){let r=e[n];if(r.length<.5)continue;let i=-1,a=-1;if(!r.checkResults||r.checkResults.length<=0)continue;if(r.checkResults[0].originalVertices&&r.checkResults[0].originalVertices.length>0){let[e,t,n,o]=r.checkResults[0].originalVertices[0];e=new c.Vector3(e.x,e.y,e.z),t=new c.Vector3(t.x,t.y,t.z),n=new c.Vector3(n.x,n.y,n.z),o=new c.Vector3(o.x,o.y,o.z),i=e.distanceTo(t),a=e.distanceTo(o);let s=h(new c.Vector3().subVectors(e,t).normalize(),new c.Vector3(0,0,1));s&&s.parallel&&([a,i]=[i,a])}if(i<0||a<0)continue;let o=new c.Vector3(r.start.x,r.start.y,r.start.z),s=new c.Vector3(r.end.x,r.end.y,r.end.z),l=o.distanceTo(s),u=r.rooftopPz-r.start.z,d=Math.ceil(l/i),f=Math.ceil(u/a),p=o.clone(),m=s.clone(),g=m.clone().add(new c.Vector3(0,0,u)),_=p.clone().add(new c.Vector3(0,0,u)),v=new c.Vector3().subVectors(m,p).normalize(),y=v.clone().negate();new c.Vector3().subVectors(g,p).normalize().clone().negate();let b=i/2;a/2;let x=l*u/(i*a);t(p,_,v,i,b,f,d,x,r),t(m,g,y,i,b,f,d,x,r),r.start=p,r.end=m;let{newOriginalPoints:S}=k(p,m,r.originalPoints);r.originalPoints=S}},j=(e,t)=>new c.Vector3(e.x,e.y,0).distanceTo(new c.Vector3(t.x,t.y,0)),M=class{results;clusterResults;constructor(){this.results=null}initLimits1(e,t){let n,r;return n=e<.3?5:e<.5?8:e<.8?10:e<1?15:e<2?30:e<3?35:e<4?40:e<5?45:e<6?50:e<7?60:e<8.5?70:e<10?80:e<15?100:e<20?120:e<25?140:e<30?160:e<35?180:e<40?200:e<45?230:e<50?250:350,t<1?r=30:t<2?r=40:t<2.5?r=45:t<3?r=60:t<3.6?r=70:t<4&&(r=80),{horizontalSubdivisions:n,verticalSubdivisions:r}}initLimits2(e,t){return{horizontalSubdivisions:Math.round(Math.max(5,Math.min(350,e/.08))),verticalSubdivisions:Math.round(Math.max(30,Math.min(80,t/.06)))}}initLimits4(e,t){return{horizontalSubdivisions:Math.round(Math.max(5,Math.min(350,65*Math.sqrt(e))))*1,verticalSubdivisions:Math.round(Math.max(30,Math.min(80,15*Math.sqrt(t))))*1}}initLimits3(e,t){let n;return n=e<2?.06:e<10?.1:.15,{horizontalSubdivisions:Math.round(Math.max(5,Math.min(350,e/n))),verticalSubdivisions:Math.round(Math.max(30,Math.min(80,t/.06)))}}createWallPlaneMeshFromFourPoints(e,t,n,r,i,a,o=16777215){let s=[e,t,n,r];for(let e=0;e<s.length;e++)this.isValidVector3(s[e])||(console.error(`点${e+1}包含无效数据:`,s[e]),s[e]=new c.Vector3(e,0,e));let l=new c.BufferGeometry,u=new Float32Array([s[0].x,s[0].y,s[0].z,s[1].x,s[1].y,s[1].z,s[2].x,s[2].y,s[2].z,s[3].x,s[3].y,s[3].z]),d=[0,1,2,0,2,3],f=new Float32Array([0,0,1,0,1,1,0,1]);l.setAttribute(`position`,new c.BufferAttribute(u,3)),l.setAttribute(`uv`,new c.BufferAttribute(f,2)),l.setIndex(d),l.computeVertexNormals();let p=new c.Mesh(l,i),m=new c.EdgesGeometry(l),h=new c.LineSegments(m,new c.LineBasicMaterial({color:o}));return a.add(h),p}createSubdividedWallPlane(e,t,n,r,i){let a=[];if(!e||!e.start||!e.end||isNaN(e.rooftopPz))return a;let o=new c.Vector3(e.start.x,e.start.y,e.start.z),s=new c.Vector3(e.end.x,e.end.y,e.end.z),l=[o,s,new c.Vector3(e.end.x,e.end.y,e.rooftopPz),new c.Vector3(e.start.x,e.start.y,e.rooftopPz)];for(let e=0;e<l.length;e++)if(!this.isValidVector3(l[e]))return console.error(`墙体基础点${e}无效:`,l[e]),a;let u=o.distanceTo(s),d=Math.abs(e.rooftopPz-e.start.z);if(u===0||d===0)return console.error(`墙体宽度或高度为0`),console.log(`line.rooftopPz`,e.rooftopPz),console.log(`bottomLeft`,o),console.log(`bottomRight`,s),console.log(e.length),a;let f=u/n,p=d/t,m=new c.Vector3().subVectors(s,o).normalize(),h=new c.Vector3(0,0,1);for(let e=0;e<t;e++)for(let t=0;t<n;t++)try{let n={points:[new c.Vector3().copy(o).add(m.clone().multiplyScalar(t*f)).add(h.clone().multiplyScalar(e*p)),new c.Vector3().copy(o).add(m.clone().multiplyScalar((t+1)*f)).add(h.clone().multiplyScalar(e*p)),new c.Vector3().copy(o).add(m.clone().multiplyScalar((t+1)*f)).add(h.clone().multiplyScalar((e+1)*p)),new c.Vector3().copy(o).add(m.clone().multiplyScalar(t*f)).add(h.clone().multiplyScalar((e+1)*p))],userData:{}};n.userData={type:`wallSegment`,row:e,col:t,originalMaterial:r,width:f,height:p},a.push(n)}catch(n){console.error(`创建墙体小平面(${e}, ${t})时出错:`,n)}return a}calculateMinRequiredPoints(e){let t=Math.floor(Math.log1p(e)*.5);return Math.max(1,t)}calculateExpectedDensity(e){return e<1?10:e<4?5:2}evaluateSegmentByArea(e,t,n,r){if(t===0)return!1;let i=e.area,a=e.minRequiredPoints,o=e.expectedDensity,s=t/Math.max(i,.001),c=t>=a,l=s>=o*.3,u=!0;return i<.5&&t>0&&(u=this.checkSmallAreaDistribution(e,n,r)),c&&l&&u}checkSmallAreaDistribution(e,t,n){if(t.length<2)return!0;let r=t.map(e=>n[e]);e.bounds;let i=e.center,a=0,o=Math.sqrt(e.area)*.3;return r.forEach(e=>{e.distanceTo(i)<=o&&a++}),a/r.length>=.3}getFailureReason(e,t){let n=e.area,r=e.minRequiredPoints,i=e.expectedDensity,a=t/Math.max(n,.001),o=[];return t===0?o.push(`无点`):(t<r&&o.push(`点数不足: ${t}/${r}`),a<i*.3&&o.push(`密度不足: ${a.toFixed(2)}/${i.toFixed(2)}`),n<.5&&t>0&&(this.checkSmallAreaDistribution(e,[],[])||o.push(`小平面分布不均`))),o.length>0?o.join(`, `):`未知原因`}calculateMedian(e){let t=[...e].sort((e,t)=>e-t),n=Math.floor(t.length/2);return t.length%2==0?(t[n-1]+t[n])/2:t[n]}isPointInBoundingBox(e,t,n=.1){return e.x>=t.min[0]-n&&e.x<=t.max[0]+n&&e.y>=t.min[1]-n&&e.y<=t.max[1]+n&&e.z>=t.min[2]-n&&e.z<=t.max[2]+n}getSegmentPlane(e){let t=e.geometry;if(!t||!t.attributes.position)return null;let n=t.attributes.position;if(n.count<3)return null;let r=new c.Vector3().fromBufferAttribute(n,0),i=new c.Vector3().fromBufferAttribute(n,1),a=new c.Vector3().fromBufferAttribute(n,2);r.applyMatrix4(e.matrixWorld),i.applyMatrix4(e.matrixWorld),a.applyMatrix4(e.matrixWorld);let o=new c.Plane;return o.setFromCoplanarPoints(r,i,a),o}analyzePointDistribution(e,t,n=.01){this.results={totalSegments:0,segmentsWithPoints:0,segmentsWithoutPoints:0,missingSegments:[],pointsInSegments:new Map,segmentsByPoint:new Map,segmentDetails:new Map,tolerance:n,analysisTime:null};let r=performance.now(),i=[];e.forEach(e=>{if(e.userData&&e.userData.type===`wallSegment`){let t=this.getWorldVerticesNew(e.points);if(t.length<4)return;let n=new c.Vector3().add(t[0]).add(t[1]).add(t[2]).add(t[3]).multiplyScalar(.25);e.userData.center=n,e.vertices=t,i.push(e)}}),this.results.totalSegments=i.length;for(let e of i){if(e.vertices.length<4)continue;let r=e.userData.center,i=`row${e.userData.row}_col${e.userData.col}`,[a,o,s,l]=e.vertices,u=new c.Vector3().subVectors(o,a).normalize(),d=new c.Vector3().subVectors(l,a).normalize(),f=a.distanceTo(o),p=a.distanceTo(l),m=new c.Plane;m.setFromCoplanarPoints(a,o,l);let h=f/2,g=p/2;for(let e=0;e<t.length;e++)this.isPointInWallSegmentOptimized(t[e],r,m,a,u,d,f,p,h,g,n)&&(this.results.pointsInSegments.has(i)||this.results.pointsInSegments.set(i,[]),this.results.pointsInSegments.get(i).push(e))}return i.forEach(e=>{let t=`row${e.userData.row}_col${e.userData.col}`,n=this.results.pointsInSegments.has(t)?this.results.pointsInSegments.get(t).length:0;n>=1?this.results.segmentsWithPoints++:(this.results.segmentsWithoutPoints++,this.results.missingSegments.push({row:e.userData.row,col:e.userData.col,width:e.userData.width,height:e.userData.height,segmentKey:t,vertices:e.vertices,center:e.userData.center,mesh:e,pointCount:n,hasPointsButNotEnough:n>0&&n<1}))}),this.results.analysisTime=performance.now()-r,this.printAnalysisResults(),this.results}isPointInWallSegmentOptimized(e,t,n,r,i,a,o,s,l,u,d){if(Math.abs(n.distanceToPoint(e))>d)return!1;let f=new c.Vector3().subVectors(e,t),p=Math.abs(f.dot(i)),m=Math.abs(f.dot(a));if(p>l&&m>u)return!1;let h=new c.Vector3().subVectors(e,r),g=h.dot(i),_=h.dot(a),v=g>=-d&&g<=o+d,y=_>=-d&&_<=s+d;return v&&y}isPointInWallSegment(e,t,n,r=.1){if(!e.geometry)return!1;try{let i=this.getWorldVertices(e);return i.length<4?!1:this.isPointInRotatedWall(i,t,n,r)}catch(e){return console.error(`判断点是否在墙体平面内时出错:`,e),!1}}isPointInRotatedWall(e,t,n,r=.1){let[i,a,o,s]=e,l=new c.Plane;if(l.setFromCoplanarPoints(i,a,s),Math.abs(l.distanceToPoint(t))>r)return!1;let u=new c.Vector3().subVectors(a,i).normalize(),d=new c.Vector3().subVectors(s,i).normalize();new c.Vector3().crossVectors(u,d).normalize();let f=new c.Vector3().subVectors(t,i),p=f.dot(u),m=f.dot(d),h=i.distanceTo(a),g=i.distanceTo(s),_=p>=-r&&p<=h+r,v=m>=-r&&m<=g+r;return _&&v&&n.distanceTo(t)<h/2&&n.distanceTo(t)<g/2}getWorldVertices(e){if(e._cachedWorldVertices&&e.matrixWorld.equals(e._cachedMatrixWorld))return e._cachedWorldVertices;let t=e.geometry.getAttribute(`position`),n=[];for(let r=0;r<Math.min(4,t.count);r++){let i=new c.Vector3;i.fromBufferAttribute(t,r),i.applyMatrix4(e.matrixWorld),n.push(i)}return e._cachedWorldVertices=n,e._cachedMatrixWorld=e.matrixWorld.clone(),n}getWorldVerticesNew(e){let t=[];for(let n=0;n<e.length;n++){let r=e[n].clone();t.push(r)}return t}getWallSegmentBounds(e){if(!e.geometry)return null;let t=e.geometry.getAttribute(`position`),n=new c.Box3;for(let r=0;r<t.count;r++){let i=new c.Vector3;i.fromBufferAttribute(t,r),i.applyMatrix4(e.matrixWorld),n.expandByPoint(i)}return{min:n.min.toArray(),max:n.max.toArray(),center:n.getCenter(new c.Vector3).toArray(),size:n.getSize(new c.Vector3).toArray()}}highlightEmptySegments(e,t,n,r=1){if(!this.results){console.warn(`请先运行分析函数`);return}let i=new Map,a=new Map;this.results.missingSegments.forEach(e=>{i.set(e.segmentKey,e),a.has(e.row)||a.set(e.row,new Map),a.get(e.row).set(e.col,{segmentKey:e.segmentKey,isEmpty:!0,visited:!1})}),e.forEach(e=>{if(e.userData&&e.userData.type===`wallSegment`){let t=`row${e.userData.row}_col${e.userData.col}`,n=i.has(t);a.has(e.userData.row)||a.set(e.userData.row,new Map),a.get(e.userData.row).has(e.userData.col)||a.get(e.userData.row).set(e.userData.col,{segmentKey:t,isEmpty:n,visited:!1})}});let o=[],s=[[0,1],[1,0],[0,-1],[-1,0]];for(let[e,t]of i){let e=a.get(t.row)?.get(t.col);if(e&&!e.visited){let n=[],c=[[t.row,t.col]];for(e.visited=!0;c.length>0;){let[e,t]=c.shift(),r=`row${e}_col${t}`;i.has(r)&&n.push(i.get(r));for(let[n,r]of s){let i=e+n,o=t+r;if(a.has(i)&&a.get(i).has(o)){let e=a.get(i).get(o);!e.visited&&e.isEmpty&&(e.visited=!0,c.push([i,o]))}}}n.length>=r&&o.push(n)}}let c=this.createClusterMaterials(o.length),l=[];return this.clusterResults={clusters:o.sort((e,t)=>t.length-e.length),highlightedMeshes:l,clusterMaterials:c,totalClusters:o.length,totalHighlighted:l.length},this.clusterResults}calculateBoundaryRegularity(e,t,n,r,i){let a=0,o=0,s=t,c=!0;for(let t=r;t<=i;t++)e.has(s)&&e.get(s).has(t)?o++:c=!1;c&&a++;let l=n,u=!0;for(let t=r;t<=i;t++)e.has(l)&&e.get(l).has(t)?o++:u=!1;u&&a++;let d=r,f=!0;for(let r=t;r<=n;r++)e.has(r)&&e.get(r).has(d)?o++:f=!1;f&&a++;let p=i,m=!0;for(let r=t;r<=n;r++)e.has(r)&&e.get(r).has(p)?o++:m=!1;return m&&a++,a/4}calculateClusterQualityScore(e){let t=1/Math.max(1,e.aspectRatio);return e.compactness*.4+e.boundaryRegularity*.4+t*.2}calculateStrictClusterBounds(e){if(e.length===0)return{width:0,height:0};let t=[];if(e.forEach(e=>{e.vertices&&e.vertices.length>0&&t.push(...e.vertices)}),t.length===0){let t=[...new Set(e.map(e=>e.row))];return{width:[...new Set(e.map(e=>e.col))].length*(e[0].bounds?.width||1),height:t.length*(e[0].bounds?.height||1)}}let n=Math.min(...t.map(e=>e.x)),r=Math.max(...t.map(e=>e.x)),i=Math.min(...t.map(e=>e.y)),a=Math.max(...t.map(e=>e.y)),o=Math.min(...t.map(e=>e.z)),s=Math.max(...t.map(e=>e.z));return{width:Math.abs(r-n),height:Math.max(Math.abs(a-i),Math.abs(s-o))}}createClusterMaterials(e){let t=[],n=[16711935,65535,16753920,16738740,9055202,3329330,`lavenderblush`,`lawngreen`,`lemonchiffon`,`lightblue`,`lightcoral`];for(let r=0;r<e;r++){let e=n[r%n.length],i=new c.MeshBasicMaterial({color:e,transparent:!0,opacity:.6,side:c.DoubleSide});t.push(i)}return t}clearHighlights(e){this.clusterResults&&this.clusterResults.highlightedMeshes&&(this.clusterResults.highlightedMeshes.forEach(t=>{e.remove(t),t.material&&t.material.dispose()}),this.results&&this.results.segmentDetails&&this.results.segmentDetails.forEach((e,t)=>{e.mesh&&e.mesh.userData.originalMaterial&&(e.mesh.material=e.mesh.userData.originalMaterial)}),this.clusterResults.highlightedMeshes=[])}getClusterInfo(){return this.clusterResults?{totalClusters:this.clusterResults.totalClusters,clusterSizes:this.clusterResults.clusters.map((e,t)=>({clusterIndex:t,size:e.length,segments:e.map(e=>e.segmentKey)})),largestCluster:Math.max(...this.clusterResults.clusters.map(e=>e.length)),smallestCluster:Math.min(...this.clusterResults.clusters.map(e=>e.length))}:null}restoreOriginalMaterials(e){let t=0;e.traverse(e=>{e.userData&&e.userData.originalMaterial&&(e.material=e.userData.originalMaterial,delete e.userData.originalMaterial,t++)}),console.log(`恢复了 ${t} 个平面的原始材质`)}createDebugHelpers(e){let t=new c.Group;return e.traverse(e=>{if(e.userData&&e.userData.type===`wallSegment`){let n=this.getWorldVertices(e);if(n.length<4)return;let r=new c.Vector3().add(n[0]).add(n[1]).add(n[2]).add(n[3]).multiplyScalar(.25),i=new c.Vector3().subVectors(n[1],n[0]).normalize(),a=new c.Vector3().subVectors(n[3],n[0]).normalize(),o=new c.Vector3().crossVectors(i,a).normalize(),s=new c.ArrowHelper(o,r,.3,16711680),l=new c.ArrowHelper(i,r,.2,65280),u=new c.ArrowHelper(a,r,.2,255);t.add(s),t.add(l),t.add(u)}}),t}printAnalysisResults(){if(!this.results){console.warn(`没有可用的分析结果`);return}this.results.missingSegments.length>0&&this.results.missingSegments.forEach(e=>{}),this.results.pointsInSegments.forEach((e,t)=>{})}isValidVector3(e){return e&&!isNaN(e.x)&&!isNaN(e.y)&&!isNaN(e.z)&&isFinite(e.x)&&isFinite(e.y)&&isFinite(e.z)}mergeWallGroupMeshes(e){let t=[];if(e.traverse(e=>{e.isMesh&&e.geometry&&t.push(e)}),t.length===0)return null;try{let e=t.map(e=>e.geometry);e.forEach((e,n)=>{e.applyMatrix4(t[n].matrixWorld)});let n=l(e),r=t[0].material.clone(),i=new c.Mesh(n,r);return i.name=`mergedWalls`,i}catch(e){return console.error(`合并网格时出错:`,e),null}}},N=e=>{let t={},n=[],r=[],i=[];return e[0].width,e[0].height,e.forEach(e=>{r.push(...e.vertices),i.push(e.center),n.push(e.vertices)}),t.originalVertices=n,t.allCenterPoints=i,t},P=e=>{let t=new c.Vector3().copy(e.start).add(e.end).multiplyScalar(.5),n=new c.Vector3().subVectors(new c.Vector3(t.x,t.y,t.z+1),t).normalize(),r=new c.Vector3(e.direction.x,e.direction.y,0);r.applyAxisAngle(n,Math.PI/2);let i=new c.Vector3().copy(t).add(r.clone().multiplyScalar(e.length)),a=new c.Vector3().copy(t).sub(r.clone().multiplyScalar(e.length));return t.z=e.start.z,i.z=e.start.z,a.z=e.start.z,{start:i,end:a,center:t,direction:r}},F=(e,t,n=!0)=>{let r=[];if(!e)return[];for(let i=0;i<e.length;i++){let a=e[i];if(!a)continue;let o=new c.Vector3(a.start.x,a.start.y,a.start.z);o.applyEuler(new c.Euler(-Math.PI/2,0,0));let s=new c.Vector3(a.end.x,a.end.y,a.end.z);if(s.applyEuler(new c.Euler(-Math.PI/2,0,0)),t.add(p(String(i),new c.Vector3((o.x+s.x)/2,(o.y+s.y)/2,(o.z+s.z)/2))),(a.isBayWindow||a.isWindow)&&n)continue;let l=new c.Plane().setFromCoplanarPoints(new c.Vector3(a.start.x,a.start.y,a.start.z),new c.Vector3(a.end.x,a.end.y,a.end.z),new c.Vector3(a.start.x,a.start.y,a.rooftopPz)),u=[];a.originalPoints.forEach((e,t)=>{let n=new c.Vector3;l.projectPoint(e,n),u.push(n)});let d=new c.Vector3(a.start.x,a.start.y,a.start.z).distanceTo(new c.Vector3(a.end.x,a.end.y,a.end.z)),f=Math.abs(a.rooftopPz-a.start.z),m=new c.MeshBasicMaterial({color:`red`,side:c.DoubleSide,wireframe:!1,transparent:!0,opacity:.5}),h=new c.Group,g=new M,{horizontalSubdivisions:_,verticalSubdivisions:v}=g.initLimits4(d,f),y=g.createSubdividedWallPlane(a,v,_,m,h);g.analyzePointDistribution(y,u,.01);let b=new c.MeshBasicMaterial({color:`black`,transparent:!0,opacity:.8,side:c.DoubleSide}),x=g.highlightEmptySegments(y,b,h),S=[];for(let e of x.clusters){let t=N(e);t&&S.push(t)}let C=P(a).direction;a.checkResults=S,a.verticalDirection=C,a.originaIndex=i,r.push(a)}return r},I=async(e,t,n,r)=>{let i=await fetch(e);if(!i)return;let a=await i.json();for(let e=0;e<a.length;e++){let t=a[e];if(!t)continue;let n=new c.Vector3(t.start.x,t.start.y,t.start.z);n.applyEuler(new c.Euler(-Math.PI/2,0,0));let i=new c.Vector3(t.end.x,t.end.y,t.end.z);i.applyEuler(new c.Euler(-Math.PI/2,0,0)),r.add(p(String(e),new c.Vector3((n.x+i.x)/2,(n.y+i.y)/2,(n.z+i.z)/2)))}a.sort((e,t)=>e.originalIndex-t.originalIndex),console.log(`lins`,a);let o=await fetch(t);if(!o)return;let s=await o.json(),u=Array.isArray(s)?s:Object.values(s),d=await fetch(n);if(console.log(`respcd`,d),!d)return;let f=await d.json();A(a,r),console.time(),a=await L(a,u,f,r),console.timeEnd();let m=new c.MeshBasicMaterial({color:`#cffd00`,transparent:!0,opacity:.4,side:c.DoubleSide}),h=new c.MeshBasicMaterial({color:`#ff0000`,transparent:!0,opacity:.4,side:c.DoubleSide});new c.LineBasicMaterial({color:`#cffd00`});let g=new c.PointsMaterial({color:65535,size:.01}),_=[],y=[];for(let e of a){if(!e.originalPoints||e.originalPoints.length===0)continue;let t=new Float32Array(e.originalPoints.length*3);e.originalPoints.forEach((e,n)=>{t[n*3]=e.x,t[n*3+1]=e.y,t[n*3+2]=e.z});let n=new c.BufferGeometry;n.setAttribute(`position`,new c.BufferAttribute(t,3)),_.push(n);let r=[];r.push(new c.Vector3(e.start.x,e.start.y,e.start.z)),r.push(new c.Vector3(e.end.x,e.end.y,e.end.z));let i=new c.BufferGeometry;i.setFromPoints(r),y.push(i)}if(_.length>0){let e=l(_,!1);e.rotateX(-Math.PI/2),r.add(new c.Points(e,g)),_.forEach(e=>e.dispose())}let b=[],x=[],S=new c.LineBasicMaterial({color:`#f30606`}),C=new c.MeshBasicMaterial({color:`#071ac4`});for(let e=0;e<u.length;e++){let t=new c.SphereGeometry(.05);t.translate(u[e].x,u[e].y,u[e].z);let n=[],r=new c.Vector3(u[e].x,u[e].y,u[e].z),i=new c.Quaternion(u[e].qx,u[e].qy,u[e].qz,u[e].qw),a=new c.Vector3(0,0,-1).applyQuaternion(i),o=r.clone().addScaledVector(a,.001);n.push(r),n.push(o);let s=new c.BufferGeometry;s.setFromPoints(n),x.push(s),b.push(t)}if(b.length>0){let e=l(b,!1);e.rotateX(-Math.PI/2),r.add(new c.Mesh(e,C)),b.forEach(e=>e.dispose())}if(x.length>0){let e=l(x,!1);e.rotateX(-Math.PI/2),r.add(new c.Line(e,S)),x.forEach(e=>e.dispose())}let w=v(f);w.rotateX(-Math.PI/2),r.add(w),console.log(`jsonpcdData`,f),a.forEach(e=>{if(e.doorAndBeamData&&e.doorAndBeamData.length>0){let t=new c.Group;e.doorAndBeamData.forEach(e=>{if(e.beamStart){let n=new c.BufferGeometry,r=new Float32Array([e.beamStart.x,e.beamStart.y,e.beamStart.z,e.beamEnd.x,e.beamEnd.y,e.beamEnd.z,e.beamEnd.x,e.beamEnd.y,e.beamEnd.z+e.beamHeight,e.beamStart.x,e.beamStart.y,e.beamStart.z+e.beamHeight]);n.setAttribute(`position`,new c.BufferAttribute(r,3)),n.setIndex([0,1,2,0,2,3]),n.computeVertexNormals(),n.rotateX(-Math.PI/2),t.add(new c.Mesh(n,m))}let n=new c.BufferGeometry,r=new Float32Array([e.doorStart.x,e.doorStart.y,e.doorStart.z,e.doorEnd.x,e.doorEnd.y,e.doorEnd.z,e.doorEnd.x,e.doorEnd.y,e.doorEnd.z+e.doorHeight,e.doorStart.x,e.doorStart.y,e.doorStart.z+e.doorHeight]);n.setAttribute(`position`,new c.BufferAttribute(r,3)),n.setIndex([0,1,2,0,2,3]),n.computeVertexNormals(),n.rotateX(-Math.PI/2),t.add(new c.Mesh(n,h));let i=`id:${e.id},Nid:${e.nearId}`,a=new c.Vector3((e.doorStart.x+e.doorEnd.x)/2,(e.doorStart.y+e.doorEnd.y)/2,e.doorStart.z+e.doorHeight);a.applyEuler(new c.Euler(-Math.PI/2,0,0)),t.add(p(i,a))}),r.add(t)}})},L=async(e,t,n)=>{if(e.length<=0||t.length<=0)return e;n||=[];for(let t of e){let e=[];if(!(!t.checkResults||t.checkResults.length<=0)){for(let n=0;n<t.checkResults.length;n++){let r=t.checkResults[n];if(r.isDoor=!1,!r.originalVertices||r.originalVertices.length==0)continue;let i=[];r.originalVertices.forEach(e=>{let t=[];e.forEach(e=>{t.push(new c.Vector3(e.x,e.y,e.z))}),i.push(t)});let a=i.length,o=Array.from({length:a},(e,t)=>t);function s(e){for(;o[e]!==e;)o[e]=o[o[e]],e=o[e];return e}for(let e=0;e<a;e++)for(let t=e+1;t<a;t++){let n=i[e],r=i[t];if(n.some(e=>r.some(t=>e.equals(t)))){let n=s(e),r=s(t);n!==r&&(o[n]=r)}}let l=new Map,u=-1,d=-1;for(let e=0;e<a;e++){let t=s(e);l.has(t)||l.set(t,{points:[],area:0,originalVertices:[],index:-1,gridWidth:-1,gridHeight:-1,centerPts:[]});let r=l.get(t),[a,o,f,p]=i[e],m=new c.Vector3().add(a).add(o).add(f).add(p).multiplyScalar(.25);if(r.points.push(...i[e]),r.area+=a.distanceTo(o)*a.distanceTo(p),r.originalVertices.push(i[e]),r.index=n,r.centerPts.push(m),u==-1||d==-1){let e=new c.Vector3().subVectors(a,o).normalize();d=a.distanceTo(p),u=a.distanceTo(o);let t=h(e,new c.Vector3(0,0,1));t&&t.parallel&&([d,u]=[u,d])}r.gridWidth=u,r.gridHeight=d}e.push(...l.values())}t.mergeCheckRegion=e}}let r=new Map;for(let t of e){let e=t.rooftopPz-t.start.z,n=Math.floor(t.rooftopPz-t.start.z);r.has(n)?r.set(n,{num:r.get(n).num+1,totalHeight:r.get(n).totalHeight+e}):r.set(n,{num:1,totalHeight:e})}let i=[...r.entries()].reduce((e,t)=>t[1].num>e[1].num?t:e),a=i[0];i[1].totalHeight/i[1].num;let o=[];for(let e=0;e<n.length;e++){let t=n[e];if(t.isFindBeam=!1,t.category==`door`){let n=[];t.coordinatesByArea.coordinates.forEach(e=>{n.push(new c.Vector3(e.x,e.y,e.z))}),n[0].equals(n[n.length-1])&&n.pop();let r=n[0].distanceTo(n[1]),i=n[0].distanceTo(n[3]);i>r&&([r,i]=[i,r],n=[n[0],n[3],n[2],n[1]]);let a=S(n[0],n[3],!0),s=S(n[1],n[2],!0);o.push({doorStartPt:a,doorEndPt:s,boxPoints:n,minZ:t.coordinatesByArea.heightData.minZ,maxZ:t.coordinatesByArea.heightData.maxZ,index:e,isFind:!1})}}for(let t=0;t<o.length;t++){let n=o[t];for(let t=0;t<e.length;t++){if(e[t].length<.8)continue;let r=[],i=S(e[t].start,e[t].end,!0);r.push(i),r.push(S(i,e[t].start,!0)),r.push(S(i,e[t].end,!0));let a=!1,o=n.boxPoints.length;for(let e=0;e<r.length;e++){let t=r[e].x,i=r[e].y;for(let e=0,r=o-1;e<o;r=e++){let o=n.boxPoints[e].x,s=n.boxPoints[e].y,c=n.boxPoints[r].x,l=n.boxPoints[r].y;s>i!=l>i&&t<(c-o)*(i-s)/(l-s)+o&&(a=!a)}if(a)break}if(a){let r=new c.Vector3(e[t].start.x,e[t].start.y,0),i=new c.Vector3(e[t].end.x,e[t].end.y,0),a=n.doorStartPt.distanceTo(n.doorEndPt),o=r.distanceTo(i),s=b(n.doorStartPt,n.doorEndPt,r,i);if(Math.abs(o-a)<.5&&s&&s.isParallel){n.doorStartPt=new c.Vector3(e[t].start.x,e[t].start.y,0),n.doorEndPt=new c.Vector3(e[t].end.x,e[t].end.y,0);break}}else continue}}let s=0,l=[],u=[];for(let r=0;r<e.length;r++){let i=e[r];if(!i.mergeCheckRegion||i.mergeCheckRegion.length==0)continue;i.doorAndBeamData=[],i.completePointAreaPercentage=-1;let d=new c.Box3;d.setFromPoints(i.originalPoints);let p=d.max.z,m=d.min.z,g=new c.Vector3(i.start.x,i.start.y,i.start.z),_=new c.Vector3(i.end.x,i.end.y,i.end.z),v=i.rooftopPz-i.start.z;if(v<a-.5||v>a+2)continue;let w=0;for(let a of i.mergeCheckRegion){let{facePoints:d,boxArea:T,totalInPlaneNum:E}=y(a.points,a.gridWidth,a.gridHeight,a.centerPts);if(!d||d.length<=0)continue;w+=a.area;let D=0;if(E!=0&&(D=a.gridHeight*a.gridWidth*E),(a.area-D)/T*100>85){let y=d[0].distanceTo(d[1]),w=y/i.length*100,T=d[0].distanceTo(d[3]),E=i.rooftopPz-d[2].z;if(m>i.start.z-.15&&i.length<1.5||i.length,T<.5||p<d[2].z||E>v/3)continue;let D=!1,O=!1,A=!1,M=[],N=d[2].clone(),P=d[3].clone(),F=d[0].clone(),I=d[1].clone(),L=!1;if(Math.abs(d[0].z-i.start.z)<.25){if(y>.2&&i.length>.5){let t=j(g,F),n=j(g,I),a=j(_,F),o=j(_,I),s=.1;if(t<s||n<s){let{newOriginalPoints:n}=k(F,I,i.originalPoints),a=[],o=[];if(n.forEach(e=>{e.z>=N.z?a.push(e):e.z<=F.z&&o.push(e)}),a.length>0){let n=new c.Vector3(i.start.x,i.start.y,0),a=new c.Vector3(i.end.x,i.end.y,0),o=new c.Vector3().subVectors(n,a).normalize(),l=new c.Vector3,u=1/0,d=-1;for(let t=0;t<e.length;t++){if(t==r)continue;let i=new c.Vector3(e[t].start.x,e[t].start.y,0),s=new c.Vector3(e[t].end.x,e[t].end.y,0);if(f(n,a,i,s)!=null)continue;let p=x(n,a,i,s);if(p&&p.angle>85){let e=b(n,a,new c.Line3(i,s).closestPointToPoint(n,!0,l),n);if(!e||e.direction!=`same`)continue;let r=f(n.clone().addScaledVector(o,10),a.clone().addScaledVector(o.clone().negate(),10),i,s);if(r==null)continue;let p=n.distanceTo(r.point);p<u&&(u=p,d=t)}else if(p&&p.angle<5){let e=j(n,i),t=j(n,s),r=j(a,i),o=j(a,s);if(e<.1||t<.1||r<.1||o<.1){u=1/0,d=-1;break}}}u!=1/0&&u<2&&e[d].length>.5&&(L=!0,t<s?(P=P.addScaledVector(o,u),F=F.addScaledVector(o,u)):(N=N.addScaledVector(o,u),I=I.addScaledVector(o,u)))}}else if(a<s||o<s){let{newOriginalPoints:t}=k(F,I,i.originalPoints),n=[],o=[];if(t.forEach(e=>{e.z>=N.z?n.push(e):e.z<=F.z&&o.push(e)}),n.length>0){let t=new c.Vector3(i.start.x,i.start.y,0),n=new c.Vector3(i.end.x,i.end.y,0),o=new c.Vector3().subVectors(n,t).normalize(),l=new c.Vector3,u=1/0,d=-1;for(let i=0;i<e.length;i++){if(i==r)continue;let a=new c.Vector3(e[i].start.x,e[i].start.y,0),s=new c.Vector3(e[i].end.x,e[i].end.y,0);if(f(t,n,a,s)!=null)continue;let p=x(t,n,a,s);if(p&&p.angle>85){let e=b(n,t,new c.Line3(a,s).closestPointToPoint(n,!0,l),n);if(!e||e.direction!=`same`)continue;let r=f(t.clone().addScaledVector(o,10),n.clone().addScaledVector(o.clone().negate(),10),a,s);if(r==null)continue;let p=n.distanceTo(r.point);p<u&&(u=p,d=i)}if(p&&p.angle<5){let e=j(t,a),r=j(t,s),i=j(n,a),o=j(n,s);if(e<.1||r<.1||i<.1||o<.1){u=1/0,d=-1;break}}}u!=1/0&&u<2&&e[d].length>.5&&(L=!0,a<s?(P=P.addScaledVector(o,u),F=F.addScaledVector(o,u)):(N=N.addScaledVector(o,u),I=I.addScaledVector(o,u)))}}y=P.distanceTo(N)}if(w<85&&y>.35||y>.7&&T>v-v/3){let e=!1;for(let n=0;n<t.length;n++){let r=n+1;if(r>=t.length)continue;let i=new c.Vector3(t[n].x,t[n].y,t[n].z),a=new c.Vector3(t[r].x,t[r].y,t[r].z);if(f(new c.Vector3(i.x,i.y,0),new c.Vector3(a.x,a.y,0),new c.Vector3(d[0].x,d[0].y,0),new c.Vector3(d[1].x,d[1].y,0),.1)!=null){e=!0;break}}for(let e=0;e<o.length;e++){let t=new c.Vector3(o[e].doorStartPt.x,o[e].doorStartPt.y,0),r=new c.Vector3(o[e].doorEndPt.x,o[e].doorEndPt.y,0),a=new c.Vector3(d[0].x,d[0].y,0),s=new c.Vector3(d[1].x,d[1].y,0),l=new c.Vector3().subVectors(t,r).normalize(),u=new c.Vector3().subVectors(a,s).normalize(),p=.2,m=f(t.clone().addScaledVector(l,p),r.clone().addScaledVector(l.clone().negate(),p),a.clone().addScaledVector(u,p),s.clone().addScaledVector(u.clone().negate(),p)),h=d[2].z<o[e].maxZ||Math.abs(d[2].z-o[e].maxZ)<.3,g=b(t.clone(),r.clone(),a.clone(),s.clone());if(m&&h)n[o[e].index].isFindBeam=!0,O=!0,o[e].isFind=!0;else if(h&&g&&g.isParallel){let a=S(t,r),s=new c.Line3(new c.Vector3(i.start.x,i.start.y,0),new c.Vector3(i.end.x,i.end.y,0)),l=new c.Vector3,u=s.closestPointToPoint(a,!0,l);a.distanceTo(u)<.3&&(n[o[e].index].isFindBeam=!0)}}if(!O)for(let e=0;e<o.length;e++){let t=new c.Vector3(o[e].doorStartPt.x,o[e].doorStartPt.y,0),r=new c.Vector3(o[e].doorEndPt.x,o[e].doorEndPt.y,0),a=new c.Vector3(d[0].x,d[0].y,0),s=new c.Vector3(d[1].x,d[1].y,0),l=d[2].z<o[e].maxZ||Math.abs(d[2].z-o[e].maxZ)<.3,u=Math.abs(d[2].z-o[e].maxZ)>.2;if(l&&!u){let l=C(t.clone(),r.clone(),a.clone(),s.clone(),{distanceThreshold:2});if(l&&l.type==`collinear_gap`){let u=S(t,r),f=new c.Line3(new c.Vector3(i.start.x,i.start.y,0),new c.Vector3(i.end.x,i.end.y,0)),p=new c.Vector3,m=f.closestPointToPoint(u,!0,p);if(u.distanceTo(m)<.15){let i=a.distanceTo(t),u=s.distanceTo(t);if(i>.25&&u>.25){let e=a.distanceTo(r),t=s.distanceTo(r);if(e>.2&&t>.2)continue}O=!0,A=!0,M.push(o[e].index),o[e].isFind=!0,n[o[e].index].isFindBeam=!0;let f=new c.Vector3().subVectors(a,s).normalize(),p=t.distanceTo(r)-l.overlap;u<i&&f.negate(),i>u?(N=d[2].clone().addScaledVector(f,p),P=d[3].clone(),F=N.clone(),I=P.clone(),F.z=I.z=d[0].z):(N=d[2].clone(),P=d[3].clone().addScaledVector(f,p),F=N.clone(),I=P.clone(),F.z=I.z=d[0].z)}}else if(l&&l.type==`collinear_overlap`&&l.maxPerpendicularDistance<.2){let i=a.distanceTo(s),c=r.distanceTo(t);(Math.abs(i-l.overlap)<.25||Math.abs(c-l.overlap)<.25)&&(O=!0,A=!0,M.push(o[e].index),o[e].isFind=!0,n[o[e].index].isFindBeam=!0)}}}(O||e||T>v/2&&y>.8&&Math.abs(m-i.start.z)<.4)&&(D=!0)}}if(D){if(!A){new c.PointsMaterial({color:65535,size:.01});let e=new c.Vector3().subVectors(F,I),t=new c.Vector3().subVectors(I,F),n=1/0,r=1/0,a=F.clone(),o=I.clone(),s=[];if(i.originalPoints.forEach(e=>{e.z-F.z<T/2+.1&&e.z-F.z>T/2-.1&&s.push(e)}),s.length>0){let{newProjectPoints:i}=k(F,I,s);for(let e of i)e.z=F.z;i.forEach(i=>{let s=new c.Vector3(i.x,i.y,F.z),l=new c.Vector3().subVectors(s,I),u=h(e,new c.Vector3().subVectors(s,F)),d=h(t,l);if(u&&u.parallel&&u.sameDirection){let e=s.distanceTo(F);e<n&&(n=e,a=s)}else if(d&&d.parallel&&d.sameDirection){let e=s.distanceTo(I);e<r&&(r=e,o=s)}})}a.equals(F)||(F=a,P=new c.Vector3(F.x,F.y,P.z)),o.equals(I)||(I=o,N=new c.Vector3(I.x,I.y,N.z))}u.includes(r)||u.push(r),i.checkResults[a.index].isDoor=!0;let t={id:s,beamStart:N,beamEnd:P,beamHeight:i.rooftopPz-d[2].z,doorStart:F,doorEnd:I,doorHeight:T,nearId:-1,type:L?`extendBeam`:`doorBeam`,isDoor:O,isPullOutDoor:A,pcbDoorIndexs:M};for(let n=0;n<l.length;n++){let r=e[l[n].linesIndex].doorAndBeamData[l[n].doorIndex];if(r.nearId!=-1)continue;let i=b(N.clone(),P.clone(),r.beamStart.clone(),r.beamEnd.clone());if(i&&i.isParallel){let e=N.distanceTo(P),n=r.beamStart.distanceTo(r.beamEnd),i=new c.Line3(N.clone(),P.clone());i.start.z=0,i.end.z=0;let a=S(r.beamStart,r.beamEnd,!0);n>e&&(i.start=new c.Vector3(r.beamStart.x,r.beamStart.y,0),i.end=new c.Vector3(r.beamEnd.x,r.beamEnd.y,0),a=S(N,P,!0));let o=new c.Vector3,l=i.closestPointToPoint(a,!0,o),u=l.distanceTo(a),d=new c.Vector3().subVectors(a,l).normalize(),p=d.clone().negate(),m=a.clone().addScaledVector(p,u+.1),h=f(i.start,i.end,a,m);if(u<.4&&h!=null){r.nearId=s,t.nearId=r.id;let i=.8;(r.isDoor||t.isDoor)&&(r.isDoor=!0,t.isDoor=!0);let a=M.some(e=>r.pcbDoorIndexs.includes(e));e<n?e>i&&!r.isPullOutDoor?(r.beamStart=N.clone().addScaledVector(p,u),r.beamEnd=P.clone().addScaledVector(p,u),r.doorStart=F.clone().addScaledVector(p,u),r.doorEnd=I.clone().addScaledVector(p,u),r.doorHeight=t.doorHeight,r.beamStart.z=r.beamEnd.z=r.doorStart.z+r.doorHeight):!A||A&&r.isPullOutDoor&&a?(t.beamStart=r.beamStart.clone().addScaledVector(d,u),t.beamEnd=r.beamEnd.clone().addScaledVector(d,u),t.doorStart=r.doorStart.clone().addScaledVector(d,u),t.doorEnd=r.doorEnd.clone().addScaledVector(d,u),t.doorHeight=r.doorHeight,t.beamStart.z=t.beamEnd.z=t.doorStart.z+t.doorHeight):(r.beamStart=N.clone().addScaledVector(p,u),r.beamEnd=P.clone().addScaledVector(p,u),r.doorStart=F.clone().addScaledVector(p,u),r.doorEnd=I.clone().addScaledVector(p,u),r.doorHeight=t.doorHeight,r.beamStart.z=r.beamEnd.z=r.doorStart.z+r.doorHeight):e>n&&(n>i&&!A?(t.beamStart=r.beamStart.clone().addScaledVector(p,u),t.beamEnd=r.beamEnd.clone().addScaledVector(p,u),t.doorStart=r.doorStart.clone().addScaledVector(p,u),t.doorEnd=r.doorEnd.clone().addScaledVector(p,u),t.doorHeight=r.doorHeight,t.beamStart.z=t.beamEnd.z=t.doorStart.z+t.doorHeight):!r.isPullOutDoor||A&&r.isPullOutDoor&&a?(r.beamStart=N.clone().addScaledVector(d,u),r.beamEnd=P.clone().addScaledVector(d,u),r.doorStart=F.clone().addScaledVector(d,u),r.doorEnd=I.clone().addScaledVector(d,u),r.doorHeight=t.doorHeight,r.beamStart.z=r.beamEnd.z=r.doorStart.z+r.doorHeight):(t.beamStart=r.beamStart.clone().addScaledVector(p,u),t.beamEnd=r.beamEnd.clone().addScaledVector(p,u),t.doorStart=r.doorStart.clone().addScaledVector(p,u),t.doorEnd=r.doorEnd.clone().addScaledVector(p,u),t.doorHeight=r.doorHeight,t.beamStart.z=t.beamEnd.z=t.doorStart.z+t.doorHeight));break}}}i.doorAndBeamData.push(t),l.push({linesIndex:r,doorIndex:i.doorAndBeamData.length-1}),s++}}}let T=i.length*(i.rooftopPz-i.start.z);i.completePointAreaPercentage=w/T*100}if(u.length>1)for(let t=0;t<u.length;t++){let n=[];for(let r=0;r<e[u[t]].doorAndBeamData.length;r++){let i=e[u[t]].doorAndBeamData[r].doorStart.clone(),a=e[u[t]].doorAndBeamData[r].doorEnd.clone();i.z=a.z=0;let o=i.distanceTo(a),s=!1;for(let n=0;n<u.length;n++)if(t!=n)for(let t=0;t<e[u[n]].doorAndBeamData.length;t++){let r=e[u[n]].doorAndBeamData[t].doorStart.clone(),l=e[u[n]].doorAndBeamData[t].doorEnd.clone();r.z=l.z=0,r.distanceTo(l);let d=C(i,a,r,l);if(d.type==`collinear_overlap`&&d.maxPerpendicularDistance<.2){let e=new c.Line3(r,l),t=S(i,a),n=e.closestPointToPoint(t,!0,new c.Vector3),u=t.distanceTo(n);if(Math.abs(o-d.overlap)<.1&&u<.05){s=!0;break}}}s||n.push(e[u[t]].doorAndBeamData[r])}n.length!=e[u[t]].doorAndBeamData.length&&(e[u[t]].doorAndBeamData=n)}let d=.65,p=[],m=[],g=.2;for(let t=0;t<o.length;t++){if(o[t].isFind)continue;let r=new c.Vector3(o[t].doorStartPt.x,o[t].doorStartPt.y,0),i=new c.Vector3(o[t].doorEndPt.x,o[t].doorEndPt.y,0),a=new c.Vector3().subVectors(r,i).normalize();for(let l=0;l<e.length;l++){if(e[l].length<1||e[l].completePointAreaPercentage>60)continue;let u=new c.Vector3(e[l].start.x,e[l].start.y,0),h=new c.Vector3(e[l].end.x,e[l].end.y,0),_=new c.Vector3().subVectors(u,h).normalize(),v=.3,y=r.clone().addScaledVector(a,v),S=i.clone().addScaledVector(a.clone().negate(),v),w=u.clone().addScaledVector(_,v),T=h.clone().addScaledVector(_.clone().negate(),v);f(y,S,w,T);let E=b(y,S,w,T),D=r.distanceTo(u)>i.distanceTo(u)?i:r,O=D.distanceTo(u),k=D.distanceTo(h);if((O<g||k<g)&&E&&!E.isParallel){let a=r.distanceTo(u)<r.distanceTo(h)?u:h,f=!1;for(let g=0;g<e.length;g++){if(g==l||p.includes(g)||e[g].length<1)continue;let _=new c.Vector3(e[g].start.x,e[g].start.y,0),v=new c.Vector3(e[g].end.x,e[g].end.y,0),y=x(u,h,_,v);if(y&&y.angle>85){let y=new c.Line3(u,h),x=new c.Vector3,S=y.closestPointToPoint(_,!0,x),w=y.closestPointToPoint(v,!0,x),T=S.distanceTo(_),E=w.distanceTo(v),D,O;if(T>E&&E>d)D=v,O=w;else if(T<E&&T>d)D=_,O=S;else continue;if(T=a.distanceTo(_),E=a.distanceTo(v),T>2&&E>2)continue;let k=b(D,O,_,v);if(!k||!k.isParallel)continue;let A=C(r,i,D,O);if(A&&A.type==`collinear_overlap`&&A.maxPerpendicularDistance<.2)continue;D.z=O.z=e[g].start.z;let j=!1;for(let e=0;e<m.length;e++){let t=C(m[e].start,m[e].end,D,O);if(t&&t.type==`collinear_overlap`&&t.maxPerpendicularDistance<.2){j=!0;break}}if(j)continue;m.push({start:D,end:O}),p.push(g);let M=e[g].rooftopPz-e[g].start.z,N=e[l].rooftopPz-e[l].start.z,P=M<N?M:N;e[g].doorAndBeamData||(e[g].doorAndBeamData=[]),e[g].doorAndBeamData.push({id:s,doorStart:D,doorEnd:O,doorHeight:P,nearId:-1,type:`onlyDoor`,isDoor:!0}),n[o[t].index].isFindBeam=!0,n[o[t].index].isFindOnlyDoor=!0,s++,f=!0;break}}if(f)break}}}let _=0;for(let t=0;t<e.length;t++){let n=e[t];n.mainBeamDataIds||=[],n.isMainBeam||=!1;let{minZ:r,maxZ:i}=O(n);if(r==null||i==null)continue;i-r;let o=n.rooftopPz-n.start.z;if(o<a-.5||o>a+2||Math.abs(r-n.start.z)>.3||n.length<1)continue;let s=new c.Vector3(n.start.x,n.start.y,0),l=new c.Vector3(n.end.x,n.end.y,0);for(let r=0;r<e.length;r++){if(r==t||e[r].length<1)continue;let a=e[r],o=new c.Vector3(a.start.x,a.start.y,0),u=new c.Vector3(a.end.x,a.end.y,0),d=b(s,l,o,u),f=S(s,l),p=new c.Line3(o,u).closestPointToPoint(f,!0,new c.Vector3),m=f.distanceTo(p),h=O(a);if(!(h.minZ===void 0||h.maxZ===void 0)&&!(h.maxZ-h.minZ>1)&&d&&d.isParallel&&m<.25&&Math.abs(i-h.minZ)<.1){let e=`mainBeamId${_}`;a.mainBeamDataIds||=[],a.isMainBeam||=!0,a.mainBeamDataIds.push(e),n.mainBeamDataIds.push(e),_++;break}}}return e},R=(e,t)=>{if(e.length<=0)return e;for(let t of e){let e=[];for(let n=0;n<t.checkResults.length;n++){let r=t.checkResults[n];if(r.isDoor=!1,!r.originalVertices||r.originalVertices.length==0)continue;let i=[];r.originalVertices.forEach(e=>{let t=[];e.forEach(e=>{t.push(new c.Vector3(e.x,e.y,e.z))}),i.push(t)});let a=i.length,o=Array.from({length:a},(e,t)=>t);function s(e){for(;o[e]!==e;)o[e]=o[o[e]],e=o[e];return e}for(let e=0;e<a;e++)for(let t=e+1;t<a;t++){let n=i[e],r=i[t];if(n.some(e=>r.some(t=>e.equals(t)))){let n=s(e),r=s(t);n!==r&&(o[n]=r)}}let l=new Map,u=-1,d=-1;for(let e=0;e<a;e++){let t=s(e);l.has(t)||l.set(t,{points:[],area:0,originalVertices:[],index:-1,gridWidth:-1,gridHeight:-1,centerPts:[]});let r=l.get(t),[a,o,f,p]=i[e],m=new c.Vector3().add(a).add(o).add(f).add(p).multiplyScalar(.25);if(r.points.push(...i[e]),r.area+=a.distanceTo(o)*a.distanceTo(p),r.originalVertices.push(i[e]),r.index=n,r.centerPts.push(m),u==-1||d==-1){let e=new c.Vector3().subVectors(a,o).normalize();d=a.distanceTo(p),u=a.distanceTo(o);let t=h(e,new c.Vector3(0,0,1));t&&t.parallel&&([d,u]=[u,d])}r.gridWidth=u,r.gridHeight=d}e.push(...l.values())}t.mergeCheckRegion=e}let n=new Map;for(let t of e){let e=t.rooftopPz-t.start.z,r=Math.floor(t.rooftopPz-t.start.z);n.has(r)?n.set(r,{num:n.get(r).num+1,totalHeight:n.get(r).totalHeight+e}):n.set(r,{num:1,totalHeight:e})}let r=[...n.entries()].reduce((e,t)=>t[1].num>e[1].num?t:e),i=r[0];r[1].totalHeight/r[1].num,new c.MeshBasicMaterial({color:`#cffd00`,transparent:!0,opacity:.4,side:c.DoubleSide}),new c.PointsMaterial({color:65535,size:.01});for(let n=0;n<e.length;n++){let r=e[n],{minZ:a,maxZ:o}=O(r);if(a==null||o==null)continue;r.doorAndBeamData=[],r.completePointAreaPercentage=-1,o-a;let s=r.rooftopPz-r.start.z;if(s<i-.5||s>i+2||Math.abs(a-r.start.z)>.3||r.length<1)continue;let l=new c.Vector3(r.start.x,r.start.y,0),u=new c.Vector3(r.end.x,r.end.y,0);for(let r=0;r<e.length;r++){if(r==n||e[r].length<1)continue;let i=e[r],a=new c.Vector3(i.start.x,i.start.y,0),s=new c.Vector3(i.end.x,i.end.y,0),d=b(l,u,a,s),f=S(l,u),m=new c.Line3(a,s).closestPointToPoint(f,!0,new c.Vector3),h=f.distanceTo(m),g=O(i);if(!(g.minZ===void 0||g.maxZ===void 0)&&!(g.maxZ-g.minZ>1)&&d&&d.isParallel&&h<.25&&Math.abs(o-g.minZ)<.1){let i=S(l,u),o=S(a,s);i.z=o.z=1,i.applyEuler(new c.Euler(-Math.PI/2,0,0)),o.applyEuler(new c.Euler(-Math.PI/2,0,0));let d=`w:${e[n].originaIndex},b:${e[r].originaIndex}`,f=`b:${e[r].originaIndex},w:${e[n].originaIndex}`;t.add(p(d,i,{textColor:`#fd0000`})),t.add(p(f,o,{textColor:`#fd0000`}))}}}return e};exports.classifySegments=C,exports.getAllGeometry=F,exports.getBeamLine=L,exports.getMainBeamLine=R,exports.getParallelism=b,exports.getPointCloudMinMax=O,exports.getPointCoverageOnQuad=E,exports.isParallel=h,exports.perpendicularInfo=x,exports.removeNoisePoints=A,exports.segmentsIntersect2D=f,exports.usegetBeamLine=I;
+Object.defineProperty(exports,Symbol.toStringTag,{value:`Module`});var e=Object.create,t=Object.defineProperty,n=Object.getOwnPropertyDescriptor,r=Object.getOwnPropertyNames,i=Object.getPrototypeOf,a=Object.prototype.hasOwnProperty,o=(e,i,o,s)=>{if(i&&typeof i==`object`||typeof i==`function`)for(var c=r(i),l=0,u=c.length,d;l<u;l++)d=c[l],!a.call(e,d)&&d!==o&&t(e,d,{get:(e=>i[e]).bind(null,d),enumerable:!(s=n(i,d))||s.enumerable});return e},s=(n,r,a)=>(a=n==null?{}:e(i(n)),o(r||!n||!n.__esModule?t(a,`default`,{value:n,enumerable:!0}):a,n));let c=require(`three`);c=s(c);function l(e,t=!1){let n=e[0].index!==null,r=new Set(Object.keys(e[0].attributes)),i=new Set(Object.keys(e[0].morphAttributes)),a={},o={},s=e[0].morphTargetsRelative,l=new c.BufferGeometry,d=0;for(let c=0;c<e.length;++c){let u=e[c],f=0;if(n!==(u.index!==null))return console.error(`THREE.BufferGeometryUtils: .mergeGeometries() failed with geometry at index `+c+`. All geometries must have compatible attributes; make sure index attribute exists among all geometries, or in none of them.`),null;for(let e in u.attributes){if(!r.has(e))return console.error(`THREE.BufferGeometryUtils: .mergeGeometries() failed with geometry at index `+c+`. All geometries must have compatible attributes; make sure "`+e+`" attribute exists among all geometries, or in none of them.`),null;a[e]===void 0&&(a[e]=[]),a[e].push(u.attributes[e]),f++}if(f!==r.size)return console.error(`THREE.BufferGeometryUtils: .mergeGeometries() failed with geometry at index `+c+`. Make sure all geometries have the same number of attributes.`),null;if(s!==u.morphTargetsRelative)return console.error(`THREE.BufferGeometryUtils: .mergeGeometries() failed with geometry at index `+c+`. .morphTargetsRelative must be consistent throughout all geometries.`),null;for(let e in u.morphAttributes){if(!i.has(e))return console.error(`THREE.BufferGeometryUtils: .mergeGeometries() failed with geometry at index `+c+`.  .morphAttributes must be consistent throughout all geometries.`),null;o[e]===void 0&&(o[e]=[]),o[e].push(u.morphAttributes[e])}if(t){let e;if(n)e=u.index.count;else if(u.attributes.position!==void 0)e=u.attributes.position.count;else return console.error(`THREE.BufferGeometryUtils: .mergeGeometries() failed with geometry at index `+c+`. The geometry must have either an index or a position attribute`),null;l.addGroup(d,e,c),d+=e}}if(n){let t=0,n=[];for(let r=0;r<e.length;++r){let i=e[r].index;for(let e=0;e<i.count;++e)n.push(i.getX(e)+t);t+=e[r].attributes.position.count}l.setIndex(n)}for(let e in a){let t=u(a[e]);if(!t)return console.error(`THREE.BufferGeometryUtils: .mergeGeometries() failed while trying to merge the `+e+` attribute.`),null;l.setAttribute(e,t)}for(let e in o){let t=o[e][0].length;if(t===0)break;l.morphAttributes=l.morphAttributes||{},l.morphAttributes[e]=[];for(let n=0;n<t;++n){let t=[];for(let r=0;r<o[e].length;++r)t.push(o[e][r][n]);let r=u(t);if(!r)return console.error(`THREE.BufferGeometryUtils: .mergeGeometries() failed while trying to merge the `+e+` morphAttribute.`),null;l.morphAttributes[e].push(r)}}return l}function u(e){let t,n,r,i=-1,a=0;for(let o=0;o<e.length;++o){let s=e[o];if(t===void 0&&(t=s.array.constructor),t!==s.array.constructor)return console.error(`THREE.BufferGeometryUtils: .mergeAttributes() failed. BufferAttribute.array must be of consistent array types across matching attributes.`),null;if(n===void 0&&(n=s.itemSize),n!==s.itemSize)return console.error(`THREE.BufferGeometryUtils: .mergeAttributes() failed. BufferAttribute.itemSize must be consistent across matching attributes.`),null;if(r===void 0&&(r=s.normalized),r!==s.normalized)return console.error(`THREE.BufferGeometryUtils: .mergeAttributes() failed. BufferAttribute.normalized must be consistent across matching attributes.`),null;if(i===-1&&(i=s.gpuType),i!==s.gpuType)return console.error(`THREE.BufferGeometryUtils: .mergeAttributes() failed. BufferAttribute.gpuType must be consistent across matching attributes.`),null;a+=s.count*n}let o=new t(a),s=new c.BufferAttribute(o,n,r),l=0;for(let t=0;t<e.length;++t){let r=e[t];if(r.isInterleavedBufferAttribute){let e=l/n;for(let t=0,i=r.count;t<i;t++)for(let i=0;i<n;i++){let n=r.getComponent(t,i);s.setComponent(t+e,i,n)}}else o.set(r.array,l);l+=r.count*n}return i!==void 0&&(s.gpuType=i),s}var d=class extends c.Object3D{constructor(e=document.createElement(`div`)){super(),this.isCSS2DObject=!0,this.element=e,this.element.style.position=`absolute`,this.element.style.userSelect=`none`,this.element.setAttribute(`draggable`,!1),this.center=new c.Vector2(.5,.5),this.addEventListener(`removed`,function(){this.traverse(function(e){e.element&&e.element instanceof e.element.ownerDocument.defaultView.Element&&e.element.parentNode!==null&&e.element.remove()})})}copy(e,t){return super.copy(e,t),this.element=e.element.cloneNode(!0),this.center=e.center,this}};new c.Vector3,new c.Matrix4,new c.Matrix4,new c.Vector3,new c.Vector3;var f=(e,t,n,r,i=1e-8)=>{let a=t.x-e.x,o=t.y-e.y,s=r.x-n.x,l=r.y-n.y,u=a*l-o*s;if(Math.abs(u)<i)return null;let d=((n.x-e.x)*l-(n.y-e.y)*s)/u,f=((n.x-e.x)*o-(n.y-e.y)*a)/u;return d<-i||d>1+i||f<-i||f>1+i?null:{point:new c.Vector3(e.x+d*a,e.y+d*o,e.z)}},p=(e,t,n={})=>{let{fontSize:r=36,fontFamily:i=`sans-serif`,backgroundColor:a=`rgba(0,0,0,0)`,textColor:o=`#ffffff`,padding:s=10,maxWidth:l=256}=n,u=document.createElement(`canvas`).getContext(`2d`);u.font=`bold ${r}px ${i}`;let d=u.measureText(e).width,f=Math.min(d+s*2,l),p=r+s*2,h=document.createElement(`canvas`);h.width=f,h.height=p;let g=h.getContext(`2d`);g.fillStyle=a,g.fillRect(0,0,f,p),g.fillStyle=o,g.font=`bold ${r}px ${i}`,g.textAlign=`center`,g.textBaseline=`middle`,d>l-s*2?m(g,e,f/2,p/2,l,r):g.fillText(e,f/2,p/2);let _=new c.CanvasTexture(h),v=new c.SpriteMaterial({map:_,depthTest:!1,transparent:!0}),y=new c.Sprite(v);y.position.set(t.x,t.y,t.z);let b=.005;return y.scale.set(f*b,p*b,1),y};function m(e,t,n,r,i,a){let o=t.split(``),s=``,c=a*1.2,l=r;for(let t=0;t<o.length;t++){let r=s+o[t];e.measureText(r).width>i&&t>0?(e.fillText(s,n,l),s=o[t],l+=c):s=r}e.fillText(s,n,l)}var h=(e,t,n=.01)=>{let r=Math.sqrt(e.x*e.x+e.y*e.y+e.z*e.z),i=Math.sqrt(t.x*t.x+t.y*t.y+t.z*t.z);if(r<1e-10||i<1e-10)return{parallel:!1,sameDirection:null,angle:NaN};let a={x:e.x/r,y:e.y/r,z:e.z/r},o={x:t.x/i,y:t.y/i,z:t.z/i},s=a.x*o.x+a.y*o.y+a.z*o.z;s=Math.max(-1,Math.min(1,s));let c=Math.acos(Math.abs(s)),l=c<=n;return{parallel:l,sameDirection:l?s>0:null,angle:c}};function g(e,t){let n=new c.BufferGeometry,r=[],i=[];e.forEach(e=>r.push(e.x,e.y,e.z)),t.forEach(e=>r.push(e.x,e.y,e.z)),i.push(0,2,1,0,3,2),i.push(4,5,6,4,6,7),i.push(0,1,5,0,5,4),i.push(1,2,6,1,6,5),i.push(2,3,7,2,7,6),i.push(3,0,4,3,4,7),n.setAttribute(`position`,new c.Float32BufferAttribute(r,3)),n.setIndex(i),n.computeVertexNormals();let a=new c.MeshBasicMaterial({color:`green`,transparent:!0,opacity:.5,side:c.DoubleSide});return new c.Mesh(n,a)}var _=class{constructor(e={}){this.config={value:e.value||`0.0000m`,position:e.position||new c.Vector3(0,0,0),style:{fontFamily:e.style?.fontFamily||`Arial, sans-serif`,fontSize:e.style?.fontSize||`16px`,fontWeight:e.style?.fontWeight||`normal`,fontStyle:e.style?.fontStyle||`normal`},backgroundOpacity:e.backgroundOpacity||.7},this.label=null,this.labelElement=null,this.createLabel()}createLabel(){return this.labelElement=document.createElement(`div`),this.labelElement.textContent=this.config.value,this.applyBlackWhiteStyle(),this.label=new d(this.labelElement),this.label.position.copy(this.config.position),this.label}applyBlackWhiteStyle(){let e=this.labelElement.style,t=this.config;e.fontFamily=t.style.fontFamily,e.fontSize=t.style.fontSize,e.fontWeight=t.style.fontWeight,e.fontStyle=t.style.fontStyle,e.color=`#FFFFFF`,e.textShadow=`0 1px 2px rgba(0, 0, 0, 0.5)`,e.backgroundColor=`rgba(0, 0, 0, ${t.backgroundOpacity})`,e.padding=`6px 12px`,e.borderRadius=`4px`,e.position=`absolute`,e.whiteSpace=`nowrap`,e.userSelect=`none`,e.pointerEvents=`none`,e.zIndex=`1000`,e.textAlign=`center`,e.lineHeight=`1.2`}updateValue(e){this.labelElement&&(this.labelElement.textContent=e)}updatePosition(e){this.label&&this.label.position.copy(e)}updateFontStyle(e){if(this.labelElement){let t=this.labelElement.style;e.fontFamily&&(t.fontFamily=e.fontFamily),e.fontSize&&(t.fontSize=e.fontSize),e.fontWeight&&(t.fontWeight=e.fontWeight),e.fontStyle&&(t.fontStyle=e.fontStyle)}}getLabel(){return this.label}getLabelElement(){return this.labelElement}dispose(){this.labelElement&&this.labelElement.parentNode&&this.labelElement.parentNode.removeChild(this.labelElement),this.label=null,this.labelElement=null}},v=e=>{let t=new c.Group;if(e&&e.length){let n=`green.red.pink.white.azure.yellow.bisque.grey.blanchedalmond.blue.blueviolet.brown.burlywood.cadetblue.chartreuse.chocolate.coral.cornflowerblue.cornsilk.crimson.cyan.darkblue.darkcyan.darkgoldenrod.darkgray.darkgreen.darkgrey.darkkhaki.darkmagenta.darkolivegreen.darkorange.darkorchid.darkred.darksalmon.darkseagreen.darkslateblue.darkslategray.darkslategrey.darkturquoise.darkviolet.deeppink.deepskyblue.dimgray.dimgrey.dodgerblue.firebrick.floralwhite.forestgreen.fuchsia.gainsboro.ghostwhite.gold.goldenrod.gray.aliceblue.greenyellow.black.honeydew.hotpink.indianred.indigo.ivory.khaki.lavender.lavenderblush.lawngreen.lemonchiffon.lightblue.lightcoral.lightcyan.lightgoldenrodyellow.lightgray.lightgreen.lightgrey.lightpink.lightsalmon.lightseagreen.lightskyblue.lightslategray.lightslategrey.lightsteelblue.lightyellow.lime.limegreen.linen.magenta.maroon.mediumaquamarine.mediumblue.mediumorchid.mediumpurple.mediumseagreen.mediumslateblue.mediumspringgreen.mediumturquoise.mediumvioletred.midnightblue.mintcream.mistyrose.moccasin.navajowhite.navy.oldlace.olive.olivedrab.orange.orangered.orchid.palegoldenrod.palegreen.paleturquoise.palevioletred.papayawhip.peachpuff.peru.pink.plum.powderblue.purple.rebeccapurple.yellow.rosybrown.royalblue.saddlebrown.salmon.sandybrown.seagreen.seashell.sienna.silver.skyblue.slateblue.slategray.slategrey.snow.springgreen.steelblue.tan.teal.thistle.tomato.turquoise.violet.wheat.white.whitesmoke.yellow.yellowgreen`.split(`.`);e.forEach((e,r)=>{let i=n[r],a=new c.Group,{pcdPoints:o,category:s,name:l,center:u,box:d,obbBox:f,contour:p,nline:m,coordinatesByWidth:h,coordinatesByArea:v}=e,y=o.length;if(s!=`door`)return;s===`balcony railing`&&console.log(`pcdObj`,e);let b=new c.Box3(new c.Vector3(d.min.x,d.min.y,d.min.z),new c.Vector3(d.max.x,d.max.y,d.max.z)),x=new c.Box3Helper(b,i);a.add(x);let S=new _({value:l,position:new c.Vector3(u.x,u.y,u.z)}).getLabel();a.add(S);let C=new c.InstancedMesh(new c.SphereGeometry(.007,8,8),new c.MeshBasicMaterial({color:i}),y);a.add(C);let w=new c.Matrix4,T=0;if(o.forEach(e=>{w.setPosition(e.x,e.y,e.z),C.setMatrixAt(T,w),T++}),C.instanceMatrix.needsUpdate=!0,m){let e=new c.BufferGeometry;e.setFromPoints([new c.Vector3(m.start.x,m.start.y,m.start.z),new c.Vector3(m.end.x,m.end.y,m.end.z)]);let t=new c.LineBasicMaterial({color:`red`,depthTest:!1,depthWrite:!1});new c.Line(e,t);let n=new c.Vector3().copy(new c.Vector3(m.start.x,m.start.y,m.start.z)).add(new c.Vector3(m.end.x,m.end.y,m.end.z)).multiplyScalar(.5),r=new c.Vector3(m.centerDir.x,m.centerDir.y,m.centerDir.z);new c.ArrowHelper(r,n,1,`red`)}if(p){let e=[];p.forEach((t,n)=>{e.push(new c.Vector3(t.x,t.y,t.z))});let n=new c.BufferGeometry;n.setFromPoints(e);let r=new c.LineBasicMaterial({color:`yellow`,depthTest:!1,depthWrite:!1}),i=new c.Line(n,r);t.add(i)}if(v){let e=[],n=[];v.coordinates.forEach((t,r)=>{e.push(new c.Vector3(t.x,t.y,t.z)),n.length<4&&n.push(new c.Vector3(t.x,t.y,t.z))});let r=new c.BufferGeometry;r.setFromPoints(e);let i=new c.LineBasicMaterial({color:`green`,depthTest:!1,depthWrite:!1}),a=new c.Line(r,i);t.add(a);let o=g(n.map(e=>new c.Vector3(e.x,e.y,v.heightData.minZ)),n.map(e=>new c.Vector3(e.x,e.y,e.z)));t.add(o);let[s,l,u]=e,d=s.distanceTo(l),f=new c.Vector3().copy(s).add(l).multiplyScalar(.5);new _({value:String(d.toFixed(3)),position:f}).getLabel();let p=l.distanceTo(u),m=new c.Vector3().copy(l).add(u).multiplyScalar(.5);new _({value:String(p.toFixed(3)),position:m}).getLabel()}if(f){let e=new c.BoxGeometry(f.size.x,f.size.y,f.size.z),n=new c.LineBasicMaterial({color:i,linewidth:2}),r=new c.EdgesGeometry(e),a=new c.LineSegments(r,n),o=new Vector3(f.center.x,f.center.y,f.center.z),s=new Quaternion(f.rotation.x,f.rotation.y,f.rotation.z,f.rotation.w);a.position.copy(o),a.quaternion.copy(s),t.add(a)}})}return t},y=(e,t,n,r)=>{let i=e,a=[];if(r&&r.length&&(a=r.map(e=>e.clone())),!i||i.length<4)return{};let o=new c.Vector3().subVectors(i[1],i[0]).normalize(),s=new c.Vector3().subVectors(i[3],i[0]).normalize(),l=o.clone().cross(s).normalize(),u=1/0,d=-1/0,f=1/0,p=-1/0;for(let e of i){let t=e.dot(o),n=e.dot(s);t<u&&(u=t),t>d&&(d=t),n<f&&(f=n),n>p&&(p=n)}let m=i[0].dot(l),g=o.clone().multiplyScalar(u).add(s.clone().multiplyScalar(f)).add(l.clone().multiplyScalar(m)),_=o.clone().multiplyScalar(d).add(s.clone().multiplyScalar(f)).add(l.clone().multiplyScalar(m)),v=o.clone().multiplyScalar(d).add(s.clone().multiplyScalar(p)).add(l.clone().multiplyScalar(m)),y=o.clone().multiplyScalar(u).add(s.clone().multiplyScalar(p)).add(l.clone().multiplyScalar(m)),b=[g,_,v,y],x=h(new c.Vector3().subVectors(g,_),new c.Vector3(0,0,1));x&&x.parallel&&(b=[g,y,v,_]);let S=0;if(r.length>0){let e=new c.Vector3().subVectors(b[1],b[0]).normalize(),r=e.clone().negate(),i=new c.Vector3().subVectors(b[3],b[0]).normalize(),o=i.clone().negate(),s=b[0].distanceTo(b[1]),l=b[0].distanceTo(b[3]),u=s/t,d=l/n,f=.001;for(let r=0;r<u;r++){let r=b[0],o=0;for(let s=0;s<d;s++){let s=r.clone().addScaledVector(e,t),l=s.clone().addScaledVector(i,n),u=r.clone().addScaledVector(i,n),d=new c.Vector3().add(r).add(s).add(l).add(u).multiplyScalar(.25);for(let e=0;e<a.length;e++)if(a[e].distanceTo(d)<f){o++;break}r=u}if(o>=d/2)break;S+=o,b[0]=b[0].addScaledVector(e,t),b[3]=b[3].addScaledVector(e,t)}s=b[0].distanceTo(b[1]),u=s/t;for(let e=0;e<u;e++){let e=b[1],o=0;for(let s=0;s<d;s++){let s=e.clone().addScaledVector(r,t),l=s.clone().addScaledVector(i,n),u=e.clone().addScaledVector(i,n),d=new c.Vector3().add(e).add(s).add(l).add(u).multiplyScalar(.25);for(let e=0;e<a.length;e++)if(a[e].distanceTo(d)<f){o++;break}e=u}if(o>=d/2)break;S+=o,b[1]=b[1].addScaledVector(r,t),b[2]=b[2].addScaledVector(r,t)}s=b[0].distanceTo(b[1]),u=s/t;for(let r=0;r<d;r++){let r=b[3],i=0;for(let s=0;s<u;s++){let s=r.clone().addScaledVector(e,t),l=s.clone().addScaledVector(o,n),u=r.clone().addScaledVector(o,n),d=new c.Vector3().add(r).add(s).add(l).add(u).multiplyScalar(.25);for(let e=0;e<a.length;e++)if(a[e].distanceTo(d)<f){i++;break}r=s}if(i>=u/2)break;S+=i,b[2]=b[2].addScaledVector(o,n),b[3]=b[3].addScaledVector(o,n)}}let C=b[0].distanceTo(b[1])*b[0].distanceTo(b[3]);return{facePoints:b,boxArea:C,totalInPlaneNum:S}},b=(e,t,n,r,i=c.MathUtils.degToRad(10))=>{let a=new c.Vector3,o=new c.Vector3;a.subVectors(new c.Vector3(t.x,t.y,0),new c.Vector3(e.x,e.y,0)).normalize(),o.subVectors(new c.Vector3(r.x,r.y,0),new c.Vector3(n.x,n.y,0)).normalize();let s=a.dot(o),l=Math.acos(Math.min(Math.abs(s),1)),u=1-Math.min(l/i,1),d=s>0?`same`:`opposite`;return{parallelism:u,angle:l,angleDeg:c.MathUtils.radToDeg(l),isParallel:l<i,direction:d,vectors:{v1:a,v2:o}}},x=(e,t,n,r)=>{let i=new c.Vector3().subVectors(t,e),a=new c.Vector3().subVectors(r,n),o=i.length(),s=a.length();if(o===0||s===0)return{rate:NaN,angle:NaN,isPerpendicular:!1,isParallel:!1};i.normalize(),a.normalize();let l=Math.abs(i.dot(a)),u=Math.min(1,Math.max(0,l)),d=Math.acos(u),f=c.MathUtils.radToDeg(d),p=1-u;return{rate:p,percent:`${(p*100).toFixed(2)}%`,angle:f,isPerpendicular:Math.abs(f-90)<.01,isParallel:f<.01}},S=(e,t,n=!1)=>new c.Vector3((e.x+t.x)/2,(e.y+t.y)/2,n?0:(e.z+t.z)/2),C=(e,t,n,r,i={})=>{let{parallelAngleMax:a=c.MathUtils.degToRad(10),collinearAngleMax:o=c.MathUtils.degToRad(7),distanceThreshold:s=2,overlapThreshold:l=.01}=i,u=new c.Vector3().subVectors(t,e),d=new c.Vector3().subVectors(r,n),f=u.clone().normalize(),p=d.clone().normalize(),m=c.MathUtils.clamp(f.dot(p),-1,1),h=Math.acos(Math.abs(m)),g=c.MathUtils.radToDeg(h);if(h>a)return{type:`not_parallel`,angleDeg:g,avgPerpendicularDistance:null,maxPerpendicularDistance:null,gap:null,overlap:null,closestDistance:null};m<0&&p.negate();let _=new c.Vector3().addVectors(f,p).normalize(),v=new c.Vector3().addVectors(e,t).multiplyScalar(.5),y=w(n,v,f),b=w(r,v,f),x=(y+b)/2,S=Math.max(y,b),C=e.dot(_),E=t.dot(_),D=n.dot(_),O=r.dot(_),k=Math.min(C,E),A=Math.max(C,E),j=Math.min(D,O),M=Math.max(D,O),N=A-k,P=M-j,F=Math.min(N,P),I=Math.min(A,M)-Math.max(k,j),L=I<0?-I:0,R=I>0?I:0,z=F>0?R/F:0,B=T(e,t,n,r),V=h<=o&&x<s,H;return V?H=z>.5?`collinear_overlap`:`collinear_gap`:h<=a&&(H=`parallel_offset`),{type:H,angleDeg:g,avgPerpendicularDistance:x,maxPerpendicularDistance:S,gap:L,overlap:R,closestDistance:B}};function w(e,t,n){let r=new c.Vector3().subVectors(e,t),i=r.dot(n),a=n.clone().multiplyScalar(i);return new c.Vector3().subVectors(r,a).length()}function T(e,t,n,r){let i=new c.Vector3().subVectors(t,e),a=new c.Vector3().subVectors(r,n),o=new c.Vector3().subVectors(n,e),s=i.dot(i),l=i.dot(a),u=a.dot(a),d=i.dot(o),f=a.dot(o),p=s*u-l*l,m,h;p<1e-10?(m=0,h=u===0?0:f/u):(m=(l*f-u*d)/p,h=(s*f-l*d)/p),m=c.MathUtils.clamp(m,0,1),h=u===0?0:c.MathUtils.clamp((l*m+f)/u,0,1),m=s===0?0:c.MathUtils.clamp((l*h-d)/s,0,1);let g=e.clone().add(i.clone().multiplyScalar(m)),_=n.clone().add(a.clone().multiplyScalar(h));return g.distanceTo(_)}var E=(e,t,n,r=.05,i=.02)=>{let[a,o,s,l]=e,u=new c.Vector3().subVectors(o,a),d=new c.Vector3().subVectors(l,a),f=new c.Vector3().crossVectors(u,d).normalize(),p=u.length(),m=d.length(),h=u.clone().normalize(),g=d.clone().normalize(),_=n-r,v=n+r,y=new c.Vector3,b=[],x=Math.ceil(p/i),S=Math.ceil(m/i),C=p/x,w=m/S,T=new Uint8Array(x*S);for(let e=0;e<t.length;e++){y.subVectors(t[e],a);let n=Math.abs(y.dot(f));if(n<_||n>v)continue;let r=y.dot(h),i=y.dot(g);if(r<0||r>p||i<0||i>m)continue;b.push(new c.Vector3(t[e].x,t[e].y,t[e].z));let o=Math.min(Math.floor(r/C),x-1),s=Math.min(Math.floor(i/w),S-1);T[s*x+o]=1}let E=0;for(let e=0;e<T.length;e++)T[e]&&E++;let D=x*S,O=p*m,k=E==0||D==0?0:E/D*O;return{filteredPoints:b,coveragePercent:E==0||D==0?0:E/D*100,quadArea:O,coveredArea:k}},D=(e,t)=>Math.round(Math.max(5,Math.min(350,t*Math.sqrt(e)))),O=e=>{let t=new c.Vector3(e.start.x,e.start.y,e.start.z),n=new c.Vector3(e.end.x,e.end.y,e.end.z),r=t.distanceTo(n),i=e.rooftopPz-e.start.z,a=D(r,100),o=D(i,70),s=i/o,l=new Map;for(let e=0;e<o;e++)l.set((e+1)*s,{count:0,minZ:1/0,maxZ:-1/0});let u=new c.Line3(t,n),d=e.originalPoints,f=new c.Vector3;for(let e=0;e<d.length;e++){let t=d[e],n=u.closestPointToPoint(t,!0,f).distanceTo(t),r=Math.floor(n/s);if(r<o){let e=l.get((r+1)*s);e.count++,t.z<e.minZ&&(e.minZ=t.z),t.z>e.maxZ&&(e.maxZ=t.z)}}let p,m,h=a/2;for(let[e,t]of l)if(t.count>=h){p=t.minZ;break}for(let[e,t]of[...l.entries()].reverse())if(t.count>=h){m=t.maxZ;break}return{minZ:p,maxZ:m}},k=(e,t,n)=>{let r=new c.Plane,i=new c.Vector3(e.x,e.y,e.z),a=new c.Vector3(t.x,t.y,t.z),o=new c.Vector3(t.x,t.y,t.z+5);r.setFromCoplanarPoints(i,a,o);let s=[];n.forEach(e=>{let t=new c.Vector3;r.projectPoint(e,t),t.z=i.z,s.push(t)});let l=new c.Vector3,u=new c.Vector3;r.projectPoint(i,l),r.projectPoint(a,u);let d=new c.Line3(l,u),f=[],p=[],m=new c.Vector3;return s.forEach((e,t)=>{d.closestPointToPoint(e,!0,m).distanceTo(e)<.001?f.push(n[t]):p.push(e)}),{newOriginalPoints:f,newProjectPoints:p}},A=e=>{function t(e,t,n,r,i,a,o,s,l){let u=new c.Line3(e,t),d=new Map,f=new c.Vector3,p=e.clone(),m=t.clone(),h=0;for(;;){let g=new Map,_=0;for(let e=0;e<l.checkResults.length;e++)if(!(l.checkResults[e].allCenterPoints.length/s*100<1))for(let t=0;t<l.checkResults[e].allCenterPoints.length;t++){if(d.has(e)&&d.get(e).removePtsIndex.includes(t))continue;let n=l.checkResults[e].allCenterPoints[t];n=new c.Vector3(n.x,n.y,n.z);let r=u.closestPointToPoint(n,!0,f),a=n.distanceTo(r);Math.abs(a-i)<.01&&(_++,g.has(e)?g.get(e).removePtsIndex.push(t):g.set(e,{index:e,removePtsIndex:[t]}))}if(h==o){e=p,t=m;break}if(_<a/20*19||_==0){g.clear(),o-h<5&&(e=p,t=m,d.clear());break}else{h++;for(let[e,t]of g)d.has(e)?d.get(e).removePtsIndex.push(...t.removePtsIndex):d.set(e,t);g.clear(),e=e.addScaledVector(n,r),t=t.addScaledVector(n,r),u.set(e,t)}}if(d.size>0&&h!=o)for(let[e,t]of d){let n=[],r=[];for(let i=0;i<l.checkResults[e].allCenterPoints.length;i++)t.removePtsIndex.includes(i)||(n.push(l.checkResults[e].allCenterPoints[i]),r.push(l.checkResults[e].originalVertices[i]));l.checkResults[e].allCenterPoints=n,l.checkResults[e].originalVertices=r}}for(let n=0;n<e.length;n++){let r=e[n];if(r.length<.5)continue;let i=-1,a=-1;if(!r.checkResults||r.checkResults.length<=0)continue;if(r.checkResults[0].originalVertices&&r.checkResults[0].originalVertices.length>0){let[e,t,n,o]=r.checkResults[0].originalVertices[0];e=new c.Vector3(e.x,e.y,e.z),t=new c.Vector3(t.x,t.y,t.z),n=new c.Vector3(n.x,n.y,n.z),o=new c.Vector3(o.x,o.y,o.z),i=e.distanceTo(t),a=e.distanceTo(o);let s=h(new c.Vector3().subVectors(e,t).normalize(),new c.Vector3(0,0,1));s&&s.parallel&&([a,i]=[i,a])}if(i<0||a<0)continue;let o=new c.Vector3(r.start.x,r.start.y,r.start.z),s=new c.Vector3(r.end.x,r.end.y,r.end.z),l=o.distanceTo(s),u=r.rooftopPz-r.start.z,d=Math.ceil(l/i),f=Math.ceil(u/a),p=o.clone(),m=s.clone(),g=m.clone().add(new c.Vector3(0,0,u)),_=p.clone().add(new c.Vector3(0,0,u)),v=new c.Vector3().subVectors(m,p).normalize(),y=v.clone().negate();new c.Vector3().subVectors(g,p).normalize().clone().negate();let b=i/2;a/2;let x=l*u/(i*a);t(p,_,v,i,b,f,d,x,r),t(m,g,y,i,b,f,d,x,r),r.start=p,r.end=m;let{newOriginalPoints:S}=k(p,m,r.originalPoints);r.originalPoints=S}},j=(e,t)=>new c.Vector3(e.x,e.y,0).distanceTo(new c.Vector3(t.x,t.y,0)),M=e=>{let t=0,n=[];for(let r=0;r<e.length;r++)for(let i=r+1;i<e.length;i++){let a=e[r].distanceTo(e[i]);a>t&&(t=a,n=[e[r],e[i]])}return n},N=class{results;clusterResults;constructor(){this.results=null}initLimits1(e,t){let n,r;return n=e<.3?5:e<.5?8:e<.8?10:e<1?15:e<2?30:e<3?35:e<4?40:e<5?45:e<6?50:e<7?60:e<8.5?70:e<10?80:e<15?100:e<20?120:e<25?140:e<30?160:e<35?180:e<40?200:e<45?230:e<50?250:350,t<1?r=30:t<2?r=40:t<2.5?r=45:t<3?r=60:t<3.6?r=70:t<4&&(r=80),{horizontalSubdivisions:n,verticalSubdivisions:r}}initLimits2(e,t){return{horizontalSubdivisions:Math.round(Math.max(5,Math.min(350,e/.08))),verticalSubdivisions:Math.round(Math.max(30,Math.min(80,t/.06)))}}initLimits4(e,t){return{horizontalSubdivisions:Math.round(Math.max(5,Math.min(350,65*Math.sqrt(e))))*1,verticalSubdivisions:Math.round(Math.max(30,Math.min(80,15*Math.sqrt(t))))*1}}initLimits3(e,t){let n;return n=e<2?.06:e<10?.1:.15,{horizontalSubdivisions:Math.round(Math.max(5,Math.min(350,e/n))),verticalSubdivisions:Math.round(Math.max(30,Math.min(80,t/.06)))}}createWallPlaneMeshFromFourPoints(e,t,n,r,i,a,o=16777215){let s=[e,t,n,r];for(let e=0;e<s.length;e++)this.isValidVector3(s[e])||(console.error(`点${e+1}包含无效数据:`,s[e]),s[e]=new c.Vector3(e,0,e));let l=new c.BufferGeometry,u=new Float32Array([s[0].x,s[0].y,s[0].z,s[1].x,s[1].y,s[1].z,s[2].x,s[2].y,s[2].z,s[3].x,s[3].y,s[3].z]),d=[0,1,2,0,2,3],f=new Float32Array([0,0,1,0,1,1,0,1]);l.setAttribute(`position`,new c.BufferAttribute(u,3)),l.setAttribute(`uv`,new c.BufferAttribute(f,2)),l.setIndex(d),l.computeVertexNormals();let p=new c.Mesh(l,i),m=new c.EdgesGeometry(l),h=new c.LineSegments(m,new c.LineBasicMaterial({color:o}));return a.add(h),p}createSubdividedWallPlane(e,t,n,r,i){let a=[];if(!e||!e.start||!e.end||isNaN(e.rooftopPz))return a;let o=new c.Vector3(e.start.x,e.start.y,e.start.z),s=new c.Vector3(e.end.x,e.end.y,e.end.z),l=[o,s,new c.Vector3(e.end.x,e.end.y,e.rooftopPz),new c.Vector3(e.start.x,e.start.y,e.rooftopPz)];for(let e=0;e<l.length;e++)if(!this.isValidVector3(l[e]))return console.error(`墙体基础点${e}无效:`,l[e]),a;let u=o.distanceTo(s),d=Math.abs(e.rooftopPz-e.start.z);if(u===0||d===0)return console.error(`墙体宽度或高度为0`),console.log(`line.rooftopPz`,e.rooftopPz),console.log(`bottomLeft`,o),console.log(`bottomRight`,s),console.log(e.length),a;let f=u/n,p=d/t,m=new c.Vector3().subVectors(s,o).normalize(),h=new c.Vector3(0,0,1);for(let e=0;e<t;e++)for(let t=0;t<n;t++)try{let n={points:[new c.Vector3().copy(o).add(m.clone().multiplyScalar(t*f)).add(h.clone().multiplyScalar(e*p)),new c.Vector3().copy(o).add(m.clone().multiplyScalar((t+1)*f)).add(h.clone().multiplyScalar(e*p)),new c.Vector3().copy(o).add(m.clone().multiplyScalar((t+1)*f)).add(h.clone().multiplyScalar((e+1)*p)),new c.Vector3().copy(o).add(m.clone().multiplyScalar(t*f)).add(h.clone().multiplyScalar((e+1)*p))],userData:{}};n.userData={type:`wallSegment`,row:e,col:t,originalMaterial:r,width:f,height:p},a.push(n)}catch(n){console.error(`创建墙体小平面(${e}, ${t})时出错:`,n)}return a}calculateMinRequiredPoints(e){let t=Math.floor(Math.log1p(e)*.5);return Math.max(1,t)}calculateExpectedDensity(e){return e<1?10:e<4?5:2}evaluateSegmentByArea(e,t,n,r){if(t===0)return!1;let i=e.area,a=e.minRequiredPoints,o=e.expectedDensity,s=t/Math.max(i,.001),c=t>=a,l=s>=o*.3,u=!0;return i<.5&&t>0&&(u=this.checkSmallAreaDistribution(e,n,r)),c&&l&&u}checkSmallAreaDistribution(e,t,n){if(t.length<2)return!0;let r=t.map(e=>n[e]);e.bounds;let i=e.center,a=0,o=Math.sqrt(e.area)*.3;return r.forEach(e=>{e.distanceTo(i)<=o&&a++}),a/r.length>=.3}getFailureReason(e,t){let n=e.area,r=e.minRequiredPoints,i=e.expectedDensity,a=t/Math.max(n,.001),o=[];return t===0?o.push(`无点`):(t<r&&o.push(`点数不足: ${t}/${r}`),a<i*.3&&o.push(`密度不足: ${a.toFixed(2)}/${i.toFixed(2)}`),n<.5&&t>0&&(this.checkSmallAreaDistribution(e,[],[])||o.push(`小平面分布不均`))),o.length>0?o.join(`, `):`未知原因`}calculateMedian(e){let t=[...e].sort((e,t)=>e-t),n=Math.floor(t.length/2);return t.length%2==0?(t[n-1]+t[n])/2:t[n]}isPointInBoundingBox(e,t,n=.1){return e.x>=t.min[0]-n&&e.x<=t.max[0]+n&&e.y>=t.min[1]-n&&e.y<=t.max[1]+n&&e.z>=t.min[2]-n&&e.z<=t.max[2]+n}getSegmentPlane(e){let t=e.geometry;if(!t||!t.attributes.position)return null;let n=t.attributes.position;if(n.count<3)return null;let r=new c.Vector3().fromBufferAttribute(n,0),i=new c.Vector3().fromBufferAttribute(n,1),a=new c.Vector3().fromBufferAttribute(n,2);r.applyMatrix4(e.matrixWorld),i.applyMatrix4(e.matrixWorld),a.applyMatrix4(e.matrixWorld);let o=new c.Plane;return o.setFromCoplanarPoints(r,i,a),o}analyzePointDistribution(e,t,n=.01){this.results={totalSegments:0,segmentsWithPoints:0,segmentsWithoutPoints:0,missingSegments:[],pointsInSegments:new Map,segmentsByPoint:new Map,segmentDetails:new Map,tolerance:n,analysisTime:null};let r=performance.now(),i=[];e.forEach(e=>{if(e.userData&&e.userData.type===`wallSegment`){let t=this.getWorldVerticesNew(e.points);if(t.length<4)return;let n=new c.Vector3().add(t[0]).add(t[1]).add(t[2]).add(t[3]).multiplyScalar(.25);e.userData.center=n,e.vertices=t,i.push(e)}}),this.results.totalSegments=i.length;for(let e of i){if(e.vertices.length<4)continue;let r=e.userData.center,i=`row${e.userData.row}_col${e.userData.col}`,[a,o,s,l]=e.vertices,u=new c.Vector3().subVectors(o,a).normalize(),d=new c.Vector3().subVectors(l,a).normalize(),f=a.distanceTo(o),p=a.distanceTo(l),m=new c.Plane;m.setFromCoplanarPoints(a,o,l);let h=f/2,g=p/2;for(let e=0;e<t.length;e++)this.isPointInWallSegmentOptimized(t[e],r,m,a,u,d,f,p,h,g,n)&&(this.results.pointsInSegments.has(i)||this.results.pointsInSegments.set(i,[]),this.results.pointsInSegments.get(i).push(e))}return i.forEach(e=>{let t=`row${e.userData.row}_col${e.userData.col}`,n=this.results.pointsInSegments.has(t)?this.results.pointsInSegments.get(t).length:0;n>=1?this.results.segmentsWithPoints++:(this.results.segmentsWithoutPoints++,this.results.missingSegments.push({row:e.userData.row,col:e.userData.col,width:e.userData.width,height:e.userData.height,segmentKey:t,vertices:e.vertices,center:e.userData.center,mesh:e,pointCount:n,hasPointsButNotEnough:n>0&&n<1}))}),this.results.analysisTime=performance.now()-r,this.printAnalysisResults(),this.results}isPointInWallSegmentOptimized(e,t,n,r,i,a,o,s,l,u,d){if(Math.abs(n.distanceToPoint(e))>d)return!1;let f=new c.Vector3().subVectors(e,t),p=Math.abs(f.dot(i)),m=Math.abs(f.dot(a));if(p>l&&m>u)return!1;let h=new c.Vector3().subVectors(e,r),g=h.dot(i),_=h.dot(a),v=g>=-d&&g<=o+d,y=_>=-d&&_<=s+d;return v&&y}isPointInWallSegment(e,t,n,r=.1){if(!e.geometry)return!1;try{let i=this.getWorldVertices(e);return i.length<4?!1:this.isPointInRotatedWall(i,t,n,r)}catch(e){return console.error(`判断点是否在墙体平面内时出错:`,e),!1}}isPointInRotatedWall(e,t,n,r=.1){let[i,a,o,s]=e,l=new c.Plane;if(l.setFromCoplanarPoints(i,a,s),Math.abs(l.distanceToPoint(t))>r)return!1;let u=new c.Vector3().subVectors(a,i).normalize(),d=new c.Vector3().subVectors(s,i).normalize();new c.Vector3().crossVectors(u,d).normalize();let f=new c.Vector3().subVectors(t,i),p=f.dot(u),m=f.dot(d),h=i.distanceTo(a),g=i.distanceTo(s),_=p>=-r&&p<=h+r,v=m>=-r&&m<=g+r;return _&&v&&n.distanceTo(t)<h/2&&n.distanceTo(t)<g/2}getWorldVertices(e){if(e._cachedWorldVertices&&e.matrixWorld.equals(e._cachedMatrixWorld))return e._cachedWorldVertices;let t=e.geometry.getAttribute(`position`),n=[];for(let r=0;r<Math.min(4,t.count);r++){let i=new c.Vector3;i.fromBufferAttribute(t,r),i.applyMatrix4(e.matrixWorld),n.push(i)}return e._cachedWorldVertices=n,e._cachedMatrixWorld=e.matrixWorld.clone(),n}getWorldVerticesNew(e){let t=[];for(let n=0;n<e.length;n++){let r=e[n].clone();t.push(r)}return t}getWallSegmentBounds(e){if(!e.geometry)return null;let t=e.geometry.getAttribute(`position`),n=new c.Box3;for(let r=0;r<t.count;r++){let i=new c.Vector3;i.fromBufferAttribute(t,r),i.applyMatrix4(e.matrixWorld),n.expandByPoint(i)}return{min:n.min.toArray(),max:n.max.toArray(),center:n.getCenter(new c.Vector3).toArray(),size:n.getSize(new c.Vector3).toArray()}}highlightEmptySegments(e,t,n,r=1){if(!this.results){console.warn(`请先运行分析函数`);return}let i=new Map,a=new Map;this.results.missingSegments.forEach(e=>{i.set(e.segmentKey,e),a.has(e.row)||a.set(e.row,new Map),a.get(e.row).set(e.col,{segmentKey:e.segmentKey,isEmpty:!0,visited:!1})}),e.forEach(e=>{if(e.userData&&e.userData.type===`wallSegment`){let t=`row${e.userData.row}_col${e.userData.col}`,n=i.has(t);a.has(e.userData.row)||a.set(e.userData.row,new Map),a.get(e.userData.row).has(e.userData.col)||a.get(e.userData.row).set(e.userData.col,{segmentKey:t,isEmpty:n,visited:!1})}});let o=[],s=[[0,1],[1,0],[0,-1],[-1,0]];for(let[e,t]of i){let e=a.get(t.row)?.get(t.col);if(e&&!e.visited){let n=[],c=[[t.row,t.col]];for(e.visited=!0;c.length>0;){let[e,t]=c.shift(),r=`row${e}_col${t}`;i.has(r)&&n.push(i.get(r));for(let[n,r]of s){let i=e+n,o=t+r;if(a.has(i)&&a.get(i).has(o)){let e=a.get(i).get(o);!e.visited&&e.isEmpty&&(e.visited=!0,c.push([i,o]))}}}n.length>=r&&o.push(n)}}let c=this.createClusterMaterials(o.length),l=[];return this.clusterResults={clusters:o.sort((e,t)=>t.length-e.length),highlightedMeshes:l,clusterMaterials:c,totalClusters:o.length,totalHighlighted:l.length},this.clusterResults}calculateBoundaryRegularity(e,t,n,r,i){let a=0,o=0,s=t,c=!0;for(let t=r;t<=i;t++)e.has(s)&&e.get(s).has(t)?o++:c=!1;c&&a++;let l=n,u=!0;for(let t=r;t<=i;t++)e.has(l)&&e.get(l).has(t)?o++:u=!1;u&&a++;let d=r,f=!0;for(let r=t;r<=n;r++)e.has(r)&&e.get(r).has(d)?o++:f=!1;f&&a++;let p=i,m=!0;for(let r=t;r<=n;r++)e.has(r)&&e.get(r).has(p)?o++:m=!1;return m&&a++,a/4}calculateClusterQualityScore(e){let t=1/Math.max(1,e.aspectRatio);return e.compactness*.4+e.boundaryRegularity*.4+t*.2}calculateStrictClusterBounds(e){if(e.length===0)return{width:0,height:0};let t=[];if(e.forEach(e=>{e.vertices&&e.vertices.length>0&&t.push(...e.vertices)}),t.length===0){let t=[...new Set(e.map(e=>e.row))];return{width:[...new Set(e.map(e=>e.col))].length*(e[0].bounds?.width||1),height:t.length*(e[0].bounds?.height||1)}}let n=Math.min(...t.map(e=>e.x)),r=Math.max(...t.map(e=>e.x)),i=Math.min(...t.map(e=>e.y)),a=Math.max(...t.map(e=>e.y)),o=Math.min(...t.map(e=>e.z)),s=Math.max(...t.map(e=>e.z));return{width:Math.abs(r-n),height:Math.max(Math.abs(a-i),Math.abs(s-o))}}createClusterMaterials(e){let t=[],n=[16711935,65535,16753920,16738740,9055202,3329330,`lavenderblush`,`lawngreen`,`lemonchiffon`,`lightblue`,`lightcoral`];for(let r=0;r<e;r++){let e=n[r%n.length],i=new c.MeshBasicMaterial({color:e,transparent:!0,opacity:.6,side:c.DoubleSide});t.push(i)}return t}clearHighlights(e){this.clusterResults&&this.clusterResults.highlightedMeshes&&(this.clusterResults.highlightedMeshes.forEach(t=>{e.remove(t),t.material&&t.material.dispose()}),this.results&&this.results.segmentDetails&&this.results.segmentDetails.forEach((e,t)=>{e.mesh&&e.mesh.userData.originalMaterial&&(e.mesh.material=e.mesh.userData.originalMaterial)}),this.clusterResults.highlightedMeshes=[])}getClusterInfo(){return this.clusterResults?{totalClusters:this.clusterResults.totalClusters,clusterSizes:this.clusterResults.clusters.map((e,t)=>({clusterIndex:t,size:e.length,segments:e.map(e=>e.segmentKey)})),largestCluster:Math.max(...this.clusterResults.clusters.map(e=>e.length)),smallestCluster:Math.min(...this.clusterResults.clusters.map(e=>e.length))}:null}restoreOriginalMaterials(e){let t=0;e.traverse(e=>{e.userData&&e.userData.originalMaterial&&(e.material=e.userData.originalMaterial,delete e.userData.originalMaterial,t++)}),console.log(`恢复了 ${t} 个平面的原始材质`)}createDebugHelpers(e){let t=new c.Group;return e.traverse(e=>{if(e.userData&&e.userData.type===`wallSegment`){let n=this.getWorldVertices(e);if(n.length<4)return;let r=new c.Vector3().add(n[0]).add(n[1]).add(n[2]).add(n[3]).multiplyScalar(.25),i=new c.Vector3().subVectors(n[1],n[0]).normalize(),a=new c.Vector3().subVectors(n[3],n[0]).normalize(),o=new c.Vector3().crossVectors(i,a).normalize(),s=new c.ArrowHelper(o,r,.3,16711680),l=new c.ArrowHelper(i,r,.2,65280),u=new c.ArrowHelper(a,r,.2,255);t.add(s),t.add(l),t.add(u)}}),t}printAnalysisResults(){if(!this.results){console.warn(`没有可用的分析结果`);return}this.results.missingSegments.length>0&&this.results.missingSegments.forEach(e=>{}),this.results.pointsInSegments.forEach((e,t)=>{})}isValidVector3(e){return e&&!isNaN(e.x)&&!isNaN(e.y)&&!isNaN(e.z)&&isFinite(e.x)&&isFinite(e.y)&&isFinite(e.z)}mergeWallGroupMeshes(e){let t=[];if(e.traverse(e=>{e.isMesh&&e.geometry&&t.push(e)}),t.length===0)return null;try{let e=t.map(e=>e.geometry);e.forEach((e,n)=>{e.applyMatrix4(t[n].matrixWorld)});let n=l(e),r=t[0].material.clone(),i=new c.Mesh(n,r);return i.name=`mergedWalls`,i}catch(e){return console.error(`合并网格时出错:`,e),null}}},P=e=>{let t={},n=[],r=[],i=[];return e[0].width,e[0].height,e.forEach(e=>{r.push(...e.vertices),i.push(e.center),n.push(e.vertices)}),t.originalVertices=n,t.allCenterPoints=i,t},F=e=>{let t=new c.Vector3().copy(e.start).add(e.end).multiplyScalar(.5),n=new c.Vector3().subVectors(new c.Vector3(t.x,t.y,t.z+1),t).normalize(),r=new c.Vector3(e.direction.x,e.direction.y,0);r.applyAxisAngle(n,Math.PI/2);let i=new c.Vector3().copy(t).add(r.clone().multiplyScalar(e.length)),a=new c.Vector3().copy(t).sub(r.clone().multiplyScalar(e.length));return t.z=e.start.z,i.z=e.start.z,a.z=e.start.z,{start:i,end:a,center:t,direction:r}},I=(e,t,n=!0)=>{let r=[];if(!e)return[];for(let i=0;i<e.length;i++){let a=e[i];if(!a)continue;let o=new c.Vector3(a.start.x,a.start.y,a.start.z);o.applyEuler(new c.Euler(-Math.PI/2,0,0));let s=new c.Vector3(a.end.x,a.end.y,a.end.z);if(s.applyEuler(new c.Euler(-Math.PI/2,0,0)),t.add(p(String(i),new c.Vector3((o.x+s.x)/2,(o.y+s.y)/2,(o.z+s.z)/2))),(a.isBayWindow||a.isWindow)&&n)continue;let l=new c.Plane().setFromCoplanarPoints(new c.Vector3(a.start.x,a.start.y,a.start.z),new c.Vector3(a.end.x,a.end.y,a.end.z),new c.Vector3(a.start.x,a.start.y,a.rooftopPz)),u=[];a.originalPoints.forEach((e,t)=>{let n=new c.Vector3;l.projectPoint(e,n),u.push(n)});let d=new c.Vector3(a.start.x,a.start.y,a.start.z).distanceTo(new c.Vector3(a.end.x,a.end.y,a.end.z)),f=Math.abs(a.rooftopPz-a.start.z),m=new c.MeshBasicMaterial({color:`red`,side:c.DoubleSide,wireframe:!1,transparent:!0,opacity:.5}),h=new c.Group,g=new N,{horizontalSubdivisions:_,verticalSubdivisions:v}=g.initLimits4(d,f),y=g.createSubdividedWallPlane(a,v,_,m,h);g.analyzePointDistribution(y,u,.01);let b=new c.MeshBasicMaterial({color:`black`,transparent:!0,opacity:.8,side:c.DoubleSide}),x=g.highlightEmptySegments(y,b,h),S=[];for(let e of x.clusters){let t=P(e);t&&S.push(t)}let C=F(a).direction;a.checkResults=S,a.verticalDirection=C,a.originaIndex=i,r.push(a)}return r},L=async(e,t,n,r)=>{let i=await fetch(e);if(!i)return;let a=await i.json(),o=[a[31],a[59]];for(let e=0;e<o.length;e++){let t=o[e];if(!t)continue;let n=new c.Vector3(t.start.x,t.start.y,t.start.z);n.applyEuler(new c.Euler(-Math.PI/2,0,0));let i=new c.Vector3(t.end.x,t.end.y,t.end.z);i.applyEuler(new c.Euler(-Math.PI/2,0,0)),r.add(p(String(e),new c.Vector3((n.x+i.x)/2,(n.y+i.y)/2,(n.z+i.z)/2)))}o.sort((e,t)=>e.originalIndex-t.originalIndex),console.log(`lins`,o);let s=await fetch(t);if(!s)return;let u=await s.json(),d=Array.isArray(u)?u:Object.values(u),f=await fetch(n);if(console.log(`respcd`,f),!f)return;let m=await f.json();A(o,r),console.time(),o=await R(o,d,m,r),console.timeEnd();let h=new c.MeshBasicMaterial({color:`#cffd00`,transparent:!0,opacity:.4,side:c.DoubleSide}),g=new c.MeshBasicMaterial({color:`#ff0000`,transparent:!0,opacity:.4,side:c.DoubleSide});new c.LineBasicMaterial({color:`#cffd00`});let _=new c.PointsMaterial({color:65535,size:.01}),y=[],b=[];for(let e of o){if(!e.originalPoints||e.originalPoints.length===0)continue;let t=new Float32Array(e.originalPoints.length*3);e.originalPoints.forEach((e,n)=>{t[n*3]=e.x,t[n*3+1]=e.y,t[n*3+2]=e.z});let n=new c.BufferGeometry;n.setAttribute(`position`,new c.BufferAttribute(t,3)),y.push(n);let r=[];r.push(new c.Vector3(e.start.x,e.start.y,e.start.z)),r.push(new c.Vector3(e.end.x,e.end.y,e.end.z));let i=new c.BufferGeometry;i.setFromPoints(r),b.push(i)}if(y.length>0){let e=l(y,!1);e.rotateX(-Math.PI/2),r.add(new c.Points(e,_)),y.forEach(e=>e.dispose())}let x=[],S=[],C=new c.LineBasicMaterial({color:`#f30606`}),w=new c.MeshBasicMaterial({color:`#071ac4`});for(let e=0;e<d.length;e++){let t=new c.SphereGeometry(.05);t.translate(d[e].x,d[e].y,d[e].z);let n=[],r=new c.Vector3(d[e].x,d[e].y,d[e].z),i=new c.Quaternion(d[e].qx,d[e].qy,d[e].qz,d[e].qw),a=new c.Vector3(0,0,-1).applyQuaternion(i),o=r.clone().addScaledVector(a,.001);n.push(r),n.push(o);let s=new c.BufferGeometry;s.setFromPoints(n),S.push(s),x.push(t)}if(x.length>0){let e=l(x,!1);e.rotateX(-Math.PI/2),r.add(new c.Mesh(e,w)),x.forEach(e=>e.dispose())}if(S.length>0){let e=l(S,!1);e.rotateX(-Math.PI/2),r.add(new c.Line(e,C)),S.forEach(e=>e.dispose())}let T=v(m);T.rotateX(-Math.PI/2),r.add(T),console.log(`jsonpcdData`,m),o.forEach(e=>{if(e.doorAndBeamData&&e.doorAndBeamData.length>0){let t=new c.Group;e.doorAndBeamData.forEach(e=>{if(e.beamStart){let n=new c.BufferGeometry,r=new Float32Array([e.beamStart.x,e.beamStart.y,e.beamStart.z,e.beamEnd.x,e.beamEnd.y,e.beamEnd.z,e.beamEnd.x,e.beamEnd.y,e.beamEnd.z+e.beamHeight,e.beamStart.x,e.beamStart.y,e.beamStart.z+e.beamHeight]);n.setAttribute(`position`,new c.BufferAttribute(r,3)),n.setIndex([0,1,2,0,2,3]),n.computeVertexNormals(),n.rotateX(-Math.PI/2),t.add(new c.Mesh(n,h))}let n=new c.BufferGeometry,r=new Float32Array([e.doorStart.x,e.doorStart.y,e.doorStart.z,e.doorEnd.x,e.doorEnd.y,e.doorEnd.z,e.doorEnd.x,e.doorEnd.y,e.doorEnd.z+e.doorHeight,e.doorStart.x,e.doorStart.y,e.doorStart.z+e.doorHeight]);n.setAttribute(`position`,new c.BufferAttribute(r,3)),n.setIndex([0,1,2,0,2,3]),n.computeVertexNormals(),n.rotateX(-Math.PI/2),t.add(new c.Mesh(n,g));let i=`id:${e.id},Nid:${e.nearId}`,a=new c.Vector3((e.doorStart.x+e.doorEnd.x)/2,(e.doorStart.y+e.doorEnd.y)/2,e.doorStart.z+e.doorHeight);a.applyEuler(new c.Euler(-Math.PI/2,0,0)),t.add(p(i,a))}),r.add(t)}})},R=async(e,t,n)=>{if(e.length<=0||t.length<=0)return e;n||=[];for(let t of e){let e=[];if(!(!t.checkResults||t.checkResults.length<=0)){for(let n=0;n<t.checkResults.length;n++){let r=t.checkResults[n];if(r.isDoor=!1,!r.originalVertices||r.originalVertices.length==0)continue;let i=[];r.originalVertices.forEach(e=>{let t=[];e.forEach(e=>{t.push(new c.Vector3(e.x,e.y,e.z))}),i.push(t)});let a=i.length,o=Array.from({length:a},(e,t)=>t);function s(e){for(;o[e]!==e;)o[e]=o[o[e]],e=o[e];return e}for(let e=0;e<a;e++)for(let t=e+1;t<a;t++){let n=i[e],r=i[t];if(n.some(e=>r.some(t=>e.equals(t)))){let n=s(e),r=s(t);n!==r&&(o[n]=r)}}let l=new Map,u=-1,d=-1;for(let e=0;e<a;e++){let t=s(e);l.has(t)||l.set(t,{points:[],area:0,originalVertices:[],index:-1,gridWidth:-1,gridHeight:-1,centerPts:[]});let r=l.get(t),[a,o,f,p]=i[e],m=new c.Vector3().add(a).add(o).add(f).add(p).multiplyScalar(.25);if(r.points.push(...i[e]),r.area+=a.distanceTo(o)*a.distanceTo(p),r.originalVertices.push(i[e]),r.index=n,r.centerPts.push(m),u==-1||d==-1){let e=new c.Vector3().subVectors(a,o).normalize();d=a.distanceTo(p),u=a.distanceTo(o);let t=h(e,new c.Vector3(0,0,1));t&&t.parallel&&([d,u]=[u,d])}r.gridWidth=u,r.gridHeight=d}e.push(...l.values())}t.mergeCheckRegion=e}}let r=new Map;for(let t of e){let e=t.rooftopPz-t.start.z,n=Math.floor(t.rooftopPz-t.start.z);r.has(n)?r.set(n,{num:r.get(n).num+1,totalHeight:r.get(n).totalHeight+e}):r.set(n,{num:1,totalHeight:e})}let i=[...r.entries()].reduce((e,t)=>t[1].num>e[1].num?t:e),a=i[0];i[1].totalHeight/i[1].num;let o=[];for(let e=0;e<n.length;e++){let t=n[e];if(t.isFindBeam=!1,t.category==`door`){let n=[];t.coordinatesByArea.coordinates.forEach(e=>{n.push(new c.Vector3(e.x,e.y,e.z))}),n[0].equals(n[n.length-1])&&n.pop();let r=n[0].distanceTo(n[1]),i=n[0].distanceTo(n[3]);i>r&&([r,i]=[i,r],n=[n[0],n[3],n[2],n[1]]);let a=S(n[0],n[3],!0),s=S(n[1],n[2],!0);o.push({doorStartPt:a,doorEndPt:s,boxPoints:n,minZ:t.coordinatesByArea.heightData.minZ,maxZ:t.coordinatesByArea.heightData.maxZ,index:e,isFind:!1})}}for(let t=0;t<o.length;t++){let n=o[t];for(let t=0;t<e.length;t++){if(e[t].length<.8)continue;let r=[],i=S(e[t].start,e[t].end,!0);r.push(i),r.push(S(i,e[t].start,!0)),r.push(S(i,e[t].end,!0));let a=!1,o=n.boxPoints.length;for(let e=0;e<r.length;e++){let t=r[e].x,i=r[e].y;for(let e=0,r=o-1;e<o;r=e++){let o=n.boxPoints[e].x,s=n.boxPoints[e].y,c=n.boxPoints[r].x,l=n.boxPoints[r].y;s>i!=l>i&&t<(c-o)*(i-s)/(l-s)+o&&(a=!a)}if(a)break}if(a){let r=new c.Vector3(e[t].start.x,e[t].start.y,0),i=new c.Vector3(e[t].end.x,e[t].end.y,0),a=n.doorStartPt.distanceTo(n.doorEndPt),o=r.distanceTo(i),s=b(n.doorStartPt,n.doorEndPt,r,i);if(Math.abs(o-a)<.5&&s&&s.isParallel){n.doorStartPt=new c.Vector3(e[t].start.x,e[t].start.y,0),n.doorEndPt=new c.Vector3(e[t].end.x,e[t].end.y,0);break}}else continue}}let s=0,l=[];new c.PointsMaterial({color:65535,size:.01});let u=[],d=[];for(let r=0;r<e.length;r++){let i=e[r];if(!i.mergeCheckRegion||i.mergeCheckRegion.length==0)continue;i.doorAndBeamData=[],i.completePointAreaPercentage=-1;let p=new c.Box3;p.setFromPoints(i.originalPoints);let m=p.max.z,g=p.min.z,_=new c.Vector3(i.start.x,i.start.y,i.start.z),v=new c.Vector3(i.end.x,i.end.y,i.end.z),w=i.rooftopPz-i.start.z;if(w<a-.5||w>a+2)continue;let T=0;for(let a of i.mergeCheckRegion){let{facePoints:p,boxArea:E,totalInPlaneNum:D}=y(a.points,a.gridWidth,a.gridHeight,a.centerPts);if(!p||p.length<=0)continue;p[0].distanceTo(p[1]);let O=p[0].distanceTo(p[3]),A=a.gridHeight*(O/a.gridHeight/2)+a.gridHeight/2,N=new c.Vector3(p[0].x,p[0].y,p[0].z+A),P=new c.Vector3(p[1].x,p[1].y,p[1].z+A),F=new c.Line3(N,P),I=[],L=new c.Vector3;a.points.forEach(e=>{F.closestPointToPoint(e,!0,L).distanceTo(e)<a.gridHeight+.01&&I.push(e)});let R=[],z=[];for(;;){let e=[];for(let t=0;t<I.length;t++)z.includes(t)||e.length==0&&(e.push(I[t]),z.push(t));for(let t=0;t<I.length;t++)if(!z.includes(t)){for(let n=0;n<e.length;n++)if(I[t].distanceTo(e[n])<a.gridWidth+.01){e.push(I[t]),z.push(t),t=-1;break}}if(e.length==0)break;R.push(e)}let B=[];for(let e=0;e<R.length;e++){let t={lineSt:new c.Vector3,lineEd:new c.Vector3,facePoints:[]},n=M(R[e]);Array.isArray(n)&&n.length==2&&(t.lineSt=n[0],t.lineEd=n[1],t.facePoints.push(new c.Vector3(n[0].x,n[0].y,p[0].z)),t.facePoints.push(new c.Vector3(n[1].x,n[1].y,p[0].z)),t.facePoints.push(new c.Vector3(n[1].x,n[1].y,p[2].z)),t.facePoints.push(new c.Vector3(n[0].x,n[0].y,p[2].z)),B.push(t))}T+=a.area;let V=0;if(D!=0&&(V=a.gridHeight*a.gridWidth*D),(a.area-V)/E*100>85&&B.length==1||B.length>1)for(let p=0;p<B.length;p++){let y=B[p],T=y.facePoints[0].distanceTo(y.facePoints[1]),E=y.facePoints[0].distanceTo(y.facePoints[3]),D=T/i.length*100,O=i.rooftopPz-y.facePoints[2].z;if(g>i.start.z-.15&&i.length<1.5||i.length,E<.5||m<y.facePoints[2].z||O>w/3*2)continue;let A=!1,M=!1,N=!1,P=[],F=y.facePoints[2].clone(),I=y.facePoints[3].clone(),L=y.facePoints[0].clone(),R=y.facePoints[1].clone(),z=!1;if(Math.abs(y.facePoints[0].z-i.start.z)<.25){if(T>.2&&i.length>.5){let t=j(_,L),n=j(_,R),a=j(v,L),o=j(v,R),s=.1;if(t<s||n<s){let{newOriginalPoints:n}=k(L,R,i.originalPoints),a=[],o=[];if(n.forEach(e=>{e.z>=F.z?a.push(e):e.z<=L.z&&o.push(e)}),a.length>0){let n=new c.Vector3(i.start.x,i.start.y,0),a=new c.Vector3(i.end.x,i.end.y,0),o=new c.Vector3().subVectors(n,a).normalize(),l=new c.Vector3,u=1/0,d=-1;for(let t=0;t<e.length;t++){if(t==r)continue;let i=new c.Vector3(e[t].start.x,e[t].start.y,0),s=new c.Vector3(e[t].end.x,e[t].end.y,0);if(f(n,a,i,s)!=null)continue;let p=x(n,a,i,s);if(p&&p.angle>85){let e=b(n,a,new c.Line3(i,s).closestPointToPoint(n,!0,l),n);if(!e||e.direction!=`same`)continue;let r=n.clone().addScaledVector(o,10),p=a.clone().addScaledVector(o.clone().negate(),10),m=new c.Vector3().subVectors(s,i).normalize(),h=f(r,p,i.clone().addScaledVector(m.clone().negate(),.15),s.clone().addScaledVector(m,.15));if(h==null)continue;let g=n.distanceTo(h.point);g<u&&(u=g,d=t)}else if(p&&p.angle<5){let e=j(n,i),t=j(n,s),r=j(a,i),o=j(a,s);if(e<.1||t<.1||r<.1||o<.1){u=1/0,d=-1;break}}}u!=1/0&&u<2&&e[d].length>.5&&(z=!0,t<s?(I=I.addScaledVector(o,u),L=L.addScaledVector(o,u)):(F=F.addScaledVector(o,u),R=R.addScaledVector(o,u)))}}else if(a<s||o<s){let{newOriginalPoints:t}=k(L,R,i.originalPoints),n=[],o=[];if(t.forEach(e=>{e.z>=F.z?n.push(e):e.z<=L.z&&o.push(e)}),n.length>0){let t=new c.Vector3(i.start.x,i.start.y,0),n=new c.Vector3(i.end.x,i.end.y,0),o=new c.Vector3().subVectors(n,t).normalize(),l=new c.Vector3,u=1/0,d=-1;for(let i=0;i<e.length;i++){if(i==r)continue;let a=new c.Vector3(e[i].start.x,e[i].start.y,0),s=new c.Vector3(e[i].end.x,e[i].end.y,0);if(f(t,n,a,s)!=null)continue;let p=x(t,n,a,s);if(p&&p.angle>85){let e=b(n,t,new c.Line3(a,s).closestPointToPoint(n,!0,l),n);if(!e||e.direction!=`same`)continue;let r=t.clone().addScaledVector(o,10),p=n.clone().addScaledVector(o.clone().negate(),10),m=new c.Vector3().subVectors(s,a).normalize(),h=f(r,p,a.clone().addScaledVector(m.clone().negate(),.15),s.clone().addScaledVector(m,.15));if(h==null)continue;let g=n.distanceTo(h.point);g<u&&(u=g,d=i)}if(p&&p.angle<5){let e=j(t,a),r=j(t,s),i=j(n,a),o=j(n,s);if(e<.1||r<.1||i<.1||o<.1){u=1/0,d=-1;break}}}u!=1/0&&u<2&&e[d].length>.5&&(z=!0,a<s?(I=I.addScaledVector(o,u),L=L.addScaledVector(o,u)):(F=F.addScaledVector(o,u),R=R.addScaledVector(o,u)))}}T=I.distanceTo(F)}if(D<85&&T>.35||T>.7&&E>w-w/3){let e=!1;for(let n=0;n<t.length;n++){let r=n+1;if(r>=t.length)continue;let i=new c.Vector3(t[n].x,t[n].y,t[n].z),a=new c.Vector3(t[r].x,t[r].y,t[r].z);if(f(new c.Vector3(i.x,i.y,0),new c.Vector3(a.x,a.y,0),new c.Vector3(y.facePoints[0].x,y.facePoints[0].y,0),new c.Vector3(y.facePoints[1].x,y.facePoints[1].y,0),.1)!=null){e=!0;break}}for(let e=0;e<o.length;e++){let t=new c.Vector3(o[e].doorStartPt.x,o[e].doorStartPt.y,0),r=new c.Vector3(o[e].doorEndPt.x,o[e].doorEndPt.y,0),a=new c.Vector3(y.facePoints[0].x,y.facePoints[0].y,0),s=new c.Vector3(y.facePoints[1].x,y.facePoints[1].y,0),l=new c.Vector3().subVectors(t,r).normalize(),u=new c.Vector3().subVectors(a,s).normalize(),d=.2,p=f(t.clone().addScaledVector(l,d),r.clone().addScaledVector(l.clone().negate(),d),a.clone().addScaledVector(u,d),s.clone().addScaledVector(u.clone().negate(),d)),m=y.facePoints[2].z<o[e].maxZ||Math.abs(y.facePoints[2].z-o[e].maxZ)<.3,h=b(t.clone(),r.clone(),a.clone(),s.clone());if(p&&m)n[o[e].index].isFindBeam=!0,M=!0,o[e].isFind=!0;else if(m&&h&&h.isParallel){let a=S(t,r),s=new c.Line3(new c.Vector3(i.start.x,i.start.y,0),new c.Vector3(i.end.x,i.end.y,0)),l=new c.Vector3,u=s.closestPointToPoint(a,!0,l);a.distanceTo(u)<.3&&(n[o[e].index].isFindBeam=!0)}}for(let e=0;e<o.length;e++){let t=new c.Vector3(o[e].doorStartPt.x,o[e].doorStartPt.y,0),r=new c.Vector3(o[e].doorEndPt.x,o[e].doorEndPt.y,0),a=new c.Vector3(y.facePoints[0].x,y.facePoints[0].y,0),s=new c.Vector3(y.facePoints[1].x,y.facePoints[1].y,0),l=y.facePoints[2].z<o[e].maxZ||Math.abs(y.facePoints[2].z-o[e].maxZ)<.3,u=Math.abs(y.facePoints[2].z-o[e].maxZ)>.2;if(l&&!u){let l=C(t.clone(),r.clone(),a.clone(),s.clone(),{distanceThreshold:2});if(l&&l.type==`collinear_gap`){let u=S(t,r),d=new c.Line3(new c.Vector3(i.start.x,i.start.y,0),new c.Vector3(i.end.x,i.end.y,0)),f=new c.Vector3,p=d.closestPointToPoint(u,!0,f);if(u.distanceTo(p)<.15){let i=a.distanceTo(t),u=s.distanceTo(t);if(i>.25&&u>.25){let e=a.distanceTo(r),t=s.distanceTo(r);if(e>.2&&t>.2)continue}M=!0,N=!0,P.push(o[e].index),o[e].isFind=!0,n[o[e].index].isFindBeam=!0;let d=new c.Vector3().subVectors(a,s).normalize(),f=t.distanceTo(r)-l.overlap;u<i&&d.negate(),i>u?(F=y.facePoints[2].clone().addScaledVector(d,f),I=y.facePoints[3].clone(),L=F.clone(),R=I.clone(),L.z=R.z=y.facePoints[0].z):(F=y.facePoints[2].clone(),I=y.facePoints[3].clone().addScaledVector(d,f),L=F.clone(),R=I.clone(),L.z=R.z=y.facePoints[0].z)}}else if(l&&l.type==`collinear_overlap`&&l.maxPerpendicularDistance<.35){let i=a.distanceTo(s),c=r.distanceTo(t);(Math.abs(i-l.overlap)<.25||Math.abs(c-l.overlap)<.25)&&(M=!0,N=!0,P.push(o[e].index),o[e].isFind=!0,n[o[e].index].isFindBeam=!0)}}}(M||e||E>w/2&&T>.8&&Math.abs(g-i.start.z)<.4)&&(A=!0)}}if(A){let t=!1;for(let e=0;e<d.length;e++){let n=C(d[e].start,d[e].end,L,R);if(n&&n.type==`collinear_overlap`&&n.maxPerpendicularDistance<.01){t=!0;break}}if(t)continue;if(d.push({start:L,end:R}),!N){let e=new c.Vector3().subVectors(L,R),t=new c.Vector3().subVectors(R,L),n=1/0,r=1/0,a=L.clone(),o=R.clone(),s=[];if(i.originalPoints.forEach(e=>{e.z-L.z<E/2+.1&&e.z-L.z>E/2-.1&&s.push(e)}),s.length>0){let{newProjectPoints:i}=k(L,R,s);for(let e of i)e.z=L.z;i.forEach(i=>{let s=new c.Vector3(i.x,i.y,L.z),l=new c.Vector3().subVectors(s,R),u=h(e,new c.Vector3().subVectors(s,L)),d=h(t,l);if(u&&u.parallel&&u.sameDirection){let e=s.distanceTo(L);e<n&&(n=e,a=s)}else if(d&&d.parallel&&d.sameDirection){let e=s.distanceTo(R);e<r&&(r=e,o=s)}})}a.equals(L)||(L=a,I=new c.Vector3(L.x,L.y,I.z)),o.equals(R)||(R=o,F=new c.Vector3(R.x,R.y,F.z))}u.includes(r)||u.push(r),i.checkResults[a.index].isDoor=!0;let n={id:s,beamStart:F,beamEnd:I,beamHeight:i.rooftopPz-y.facePoints[2].z,doorStart:L,doorEnd:R,doorHeight:E,nearId:-1,type:z?`extendBeam`:`doorBeam`,isDoor:M,isPullOutDoor:N,pcbDoorIndexs:P};for(let t=0;t<l.length;t++){let r=e[l[t].linesIndex].doorAndBeamData[l[t].doorIndex];if(r.nearId!=-1)continue;let i=b(F.clone(),I.clone(),r.beamStart.clone(),r.beamEnd.clone());if(i&&i.isParallel){let e=F.distanceTo(I),t=r.beamStart.distanceTo(r.beamEnd),i=new c.Line3(F.clone(),I.clone());i.start.z=0,i.end.z=0;let a=S(r.beamStart,r.beamEnd,!0);t>e&&(i.start=new c.Vector3(r.beamStart.x,r.beamStart.y,0),i.end=new c.Vector3(r.beamEnd.x,r.beamEnd.y,0),a=S(F,I,!0));let o=new c.Vector3,l=i.closestPointToPoint(a,!0,o),u=l.distanceTo(a),d=new c.Vector3().subVectors(a,l).normalize(),p=d.clone().negate(),m=a.clone().addScaledVector(p,u+.1),h=f(i.start,i.end,a,m);if(u<.4&&h!=null){r.nearId=s,n.nearId=r.id;let i=.8;(r.isDoor||n.isDoor)&&(r.isDoor=!0,n.isDoor=!0);let a=P.some(e=>r.pcbDoorIndexs.includes(e));e<t?e>i&&!r.isPullOutDoor?(r.beamStart=F.clone().addScaledVector(p,u),r.beamEnd=I.clone().addScaledVector(p,u),r.doorStart=L.clone().addScaledVector(p,u),r.doorEnd=R.clone().addScaledVector(p,u),r.doorHeight=n.doorHeight,r.beamStart.z=r.beamEnd.z=r.doorStart.z+r.doorHeight):!N||N&&r.isPullOutDoor&&a?(n.beamStart=r.beamStart.clone().addScaledVector(d,u),n.beamEnd=r.beamEnd.clone().addScaledVector(d,u),n.doorStart=r.doorStart.clone().addScaledVector(d,u),n.doorEnd=r.doorEnd.clone().addScaledVector(d,u),n.doorHeight=r.doorHeight,n.beamStart.z=n.beamEnd.z=n.doorStart.z+n.doorHeight):(r.beamStart=F.clone().addScaledVector(p,u),r.beamEnd=I.clone().addScaledVector(p,u),r.doorStart=L.clone().addScaledVector(p,u),r.doorEnd=R.clone().addScaledVector(p,u),r.doorHeight=n.doorHeight,r.beamStart.z=r.beamEnd.z=r.doorStart.z+r.doorHeight):e>t&&(t>i&&!N?(n.beamStart=r.beamStart.clone().addScaledVector(p,u),n.beamEnd=r.beamEnd.clone().addScaledVector(p,u),n.doorStart=r.doorStart.clone().addScaledVector(p,u),n.doorEnd=r.doorEnd.clone().addScaledVector(p,u),n.doorHeight=r.doorHeight,n.beamStart.z=n.beamEnd.z=n.doorStart.z+n.doorHeight):!r.isPullOutDoor||N&&r.isPullOutDoor&&a?(r.beamStart=F.clone().addScaledVector(d,u),r.beamEnd=I.clone().addScaledVector(d,u),r.doorStart=L.clone().addScaledVector(d,u),r.doorEnd=R.clone().addScaledVector(d,u),r.doorHeight=n.doorHeight,r.beamStart.z=r.beamEnd.z=r.doorStart.z+r.doorHeight):(n.beamStart=r.beamStart.clone().addScaledVector(p,u),n.beamEnd=r.beamEnd.clone().addScaledVector(p,u),n.doorStart=r.doorStart.clone().addScaledVector(p,u),n.doorEnd=r.doorEnd.clone().addScaledVector(p,u),n.doorHeight=r.doorHeight,n.beamStart.z=n.beamEnd.z=n.doorStart.z+n.doorHeight));break}}}if(i.doorAndBeamData.push(n),l.push({linesIndex:r,doorIndex:i.doorAndBeamData.length-1}),s++,N&&T<2)break}}}let E=i.length*(i.rooftopPz-i.start.z);i.completePointAreaPercentage=T/E*100}if(u.length>1)for(let t=0;t<u.length;t++){let n=[];for(let r=0;r<e[u[t]].doorAndBeamData.length;r++){let i=e[u[t]].doorAndBeamData[r].doorStart.clone(),a=e[u[t]].doorAndBeamData[r].doorEnd.clone();i.z=a.z=0;let o=i.distanceTo(a),s=!1;for(let n=0;n<u.length;n++)if(t!=n)for(let t=0;t<e[u[n]].doorAndBeamData.length;t++){let r=e[u[n]].doorAndBeamData[t].doorStart.clone(),l=e[u[n]].doorAndBeamData[t].doorEnd.clone();r.z=l.z=0,r.distanceTo(l);let d=C(i,a,r,l);if(d.type==`collinear_overlap`&&d.maxPerpendicularDistance<.2){let e=new c.Line3(r,l),t=S(i,a),n=e.closestPointToPoint(t,!0,new c.Vector3),u=t.distanceTo(n);if(Math.abs(o-d.overlap)<.1&&u<.05){s=!0;break}}}s||n.push(e[u[t]].doorAndBeamData[r])}n.length!=e[u[t]].doorAndBeamData.length&&(e[u[t]].doorAndBeamData=n)}let p=.65,m=[],g=[],_=.2;for(let t=0;t<o.length;t++){if(o[t].isFind)continue;let r=new c.Vector3(o[t].doorStartPt.x,o[t].doorStartPt.y,0),i=new c.Vector3(o[t].doorEndPt.x,o[t].doorEndPt.y,0),a=new c.Vector3().subVectors(r,i).normalize();for(let l=0;l<e.length;l++){if(e[l].length<1||e[l].completePointAreaPercentage>60)continue;let u=new c.Vector3(e[l].start.x,e[l].start.y,0),d=new c.Vector3(e[l].end.x,e[l].end.y,0),h=new c.Vector3().subVectors(u,d).normalize(),v=.3,y=r.clone().addScaledVector(a,v),S=i.clone().addScaledVector(a.clone().negate(),v),w=u.clone().addScaledVector(h,v),T=d.clone().addScaledVector(h.clone().negate(),v);f(y,S,w,T);let E=b(y,S,w,T),D=r.distanceTo(u)>i.distanceTo(u)?i:r,O=D.distanceTo(u),k=D.distanceTo(d);if((O<_||k<_)&&E&&!E.isParallel){let a=r.distanceTo(u)<r.distanceTo(d)?u:d,f=!1;for(let h=0;h<e.length;h++){if(h==l||m.includes(h)||e[h].length<1)continue;let _=new c.Vector3(e[h].start.x,e[h].start.y,0),v=new c.Vector3(e[h].end.x,e[h].end.y,0),y=x(u,d,_,v);if(y&&y.angle>85){let y=new c.Line3(u,d),x=new c.Vector3,S=y.closestPointToPoint(_,!0,x),w=y.closestPointToPoint(v,!0,x),T=S.distanceTo(_),E=w.distanceTo(v),D,O;if(T>E&&E>p)D=v,O=w;else if(T<E&&T>p)D=_,O=S;else continue;if(T=a.distanceTo(_),E=a.distanceTo(v),T>2&&E>2)continue;let k=b(D,O,_,v);if(!k||!k.isParallel)continue;let A=C(r,i,D,O);if(A&&A.type==`collinear_overlap`&&A.maxPerpendicularDistance<.2)continue;D.z=O.z=e[h].start.z;let j=!1;for(let e=0;e<g.length;e++){let t=C(g[e].start,g[e].end,D,O);if(t&&t.type==`collinear_overlap`&&t.maxPerpendicularDistance<.2){j=!0;break}}if(j)continue;g.push({start:D,end:O}),m.push(h);let M=e[h].rooftopPz-e[h].start.z,N=e[l].rooftopPz-e[l].start.z,P=M<N?M:N;e[h].doorAndBeamData||(e[h].doorAndBeamData=[]),e[h].doorAndBeamData.push({id:s,beamStart:D,beamEnd:D,beamHeight:P,doorStart:D,doorEnd:O,doorHeight:P,nearId:-1,type:`onlyDoor`,isDoor:!0}),n[o[t].index].isFindBeam=!0,n[o[t].index].isFindOnlyDoor=!0,s++,f=!0;break}}if(f)break}}}let v=0;for(let t=0;t<e.length;t++){let n=e[t];n.mainBeamDataIds||=[],n.isMainBeam||=!1;let{minZ:r,maxZ:i}=O(n);if(r==null||i==null)continue;i-r;let o=n.rooftopPz-n.start.z;if(o<a-.5||o>a+2||Math.abs(r-n.start.z)>.3||n.length<1)continue;let s=new c.Vector3(n.start.x,n.start.y,0),l=new c.Vector3(n.end.x,n.end.y,0);for(let r=0;r<e.length;r++){if(r==t||e[r].length<1)continue;let a=e[r],o=new c.Vector3(a.start.x,a.start.y,0),u=new c.Vector3(a.end.x,a.end.y,0),d=b(s,l,o,u),f=S(s,l),p=new c.Line3(o,u).closestPointToPoint(f,!0,new c.Vector3),m=f.distanceTo(p),h=O(a);if(!(h.minZ===void 0||h.maxZ===void 0)&&!(h.maxZ-h.minZ>1)&&d&&d.isParallel&&m<.25&&Math.abs(i-h.minZ)<.1){let e=`mainBeamId${v}`;a.mainBeamDataIds||=[],a.isMainBeam||=!0,a.mainBeamDataIds.push(e),n.mainBeamDataIds.push(e),v++;break}}}return e},z=(e,t)=>{if(e.length<=0)return e;for(let t of e){let e=[];for(let n=0;n<t.checkResults.length;n++){let r=t.checkResults[n];if(r.isDoor=!1,!r.originalVertices||r.originalVertices.length==0)continue;let i=[];r.originalVertices.forEach(e=>{let t=[];e.forEach(e=>{t.push(new c.Vector3(e.x,e.y,e.z))}),i.push(t)});let a=i.length,o=Array.from({length:a},(e,t)=>t);function s(e){for(;o[e]!==e;)o[e]=o[o[e]],e=o[e];return e}for(let e=0;e<a;e++)for(let t=e+1;t<a;t++){let n=i[e],r=i[t];if(n.some(e=>r.some(t=>e.equals(t)))){let n=s(e),r=s(t);n!==r&&(o[n]=r)}}let l=new Map,u=-1,d=-1;for(let e=0;e<a;e++){let t=s(e);l.has(t)||l.set(t,{points:[],area:0,originalVertices:[],index:-1,gridWidth:-1,gridHeight:-1,centerPts:[]});let r=l.get(t),[a,o,f,p]=i[e],m=new c.Vector3().add(a).add(o).add(f).add(p).multiplyScalar(.25);if(r.points.push(...i[e]),r.area+=a.distanceTo(o)*a.distanceTo(p),r.originalVertices.push(i[e]),r.index=n,r.centerPts.push(m),u==-1||d==-1){let e=new c.Vector3().subVectors(a,o).normalize();d=a.distanceTo(p),u=a.distanceTo(o);let t=h(e,new c.Vector3(0,0,1));t&&t.parallel&&([d,u]=[u,d])}r.gridWidth=u,r.gridHeight=d}e.push(...l.values())}t.mergeCheckRegion=e}let n=new Map;for(let t of e){let e=t.rooftopPz-t.start.z,r=Math.floor(t.rooftopPz-t.start.z);n.has(r)?n.set(r,{num:n.get(r).num+1,totalHeight:n.get(r).totalHeight+e}):n.set(r,{num:1,totalHeight:e})}let r=[...n.entries()].reduce((e,t)=>t[1].num>e[1].num?t:e),i=r[0];r[1].totalHeight/r[1].num,new c.MeshBasicMaterial({color:`#cffd00`,transparent:!0,opacity:.4,side:c.DoubleSide}),new c.PointsMaterial({color:65535,size:.01});for(let n=0;n<e.length;n++){let r=e[n],{minZ:a,maxZ:o}=O(r);if(a==null||o==null)continue;r.doorAndBeamData=[],r.completePointAreaPercentage=-1,o-a;let s=r.rooftopPz-r.start.z;if(s<i-.5||s>i+2||Math.abs(a-r.start.z)>.3||r.length<1)continue;let l=new c.Vector3(r.start.x,r.start.y,0),u=new c.Vector3(r.end.x,r.end.y,0);for(let r=0;r<e.length;r++){if(r==n||e[r].length<1)continue;let i=e[r],a=new c.Vector3(i.start.x,i.start.y,0),s=new c.Vector3(i.end.x,i.end.y,0),d=b(l,u,a,s),f=S(l,u),m=new c.Line3(a,s).closestPointToPoint(f,!0,new c.Vector3),h=f.distanceTo(m),g=O(i);if(!(g.minZ===void 0||g.maxZ===void 0)&&!(g.maxZ-g.minZ>1)&&d&&d.isParallel&&h<.25&&Math.abs(o-g.minZ)<.1){let i=S(l,u),o=S(a,s);i.z=o.z=1,i.applyEuler(new c.Euler(-Math.PI/2,0,0)),o.applyEuler(new c.Euler(-Math.PI/2,0,0));let d=`w:${e[n].originaIndex},b:${e[r].originaIndex}`,f=`b:${e[r].originaIndex},w:${e[n].originaIndex}`;t.add(p(d,i,{textColor:`#fd0000`})),t.add(p(f,o,{textColor:`#fd0000`}))}}}return e};exports.classifySegments=C,exports.getAllGeometry=I,exports.getBeamLine=R,exports.getMainBeamLine=z,exports.getParallelism=b,exports.getPointCloudMinMax=O,exports.getPointCoverageOnQuad=E,exports.isParallel=h,exports.perpendicularInfo=x,exports.removeNoisePoints=A,exports.segmentsIntersect2D=f,exports.usegetBeamLine=L;

package/dist/index.js

@@ -538,7 +538,14 @@ var T = (t, n, r, i = .05, a = .02) => {
 		let { newOriginalPoints: S } = O(p, h, i.originalPoints);
 		i.originalPoints = S;
 	}
-}, A = (t, n) => new e.Vector3(t.x, t.y, 0).distanceTo(new e.Vector3(n.x, n.y, 0)), j = class {
+}, A = (t, n) => new e.Vector3(t.x, t.y, 0).distanceTo(new e.Vector3(n.x, n.y, 0)), j = (e) => {
+	let t = 0, n = [];
+	for (let r = 0; r < e.length; r++) for (let i = r + 1; i < e.length; i++) {
+		let a = e[r].distanceTo(e[i]);
+		a > t && (t = a, n = [e[r], e[i]]);
+	}
+	return n;
+}, M = class {
 	results;
 	clusterResults;
 	constructor() {
@@ -970,12 +977,12 @@ var T = (t, n, r, i = .05, a = .02) => {
 			return console.error("合并网格时出错:", e), null;
 		}
 	}
-}, M = (e) => {
+}, N = (e) => {
 	let t = {}, n = [], r = [], i = [];
 	return e[0].width, e[0].height, e.forEach((e) => {
 		r.push(...e.vertices), i.push(e.center), n.push(e.vertices);
 	}), t.originalVertices = n, t.allCenterPoints = i, t;
-}, N = (e) => {
+}, P = (e) => {
 	let t = new s().copy(e.start).add(e.end).multiplyScalar(.5), n = new s().subVectors(new s(t.x, t.y, t.z + 1), t).normalize(), r = new s(e.direction.x, e.direction.y, 0);
 	r.applyAxisAngle(n, Math.PI / 2);
 	let i = new s().copy(t).add(r.clone().multiplyScalar(e.length)), a = new s().copy(t).sub(r.clone().multiplyScalar(e.length));
@@ -985,7 +992,7 @@ var T = (t, n, r, i = .05, a = .02) => {
 		center: t,
 		direction: r
 	};
-}, P = (t, n, r = !0) => {
+}, F = (t, n, r = !0) => {
 	let i = [];
 	if (!t) return [];
 	for (let o = 0; o < t.length; o++) {
@@ -1006,7 +1013,7 @@ var T = (t, n, r, i = .05, a = .02) => {
 			wireframe: !1,
 			transparent: !0,
 			opacity: .5
-		}), _ = new e.Group(), v = new j(), { horizontalSubdivisions: y, verticalSubdivisions: b } = v.initLimits4(m, h), x = v.createSubdividedWallPlane(c, b, y, g, _);
+		}), _ = new e.Group(), v = new M(), { horizontalSubdivisions: y, verticalSubdivisions: b } = v.initLimits4(m, h), x = v.createSubdividedWallPlane(c, b, y, g, _);
 		v.analyzePointDistribution(x, p, .01);
 		let S = new e.MeshBasicMaterial({
 			color: "black",
@@ -1015,84 +1022,84 @@ var T = (t, n, r, i = .05, a = .02) => {
 			side: e.DoubleSide
 		}), C = v.highlightEmptySegments(x, S, _), w = [];
 		for (let e of C.clusters) {
-			let t = M(e);
+			let t = N(e);
 			t && w.push(t);
 		}
-		let T = N(c).direction;
+		let T = P(c).direction;
 		c.checkResults = w, c.verticalDirection = T, c.originaIndex = o, i.push(c);
 	}
 	return i;
-}, F = async (t, n, r, i) => {
+}, I = async (t, n, r, i) => {
 	let a = await fetch(t);
 	if (!a) return;
-	let o = await a.json();
-	for (let t = 0; t < o.length; t++) {
-		let n = o[t];
+	let o = await a.json(), s = [o[31], o[59]];
+	for (let t = 0; t < s.length; t++) {
+		let n = s[t];
 		if (!n) continue;
 		let r = new e.Vector3(n.start.x, n.start.y, n.start.z);
 		r.applyEuler(new e.Euler(-Math.PI / 2, 0, 0));
 		let a = new e.Vector3(n.end.x, n.end.y, n.end.z);
 		a.applyEuler(new e.Euler(-Math.PI / 2, 0, 0)), i.add(f(String(t), new e.Vector3((r.x + a.x) / 2, (r.y + a.y) / 2, (r.z + a.z) / 2)));
 	}
-	o.sort((e, t) => e.originalIndex - t.originalIndex), console.log("lins", o);
-	let s = await fetch(n);
-	if (!s) return;
-	let l = await s.json(), u = Array.isArray(l) ? l : Object.values(l), d = await fetch(r);
-	if (console.log("respcd", d), !d) return;
-	let p = await d.json();
-	k(o, i), console.time(), o = await I(o, u, p, i), console.timeEnd();
-	let m = new e.MeshBasicMaterial({
+	s.sort((e, t) => e.originalIndex - t.originalIndex), console.log("lins", s);
+	let l = await fetch(n);
+	if (!l) return;
+	let u = await l.json(), d = Array.isArray(u) ? u : Object.values(u), p = await fetch(r);
+	if (console.log("respcd", p), !p) return;
+	let m = await p.json();
+	k(s, i), console.time(), s = await L(s, d, m, i), console.timeEnd();
+	let h = new e.MeshBasicMaterial({
 		color: "#cffd00",
 		transparent: !0,
 		opacity: .4,
 		side: e.DoubleSide
-	}), h = new e.MeshBasicMaterial({
+	}), g = new e.MeshBasicMaterial({
 		color: "#ff0000",
 		transparent: !0,
 		opacity: .4,
 		side: e.DoubleSide
 	});
 	new e.LineBasicMaterial({ color: "#cffd00" });
-	let g = new e.PointsMaterial({
+	let v = new e.PointsMaterial({
 		color: 65535,
 		size: .01
-	}), v = [], y = [];
-	for (let t of o) {
+	}), y = [], b = [];
+	for (let t of s) {
 		if (!t.originalPoints || t.originalPoints.length === 0) continue;
 		let n = new Float32Array(t.originalPoints.length * 3);
 		t.originalPoints.forEach((e, t) => {
 			n[t * 3] = e.x, n[t * 3 + 1] = e.y, n[t * 3 + 2] = e.z;
 		});
 		let r = new e.BufferGeometry();
-		r.setAttribute("position", new e.BufferAttribute(n, 3)), v.push(r);
+		r.setAttribute("position", new e.BufferAttribute(n, 3)), y.push(r);
 		let i = [];
 		i.push(new e.Vector3(t.start.x, t.start.y, t.start.z)), i.push(new e.Vector3(t.end.x, t.end.y, t.end.z));
 		let a = new e.BufferGeometry();
-		a.setFromPoints(i), y.push(a);
+		a.setFromPoints(i), b.push(a);
 	}
-	if (v.length > 0) {
-		let t = c(v, !1);
-		t.rotateX(-Math.PI / 2), i.add(new e.Points(t, g)), v.forEach((e) => e.dispose());
+	if (y.length > 0) {
+		let t = c(y, !1);
+		t.rotateX(-Math.PI / 2), i.add(new e.Points(t, v)), y.forEach((e) => e.dispose());
 	}
-	let b = [], x = [], S = new e.LineBasicMaterial({ color: "#f30606" }), C = new e.MeshBasicMaterial({ color: "#071ac4" });
-	for (let t = 0; t < u.length; t++) {
+	let x = [], S = [], C = new e.LineBasicMaterial({ color: "#f30606" }), w = new e.MeshBasicMaterial({ color: "#071ac4" });
+	for (let t = 0; t < d.length; t++) {
 		let n = new e.SphereGeometry(.05);
-		n.translate(u[t].x, u[t].y, u[t].z);
-		let r = [], i = new e.Vector3(u[t].x, u[t].y, u[t].z), a = new e.Quaternion(u[t].qx, u[t].qy, u[t].qz, u[t].qw), o = new e.Vector3(0, 0, -1).applyQuaternion(a), s = i.clone().addScaledVector(o, .001);
+		n.translate(d[t].x, d[t].y, d[t].z);
+		let r = [], i = new e.Vector3(d[t].x, d[t].y, d[t].z), a = new e.Quaternion(d[t].qx, d[t].qy, d[t].qz, d[t].qw), o = new e.Vector3(0, 0, -1).applyQuaternion(a), s = i.clone().addScaledVector(o, .001);
 		r.push(i), r.push(s);
 		let c = new e.BufferGeometry();
-		c.setFromPoints(r), x.push(c), b.push(n);
-	}
-	if (b.length > 0) {
-		let t = c(b, !1);
-		t.rotateX(-Math.PI / 2), i.add(new e.Mesh(t, C)), b.forEach((e) => e.dispose());
+		c.setFromPoints(r), S.push(c), x.push(n);
 	}
 	if (x.length > 0) {
 		let t = c(x, !1);
-		t.rotateX(-Math.PI / 2), i.add(new e.Line(t, S)), x.forEach((e) => e.dispose());
+		t.rotateX(-Math.PI / 2), i.add(new e.Mesh(t, w)), x.forEach((e) => e.dispose());
+	}
+	if (S.length > 0) {
+		let t = c(S, !1);
+		t.rotateX(-Math.PI / 2), i.add(new e.Line(t, C)), S.forEach((e) => e.dispose());
 	}
-	let w = _(p);
-	w.rotateX(-Math.PI / 2), i.add(w), console.log("jsonpcdData", p), o.forEach((t) => {
+	let T = _(m);
+	T.rotateX(-Math.PI / 2), i.add(T), console.log("jsonpcdData", m), s.forEach((t) => {
 		if (t.doorAndBeamData && t.doorAndBeamData.length > 0) {
 			let n = new e.Group();
 			t.doorAndBeamData.forEach((t) => {
@@ -1118,7 +1125,7 @@ var T = (t, n, r, i = .05, a = .02) => {
 						0,
 						2,
 						3
-					]), r.computeVertexNormals(), r.rotateX(-Math.PI / 2), n.add(new e.Mesh(r, m));
+					]), r.computeVertexNormals(), r.rotateX(-Math.PI / 2), n.add(new e.Mesh(r, h));
 				}
 				let r = new e.BufferGeometry(), i = new Float32Array([
 					t.doorStart.x,
@@ -1141,13 +1148,13 @@ var T = (t, n, r, i = .05, a = .02) => {
 					0,
 					2,
 					3
-				]), r.computeVertexNormals(), r.rotateX(-Math.PI / 2), n.add(new e.Mesh(r, h));
+				]), r.computeVertexNormals(), r.rotateX(-Math.PI / 2), n.add(new e.Mesh(r, g));
 				let a = `id:${t.id},Nid:${t.nearId}`, o = new e.Vector3((t.doorStart.x + t.doorEnd.x) / 2, (t.doorStart.y + t.doorEnd.y) / 2, t.doorStart.z + t.doorHeight);
 				o.applyEuler(new e.Euler(-Math.PI / 2, 0, 0)), n.add(f(a, o));
 			}), i.add(n);
 		}
 	});
-}, I = async (t, n, r) => {
+}, L = async (t, n, r) => {
 	if (t.length <= 0 || n.length <= 0) return t;
 	r ||= [];
 	for (let n of t) {
@@ -1265,32 +1272,64 @@ var T = (t, n, r, i = .05, a = .02) => {
 			} else continue;
 		}
 	}
-	let c = 0, l = [], u = [];
+	let c = 0, l = [];
+	new e.PointsMaterial({
+		color: 65535,
+		size: .01
+	});
+	let u = [], f = [];
 	for (let i = 0; i < t.length; i++) {
 		let a = t[i];
 		if (!a.mergeCheckRegion || a.mergeCheckRegion.length == 0) continue;
 		a.doorAndBeamData = [], a.completePointAreaPercentage = -1;
-		let f = new e.Box3();
-		f.setFromPoints(a.originalPoints);
-		let p = f.max.z, h = f.min.z, g = new e.Vector3(a.start.x, a.start.y, a.start.z), _ = new e.Vector3(a.end.x, a.end.y, a.end.z), C = a.rooftopPz - a.start.z;
-		if (C < o - .5 || C > o + 2) continue;
-		let w = 0;
+		let p = new e.Box3();
+		p.setFromPoints(a.originalPoints);
+		let h = p.max.z, g = p.min.z, _ = new e.Vector3(a.start.x, a.start.y, a.start.z), C = new e.Vector3(a.end.x, a.end.y, a.end.z), w = a.rooftopPz - a.start.z;
+		if (w < o - .5 || w > o + 2) continue;
+		let T = 0;
 		for (let o of a.mergeCheckRegion) {
-			let { facePoints: f, boxArea: T, totalInPlaneNum: E } = v(o.points, o.gridWidth, o.gridHeight, o.centerPts);
-			if (!f || f.length <= 0) continue;
-			w += o.area;
-			let D = 0;
-			if (E != 0 && (D = o.gridHeight * o.gridWidth * E), (o.area - D) / T * 100 > 85) {
-				let v = f[0].distanceTo(f[1]), w = v / a.length * 100, T = f[0].distanceTo(f[3]), E = a.rooftopPz - f[2].z;
-				if (h > a.start.z - .15 && a.length < 1.5 || a.length, T < .5 || p < f[2].z || E > C / 3) continue;
-				let D = !1, k = !1, j = !1, M = [], N = f[2].clone(), P = f[3].clone(), F = f[0].clone(), I = f[1].clone(), L = !1;
-				if (Math.abs(f[0].z - a.start.z) < .25) {
-					if (v > .2 && a.length > .5) {
-						let n = A(g, F), r = A(g, I), o = A(_, F), s = A(_, I), c = .1;
+			let { facePoints: p, boxArea: E, totalInPlaneNum: D } = v(o.points, o.gridWidth, o.gridHeight, o.centerPts);
+			if (!p || p.length <= 0) continue;
+			p[0].distanceTo(p[1]);
+			let k = p[0].distanceTo(p[3]), M = o.gridHeight * (k / o.gridHeight / 2) + o.gridHeight / 2, N = new e.Vector3(p[0].x, p[0].y, p[0].z + M), P = new e.Vector3(p[1].x, p[1].y, p[1].z + M), F = new e.Line3(N, P), I = [], L = new e.Vector3();
+			o.points.forEach((e) => {
+				F.closestPointToPoint(e, !0, L).distanceTo(e) < o.gridHeight + .01 && I.push(e);
+			});
+			let R = [], z = [];
+			for (;;) {
+				let e = [];
+				for (let t = 0; t < I.length; t++) z.includes(t) || e.length == 0 && (e.push(I[t]), z.push(t));
+				for (let t = 0; t < I.length; t++) if (!z.includes(t)) {
+					for (let n = 0; n < e.length; n++) if (I[t].distanceTo(e[n]) < o.gridWidth + .01) {
+						e.push(I[t]), z.push(t), t = -1;
+						break;
+					}
+				}
+				if (e.length == 0) break;
+				R.push(e);
+			}
+			let B = [];
+			for (let t = 0; t < R.length; t++) {
+				let n = {
+					lineSt: new e.Vector3(),
+					lineEd: new e.Vector3(),
+					facePoints: []
+				}, r = j(R[t]);
+				Array.isArray(r) && r.length == 2 && (n.lineSt = r[0], n.lineEd = r[1], n.facePoints.push(new e.Vector3(r[0].x, r[0].y, p[0].z)), n.facePoints.push(new e.Vector3(r[1].x, r[1].y, p[0].z)), n.facePoints.push(new e.Vector3(r[1].x, r[1].y, p[2].z)), n.facePoints.push(new e.Vector3(r[0].x, r[0].y, p[2].z)), B.push(n));
+			}
+			T += o.area;
+			let V = 0;
+			if (D != 0 && (V = o.gridHeight * o.gridWidth * D), (o.area - V) / E * 100 > 85 && B.length == 1 || B.length > 1) for (let p = 0; p < B.length; p++) {
+				let v = B[p], T = v.facePoints[0].distanceTo(v.facePoints[1]), E = v.facePoints[0].distanceTo(v.facePoints[3]), D = T / a.length * 100, k = a.rooftopPz - v.facePoints[2].z;
+				if (g > a.start.z - .15 && a.length < 1.5 || a.length, E < .5 || h < v.facePoints[2].z || k > w / 3 * 2) continue;
+				let j = !1, M = !1, N = !1, P = [], F = v.facePoints[2].clone(), I = v.facePoints[3].clone(), L = v.facePoints[0].clone(), R = v.facePoints[1].clone(), z = !1;
+				if (Math.abs(v.facePoints[0].z - a.start.z) < .25) {
+					if (T > .2 && a.length > .5) {
+						let n = A(_, L), r = A(_, R), o = A(C, L), s = A(C, R), c = .1;
 						if (n < c || r < c) {
-							let { newOriginalPoints: r } = O(F, I, a.originalPoints), o = [], s = [];
+							let { newOriginalPoints: r } = O(L, R, a.originalPoints), o = [], s = [];
 							if (r.forEach((e) => {
-								e.z >= N.z ? o.push(e) : e.z <= F.z && s.push(e);
+								e.z >= F.z ? o.push(e) : e.z <= L.z && s.push(e);
 							}), o.length > 0) {
 								let r = new e.Vector3(a.start.x, a.start.y, 0), o = new e.Vector3(a.end.x, a.end.y, 0), s = new e.Vector3().subVectors(r, o).normalize(), l = new e.Vector3(), u = Infinity, f = -1;
 								for (let n = 0; n < t.length; n++) {
@@ -1301,10 +1340,10 @@ var T = (t, n, r, i = .05, a = .02) => {
 									if (p && p.angle > 85) {
 										let t = y(r, o, new e.Line3(a, c).closestPointToPoint(r, !0, l), r);
 										if (!t || t.direction != "same") continue;
-										let i = d(r.clone().addScaledVector(s, 10), o.clone().addScaledVector(s.clone().negate(), 10), a, c);
-										if (i == null) continue;
-										let p = r.distanceTo(i.point);
-										p < u && (u = p, f = n);
+										let i = r.clone().addScaledVector(s, 10), p = o.clone().addScaledVector(s.clone().negate(), 10), m = new e.Vector3().subVectors(c, a).normalize(), h = d(i, p, a.clone().addScaledVector(m.clone().negate(), .15), c.clone().addScaledVector(m, .15));
+										if (h == null) continue;
+										let g = r.distanceTo(h.point);
+										g < u && (u = g, f = n);
 									} else if (p && p.angle < 5) {
 										let e = A(r, a), t = A(r, c), n = A(o, a), i = A(o, c);
 										if (e < .1 || t < .1 || n < .1 || i < .1) {
@@ -1313,12 +1352,12 @@ var T = (t, n, r, i = .05, a = .02) => {
 										}
 									}
 								}
-								u != Infinity && u < 2 && t[f].length > .5 && (L = !0, n < c ? (P = P.addScaledVector(s, u), F = F.addScaledVector(s, u)) : (N = N.addScaledVector(s, u), I = I.addScaledVector(s, u)));
+								u != Infinity && u < 2 && t[f].length > .5 && (z = !0, n < c ? (I = I.addScaledVector(s, u), L = L.addScaledVector(s, u)) : (F = F.addScaledVector(s, u), R = R.addScaledVector(s, u)));
 							}
 						} else if (o < c || s < c) {
-							let { newOriginalPoints: n } = O(F, I, a.originalPoints), r = [], s = [];
+							let { newOriginalPoints: n } = O(L, R, a.originalPoints), r = [], s = [];
 							if (n.forEach((e) => {
-								e.z >= N.z ? r.push(e) : e.z <= F.z && s.push(e);
+								e.z >= F.z ? r.push(e) : e.z <= L.z && s.push(e);
 							}), r.length > 0) {
 								let n = new e.Vector3(a.start.x, a.start.y, 0), r = new e.Vector3(a.end.x, a.end.y, 0), s = new e.Vector3().subVectors(r, n).normalize(), l = new e.Vector3(), u = Infinity, f = -1;
 								for (let a = 0; a < t.length; a++) {
@@ -1329,10 +1368,10 @@ var T = (t, n, r, i = .05, a = .02) => {
 									if (p && p.angle > 85) {
 										let t = y(r, n, new e.Line3(o, c).closestPointToPoint(r, !0, l), r);
 										if (!t || t.direction != "same") continue;
-										let i = d(n.clone().addScaledVector(s, 10), r.clone().addScaledVector(s.clone().negate(), 10), o, c);
-										if (i == null) continue;
-										let p = r.distanceTo(i.point);
-										p < u && (u = p, f = a);
+										let i = n.clone().addScaledVector(s, 10), p = r.clone().addScaledVector(s.clone().negate(), 10), m = new e.Vector3().subVectors(c, o).normalize(), h = d(i, p, o.clone().addScaledVector(m.clone().negate(), .15), c.clone().addScaledVector(m, .15));
+										if (h == null) continue;
+										let g = r.distanceTo(h.point);
+										g < u && (u = g, f = a);
 									}
 									if (p && p.angle < 5) {
 										let e = A(n, o), t = A(n, c), i = A(r, o), a = A(r, c);
@@ -1342,124 +1381,132 @@ var T = (t, n, r, i = .05, a = .02) => {
 										}
 									}
 								}
-								u != Infinity && u < 2 && t[f].length > .5 && (L = !0, o < c ? (P = P.addScaledVector(s, u), F = F.addScaledVector(s, u)) : (N = N.addScaledVector(s, u), I = I.addScaledVector(s, u)));
+								u != Infinity && u < 2 && t[f].length > .5 && (z = !0, o < c ? (I = I.addScaledVector(s, u), L = L.addScaledVector(s, u)) : (F = F.addScaledVector(s, u), R = R.addScaledVector(s, u)));
 							}
 						}
-						v = P.distanceTo(N);
+						T = I.distanceTo(F);
 					}
-					if (w < 85 && v > .35 || v > .7 && T > C - C / 3) {
+					if (D < 85 && T > .35 || T > .7 && E > w - w / 3) {
 						let t = !1;
 						for (let r = 0; r < n.length; r++) {
 							let i = r + 1;
 							if (i >= n.length) continue;
 							let a = new e.Vector3(n[r].x, n[r].y, n[r].z), o = new e.Vector3(n[i].x, n[i].y, n[i].z);
-							if (d(new e.Vector3(a.x, a.y, 0), new e.Vector3(o.x, o.y, 0), new e.Vector3(f[0].x, f[0].y, 0), new e.Vector3(f[1].x, f[1].y, 0), .1) != null) {
+							if (d(new e.Vector3(a.x, a.y, 0), new e.Vector3(o.x, o.y, 0), new e.Vector3(v.facePoints[0].x, v.facePoints[0].y, 0), new e.Vector3(v.facePoints[1].x, v.facePoints[1].y, 0), .1) != null) {
 								t = !0;
 								break;
 							}
 						}
 						for (let t = 0; t < s.length; t++) {
-							let n = new e.Vector3(s[t].doorStartPt.x, s[t].doorStartPt.y, 0), i = new e.Vector3(s[t].doorEndPt.x, s[t].doorEndPt.y, 0), o = new e.Vector3(f[0].x, f[0].y, 0), c = new e.Vector3(f[1].x, f[1].y, 0), l = new e.Vector3().subVectors(n, i).normalize(), u = new e.Vector3().subVectors(o, c).normalize(), p = .2, m = d(n.clone().addScaledVector(l, p), i.clone().addScaledVector(l.clone().negate(), p), o.clone().addScaledVector(u, p), c.clone().addScaledVector(u.clone().negate(), p)), h = f[2].z < s[t].maxZ || Math.abs(f[2].z - s[t].maxZ) < .3, g = y(n.clone(), i.clone(), o.clone(), c.clone());
-							if (m && h) r[s[t].index].isFindBeam = !0, k = !0, s[t].isFind = !0;
-							else if (h && g && g.isParallel) {
+							let n = new e.Vector3(s[t].doorStartPt.x, s[t].doorStartPt.y, 0), i = new e.Vector3(s[t].doorEndPt.x, s[t].doorEndPt.y, 0), o = new e.Vector3(v.facePoints[0].x, v.facePoints[0].y, 0), c = new e.Vector3(v.facePoints[1].x, v.facePoints[1].y, 0), l = new e.Vector3().subVectors(n, i).normalize(), u = new e.Vector3().subVectors(o, c).normalize(), f = .2, p = d(n.clone().addScaledVector(l, f), i.clone().addScaledVector(l.clone().negate(), f), o.clone().addScaledVector(u, f), c.clone().addScaledVector(u.clone().negate(), f)), m = v.facePoints[2].z < s[t].maxZ || Math.abs(v.facePoints[2].z - s[t].maxZ) < .3, h = y(n.clone(), i.clone(), o.clone(), c.clone());
+							if (p && m) r[s[t].index].isFindBeam = !0, M = !0, s[t].isFind = !0;
+							else if (m && h && h.isParallel) {
 								let o = x(n, i), c = new e.Line3(new e.Vector3(a.start.x, a.start.y, 0), new e.Vector3(a.end.x, a.end.y, 0)), l = new e.Vector3(), u = c.closestPointToPoint(o, !0, l);
 								o.distanceTo(u) < .3 && (r[s[t].index].isFindBeam = !0);
 							}
 						}
-						if (!k) for (let t = 0; t < s.length; t++) {
-							let n = new e.Vector3(s[t].doorStartPt.x, s[t].doorStartPt.y, 0), i = new e.Vector3(s[t].doorEndPt.x, s[t].doorEndPt.y, 0), o = new e.Vector3(f[0].x, f[0].y, 0), c = new e.Vector3(f[1].x, f[1].y, 0), l = f[2].z < s[t].maxZ || Math.abs(f[2].z - s[t].maxZ) < .3, u = Math.abs(f[2].z - s[t].maxZ) > .2;
+						for (let t = 0; t < s.length; t++) {
+							let n = new e.Vector3(s[t].doorStartPt.x, s[t].doorStartPt.y, 0), i = new e.Vector3(s[t].doorEndPt.x, s[t].doorEndPt.y, 0), o = new e.Vector3(v.facePoints[0].x, v.facePoints[0].y, 0), c = new e.Vector3(v.facePoints[1].x, v.facePoints[1].y, 0), l = v.facePoints[2].z < s[t].maxZ || Math.abs(v.facePoints[2].z - s[t].maxZ) < .3, u = Math.abs(v.facePoints[2].z - s[t].maxZ) > .2;
 							if (l && !u) {
 								let l = S(n.clone(), i.clone(), o.clone(), c.clone(), { distanceThreshold: 2 });
 								if (l && l.type == "collinear_gap") {
-									let u = x(n, i), d = new e.Line3(new e.Vector3(a.start.x, a.start.y, 0), new e.Vector3(a.end.x, a.end.y, 0)), p = new e.Vector3(), m = d.closestPointToPoint(u, !0, p);
-									if (u.distanceTo(m) < .15) {
+									let u = x(n, i), d = new e.Line3(new e.Vector3(a.start.x, a.start.y, 0), new e.Vector3(a.end.x, a.end.y, 0)), f = new e.Vector3(), p = d.closestPointToPoint(u, !0, f);
+									if (u.distanceTo(p) < .15) {
 										let a = o.distanceTo(n), u = c.distanceTo(n);
 										if (a > .25 && u > .25) {
 											let e = o.distanceTo(i), t = c.distanceTo(i);
 											if (e > .2 && t > .2) continue;
 										}
-										k = !0, j = !0, M.push(s[t].index), s[t].isFind = !0, r[s[t].index].isFindBeam = !0;
-										let d = new e.Vector3().subVectors(o, c).normalize(), p = n.distanceTo(i) - l.overlap;
-										u < a && d.negate(), a > u ? (N = f[2].clone().addScaledVector(d, p), P = f[3].clone(), F = N.clone(), I = P.clone(), F.z = I.z = f[0].z) : (N = f[2].clone(), P = f[3].clone().addScaledVector(d, p), F = N.clone(), I = P.clone(), F.z = I.z = f[0].z);
+										M = !0, N = !0, P.push(s[t].index), s[t].isFind = !0, r[s[t].index].isFindBeam = !0;
+										let d = new e.Vector3().subVectors(o, c).normalize(), f = n.distanceTo(i) - l.overlap;
+										u < a && d.negate(), a > u ? (F = v.facePoints[2].clone().addScaledVector(d, f), I = v.facePoints[3].clone(), L = F.clone(), R = I.clone(), L.z = R.z = v.facePoints[0].z) : (F = v.facePoints[2].clone(), I = v.facePoints[3].clone().addScaledVector(d, f), L = F.clone(), R = I.clone(), L.z = R.z = v.facePoints[0].z);
 									}
-								} else if (l && l.type == "collinear_overlap" && l.maxPerpendicularDistance < .2) {
+								} else if (l && l.type == "collinear_overlap" && l.maxPerpendicularDistance < .35) {
 									let e = o.distanceTo(c), a = i.distanceTo(n);
-									(Math.abs(e - l.overlap) < .25 || Math.abs(a - l.overlap) < .25) && (k = !0, j = !0, M.push(s[t].index), s[t].isFind = !0, r[s[t].index].isFindBeam = !0);
+									(Math.abs(e - l.overlap) < .25 || Math.abs(a - l.overlap) < .25) && (M = !0, N = !0, P.push(s[t].index), s[t].isFind = !0, r[s[t].index].isFindBeam = !0);
 								}
 							}
 						}
-						(k || t || T > C / 2 && v > .8 && Math.abs(h - a.start.z) < .4) && (D = !0);
+						(M || t || E > w / 2 && T > .8 && Math.abs(g - a.start.z) < .4) && (j = !0);
 					}
 				}
-				if (D) {
-					if (!j) {
-						new e.PointsMaterial({
-							color: 65535,
-							size: .01
-						});
-						let t = new e.Vector3().subVectors(F, I), n = new e.Vector3().subVectors(I, F), r = Infinity, i = Infinity, o = F.clone(), s = I.clone(), c = [];
+				if (j) {
+					let n = !1;
+					for (let e = 0; e < f.length; e++) {
+						let t = S(f[e].start, f[e].end, L, R);
+						if (t && t.type == "collinear_overlap" && t.maxPerpendicularDistance < .01) {
+							n = !0;
+							break;
+						}
+					}
+					if (n) continue;
+					if (f.push({
+						start: L,
+						end: R
+					}), !N) {
+						let t = new e.Vector3().subVectors(L, R), n = new e.Vector3().subVectors(R, L), r = Infinity, i = Infinity, o = L.clone(), s = R.clone(), c = [];
 						if (a.originalPoints.forEach((e) => {
-							e.z - F.z < T / 2 + .1 && e.z - F.z > T / 2 - .1 && c.push(e);
+							e.z - L.z < E / 2 + .1 && e.z - L.z > E / 2 - .1 && c.push(e);
 						}), c.length > 0) {
-							let { newProjectPoints: a } = O(F, I, c);
-							for (let e of a) e.z = F.z;
+							let { newProjectPoints: a } = O(L, R, c);
+							for (let e of a) e.z = L.z;
 							a.forEach((a) => {
-								let c = new e.Vector3(a.x, a.y, F.z), l = new e.Vector3().subVectors(c, I), u = m(t, new e.Vector3().subVectors(c, F)), d = m(n, l);
+								let c = new e.Vector3(a.x, a.y, L.z), l = new e.Vector3().subVectors(c, R), u = m(t, new e.Vector3().subVectors(c, L)), d = m(n, l);
 								if (u && u.parallel && u.sameDirection) {
-									let e = c.distanceTo(F);
+									let e = c.distanceTo(L);
 									e < r && (r = e, o = c);
 								} else if (d && d.parallel && d.sameDirection) {
-									let e = c.distanceTo(I);
+									let e = c.distanceTo(R);
 									e < i && (i = e, s = c);
 								}
 							});
 						}
-						o.equals(F) || (F = o, P = new e.Vector3(F.x, F.y, P.z)), s.equals(I) || (I = s, N = new e.Vector3(I.x, I.y, N.z));
+						o.equals(L) || (L = o, I = new e.Vector3(L.x, L.y, I.z)), s.equals(R) || (R = s, F = new e.Vector3(R.x, R.y, F.z));
 					}
 					u.includes(i) || u.push(i), a.checkResults[o.index].isDoor = !0;
-					let n = {
+					let r = {
 						id: c,
-						beamStart: N,
-						beamEnd: P,
-						beamHeight: a.rooftopPz - f[2].z,
-						doorStart: F,
-						doorEnd: I,
-						doorHeight: T,
+						beamStart: F,
+						beamEnd: I,
+						beamHeight: a.rooftopPz - v.facePoints[2].z,
+						doorStart: L,
+						doorEnd: R,
+						doorHeight: E,
 						nearId: -1,
-						type: L ? "extendBeam" : "doorBeam",
-						isDoor: k,
-						isPullOutDoor: j,
-						pcbDoorIndexs: M
+						type: z ? "extendBeam" : "doorBeam",
+						isDoor: M,
+						isPullOutDoor: N,
+						pcbDoorIndexs: P
 					};
-					for (let r = 0; r < l.length; r++) {
-						let i = t[l[r].linesIndex].doorAndBeamData[l[r].doorIndex];
+					for (let n = 0; n < l.length; n++) {
+						let i = t[l[n].linesIndex].doorAndBeamData[l[n].doorIndex];
 						if (i.nearId != -1) continue;
-						let a = y(N.clone(), P.clone(), i.beamStart.clone(), i.beamEnd.clone());
+						let a = y(F.clone(), I.clone(), i.beamStart.clone(), i.beamEnd.clone());
 						if (a && a.isParallel) {
-							let t = N.distanceTo(P), r = i.beamStart.distanceTo(i.beamEnd), a = new e.Line3(N.clone(), P.clone());
+							let t = F.distanceTo(I), n = i.beamStart.distanceTo(i.beamEnd), a = new e.Line3(F.clone(), I.clone());
 							a.start.z = 0, a.end.z = 0;
 							let o = x(i.beamStart, i.beamEnd, !0);
-							r > t && (a.start = new e.Vector3(i.beamStart.x, i.beamStart.y, 0), a.end = new e.Vector3(i.beamEnd.x, i.beamEnd.y, 0), o = x(N, P, !0));
+							n > t && (a.start = new e.Vector3(i.beamStart.x, i.beamStart.y, 0), a.end = new e.Vector3(i.beamEnd.x, i.beamEnd.y, 0), o = x(F, I, !0));
 							let s = new e.Vector3(), l = a.closestPointToPoint(o, !0, s), u = l.distanceTo(o), f = new e.Vector3().subVectors(o, l).normalize(), p = f.clone().negate(), m = o.clone().addScaledVector(p, u + .1), h = d(a.start, a.end, o, m);
 							if (u < .4 && h != null) {
-								i.nearId = c, n.nearId = i.id;
+								i.nearId = c, r.nearId = i.id;
 								let e = .8;
-								(i.isDoor || n.isDoor) && (i.isDoor = !0, n.isDoor = !0);
-								let a = M.some((e) => i.pcbDoorIndexs.includes(e));
-								t < r ? t > e && !i.isPullOutDoor ? (i.beamStart = N.clone().addScaledVector(p, u), i.beamEnd = P.clone().addScaledVector(p, u), i.doorStart = F.clone().addScaledVector(p, u), i.doorEnd = I.clone().addScaledVector(p, u), i.doorHeight = n.doorHeight, i.beamStart.z = i.beamEnd.z = i.doorStart.z + i.doorHeight) : !j || j && i.isPullOutDoor && a ? (n.beamStart = i.beamStart.clone().addScaledVector(f, u), n.beamEnd = i.beamEnd.clone().addScaledVector(f, u), n.doorStart = i.doorStart.clone().addScaledVector(f, u), n.doorEnd = i.doorEnd.clone().addScaledVector(f, u), n.doorHeight = i.doorHeight, n.beamStart.z = n.beamEnd.z = n.doorStart.z + n.doorHeight) : (i.beamStart = N.clone().addScaledVector(p, u), i.beamEnd = P.clone().addScaledVector(p, u), i.doorStart = F.clone().addScaledVector(p, u), i.doorEnd = I.clone().addScaledVector(p, u), i.doorHeight = n.doorHeight, i.beamStart.z = i.beamEnd.z = i.doorStart.z + i.doorHeight) : t > r && (r > e && !j ? (n.beamStart = i.beamStart.clone().addScaledVector(p, u), n.beamEnd = i.beamEnd.clone().addScaledVector(p, u), n.doorStart = i.doorStart.clone().addScaledVector(p, u), n.doorEnd = i.doorEnd.clone().addScaledVector(p, u), n.doorHeight = i.doorHeight, n.beamStart.z = n.beamEnd.z = n.doorStart.z + n.doorHeight) : !i.isPullOutDoor || j && i.isPullOutDoor && a ? (i.beamStart = N.clone().addScaledVector(f, u), i.beamEnd = P.clone().addScaledVector(f, u), i.doorStart = F.clone().addScaledVector(f, u), i.doorEnd = I.clone().addScaledVector(f, u), i.doorHeight = n.doorHeight, i.beamStart.z = i.beamEnd.z = i.doorStart.z + i.doorHeight) : (n.beamStart = i.beamStart.clone().addScaledVector(p, u), n.beamEnd = i.beamEnd.clone().addScaledVector(p, u), n.doorStart = i.doorStart.clone().addScaledVector(p, u), n.doorEnd = i.doorEnd.clone().addScaledVector(p, u), n.doorHeight = i.doorHeight, n.beamStart.z = n.beamEnd.z = n.doorStart.z + n.doorHeight));
+								(i.isDoor || r.isDoor) && (i.isDoor = !0, r.isDoor = !0);
+								let a = P.some((e) => i.pcbDoorIndexs.includes(e));
+								t < n ? t > e && !i.isPullOutDoor ? (i.beamStart = F.clone().addScaledVector(p, u), i.beamEnd = I.clone().addScaledVector(p, u), i.doorStart = L.clone().addScaledVector(p, u), i.doorEnd = R.clone().addScaledVector(p, u), i.doorHeight = r.doorHeight, i.beamStart.z = i.beamEnd.z = i.doorStart.z + i.doorHeight) : !N || N && i.isPullOutDoor && a ? (r.beamStart = i.beamStart.clone().addScaledVector(f, u), r.beamEnd = i.beamEnd.clone().addScaledVector(f, u), r.doorStart = i.doorStart.clone().addScaledVector(f, u), r.doorEnd = i.doorEnd.clone().addScaledVector(f, u), r.doorHeight = i.doorHeight, r.beamStart.z = r.beamEnd.z = r.doorStart.z + r.doorHeight) : (i.beamStart = F.clone().addScaledVector(p, u), i.beamEnd = I.clone().addScaledVector(p, u), i.doorStart = L.clone().addScaledVector(p, u), i.doorEnd = R.clone().addScaledVector(p, u), i.doorHeight = r.doorHeight, i.beamStart.z = i.beamEnd.z = i.doorStart.z + i.doorHeight) : t > n && (n > e && !N ? (r.beamStart = i.beamStart.clone().addScaledVector(p, u), r.beamEnd = i.beamEnd.clone().addScaledVector(p, u), r.doorStart = i.doorStart.clone().addScaledVector(p, u), r.doorEnd = i.doorEnd.clone().addScaledVector(p, u), r.doorHeight = i.doorHeight, r.beamStart.z = r.beamEnd.z = r.doorStart.z + r.doorHeight) : !i.isPullOutDoor || N && i.isPullOutDoor && a ? (i.beamStart = F.clone().addScaledVector(f, u), i.beamEnd = I.clone().addScaledVector(f, u), i.doorStart = L.clone().addScaledVector(f, u), i.doorEnd = R.clone().addScaledVector(f, u), i.doorHeight = r.doorHeight, i.beamStart.z = i.beamEnd.z = i.doorStart.z + i.doorHeight) : (r.beamStart = i.beamStart.clone().addScaledVector(p, u), r.beamEnd = i.beamEnd.clone().addScaledVector(p, u), r.doorStart = i.doorStart.clone().addScaledVector(p, u), r.doorEnd = i.doorEnd.clone().addScaledVector(p, u), r.doorHeight = i.doorHeight, r.beamStart.z = r.beamEnd.z = r.doorStart.z + r.doorHeight));
 								break;
 							}
 						}
 					}
-					a.doorAndBeamData.push(n), l.push({
+					if (a.doorAndBeamData.push(r), l.push({
 						linesIndex: i,
 						doorIndex: a.doorAndBeamData.length - 1
-					}), c++;
+					}), c++, N && T < 2) break;
 				}
 			}
 		}
-		let T = a.length * (a.rooftopPz - a.start.z);
-		a.completePointAreaPercentage = w / T * 100;
+		let E = a.length * (a.rooftopPz - a.start.z);
+		a.completePointAreaPercentage = T / E * 100;
 	}
 	if (u.length > 1) for (let n = 0; n < u.length; n++) {
 		let r = [];
@@ -1483,47 +1530,50 @@ var T = (t, n, r, i = .05, a = .02) => {
 		}
 		r.length != t[u[n]].doorAndBeamData.length && (t[u[n]].doorAndBeamData = r);
 	}
-	let f = .65, p = [], h = [], g = .2;
+	let p = .65, h = [], g = [], _ = .2;
 	for (let n = 0; n < s.length; n++) {
 		if (s[n].isFind) continue;
 		let i = new e.Vector3(s[n].doorStartPt.x, s[n].doorStartPt.y, 0), a = new e.Vector3(s[n].doorEndPt.x, s[n].doorEndPt.y, 0), o = new e.Vector3().subVectors(i, a).normalize();
 		for (let l = 0; l < t.length; l++) {
 			if (t[l].length < 1 || t[l].completePointAreaPercentage > 60) continue;
-			let u = new e.Vector3(t[l].start.x, t[l].start.y, 0), m = new e.Vector3(t[l].end.x, t[l].end.y, 0), _ = new e.Vector3().subVectors(u, m).normalize(), v = .3, x = i.clone().addScaledVector(o, v), C = a.clone().addScaledVector(o.clone().negate(), v), w = u.clone().addScaledVector(_, v), T = m.clone().addScaledVector(_.clone().negate(), v);
+			let u = new e.Vector3(t[l].start.x, t[l].start.y, 0), f = new e.Vector3(t[l].end.x, t[l].end.y, 0), m = new e.Vector3().subVectors(u, f).normalize(), v = .3, x = i.clone().addScaledVector(o, v), C = a.clone().addScaledVector(o.clone().negate(), v), w = u.clone().addScaledVector(m, v), T = f.clone().addScaledVector(m.clone().negate(), v);
 			d(x, C, w, T);
-			let E = y(x, C, w, T), D = i.distanceTo(u) > a.distanceTo(u) ? a : i, O = D.distanceTo(u), k = D.distanceTo(m);
-			if ((O < g || k < g) && E && !E.isParallel) {
-				let o = i.distanceTo(u) < i.distanceTo(m) ? u : m, d = !1;
-				for (let g = 0; g < t.length; g++) {
-					if (g == l || p.includes(g) || t[g].length < 1) continue;
-					let _ = new e.Vector3(t[g].start.x, t[g].start.y, 0), v = new e.Vector3(t[g].end.x, t[g].end.y, 0), x = b(u, m, _, v);
+			let E = y(x, C, w, T), D = i.distanceTo(u) > a.distanceTo(u) ? a : i, O = D.distanceTo(u), k = D.distanceTo(f);
+			if ((O < _ || k < _) && E && !E.isParallel) {
+				let o = i.distanceTo(u) < i.distanceTo(f) ? u : f, d = !1;
+				for (let m = 0; m < t.length; m++) {
+					if (m == l || h.includes(m) || t[m].length < 1) continue;
+					let _ = new e.Vector3(t[m].start.x, t[m].start.y, 0), v = new e.Vector3(t[m].end.x, t[m].end.y, 0), x = b(u, f, _, v);
 					if (x && x.angle > 85) {
-						let b = new e.Line3(u, m), x = new e.Vector3(), C = b.closestPointToPoint(_, !0, x), w = b.closestPointToPoint(v, !0, x), T = C.distanceTo(_), E = w.distanceTo(v), D, O;
-						if (T > E && E > f) D = v, O = w;
-						else if (T < E && T > f) D = _, O = C;
+						let b = new e.Line3(u, f), x = new e.Vector3(), C = b.closestPointToPoint(_, !0, x), w = b.closestPointToPoint(v, !0, x), T = C.distanceTo(_), E = w.distanceTo(v), D, O;
+						if (T > E && E > p) D = v, O = w;
+						else if (T < E && T > p) D = _, O = C;
 						else continue;
 						if (T = o.distanceTo(_), E = o.distanceTo(v), T > 2 && E > 2) continue;
 						let k = y(D, O, _, v);
 						if (!k || !k.isParallel) continue;
 						let A = S(i, a, D, O);
 						if (A && A.type == "collinear_overlap" && A.maxPerpendicularDistance < .2) continue;
-						D.z = O.z = t[g].start.z;
+						D.z = O.z = t[m].start.z;
 						let j = !1;
-						for (let e = 0; e < h.length; e++) {
-							let t = S(h[e].start, h[e].end, D, O);
+						for (let e = 0; e < g.length; e++) {
+							let t = S(g[e].start, g[e].end, D, O);
 							if (t && t.type == "collinear_overlap" && t.maxPerpendicularDistance < .2) {
 								j = !0;
 								break;
 							}
 						}
 						if (j) continue;
-						h.push({
+						g.push({
 							start: D,
 							end: O
-						}), p.push(g);
-						let M = t[g].rooftopPz - t[g].start.z, N = t[l].rooftopPz - t[l].start.z, P = M < N ? M : N;
-						t[g].doorAndBeamData || (t[g].doorAndBeamData = []), t[g].doorAndBeamData.push({
+						}), h.push(m);
+						let M = t[m].rooftopPz - t[m].start.z, N = t[l].rooftopPz - t[l].start.z, P = M < N ? M : N;
+						t[m].doorAndBeamData || (t[m].doorAndBeamData = []), t[m].doorAndBeamData.push({
 							id: c,
+							beamStart: D,
+							beamEnd: D,
+							beamHeight: P,
 							doorStart: D,
 							doorEnd: O,
 							doorHeight: P,
@@ -1538,7 +1588,7 @@ var T = (t, n, r, i = .05, a = .02) => {
 			}
 		}
 	}
-	let _ = 0;
+	let C = 0;
 	for (let n = 0; n < t.length; n++) {
 		let r = t[n];
 		r.mainBeamDataIds ||= [], r.isMainBeam ||= !1;
@@ -1552,14 +1602,14 @@ var T = (t, n, r, i = .05, a = .02) => {
 			if (i == n || t[i].length < 1) continue;
 			let o = t[i], s = new e.Vector3(o.start.x, o.start.y, 0), u = new e.Vector3(o.end.x, o.end.y, 0), d = y(c, l, s, u), f = x(c, l), p = new e.Line3(s, u).closestPointToPoint(f, !0, new e.Vector3()), m = f.distanceTo(p), h = D(o);
 			if (!(h.minZ === void 0 || h.maxZ === void 0) && !(h.maxZ - h.minZ > 1) && d && d.isParallel && m < .25 && Math.abs(a - h.minZ) < .1) {
-				let e = `mainBeamId${_}`;
-				o.mainBeamDataIds ||= [], o.isMainBeam ||= !0, o.mainBeamDataIds.push(e), r.mainBeamDataIds.push(e), _++;
+				let e = `mainBeamId${C}`;
+				o.mainBeamDataIds ||= [], o.isMainBeam ||= !0, o.mainBeamDataIds.push(e), r.mainBeamDataIds.push(e), C++;
 				break;
 			}
 		}
 	}
 	return t;
-}, L = (t, n) => {
+}, R = (t, n) => {
 	if (t.length <= 0) return t;
 	for (let n of t) {
 		let t = [];
@@ -1652,4 +1702,4 @@ var T = (t, n, r, i = .05, a = .02) => {
 	return t;
 };
 //#endregion
-export { S as classifySegments, P as getAllGeometry, I as getBeamLine, L as getMainBeamLine, y as getParallelism, D as getPointCloudMinMax, T as getPointCoverageOnQuad, m as isParallel, b as perpendicularInfo, k as removeNoisePoints, d as segmentsIntersect2D, F as usegetBeamLine };
+export { S as classifySegments, F as getAllGeometry, L as getBeamLine, R as getMainBeamLine, y as getParallelism, D as getPointCloudMinMax, T as getPointCoverageOnQuad, m as isParallel, b as perpendicularInfo, k as removeNoisePoints, d as segmentsIntersect2D, I as usegetBeamLine };

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "rm-graphical-computing",
-  "version": "1.0.6",
+  "version": "1.0.8",
   "description": "Three.js wall segment processing utilities.",
   "main": "dist/index.cjs",
   "module": "dist/index.js",