rm-graphical-computing 1.0.10 → 1.0.12

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
package/dist/index.cjs CHANGED
@@ -1 +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=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();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 R(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)}})},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/3*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;if(Math.abs(y.facePoints[2].z-o[e].maxZ),l){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;d.push({start:L,end:R});{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<.7||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,h=!1;for(let v=0;v<e.length;v++){if(v==l||m.includes(v)||e[v].length<.7)continue;let y=new c.Vector3(e[v].start.x,e[v].start.y,0),S=new c.Vector3(e[v].end.x,e[v].end.y,0),w=new c.Vector3().subVectors(y,S).normalize(),T=w.clone().negate(),E=y.clone().addScaledVector(w,5),D=S.clone().addScaledVector(T,5),A=x(u,d,y,S),j=b(u,d,y,S),M=f(u,d,E,D),N=!0;if(M!=null&&(N=M.point.distanceTo(u)<_||M.point.distanceTo(d)<_),(A&&A.angle>85||j&&j.isParallel)&&N){let f=new c.Line3(u,d),_=new c.Vector3,x=f.closestPointToPoint(y,!0,_),w=f.closestPointToPoint(S,!0,_),T=x.distanceTo(y),E=w.distanceTo(S),D,A;if(T>E&&E>p)D=S,A=w;else if(T<E&&T>p)D=y,A=x;else continue;if(T=a.distanceTo(y),E=a.distanceTo(S),T>2&&E>2)continue;let j=b(D,A,y,S);if(!j||!j.isParallel)continue;let M=C(r,i,D,A);if(M&&M.type==`collinear_overlap`&&M.maxPerpendicularDistance<.2)continue;D.z=A.z=e[v].start.z;let N=!1;for(let e=0;e<g.length;e++){let t=C(g[e].start,g[e].end,D,A);if(t&&t.type==`collinear_overlap`&&t.maxPerpendicularDistance<.2){N=!0;break}}if(N)continue;g.push({start:D,end:A}),m.push(v),console.log(`doorDatas[n]`,o[t],l,v,O,k);let P=e[v].rooftopPz-e[v].start.z,F=e[l].rooftopPz-e[l].start.z,I=P<F?P:F;e[v].doorAndBeamData||(e[v].doorAndBeamData=[]),e[v].doorAndBeamData.push({id:s,beamStart:D,beamEnd:D,beamHeight:I,doorStart:D,doorEnd:A,doorHeight:I,nearId:-1,type:`onlyDoor`,isDoor:!0}),n[o[t].index].isFindBeam=!0,n[o[t].index].isFindOnlyDoor=!0,s++,h=!0;break}}if(h)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;
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=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();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 R(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)}})},R=async(e,t,n,r)=>{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 i=new Map;for(let t of e){let e=t.rooftopPz-t.start.z,n=Math.floor(t.rooftopPz-t.start.z);i.has(n)?i.set(n,{num:i.get(n).num+1,totalHeight:i.get(n).totalHeight+e}):i.set(n,{num:1,totalHeight:e})}let a=[...i.entries()].reduce((e,t)=>t[1].num>e[1].num?t:e),o=a[0];a[1].totalHeight/a[1].num;let s=[];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),o=S(n[1],n[2],!0);s.push({doorStartPt:a,doorEndPt:o,boxPoints:n,minZ:t.coordinatesByArea.heightData.minZ,maxZ:t.coordinatesByArea.heightData.maxZ,index:e,isFind:!1})}}for(let t=0;t<s.length;t++){let n=s[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 l=0,u=[],d=[],p=[];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 a=new c.Box3;a.setFromPoints(i.originalPoints);let m=a.max.z,g=a.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<o-.5||w>o+2)continue;let T=0;for(let a of i.mergeCheckRegion){let{facePoints:o,boxArea:E,totalInPlaneNum:D}=y(a.points,a.gridWidth,a.gridHeight,a.centerPts);if(!o||o.length<=0)continue;o[0].distanceTo(o[1]);let O=o[0].distanceTo(o[3]),A=a.gridHeight*(O/a.gridHeight/3*2)+a.gridHeight/2,N=new c.Vector3(o[0].x,o[0].y,o[0].z+A),P=new c.Vector3(o[1].x,o[1].y,o[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,o[0].z)),t.facePoints.push(new c.Vector3(n[1].x,n[1].y,o[0].z)),t.facePoints.push(new c.Vector3(n[1].x,n[1].y,o[2].z)),t.facePoints.push(new c.Vector3(n[0].x,n[0].y,o[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>70&&B.length==1||B.length>1)for(let o=0;o<B.length;o++){let y=B[o],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<1.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}}}if(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{let t=.3,n=_.clone().addScaledVector(o,t),i=_.clone().addScaledVector(o.clone().negate(),t);n.z=i.z=0;let a=!1;for(let o=0;o<e.length;o++){if(o==r)continue;let s=new c.Vector3(e[o].start.x,e[o].start.y,0),l=new c.Vector3(e[o].end.x,e[o].end.y,0),u=new c.Vector3().subVectors(l,s).normalize();if(s=s.clone().addScaledVector(u.clone().negate(),t),l=l.clone().addScaledVector(u,t),a=f(n,i,s,l)!=null,a)break}if(!a)continue}}}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.clone().negate(),10),p=n.clone().addScaledVector(o,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}}}if(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)));else{let t=.3,n=v.clone().addScaledVector(o,t),i=v.clone().addScaledVector(o.clone().negate(),t);n.z=i.z=0;let a=!1;for(let o=0;o<e.length;o++){if(o==r)continue;let s=new c.Vector3(e[o].start.x,e[o].start.y,0),l=new c.Vector3(e[o].end.x,e[o].end.y,0),u=new c.Vector3().subVectors(l,s).normalize();if(s=s.clone().addScaledVector(u.clone().negate(),t),l=l.clone().addScaledVector(u,t),a=f(n,i,s,l)!=null,a)break}if(!a)continue}}}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<s.length;e++){let t=new c.Vector3(s[e].doorStartPt.x,s[e].doorStartPt.y,0),r=new c.Vector3(s[e].doorEndPt.x,s[e].doorEndPt.y,0),i=new c.Vector3(y.facePoints[0].x,y.facePoints[0].y,0),a=new c.Vector3(y.facePoints[1].x,y.facePoints[1].y,0),o=new c.Vector3().subVectors(t,r).normalize(),l=new c.Vector3().subVectors(i,a).normalize(),u=.2,d=f(t.clone().addScaledVector(o,u),r.clone().addScaledVector(o.clone().negate(),u),i.clone().addScaledVector(l,u),a.clone().addScaledVector(l.clone().negate(),u)),p=y.facePoints[2].z<s[e].maxZ||Math.abs(y.facePoints[2].z-s[e].maxZ)<.3;b(t.clone(),r.clone(),i.clone(),a.clone()),d&&p&&(n[s[e].index].isFindBeam=!0,M=!0,s[e].isFind=!0)}for(let e=0;e<s.length;e++){let t=new c.Vector3(s[e].doorStartPt.x,s[e].doorStartPt.y,0),r=new c.Vector3(s[e].doorEndPt.x,s[e].doorEndPt.y,0),a=new c.Vector3(y.facePoints[0].x,y.facePoints[0].y,0),o=new c.Vector3(y.facePoints[1].x,y.facePoints[1].y,0),l=y.facePoints[2].z<s[e].maxZ||Math.abs(y.facePoints[2].z-s[e].maxZ)<.3;if(Math.abs(y.facePoints[2].z-s[e].maxZ),l){let l=C(t.clone(),r.clone(),a.clone(),o.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=o.distanceTo(t);if(i>.25&&u>.25){let e=a.distanceTo(r),t=o.distanceTo(r);if(e>.2&&t>.2)continue}M=!0,N=!0,P.push(s[e].index),s[e].isFind=!0,n[s[e].index].isFindBeam=!0;let d=new c.Vector3().subVectors(a,o).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=I.clone(),R=F.clone(),L.z=R.z=y.facePoints[0].z):(F=y.facePoints[2].clone(),I=y.facePoints[3].clone().addScaledVector(d,f),L=I.clone(),R=F.clone(),L.z=R.z=y.facePoints[0].z)}}else if(l&&l.type==`collinear_overlap`&&l.maxPerpendicularDistance<.35){let i=a.distanceTo(o),c=r.distanceTo(t);(Math.abs(i-l.overlap)<.25||Math.abs(c-l.overlap)<.25)&&(M=!0,N=!0,P.push(s[e].index),s[e].isFind=!0,n[s[e].index].isFindBeam=!0)}}}console.log(`isDoor`,M),console.log(`isInsert`,e),(M||e||E>w/2&&T>.5&&Math.abs(g-i.start.z)<.4)&&(A=!0)}}if(A){let t=!1;for(let e=0;e<p.length;e++){let n=C(p[e].start,p[e].end,L,R);if(n&&n.type==`collinear_overlap`&&n.maxPerpendicularDistance<.01){t=!0;break}}if(t)continue;p.push({start:L,end:R});{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))}console.log(`id:`,l,r),d.includes(r)||d.push(r),i.checkResults[a.index].isDoor=!0;let n={id:l,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<u.length;t++){let r=e[u[t].linesIndex].doorAndBeamData[u[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,s=i.closestPointToPoint(a,!0,o),u=s.distanceTo(a),d=new c.Vector3().subVectors(a,s).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=l,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),u.push({linesIndex:r,doorIndex:i.doorAndBeamData.length-1}),l++,N&&T<2)break}}}let E=i.length*(i.rooftopPz-i.start.z);i.completePointAreaPercentage=T/E*100}if(d.length>1)for(let t=0;t<d.length;t++){let n=[];for(let r=0;r<e[d[t]].doorAndBeamData.length;r++){let i=e[d[t]].doorAndBeamData[r].doorStart.clone(),a=e[d[t]].doorAndBeamData[r].doorEnd.clone();i.z=a.z=0;let o=i.distanceTo(a),s=!1;for(let n=0;n<d.length;n++)if(t!=n)for(let t=0;t<e[d[n]].doorAndBeamData.length;t++){let r=e[d[n]].doorAndBeamData[t].doorStart.clone(),l=e[d[n]].doorAndBeamData[t].doorEnd.clone();r.z=l.z=0,r.distanceTo(l);let u=C(i,a,r,l);if(u.type==`collinear_overlap`&&u.maxPerpendicularDistance<.2){let e=new c.Line3(r,l),t=S(i,a),n=e.closestPointToPoint(t,!0,new c.Vector3),d=t.distanceTo(n);if(Math.abs(o-u.overlap)<.1&&d<.05){s=!0;break}}}s||n.push(e[d[t]].doorAndBeamData[r])}n.length!=e[d[t]].doorAndBeamData.length&&(e[d[t]].doorAndBeamData=n)}let m=.65,g=[],_=[],v=.2;for(let t=0;t<s.length;t++){if(s[t].isFind)continue;let r=new c.Vector3(s[t].doorStartPt.x,s[t].doorStartPt.y,0),i=new c.Vector3(s[t].doorEndPt.x,s[t].doorEndPt.y,0),a=new c.Vector3().subVectors(r,i).normalize();for(let o=0;o<e.length;o++){if(e[o].length<.7||e[o].completePointAreaPercentage>60)continue;let u=new c.Vector3(e[o].start.x,e[o].start.y,0),d=new c.Vector3(e[o].end.x,e[o].end.y,0),p=new c.Vector3().subVectors(u,d).normalize(),h=.3,y=r.clone().addScaledVector(a,h),S=i.clone().addScaledVector(a.clone().negate(),h),w=u.clone().addScaledVector(p,h),T=d.clone().addScaledVector(p.clone().negate(),h);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<v||k<v)&&E&&!E.isParallel){let a=r.distanceTo(u)<r.distanceTo(d)?u:d,p=!1;for(let h=0;h<e.length;h++){if(h==o||g.includes(h)||e[h].length<.7)continue;let y=new c.Vector3(e[h].start.x,e[h].start.y,0),S=new c.Vector3(e[h].end.x,e[h].end.y,0),w=new c.Vector3().subVectors(y,S).normalize(),T=w.clone().negate(),E=y.clone().addScaledVector(w,5),D=S.clone().addScaledVector(T,5),O=x(u,d,y,S),k=b(u,d,y,S),A=f(u,d,E,D),j=!0;if(A!=null&&(j=A.point.distanceTo(u)<v||A.point.distanceTo(d)<v),(O&&O.angle>85||k&&k.isParallel)&&j){let f=new c.Line3(u,d),v=new c.Vector3,x=f.closestPointToPoint(y,!0,v),w=f.closestPointToPoint(S,!0,v),T=x.distanceTo(y),E=w.distanceTo(S),D,O;if(T>E&&E>m)D=S,O=w;else if(T<E&&T>m)D=y,O=x;else continue;if(T=a.distanceTo(y),E=a.distanceTo(S),T>2&&E>2)continue;let k=b(D,O,y,S);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<_.length;e++){let t=C(_[e].start,_[e].end,D,O);if(t&&t.type==`collinear_overlap`&&t.maxPerpendicularDistance<.2){j=!0;break}}if(j)continue;_.push({start:D,end:O}),g.push(h);let M=e[h].rooftopPz-e[h].start.z,N=e[o].rooftopPz-e[o].start.z,P=M<N?M:N;e[h].doorAndBeamData||(e[h].doorAndBeamData=[]),e[h].doorAndBeamData.push({id:l,beamStart:D,beamEnd:D,beamHeight:P,doorStart:D,doorEnd:O,doorHeight:P,nearId:-1,type:`onlyDoor`,isDoor:!0}),n[s[t].index].isFindBeam=!0,n[s[t].index].isFindOnlyDoor=!0,l++,p=!0;break}}if(p)break}}}let w=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 a=n.rooftopPz-n.start.z;if(a<o-.5||a>o+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${w}`;a.mainBeamDataIds||=[],a.isMainBeam||=!0,a.mainBeamDataIds.push(e),n.mainBeamDataIds.push(e),w++;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 CHANGED
@@ -1154,7 +1154,7 @@ var T = (t, n, r, i = .05, a = .02) => {
1154
1154
  }), i.add(n);
1155
1155
  }
1156
1156
  });
1157
- }, L = async (t, n, r) => {
1157
+ }, L = async (t, n, r, i) => {
1158
1158
  if (t.length <= 0 || n.length <= 0) return t;
1159
1159
  r ||= [];
1160
1160
  for (let n of t) {
@@ -1208,20 +1208,20 @@ var T = (t, n, r, i = .05, a = .02) => {
1208
1208
  n.mergeCheckRegion = t;
1209
1209
  }
1210
1210
  }
1211
- let i = /* @__PURE__ */ new Map();
1211
+ let a = /* @__PURE__ */ new Map();
1212
1212
  for (let e of t) {
1213
1213
  let t = e.rooftopPz - e.start.z, n = Math.floor(e.rooftopPz - e.start.z);
1214
- i.has(n) ? i.set(n, {
1215
- num: i.get(n).num + 1,
1216
- totalHeight: i.get(n).totalHeight + t
1217
- }) : i.set(n, {
1214
+ a.has(n) ? a.set(n, {
1215
+ num: a.get(n).num + 1,
1216
+ totalHeight: a.get(n).totalHeight + t
1217
+ }) : a.set(n, {
1218
1218
  num: 1,
1219
1219
  totalHeight: t
1220
1220
  });
1221
1221
  }
1222
- let a = [...i.entries()].reduce((e, t) => t[1].num > e[1].num ? t : e), o = a[0];
1223
- a[1].totalHeight / a[1].num;
1224
- let s = [];
1222
+ let o = [...a.entries()].reduce((e, t) => t[1].num > e[1].num ? t : e), s = o[0];
1223
+ o[1].totalHeight / o[1].num;
1224
+ let c = [];
1225
1225
  for (let t = 0; t < r.length; t++) {
1226
1226
  let n = r[t];
1227
1227
  if (n.isFindBeam = !1, n.category == "door") {
@@ -1236,10 +1236,10 @@ var T = (t, n, r, i = .05, a = .02) => {
1236
1236
  r[2],
1237
1237
  r[1]
1238
1238
  ]);
1239
- let o = x(r[0], r[3], !0), c = x(r[1], r[2], !0);
1240
- s.push({
1239
+ let o = x(r[0], r[3], !0), s = x(r[1], r[2], !0);
1240
+ c.push({
1241
1241
  doorStartPt: o,
1242
- doorEndPt: c,
1242
+ doorEndPt: s,
1243
1243
  boxPoints: r,
1244
1244
  minZ: n.coordinatesByArea.heightData.minZ,
1245
1245
  maxZ: n.coordinatesByArea.heightData.maxZ,
@@ -1248,8 +1248,8 @@ var T = (t, n, r, i = .05, a = .02) => {
1248
1248
  });
1249
1249
  }
1250
1250
  }
1251
- for (let n = 0; n < s.length; n++) {
1252
- let r = s[n];
1251
+ for (let n = 0; n < c.length; n++) {
1252
+ let r = c[n];
1253
1253
  for (let n = 0; n < t.length; n++) {
1254
1254
  if (t[n].length < .8) continue;
1255
1255
  let i = [], a = x(t[n].start, t[n].end, !0);
@@ -1272,26 +1272,21 @@ var T = (t, n, r, i = .05, a = .02) => {
1272
1272
  } else continue;
1273
1273
  }
1274
1274
  }
1275
- let c = 0, l = [];
1276
- new e.PointsMaterial({
1277
- color: 65535,
1278
- size: .01
1279
- });
1280
- let u = [], f = [];
1275
+ let l = 0, u = [], f = [], p = [];
1281
1276
  for (let i = 0; i < t.length; i++) {
1282
1277
  let a = t[i];
1283
1278
  if (!a.mergeCheckRegion || a.mergeCheckRegion.length == 0) continue;
1284
1279
  a.doorAndBeamData = [], a.completePointAreaPercentage = -1;
1285
- let p = new e.Box3();
1286
- p.setFromPoints(a.originalPoints);
1287
- 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;
1288
- if (w < o - .5 || w > o + 2) continue;
1280
+ let o = new e.Box3();
1281
+ o.setFromPoints(a.originalPoints);
1282
+ let h = o.max.z, g = o.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;
1283
+ if (w < s - .5 || w > s + 2) continue;
1289
1284
  let T = 0;
1290
1285
  for (let o of a.mergeCheckRegion) {
1291
- let { facePoints: p, boxArea: E, totalInPlaneNum: D } = v(o.points, o.gridWidth, o.gridHeight, o.centerPts);
1292
- if (!p || p.length <= 0) continue;
1293
- p[0].distanceTo(p[1]);
1294
- let k = p[0].distanceTo(p[3]), M = o.gridHeight * (k / o.gridHeight / 3 * 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();
1286
+ let { facePoints: s, boxArea: E, totalInPlaneNum: D } = v(o.points, o.gridWidth, o.gridHeight, o.centerPts);
1287
+ if (!s || s.length <= 0) continue;
1288
+ s[0].distanceTo(s[1]);
1289
+ let k = s[0].distanceTo(s[3]), M = o.gridHeight * (k / o.gridHeight / 3 * 2) + o.gridHeight / 2, N = new e.Vector3(s[0].x, s[0].y, s[0].z + M), P = new e.Vector3(s[1].x, s[1].y, s[1].z + M), F = new e.Line3(N, P), I = [], L = new e.Vector3();
1295
1290
  o.points.forEach((e) => {
1296
1291
  F.closestPointToPoint(e, !0, L).distanceTo(e) < o.gridHeight + .01 && I.push(e);
1297
1292
  });
@@ -1315,13 +1310,13 @@ var T = (t, n, r, i = .05, a = .02) => {
1315
1310
  lineEd: new e.Vector3(),
1316
1311
  facePoints: []
1317
1312
  }, r = j(R[t]);
1318
- 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));
1313
+ 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, s[0].z)), n.facePoints.push(new e.Vector3(r[1].x, r[1].y, s[0].z)), n.facePoints.push(new e.Vector3(r[1].x, r[1].y, s[2].z)), n.facePoints.push(new e.Vector3(r[0].x, r[0].y, s[2].z)), B.push(n));
1319
1314
  }
1320
1315
  T += o.area;
1321
1316
  let V = 0;
1322
- 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++) {
1323
- 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;
1324
- if (g > a.start.z - .15 && a.length < 1.5 || a.length, E < .5 || h < v.facePoints[2].z || k > w / 3 * 2) continue;
1317
+ if (D != 0 && (V = o.gridHeight * o.gridWidth * D), (o.area - V) / E * 100 > 70 && B.length == 1 || B.length > 1) for (let s = 0; s < B.length; s++) {
1318
+ let v = B[s], 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;
1319
+ if (g > a.start.z - .15 && a.length < 1.5 || a.length, E < 1.5 || h < v.facePoints[2].z || k > w / 3 * 2) continue;
1325
1320
  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;
1326
1321
  if (Math.abs(v.facePoints[0].z - a.start.z) < .25) {
1327
1322
  if (T > .2 && a.length > .5) {
@@ -1352,7 +1347,18 @@ var T = (t, n, r, i = .05, a = .02) => {
1352
1347
  }
1353
1348
  }
1354
1349
  }
1355
- 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)));
1350
+ if (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)));
1351
+ else {
1352
+ let n = .3, r = _.clone().addScaledVector(s, n), a = _.clone().addScaledVector(s.clone().negate(), n);
1353
+ r.z = a.z = 0;
1354
+ let o = !1;
1355
+ for (let s = 0; s < t.length; s++) {
1356
+ if (s == i) continue;
1357
+ let c = new e.Vector3(t[s].start.x, t[s].start.y, 0), l = new e.Vector3(t[s].end.x, t[s].end.y, 0), u = new e.Vector3().subVectors(l, c).normalize();
1358
+ if (c = c.clone().addScaledVector(u.clone().negate(), n), l = l.clone().addScaledVector(u, n), o = d(r, a, c, l) != null, o) break;
1359
+ }
1360
+ if (!o) continue;
1361
+ }
1356
1362
  }
1357
1363
  } else if (o < c || s < c) {
1358
1364
  let { newOriginalPoints: n } = O(L, R, a.originalPoints), r = [], s = [];
@@ -1368,7 +1374,7 @@ var T = (t, n, r, i = .05, a = .02) => {
1368
1374
  if (p && p.angle > 85) {
1369
1375
  let t = y(r, n, new e.Line3(o, c).closestPointToPoint(r, !0, l), r);
1370
1376
  if (!t || t.direction != "same") continue;
1371
- 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));
1377
+ let i = n.clone().addScaledVector(s.clone().negate(), 10), p = r.clone().addScaledVector(s, 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));
1372
1378
  if (h == null) continue;
1373
1379
  let g = r.distanceTo(h.point);
1374
1380
  g < u && (u = g, f = a);
@@ -1381,7 +1387,18 @@ var T = (t, n, r, i = .05, a = .02) => {
1381
1387
  }
1382
1388
  }
1383
1389
  }
1384
- 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)));
1390
+ if (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)));
1391
+ else {
1392
+ let n = .3, r = C.clone().addScaledVector(s, n), a = C.clone().addScaledVector(s.clone().negate(), n);
1393
+ r.z = a.z = 0;
1394
+ let o = !1;
1395
+ for (let s = 0; s < t.length; s++) {
1396
+ if (s == i) continue;
1397
+ let c = new e.Vector3(t[s].start.x, t[s].start.y, 0), l = new e.Vector3(t[s].end.x, t[s].end.y, 0), u = new e.Vector3().subVectors(l, c).normalize();
1398
+ if (c = c.clone().addScaledVector(u.clone().negate(), n), l = l.clone().addScaledVector(u, n), o = d(r, a, c, l) != null, o) break;
1399
+ }
1400
+ if (!o) continue;
1401
+ }
1385
1402
  }
1386
1403
  }
1387
1404
  T = I.distanceTo(F);
@@ -1397,50 +1414,46 @@ var T = (t, n, r, i = .05, a = .02) => {
1397
1414
  break;
1398
1415
  }
1399
1416
  }
1400
- for (let t = 0; t < s.length; t++) {
1401
- 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());
1402
- if (p && m) r[s[t].index].isFindBeam = !0, M = !0, s[t].isFind = !0;
1403
- else if (m && h && h.isParallel) {
1404
- 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);
1405
- o.distanceTo(u) < .3 && (r[s[t].index].isFindBeam = !0);
1406
- }
1417
+ for (let t = 0; t < c.length; t++) {
1418
+ let n = new e.Vector3(c[t].doorStartPt.x, c[t].doorStartPt.y, 0), i = new e.Vector3(c[t].doorEndPt.x, c[t].doorEndPt.y, 0), a = new e.Vector3(v.facePoints[0].x, v.facePoints[0].y, 0), o = new e.Vector3(v.facePoints[1].x, v.facePoints[1].y, 0), s = new e.Vector3().subVectors(n, i).normalize(), l = new e.Vector3().subVectors(a, o).normalize(), u = .2, f = d(n.clone().addScaledVector(s, u), i.clone().addScaledVector(s.clone().negate(), u), a.clone().addScaledVector(l, u), o.clone().addScaledVector(l.clone().negate(), u)), p = v.facePoints[2].z < c[t].maxZ || Math.abs(v.facePoints[2].z - c[t].maxZ) < .3;
1419
+ y(n.clone(), i.clone(), a.clone(), o.clone()), f && p && (r[c[t].index].isFindBeam = !0, M = !0, c[t].isFind = !0);
1407
1420
  }
1408
- for (let t = 0; t < s.length; t++) {
1409
- 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;
1410
- if (Math.abs(v.facePoints[2].z - s[t].maxZ), l) {
1411
- let l = S(n.clone(), i.clone(), o.clone(), c.clone(), { distanceThreshold: 2 });
1421
+ for (let t = 0; t < c.length; t++) {
1422
+ let n = new e.Vector3(c[t].doorStartPt.x, c[t].doorStartPt.y, 0), i = new e.Vector3(c[t].doorEndPt.x, c[t].doorEndPt.y, 0), o = new e.Vector3(v.facePoints[0].x, v.facePoints[0].y, 0), s = new e.Vector3(v.facePoints[1].x, v.facePoints[1].y, 0), l = v.facePoints[2].z < c[t].maxZ || Math.abs(v.facePoints[2].z - c[t].maxZ) < .3;
1423
+ if (Math.abs(v.facePoints[2].z - c[t].maxZ), l) {
1424
+ let l = S(n.clone(), i.clone(), o.clone(), s.clone(), { distanceThreshold: 2 });
1412
1425
  if (l && l.type == "collinear_gap") {
1413
1426
  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);
1414
1427
  if (u.distanceTo(p) < .15) {
1415
- let a = o.distanceTo(n), u = c.distanceTo(n);
1428
+ let a = o.distanceTo(n), u = s.distanceTo(n);
1416
1429
  if (a > .25 && u > .25) {
1417
- let e = o.distanceTo(i), t = c.distanceTo(i);
1430
+ let e = o.distanceTo(i), t = s.distanceTo(i);
1418
1431
  if (e > .2 && t > .2) continue;
1419
1432
  }
1420
- M = !0, N = !0, P.push(s[t].index), s[t].isFind = !0, r[s[t].index].isFindBeam = !0;
1421
- let d = new e.Vector3().subVectors(o, c).normalize(), f = n.distanceTo(i) - l.overlap;
1422
- 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);
1433
+ M = !0, N = !0, P.push(c[t].index), c[t].isFind = !0, r[c[t].index].isFindBeam = !0;
1434
+ let d = new e.Vector3().subVectors(o, s).normalize(), f = n.distanceTo(i) - l.overlap;
1435
+ u < a && d.negate(), a > u ? (F = v.facePoints[2].clone().addScaledVector(d, f), I = v.facePoints[3].clone(), L = I.clone(), R = F.clone(), L.z = R.z = v.facePoints[0].z) : (F = v.facePoints[2].clone(), I = v.facePoints[3].clone().addScaledVector(d, f), L = I.clone(), R = F.clone(), L.z = R.z = v.facePoints[0].z);
1423
1436
  }
1424
1437
  } else if (l && l.type == "collinear_overlap" && l.maxPerpendicularDistance < .35) {
1425
- let e = o.distanceTo(c), a = i.distanceTo(n);
1426
- (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);
1438
+ let e = o.distanceTo(s), a = i.distanceTo(n);
1439
+ (Math.abs(e - l.overlap) < .25 || Math.abs(a - l.overlap) < .25) && (M = !0, N = !0, P.push(c[t].index), c[t].isFind = !0, r[c[t].index].isFindBeam = !0);
1427
1440
  }
1428
1441
  }
1429
1442
  }
1430
- (M || t || E > w / 2 && T > .8 && Math.abs(g - a.start.z) < .4) && (j = !0);
1443
+ console.log("isDoor", M), console.log("isInsert", t), (M || t || E > w / 2 && T > .5 && Math.abs(g - a.start.z) < .4) && (j = !0);
1431
1444
  }
1432
1445
  }
1433
1446
  if (j) {
1434
1447
  let n = !1;
1435
- for (let e = 0; e < f.length; e++) {
1436
- let t = S(f[e].start, f[e].end, L, R);
1448
+ for (let e = 0; e < p.length; e++) {
1449
+ let t = S(p[e].start, p[e].end, L, R);
1437
1450
  if (t && t.type == "collinear_overlap" && t.maxPerpendicularDistance < .01) {
1438
1451
  n = !0;
1439
1452
  break;
1440
1453
  }
1441
1454
  }
1442
1455
  if (n) continue;
1443
- f.push({
1456
+ p.push({
1444
1457
  start: L,
1445
1458
  end: R
1446
1459
  });
@@ -1464,9 +1477,9 @@ var T = (t, n, r, i = .05, a = .02) => {
1464
1477
  }
1465
1478
  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));
1466
1479
  }
1467
- u.includes(i) || u.push(i), a.checkResults[o.index].isDoor = !0;
1480
+ console.log("id:", l, i), f.includes(i) || f.push(i), a.checkResults[o.index].isDoor = !0;
1468
1481
  let r = {
1469
- id: c,
1482
+ id: l,
1470
1483
  beamStart: F,
1471
1484
  beamEnd: I,
1472
1485
  beamHeight: a.rooftopPz - v.facePoints[2].z,
@@ -1479,8 +1492,8 @@ var T = (t, n, r, i = .05, a = .02) => {
1479
1492
  isPullOutDoor: N,
1480
1493
  pcbDoorIndexs: P
1481
1494
  };
1482
- for (let n = 0; n < l.length; n++) {
1483
- let i = t[l[n].linesIndex].doorAndBeamData[l[n].doorIndex];
1495
+ for (let n = 0; n < u.length; n++) {
1496
+ let i = t[u[n].linesIndex].doorAndBeamData[u[n].doorIndex];
1484
1497
  if (i.nearId != -1) continue;
1485
1498
  let a = y(F.clone(), I.clone(), i.beamStart.clone(), i.beamEnd.clone());
1486
1499
  if (a && a.isParallel) {
@@ -1488,9 +1501,9 @@ var T = (t, n, r, i = .05, a = .02) => {
1488
1501
  a.start.z = 0, a.end.z = 0;
1489
1502
  let o = x(i.beamStart, i.beamEnd, !0);
1490
1503
  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));
1491
- 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);
1504
+ let s = new e.Vector3(), c = a.closestPointToPoint(o, !0, s), u = c.distanceTo(o), f = new e.Vector3().subVectors(o, c).normalize(), p = f.clone().negate(), m = o.clone().addScaledVector(p, u + .1), h = d(a.start, a.end, o, m);
1492
1505
  if (u < .4 && h != null) {
1493
- i.nearId = c, r.nearId = i.id;
1506
+ i.nearId = l, r.nearId = i.id;
1494
1507
  let e = .8;
1495
1508
  (i.isDoor || r.isDoor) && (i.isDoor = !0, r.isDoor = !0);
1496
1509
  let a = P.some((e) => i.pcbDoorIndexs.includes(e));
@@ -1499,112 +1512,112 @@ var T = (t, n, r, i = .05, a = .02) => {
1499
1512
  }
1500
1513
  }
1501
1514
  }
1502
- if (a.doorAndBeamData.push(r), l.push({
1515
+ if (a.doorAndBeamData.push(r), u.push({
1503
1516
  linesIndex: i,
1504
1517
  doorIndex: a.doorAndBeamData.length - 1
1505
- }), c++, N && T < 2) break;
1518
+ }), l++, N && T < 2) break;
1506
1519
  }
1507
1520
  }
1508
1521
  }
1509
1522
  let E = a.length * (a.rooftopPz - a.start.z);
1510
1523
  a.completePointAreaPercentage = T / E * 100;
1511
1524
  }
1512
- if (u.length > 1) for (let n = 0; n < u.length; n++) {
1525
+ if (f.length > 1) for (let n = 0; n < f.length; n++) {
1513
1526
  let r = [];
1514
- for (let i = 0; i < t[u[n]].doorAndBeamData.length; i++) {
1515
- let a = t[u[n]].doorAndBeamData[i].doorStart.clone(), o = t[u[n]].doorAndBeamData[i].doorEnd.clone();
1527
+ for (let i = 0; i < t[f[n]].doorAndBeamData.length; i++) {
1528
+ let a = t[f[n]].doorAndBeamData[i].doorStart.clone(), o = t[f[n]].doorAndBeamData[i].doorEnd.clone();
1516
1529
  a.z = o.z = 0;
1517
1530
  let s = a.distanceTo(o), c = !1;
1518
- for (let r = 0; r < u.length; r++) if (n != r) for (let n = 0; n < t[u[r]].doorAndBeamData.length; n++) {
1519
- let i = t[u[r]].doorAndBeamData[n].doorStart.clone(), l = t[u[r]].doorAndBeamData[n].doorEnd.clone();
1531
+ for (let r = 0; r < f.length; r++) if (n != r) for (let n = 0; n < t[f[r]].doorAndBeamData.length; n++) {
1532
+ let i = t[f[r]].doorAndBeamData[n].doorStart.clone(), l = t[f[r]].doorAndBeamData[n].doorEnd.clone();
1520
1533
  i.z = l.z = 0, i.distanceTo(l);
1521
- let d = S(a, o, i, l);
1522
- if (d.type == "collinear_overlap" && d.maxPerpendicularDistance < .2) {
1523
- let t = new e.Line3(i, l), n = x(a, o), r = t.closestPointToPoint(n, !0, new e.Vector3()), u = n.distanceTo(r);
1524
- if (Math.abs(s - d.overlap) < .1 && u < .05) {
1534
+ let u = S(a, o, i, l);
1535
+ if (u.type == "collinear_overlap" && u.maxPerpendicularDistance < .2) {
1536
+ let t = new e.Line3(i, l), n = x(a, o), r = t.closestPointToPoint(n, !0, new e.Vector3()), d = n.distanceTo(r);
1537
+ if (Math.abs(s - u.overlap) < .1 && d < .05) {
1525
1538
  c = !0;
1526
1539
  break;
1527
1540
  }
1528
1541
  }
1529
1542
  }
1530
- c || r.push(t[u[n]].doorAndBeamData[i]);
1543
+ c || r.push(t[f[n]].doorAndBeamData[i]);
1531
1544
  }
1532
- r.length != t[u[n]].doorAndBeamData.length && (t[u[n]].doorAndBeamData = r);
1533
- }
1534
- let p = .65, h = [], g = [], _ = .2;
1535
- for (let n = 0; n < s.length; n++) {
1536
- if (s[n].isFind) continue;
1537
- 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();
1538
- for (let l = 0; l < t.length; l++) {
1539
- if (t[l].length < .7 || t[l].completePointAreaPercentage > 60) continue;
1540
- 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);
1541
- d(x, C, w, T);
1542
- let E = y(x, C, w, T), D = i.distanceTo(u) > a.distanceTo(u) ? a : i, O = D.distanceTo(u), k = D.distanceTo(f);
1543
- if ((O < _ || k < _) && E && !E.isParallel) {
1544
- let o = i.distanceTo(u) < i.distanceTo(f) ? u : f, m = !1;
1545
- for (let v = 0; v < t.length; v++) {
1546
- if (v == l || h.includes(v) || t[v].length < .7) continue;
1547
- let x = new e.Vector3(t[v].start.x, t[v].start.y, 0), C = new e.Vector3(t[v].end.x, t[v].end.y, 0), w = new e.Vector3().subVectors(x, C).normalize(), T = w.clone().negate(), E = x.clone().addScaledVector(w, 5), D = C.clone().addScaledVector(T, 5), A = b(u, f, x, C), j = y(u, f, x, C), M = d(u, f, E, D), N = !0;
1548
- if (M != null && (N = M.point.distanceTo(u) < _ || M.point.distanceTo(f) < _), (A && A.angle > 85 || j && j.isParallel) && N) {
1549
- let d = new e.Line3(u, f), _ = new e.Vector3(), b = d.closestPointToPoint(x, !0, _), w = d.closestPointToPoint(C, !0, _), T = b.distanceTo(x), E = w.distanceTo(C), D, A;
1550
- if (T > E && E > p) D = C, A = w;
1551
- else if (T < E && T > p) D = x, A = b;
1545
+ r.length != t[f[n]].doorAndBeamData.length && (t[f[n]].doorAndBeamData = r);
1546
+ }
1547
+ let h = .65, g = [], _ = [], C = .2;
1548
+ for (let n = 0; n < c.length; n++) {
1549
+ if (c[n].isFind) continue;
1550
+ let i = new e.Vector3(c[n].doorStartPt.x, c[n].doorStartPt.y, 0), a = new e.Vector3(c[n].doorEndPt.x, c[n].doorEndPt.y, 0), o = new e.Vector3().subVectors(i, a).normalize();
1551
+ for (let s = 0; s < t.length; s++) {
1552
+ if (t[s].length < .7 || t[s].completePointAreaPercentage > 60) continue;
1553
+ let u = new e.Vector3(t[s].start.x, t[s].start.y, 0), f = new e.Vector3(t[s].end.x, t[s].end.y, 0), p = new e.Vector3().subVectors(u, f).normalize(), m = .3, v = i.clone().addScaledVector(o, m), x = a.clone().addScaledVector(o.clone().negate(), m), w = u.clone().addScaledVector(p, m), T = f.clone().addScaledVector(p.clone().negate(), m);
1554
+ d(v, x, w, T);
1555
+ let E = y(v, x, w, T), D = i.distanceTo(u) > a.distanceTo(u) ? a : i, O = D.distanceTo(u), k = D.distanceTo(f);
1556
+ if ((O < C || k < C) && E && !E.isParallel) {
1557
+ let o = i.distanceTo(u) < i.distanceTo(f) ? u : f, p = !1;
1558
+ for (let m = 0; m < t.length; m++) {
1559
+ if (m == s || g.includes(m) || t[m].length < .7) continue;
1560
+ let v = new e.Vector3(t[m].start.x, t[m].start.y, 0), x = new e.Vector3(t[m].end.x, t[m].end.y, 0), w = new e.Vector3().subVectors(v, x).normalize(), T = w.clone().negate(), E = v.clone().addScaledVector(w, 5), D = x.clone().addScaledVector(T, 5), O = b(u, f, v, x), k = y(u, f, v, x), A = d(u, f, E, D), j = !0;
1561
+ if (A != null && (j = A.point.distanceTo(u) < C || A.point.distanceTo(f) < C), (O && O.angle > 85 || k && k.isParallel) && j) {
1562
+ let d = new e.Line3(u, f), b = new e.Vector3(), C = d.closestPointToPoint(v, !0, b), w = d.closestPointToPoint(x, !0, b), T = C.distanceTo(v), E = w.distanceTo(x), D, O;
1563
+ if (T > E && E > h) D = x, O = w;
1564
+ else if (T < E && T > h) D = v, O = C;
1552
1565
  else continue;
1553
- if (T = o.distanceTo(x), E = o.distanceTo(C), T > 2 && E > 2) continue;
1554
- let j = y(D, A, x, C);
1555
- if (!j || !j.isParallel) continue;
1556
- let M = S(i, a, D, A);
1557
- if (M && M.type == "collinear_overlap" && M.maxPerpendicularDistance < .2) continue;
1558
- D.z = A.z = t[v].start.z;
1559
- let N = !1;
1560
- for (let e = 0; e < g.length; e++) {
1561
- let t = S(g[e].start, g[e].end, D, A);
1566
+ if (T = o.distanceTo(v), E = o.distanceTo(x), T > 2 && E > 2) continue;
1567
+ let k = y(D, O, v, x);
1568
+ if (!k || !k.isParallel) continue;
1569
+ let A = S(i, a, D, O);
1570
+ if (A && A.type == "collinear_overlap" && A.maxPerpendicularDistance < .2) continue;
1571
+ D.z = O.z = t[m].start.z;
1572
+ let j = !1;
1573
+ for (let e = 0; e < _.length; e++) {
1574
+ let t = S(_[e].start, _[e].end, D, O);
1562
1575
  if (t && t.type == "collinear_overlap" && t.maxPerpendicularDistance < .2) {
1563
- N = !0;
1576
+ j = !0;
1564
1577
  break;
1565
1578
  }
1566
1579
  }
1567
- if (N) continue;
1568
- g.push({
1580
+ if (j) continue;
1581
+ _.push({
1569
1582
  start: D,
1570
- end: A
1571
- }), h.push(v), console.log("doorDatas[n]", s[n], l, v, O, k);
1572
- let P = t[v].rooftopPz - t[v].start.z, F = t[l].rooftopPz - t[l].start.z, I = P < F ? P : F;
1573
- t[v].doorAndBeamData || (t[v].doorAndBeamData = []), t[v].doorAndBeamData.push({
1574
- id: c,
1583
+ end: O
1584
+ }), g.push(m);
1585
+ let M = t[m].rooftopPz - t[m].start.z, N = t[s].rooftopPz - t[s].start.z, P = M < N ? M : N;
1586
+ t[m].doorAndBeamData || (t[m].doorAndBeamData = []), t[m].doorAndBeamData.push({
1587
+ id: l,
1575
1588
  beamStart: D,
1576
1589
  beamEnd: D,
1577
- beamHeight: I,
1590
+ beamHeight: P,
1578
1591
  doorStart: D,
1579
- doorEnd: A,
1580
- doorHeight: I,
1592
+ doorEnd: O,
1593
+ doorHeight: P,
1581
1594
  nearId: -1,
1582
1595
  type: "onlyDoor",
1583
1596
  isDoor: !0
1584
- }), r[s[n].index].isFindBeam = !0, r[s[n].index].isFindOnlyDoor = !0, c++, m = !0;
1597
+ }), r[c[n].index].isFindBeam = !0, r[c[n].index].isFindOnlyDoor = !0, l++, p = !0;
1585
1598
  break;
1586
1599
  }
1587
1600
  }
1588
- if (m) break;
1601
+ if (p) break;
1589
1602
  }
1590
1603
  }
1591
1604
  }
1592
- let C = 0;
1605
+ let w = 0;
1593
1606
  for (let n = 0; n < t.length; n++) {
1594
1607
  let r = t[n];
1595
1608
  r.mainBeamDataIds ||= [], r.isMainBeam ||= !1;
1596
1609
  let { minZ: i, maxZ: a } = D(r);
1597
1610
  if (i == null || a == null) continue;
1598
1611
  a - i;
1599
- let s = r.rooftopPz - r.start.z;
1600
- if (s < o - .5 || s > o + 2 || Math.abs(i - r.start.z) > .3 || r.length < 1) continue;
1612
+ let o = r.rooftopPz - r.start.z;
1613
+ if (o < s - .5 || o > s + 2 || Math.abs(i - r.start.z) > .3 || r.length < 1) continue;
1601
1614
  let c = new e.Vector3(r.start.x, r.start.y, 0), l = new e.Vector3(r.end.x, r.end.y, 0);
1602
1615
  for (let i = 0; i < t.length; i++) {
1603
1616
  if (i == n || t[i].length < 1) continue;
1604
1617
  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);
1605
1618
  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) {
1606
- let e = `mainBeamId${C}`;
1607
- o.mainBeamDataIds ||= [], o.isMainBeam ||= !0, o.mainBeamDataIds.push(e), r.mainBeamDataIds.push(e), C++;
1619
+ let e = `mainBeamId${w}`;
1620
+ o.mainBeamDataIds ||= [], o.isMainBeam ||= !0, o.mainBeamDataIds.push(e), r.mainBeamDataIds.push(e), w++;
1608
1621
  break;
1609
1622
  }
1610
1623
  }
package/package.json CHANGED
@@ -1,6 +1,6 @@
1
1
  {
2
2
  "name": "rm-graphical-computing",
3
- "version": "1.0.10",
3
+ "version": "1.0.12",
4
4
  "description": "Three.js wall segment processing utilities.",
5
5
  "main": "dist/index.cjs",
6
6
  "module": "dist/index.js",