rm-graphical-computing 1.0.7 → 1.0.9
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 +1 -1
- package/dist/index.js +69 -33
- package/package.json +1 -1
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/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)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=f(n.clone().addScaledVector(o,10),a.clone().addScaledVector(o.clone().negate(),10),i,s);if(r==null)continue;let p=n.distanceTo(r.point);p<u&&(u=p,d=t)}else if(p&&p.angle<5){let e=j(n,i),t=j(n,s),r=j(a,i),o=j(a,s);if(e<.1||t<.1||r<.1||o<.1){u=1/0,d=-1;break}}}u!=1/0&&u<2&&e[d].length>.5&&(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=f(t.clone().addScaledVector(o,10),n.clone().addScaledVector(o.clone().negate(),10),a,s);if(r==null)continue;let p=n.distanceTo(r.point);p<u&&(u=p,d=i)}if(p&&p.angle<5){let e=j(t,a),r=j(t,s),i=j(n,a),o=j(n,s);if(e<.1||r<.1||i<.1||o<.1){u=1/0,d=-1;break}}}u!=1/0&&u<2&&e[d].length>.5&&(z=!0,a<s?(I=I.addScaledVector(o,u),L=L.addScaledVector(o,u)):(F=F.addScaledVector(o,u),R=R.addScaledVector(o,u)))}}T=I.distanceTo(F)}if(D<85&&T>.35||T>.7&&E>w-w/3){let e=!1;for(let n=0;n<t.length;n++){let r=n+1;if(r>=t.length)continue;let i=new c.Vector3(t[n].x,t[n].y,t[n].z),a=new c.Vector3(t[r].x,t[r].y,t[r].z);if(f(new c.Vector3(i.x,i.y,0),new c.Vector3(a.x,a.y,0),new c.Vector3(y.facePoints[0].x,y.facePoints[0].y,0),new c.Vector3(y.facePoints[1].x,y.facePoints[1].y,0),.1)!=null){e=!0;break}}for(let e=0;e<o.length;e++){let t=new c.Vector3(o[e].doorStartPt.x,o[e].doorStartPt.y,0),r=new c.Vector3(o[e].doorEndPt.x,o[e].doorEndPt.y,0),a=new c.Vector3(y.facePoints[0].x,y.facePoints[0].y,0),s=new c.Vector3(y.facePoints[1].x,y.facePoints[1].y,0),l=new c.Vector3().subVectors(t,r).normalize(),u=new c.Vector3().subVectors(a,s).normalize(),d=.2,p=f(t.clone().addScaledVector(l,d),r.clone().addScaledVector(l.clone().negate(),d),a.clone().addScaledVector(u,d),s.clone().addScaledVector(u.clone().negate(),d)),m=y.facePoints[2].z<o[e].maxZ||Math.abs(y.facePoints[2].z-o[e].maxZ)<.3,h=b(t.clone(),r.clone(),a.clone(),s.clone());if(p&&m)n[o[e].index].isFindBeam=!0,M=!0,o[e].isFind=!0;else if(m&&h&&h.isParallel){let a=S(t,r),s=new c.Line3(new c.Vector3(i.start.x,i.start.y,0),new c.Vector3(i.end.x,i.end.y,0)),l=new c.Vector3,u=s.closestPointToPoint(a,!0,l);a.distanceTo(u)<.3&&(n[o[e].index].isFindBeam=!0)}}for(let e=0;e<o.length;e++){let t=new c.Vector3(o[e].doorStartPt.x,o[e].doorStartPt.y,0),r=new c.Vector3(o[e].doorEndPt.x,o[e].doorEndPt.y,0),a=new c.Vector3(y.facePoints[0].x,y.facePoints[0].y,0),s=new c.Vector3(y.facePoints[1].x,y.facePoints[1].y,0),l=y.facePoints[2].z<o[e].maxZ||Math.abs(y.facePoints[2].z-o[e].maxZ)<.3,u=Math.abs(y.facePoints[2].z-o[e].maxZ)>.2;if(l&&!u){let l=C(t.clone(),r.clone(),a.clone(),s.clone(),{distanceThreshold:2});if(l&&l.type==`collinear_gap`){let u=S(t,r),d=new c.Line3(new c.Vector3(i.start.x,i.start.y,0),new c.Vector3(i.end.x,i.end.y,0)),f=new c.Vector3,p=d.closestPointToPoint(u,!0,f);if(u.distanceTo(p)<.15){let i=a.distanceTo(t),u=s.distanceTo(t);if(i>.25&&u>.25){let e=a.distanceTo(r),t=s.distanceTo(r);if(e>.2&&t>.2)continue}M=!0,N=!0,P.push(o[e].index),o[e].isFind=!0,n[o[e].index].isFindBeam=!0;let d=new c.Vector3().subVectors(a,s).normalize(),f=t.distanceTo(r)-l.overlap;u<i&&d.negate(),i>u?(F=y.facePoints[2].clone().addScaledVector(d,f),I=y.facePoints[3].clone(),L=F.clone(),R=I.clone(),L.z=R.z=y.facePoints[0].z):(F=y.facePoints[2].clone(),I=y.facePoints[3].clone().addScaledVector(d,f),L=F.clone(),R=I.clone(),L.z=R.z=y.facePoints[0].z)}}else if(l&&l.type==`collinear_overlap`&&l.maxPerpendicularDistance<.35){let i=a.distanceTo(s),c=r.distanceTo(t);(Math.abs(i-l.overlap)<.25||Math.abs(c-l.overlap)<.25)&&(M=!0,N=!0,P.push(o[e].index),o[e].isFind=!0,n[o[e].index].isFindBeam=!0)}}}(M||e||E>w/2&&T>.8&&Math.abs(g-i.start.z)<.4)&&(A=!0)}}if(A){let t=!1;for(let e=0;e<d.length;e++){let n=C(d[e].start,d[e].end,L,R);if(n&&n.type==`collinear_overlap`&&n.maxPerpendicularDistance<.01){t=!0;break}}if(t)continue;if(d.push({start:L,end:R}),!N){let e=new c.Vector3().subVectors(L,R),t=new c.Vector3().subVectors(R,L),n=1/0,r=1/0,a=L.clone(),o=R.clone(),s=[];if(i.originalPoints.forEach(e=>{e.z-L.z<E/2+.1&&e.z-L.z>E/2-.1&&s.push(e)}),s.length>0){let{newProjectPoints:i}=k(L,R,s);for(let e of i)e.z=L.z;i.forEach(i=>{let s=new c.Vector3(i.x,i.y,L.z),l=new c.Vector3().subVectors(s,R),u=h(e,new c.Vector3().subVectors(s,L)),d=h(t,l);if(u&&u.parallel&&u.sameDirection){let e=s.distanceTo(L);e<n&&(n=e,a=s)}else if(d&&d.parallel&&d.sameDirection){let e=s.distanceTo(R);e<r&&(r=e,o=s)}})}a.equals(L)||(L=a,I=new c.Vector3(L.x,L.y,I.z)),o.equals(R)||(R=o,F=new c.Vector3(R.x,R.y,F.z))}u.includes(r)||u.push(r),i.checkResults[a.index].isDoor=!0;let n={id:s,beamStart:F,beamEnd:I,beamHeight:i.rooftopPz-y.facePoints[2].z,doorStart:L,doorEnd:R,doorHeight:E,nearId:-1,type:z?`extendBeam`:`doorBeam`,isDoor:M,isPullOutDoor:N,pcbDoorIndexs:P};for(let t=0;t<l.length;t++){let r=e[l[t].linesIndex].doorAndBeamData[l[t].doorIndex];if(r.nearId!=-1)continue;let i=b(F.clone(),I.clone(),r.beamStart.clone(),r.beamEnd.clone());if(i&&i.isParallel){let e=F.distanceTo(I),t=r.beamStart.distanceTo(r.beamEnd),i=new c.Line3(F.clone(),I.clone());i.start.z=0,i.end.z=0;let a=S(r.beamStart,r.beamEnd,!0);t>e&&(i.start=new c.Vector3(r.beamStart.x,r.beamStart.y,0),i.end=new c.Vector3(r.beamEnd.x,r.beamEnd.y,0),a=S(F,I,!0));let o=new c.Vector3,l=i.closestPointToPoint(a,!0,o),u=l.distanceTo(a),d=new c.Vector3().subVectors(a,l).normalize(),p=d.clone().negate(),m=a.clone().addScaledVector(p,u+.1),h=f(i.start,i.end,a,m);if(u<.4&&h!=null){r.nearId=s,n.nearId=r.id;let i=.8;(r.isDoor||n.isDoor)&&(r.isDoor=!0,n.isDoor=!0);let a=P.some(e=>r.pcbDoorIndexs.includes(e));e<t?e>i&&!r.isPullOutDoor?(r.beamStart=F.clone().addScaledVector(p,u),r.beamEnd=I.clone().addScaledVector(p,u),r.doorStart=L.clone().addScaledVector(p,u),r.doorEnd=R.clone().addScaledVector(p,u),r.doorHeight=n.doorHeight,r.beamStart.z=r.beamEnd.z=r.doorStart.z+r.doorHeight):!N||N&&r.isPullOutDoor&&a?(n.beamStart=r.beamStart.clone().addScaledVector(d,u),n.beamEnd=r.beamEnd.clone().addScaledVector(d,u),n.doorStart=r.doorStart.clone().addScaledVector(d,u),n.doorEnd=r.doorEnd.clone().addScaledVector(d,u),n.doorHeight=r.doorHeight,n.beamStart.z=n.beamEnd.z=n.doorStart.z+n.doorHeight):(r.beamStart=F.clone().addScaledVector(p,u),r.beamEnd=I.clone().addScaledVector(p,u),r.doorStart=L.clone().addScaledVector(p,u),r.doorEnd=R.clone().addScaledVector(p,u),r.doorHeight=n.doorHeight,r.beamStart.z=r.beamEnd.z=r.doorStart.z+r.doorHeight):e>t&&(t>i&&!N?(n.beamStart=r.beamStart.clone().addScaledVector(p,u),n.beamEnd=r.beamEnd.clone().addScaledVector(p,u),n.doorStart=r.doorStart.clone().addScaledVector(p,u),n.doorEnd=r.doorEnd.clone().addScaledVector(p,u),n.doorHeight=r.doorHeight,n.beamStart.z=n.beamEnd.z=n.doorStart.z+n.doorHeight):!r.isPullOutDoor||N&&r.isPullOutDoor&&a?(r.beamStart=F.clone().addScaledVector(d,u),r.beamEnd=I.clone().addScaledVector(d,u),r.doorStart=L.clone().addScaledVector(d,u),r.doorEnd=R.clone().addScaledVector(d,u),r.doorHeight=n.doorHeight,r.beamStart.z=r.beamEnd.z=r.doorStart.z+r.doorHeight):(n.beamStart=r.beamStart.clone().addScaledVector(p,u),n.beamEnd=r.beamEnd.clone().addScaledVector(p,u),n.doorStart=r.doorStart.clone().addScaledVector(p,u),n.doorEnd=r.doorEnd.clone().addScaledVector(p,u),n.doorHeight=r.doorHeight,n.beamStart.z=n.beamEnd.z=n.doorStart.z+n.doorHeight));break}}}if(i.doorAndBeamData.push(n),l.push({linesIndex:r,doorIndex:i.doorAndBeamData.length-1}),s++,N&&T<2)break}}}let E=i.length*(i.rooftopPz-i.start.z);i.completePointAreaPercentage=T/E*100}if(u.length>1)for(let t=0;t<u.length;t++){let n=[];for(let r=0;r<e[u[t]].doorAndBeamData.length;r++){let i=e[u[t]].doorAndBeamData[r].doorStart.clone(),a=e[u[t]].doorAndBeamData[r].doorEnd.clone();i.z=a.z=0;let o=i.distanceTo(a),s=!1;for(let n=0;n<u.length;n++)if(t!=n)for(let t=0;t<e[u[n]].doorAndBeamData.length;t++){let r=e[u[n]].doorAndBeamData[t].doorStart.clone(),l=e[u[n]].doorAndBeamData[t].doorEnd.clone();r.z=l.z=0,r.distanceTo(l);let d=C(i,a,r,l);if(d.type==`collinear_overlap`&&d.maxPerpendicularDistance<.2){let e=new c.Line3(r,l),t=S(i,a),n=e.closestPointToPoint(t,!0,new c.Vector3),u=t.distanceTo(n);if(Math.abs(o-d.overlap)<.1&&u<.05){s=!0;break}}}s||n.push(e[u[t]].doorAndBeamData[r])}n.length!=e[u[t]].doorAndBeamData.length&&(e[u[t]].doorAndBeamData=n)}let p=.65,m=[],g=[],_=.2;for(let t=0;t<o.length;t++){if(o[t].isFind)continue;let r=new c.Vector3(o[t].doorStartPt.x,o[t].doorStartPt.y,0),i=new c.Vector3(o[t].doorEndPt.x,o[t].doorEndPt.y,0),a=new c.Vector3().subVectors(r,i).normalize();for(let l=0;l<e.length;l++){if(e[l].length<1||e[l].completePointAreaPercentage>60)continue;let u=new c.Vector3(e[l].start.x,e[l].start.y,0),d=new c.Vector3(e[l].end.x,e[l].end.y,0),h=new c.Vector3().subVectors(u,d).normalize(),v=.3,y=r.clone().addScaledVector(a,v),S=i.clone().addScaledVector(a.clone().negate(),v),w=u.clone().addScaledVector(h,v),T=d.clone().addScaledVector(h.clone().negate(),v);f(y,S,w,T);let E=b(y,S,w,T),D=r.distanceTo(u)>i.distanceTo(u)?i:r,O=D.distanceTo(u),k=D.distanceTo(d);if((O<_||k<_)&&E&&!E.isParallel){let a=r.distanceTo(u)<r.distanceTo(d)?u:d,f=!1;for(let h=0;h<e.length;h++){if(h==l||m.includes(h)||e[h].length<1)continue;let _=new c.Vector3(e[h].start.x,e[h].start.y,0),v=new c.Vector3(e[h].end.x,e[h].end.y,0),y=x(u,d,_,v);if(y&&y.angle>85){let y=new c.Line3(u,d),x=new c.Vector3,S=y.closestPointToPoint(_,!0,x),w=y.closestPointToPoint(v,!0,x),T=S.distanceTo(_),E=w.distanceTo(v),D,O;if(T>E&&E>p)D=v,O=w;else if(T<E&&T>p)D=_,O=S;else continue;if(T=a.distanceTo(_),E=a.distanceTo(v),T>2&&E>2)continue;let k=b(D,O,_,v);if(!k||!k.isParallel)continue;let A=C(r,i,D,O);if(A&&A.type==`collinear_overlap`&&A.maxPerpendicularDistance<.2)continue;D.z=O.z=e[h].start.z;let j=!1;for(let e=0;e<g.length;e++){let t=C(g[e].start,g[e].end,D,O);if(t&&t.type==`collinear_overlap`&&t.maxPerpendicularDistance<.2){j=!0;break}}if(j)continue;g.push({start:D,end:O}),m.push(h);let M=e[h].rooftopPz-e[h].start.z,N=e[l].rooftopPz-e[l].start.z,P=M<N?M:N;e[h].doorAndBeamData||(e[h].doorAndBeamData=[]),e[h].doorAndBeamData.push({id:s,doorStart:D,doorEnd:O,doorHeight:P,nearId:-1,type:`onlyDoor`,isDoor:!0}),n[o[t].index].isFindBeam=!0,n[o[t].index].isFindOnlyDoor=!0,s++,f=!0;break}}if(f)break}}}let v=0;for(let t=0;t<e.length;t++){let n=e[t];n.mainBeamDataIds||=[],n.isMainBeam||=!1;let{minZ:r,maxZ:i}=O(n);if(r==null||i==null)continue;i-r;let o=n.rooftopPz-n.start.z;if(o<a-.5||o>a+2||Math.abs(r-n.start.z)>.3||n.length<1)continue;let s=new c.Vector3(n.start.x,n.start.y,0),l=new c.Vector3(n.end.x,n.end.y,0);for(let r=0;r<e.length;r++){if(r==t||e[r].length<1)continue;let a=e[r],o=new c.Vector3(a.start.x,a.start.y,0),u=new c.Vector3(a.end.x,a.end.y,0),d=b(s,l,o,u),f=S(s,l),p=new c.Line3(o,u).closestPointToPoint(f,!0,new c.Vector3),m=f.distanceTo(p),h=O(a);if(!(h.minZ===void 0||h.maxZ===void 0)&&!(h.maxZ-h.minZ>1)&&d&&d.isParallel&&m<.25&&Math.abs(i-h.minZ)<.1){let e=`mainBeamId${v}`;a.mainBeamDataIds||=[],a.isMainBeam||=!0,a.mainBeamDataIds.push(e),n.mainBeamDataIds.push(e),v++;break}}}return e},z=(e,t)=>{if(e.length<=0)return e;for(let t of e){let e=[];for(let n=0;n<t.checkResults.length;n++){let r=t.checkResults[n];if(r.isDoor=!1,!r.originalVertices||r.originalVertices.length==0)continue;let i=[];r.originalVertices.forEach(e=>{let t=[];e.forEach(e=>{t.push(new c.Vector3(e.x,e.y,e.z))}),i.push(t)});let a=i.length,o=Array.from({length:a},(e,t)=>t);function s(e){for(;o[e]!==e;)o[e]=o[o[e]],e=o[e];return e}for(let e=0;e<a;e++)for(let t=e+1;t<a;t++){let n=i[e],r=i[t];if(n.some(e=>r.some(t=>e.equals(t)))){let n=s(e),r=s(t);n!==r&&(o[n]=r)}}let l=new Map,u=-1,d=-1;for(let e=0;e<a;e++){let t=s(e);l.has(t)||l.set(t,{points:[],area:0,originalVertices:[],index:-1,gridWidth:-1,gridHeight:-1,centerPts:[]});let r=l.get(t),[a,o,f,p]=i[e],m=new c.Vector3().add(a).add(o).add(f).add(p).multiplyScalar(.25);if(r.points.push(...i[e]),r.area+=a.distanceTo(o)*a.distanceTo(p),r.originalVertices.push(i[e]),r.index=n,r.centerPts.push(m),u==-1||d==-1){let e=new c.Vector3().subVectors(a,o).normalize();d=a.distanceTo(p),u=a.distanceTo(o);let t=h(e,new c.Vector3(0,0,1));t&&t.parallel&&([d,u]=[u,d])}r.gridWidth=u,r.gridHeight=d}e.push(...l.values())}t.mergeCheckRegion=e}let n=new Map;for(let t of e){let e=t.rooftopPz-t.start.z,r=Math.floor(t.rooftopPz-t.start.z);n.has(r)?n.set(r,{num:n.get(r).num+1,totalHeight:n.get(r).totalHeight+e}):n.set(r,{num:1,totalHeight:e})}let r=[...n.entries()].reduce((e,t)=>t[1].num>e[1].num?t:e),i=r[0];r[1].totalHeight/r[1].num,new c.MeshBasicMaterial({color:`#cffd00`,transparent:!0,opacity:.4,side:c.DoubleSide}),new c.PointsMaterial({color:65535,size:.01});for(let n=0;n<e.length;n++){let r=e[n],{minZ:a,maxZ:o}=O(r);if(a==null||o==null)continue;r.doorAndBeamData=[],r.completePointAreaPercentage=-1,o-a;let s=r.rooftopPz-r.start.z;if(s<i-.5||s>i+2||Math.abs(a-r.start.z)>.3||r.length<1)continue;let l=new c.Vector3(r.start.x,r.start.y,0),u=new c.Vector3(r.end.x,r.end.y,0);for(let r=0;r<e.length;r++){if(r==n||e[r].length<1)continue;let i=e[r],a=new c.Vector3(i.start.x,i.start.y,0),s=new c.Vector3(i.end.x,i.end.y,0),d=b(l,u,a,s),f=S(l,u),m=new c.Line3(a,s).closestPointToPoint(f,!0,new c.Vector3),h=f.distanceTo(m),g=O(i);if(!(g.minZ===void 0||g.maxZ===void 0)&&!(g.maxZ-g.minZ>1)&&d&&d.isParallel&&h<.25&&Math.abs(o-g.minZ)<.1){let i=S(l,u),o=S(a,s);i.z=o.z=1,i.applyEuler(new c.Euler(-Math.PI/2,0,0)),o.applyEuler(new c.Euler(-Math.PI/2,0,0));let d=`w:${e[n].originaIndex},b:${e[r].originaIndex}`,f=`b:${e[r].originaIndex},w:${e[n].originaIndex}`;t.add(p(d,i,{textColor:`#fd0000`})),t.add(p(f,o,{textColor:`#fd0000`}))}}}return e};exports.classifySegments=C,exports.getAllGeometry=I,exports.getBeamLine=R,exports.getMainBeamLine=z,exports.getParallelism=b,exports.getPointCloudMinMax=O,exports.getPointCoverageOnQuad=E,exports.isParallel=h,exports.perpendicularInfo=x,exports.removeNoisePoints=A,exports.segmentsIntersect2D=f,exports.usegetBeamLine=L;
|
|
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}),g=new c.LineBasicMaterial({color:`#cffd00`}),_=[];for(let e=0;e<a.length;e++){let t=a[e];if(t.checkResults)for(let e of t.checkResults)for(let t of e.originalVertices){let[e,n,r,i]=t,a=new c.BufferGeometry,o=new Float32Array([e.x,e.y,e.z,n.x,n.y,n.z,r.x,r.y,r.z,i.x,i.y,i.z]);a.setAttribute(`position`,new c.BufferAttribute(o,3)),a.setIndex([0,1,2,0,2,3]),a.computeVertexNormals();let s=new c.EdgesGeometry(a,1);_.push(s)}}if(_.length>0){let e=l(_,!1);e.rotateX(-Math.PI/2),r.add(new c.LineSegments(e,g)),_.forEach(e=>e.dispose())}let y=new c.PointsMaterial({color:65535,size:.01}),b=[],x=[];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)),b.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),x.push(i)}if(b.length>0){let e=l(b,!1);e.rotateX(-Math.PI/2),r.add(new c.Points(e,y)),b.forEach(e=>e.dispose())}let S=[],C=[],w=new c.LineBasicMaterial({color:`#f30606`}),T=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),C.push(s),S.push(t)}if(S.length>0){let e=l(S,!1);e.rotateX(-Math.PI/2),r.add(new c.Mesh(e,T)),S.forEach(e=>e.dispose())}if(C.length>0){let e=l(C,!1);e.rotateX(-Math.PI/2),r.add(new c.Line(e,w)),C.forEach(e=>e.dispose())}let E=v(f);E.rotateX(-Math.PI/2),r.add(E),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/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)}}}console.log(`isDoor`,M),(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<1||e[l].completePointAreaPercentage>60)continue;let u=new c.Vector3(e[l].start.x,e[l].start.y,0),d=new c.Vector3(e[l].end.x,e[l].end.y,0),h=new c.Vector3().subVectors(u,d).normalize(),v=.3,y=r.clone().addScaledVector(a,v),S=i.clone().addScaledVector(a.clone().negate(),v),w=u.clone().addScaledVector(h,v),T=d.clone().addScaledVector(h.clone().negate(),v);f(y,S,w,T);let E=b(y,S,w,T),D=r.distanceTo(u)>i.distanceTo(u)?i:r,O=D.distanceTo(u),k=D.distanceTo(d);if((O<_||k<_)&&E&&!E.isParallel){let a=r.distanceTo(u)<r.distanceTo(d)?u:d,f=!1;for(let h=0;h<e.length;h++){if(h==l||m.includes(h)||e[h].length<1)continue;let _=new c.Vector3(e[h].start.x,e[h].start.y,0),v=new c.Vector3(e[h].end.x,e[h].end.y,0),y=x(u,d,_,v);if(y&&y.angle>85){let y=new c.Line3(u,d),x=new c.Vector3,S=y.closestPointToPoint(_,!0,x),w=y.closestPointToPoint(v,!0,x),T=S.distanceTo(_),E=w.distanceTo(v),D,O;if(T>E&&E>p)D=v,O=w;else if(T<E&&T>p)D=_,O=S;else continue;if(T=a.distanceTo(_),E=a.distanceTo(v),T>2&&E>2)continue;let k=b(D,O,_,v);if(!k||!k.isParallel)continue;let A=C(r,i,D,O);if(A&&A.type==`collinear_overlap`&&A.maxPerpendicularDistance<.2)continue;D.z=O.z=e[h].start.z;let j=!1;for(let e=0;e<g.length;e++){let t=C(g[e].start,g[e].end,D,O);if(t&&t.type==`collinear_overlap`&&t.maxPerpendicularDistance<.2){j=!0;break}}if(j)continue;g.push({start:D,end:O}),m.push(h);let M=e[h].rooftopPz-e[h].start.z,N=e[l].rooftopPz-e[l].start.z,P=M<N?M:N;e[h].doorAndBeamData||(e[h].doorAndBeamData=[]),e[h].doorAndBeamData.push({id:s,beamStart:D,beamEnd:D,beamHeight:P,doorStart:D,doorEnd:O,doorHeight:P,nearId:-1,type:`onlyDoor`,isDoor:!0}),n[o[t].index].isFindBeam=!0,n[o[t].index].isFindOnlyDoor=!0,s++,f=!0;break}}if(f)break}}}let v=0;for(let t=0;t<e.length;t++){let n=e[t];n.mainBeamDataIds||=[],n.isMainBeam||=!1;let{minZ:r,maxZ:i}=O(n);if(r==null||i==null)continue;i-r;let o=n.rooftopPz-n.start.z;if(o<a-.5||o>a+2||Math.abs(r-n.start.z)>.3||n.length<1)continue;let s=new c.Vector3(n.start.x,n.start.y,0),l=new c.Vector3(n.end.x,n.end.y,0);for(let r=0;r<e.length;r++){if(r==t||e[r].length<1)continue;let a=e[r],o=new c.Vector3(a.start.x,a.start.y,0),u=new c.Vector3(a.end.x,a.end.y,0),d=b(s,l,o,u),f=S(s,l),p=new c.Line3(o,u).closestPointToPoint(f,!0,new c.Vector3),m=f.distanceTo(p),h=O(a);if(!(h.minZ===void 0||h.maxZ===void 0)&&!(h.maxZ-h.minZ>1)&&d&&d.isParallel&&m<.25&&Math.abs(i-h.minZ)<.1){let e=`mainBeamId${v}`;a.mainBeamDataIds||=[],a.isMainBeam||=!0,a.mainBeamDataIds.push(e),n.mainBeamDataIds.push(e),v++;break}}}return e},z=(e,t)=>{if(e.length<=0)return e;for(let t of e){let e=[];for(let n=0;n<t.checkResults.length;n++){let r=t.checkResults[n];if(r.isDoor=!1,!r.originalVertices||r.originalVertices.length==0)continue;let i=[];r.originalVertices.forEach(e=>{let t=[];e.forEach(e=>{t.push(new c.Vector3(e.x,e.y,e.z))}),i.push(t)});let a=i.length,o=Array.from({length:a},(e,t)=>t);function s(e){for(;o[e]!==e;)o[e]=o[o[e]],e=o[e];return e}for(let e=0;e<a;e++)for(let t=e+1;t<a;t++){let n=i[e],r=i[t];if(n.some(e=>r.some(t=>e.equals(t)))){let n=s(e),r=s(t);n!==r&&(o[n]=r)}}let l=new Map,u=-1,d=-1;for(let e=0;e<a;e++){let t=s(e);l.has(t)||l.set(t,{points:[],area:0,originalVertices:[],index:-1,gridWidth:-1,gridHeight:-1,centerPts:[]});let r=l.get(t),[a,o,f,p]=i[e],m=new c.Vector3().add(a).add(o).add(f).add(p).multiplyScalar(.25);if(r.points.push(...i[e]),r.area+=a.distanceTo(o)*a.distanceTo(p),r.originalVertices.push(i[e]),r.index=n,r.centerPts.push(m),u==-1||d==-1){let e=new c.Vector3().subVectors(a,o).normalize();d=a.distanceTo(p),u=a.distanceTo(o);let t=h(e,new c.Vector3(0,0,1));t&&t.parallel&&([d,u]=[u,d])}r.gridWidth=u,r.gridHeight=d}e.push(...l.values())}t.mergeCheckRegion=e}let n=new Map;for(let t of e){let e=t.rooftopPz-t.start.z,r=Math.floor(t.rooftopPz-t.start.z);n.has(r)?n.set(r,{num:n.get(r).num+1,totalHeight:n.get(r).totalHeight+e}):n.set(r,{num:1,totalHeight:e})}let r=[...n.entries()].reduce((e,t)=>t[1].num>e[1].num?t:e),i=r[0];r[1].totalHeight/r[1].num,new c.MeshBasicMaterial({color:`#cffd00`,transparent:!0,opacity:.4,side:c.DoubleSide}),new c.PointsMaterial({color:65535,size:.01});for(let n=0;n<e.length;n++){let r=e[n],{minZ:a,maxZ:o}=O(r);if(a==null||o==null)continue;r.doorAndBeamData=[],r.completePointAreaPercentage=-1,o-a;let s=r.rooftopPz-r.start.z;if(s<i-.5||s>i+2||Math.abs(a-r.start.z)>.3||r.length<1)continue;let l=new c.Vector3(r.start.x,r.start.y,0),u=new c.Vector3(r.end.x,r.end.y,0);for(let r=0;r<e.length;r++){if(r==n||e[r].length<1)continue;let i=e[r],a=new c.Vector3(i.start.x,i.start.y,0),s=new c.Vector3(i.end.x,i.end.y,0),d=b(l,u,a,s),f=S(l,u),m=new c.Line3(a,s).closestPointToPoint(f,!0,new c.Vector3),h=f.distanceTo(m),g=O(i);if(!(g.minZ===void 0||g.maxZ===void 0)&&!(g.maxZ-g.minZ>1)&&d&&d.isParallel&&h<.25&&Math.abs(o-g.minZ)<.1){let i=S(l,u),o=S(a,s);i.z=o.z=1,i.applyEuler(new c.Euler(-Math.PI/2,0,0)),o.applyEuler(new c.Euler(-Math.PI/2,0,0));let d=`w:${e[n].originaIndex},b:${e[r].originaIndex}`,f=`b:${e[r].originaIndex},w:${e[n].originaIndex}`;t.add(p(d,i,{textColor:`#fd0000`})),t.add(p(f,o,{textColor:`#fd0000`}))}}}return e};exports.classifySegments=C,exports.getAllGeometry=I,exports.getBeamLine=R,exports.getMainBeamLine=z,exports.getParallelism=b,exports.getPointCloudMinMax=O,exports.getPointCoverageOnQuad=E,exports.isParallel=h,exports.perpendicularInfo=x,exports.removeNoisePoints=A,exports.segmentsIntersect2D=f,exports.usegetBeamLine=L;
|
package/dist/index.js
CHANGED
|
@@ -1058,12 +1058,44 @@ var T = (t, n, r, i = .05, a = .02) => {
|
|
|
1058
1058
|
transparent: !0,
|
|
1059
1059
|
opacity: .4,
|
|
1060
1060
|
side: e.DoubleSide
|
|
1061
|
-
});
|
|
1062
|
-
|
|
1063
|
-
|
|
1061
|
+
}), g = new e.LineBasicMaterial({ color: "#cffd00" }), v = [];
|
|
1062
|
+
for (let t = 0; t < o.length; t++) {
|
|
1063
|
+
let n = o[t];
|
|
1064
|
+
if (n.checkResults) for (let t of n.checkResults) for (let n of t.originalVertices) {
|
|
1065
|
+
let [t, r, i, a] = n, o = new e.BufferGeometry(), s = new Float32Array([
|
|
1066
|
+
t.x,
|
|
1067
|
+
t.y,
|
|
1068
|
+
t.z,
|
|
1069
|
+
r.x,
|
|
1070
|
+
r.y,
|
|
1071
|
+
r.z,
|
|
1072
|
+
i.x,
|
|
1073
|
+
i.y,
|
|
1074
|
+
i.z,
|
|
1075
|
+
a.x,
|
|
1076
|
+
a.y,
|
|
1077
|
+
a.z
|
|
1078
|
+
]);
|
|
1079
|
+
o.setAttribute("position", new e.BufferAttribute(s, 3)), o.setIndex([
|
|
1080
|
+
0,
|
|
1081
|
+
1,
|
|
1082
|
+
2,
|
|
1083
|
+
0,
|
|
1084
|
+
2,
|
|
1085
|
+
3
|
|
1086
|
+
]), o.computeVertexNormals();
|
|
1087
|
+
let c = new e.EdgesGeometry(o, 1);
|
|
1088
|
+
v.push(c);
|
|
1089
|
+
}
|
|
1090
|
+
}
|
|
1091
|
+
if (v.length > 0) {
|
|
1092
|
+
let t = c(v, !1);
|
|
1093
|
+
t.rotateX(-Math.PI / 2), i.add(new e.LineSegments(t, g)), v.forEach((e) => e.dispose());
|
|
1094
|
+
}
|
|
1095
|
+
let y = new e.PointsMaterial({
|
|
1064
1096
|
color: 65535,
|
|
1065
1097
|
size: .01
|
|
1066
|
-
}),
|
|
1098
|
+
}), b = [], x = [];
|
|
1067
1099
|
for (let t of o) {
|
|
1068
1100
|
if (!t.originalPoints || t.originalPoints.length === 0) continue;
|
|
1069
1101
|
let n = new Float32Array(t.originalPoints.length * 3);
|
|
@@ -1071,35 +1103,35 @@ var T = (t, n, r, i = .05, a = .02) => {
|
|
|
1071
1103
|
n[t * 3] = e.x, n[t * 3 + 1] = e.y, n[t * 3 + 2] = e.z;
|
|
1072
1104
|
});
|
|
1073
1105
|
let r = new e.BufferGeometry();
|
|
1074
|
-
r.setAttribute("position", new e.BufferAttribute(n, 3)),
|
|
1106
|
+
r.setAttribute("position", new e.BufferAttribute(n, 3)), b.push(r);
|
|
1075
1107
|
let i = [];
|
|
1076
1108
|
i.push(new e.Vector3(t.start.x, t.start.y, t.start.z)), i.push(new e.Vector3(t.end.x, t.end.y, t.end.z));
|
|
1077
1109
|
let a = new e.BufferGeometry();
|
|
1078
|
-
a.setFromPoints(i),
|
|
1110
|
+
a.setFromPoints(i), x.push(a);
|
|
1079
1111
|
}
|
|
1080
|
-
if (
|
|
1081
|
-
let t = c(
|
|
1082
|
-
t.rotateX(-Math.PI / 2), i.add(new e.Points(t,
|
|
1112
|
+
if (b.length > 0) {
|
|
1113
|
+
let t = c(b, !1);
|
|
1114
|
+
t.rotateX(-Math.PI / 2), i.add(new e.Points(t, y)), b.forEach((e) => e.dispose());
|
|
1083
1115
|
}
|
|
1084
|
-
let
|
|
1116
|
+
let S = [], C = [], w = new e.LineBasicMaterial({ color: "#f30606" }), T = new e.MeshBasicMaterial({ color: "#071ac4" });
|
|
1085
1117
|
for (let t = 0; t < u.length; t++) {
|
|
1086
1118
|
let n = new e.SphereGeometry(.05);
|
|
1087
1119
|
n.translate(u[t].x, u[t].y, u[t].z);
|
|
1088
1120
|
let r = [], i = new e.Vector3(u[t].x, u[t].y, u[t].z), a = new e.Quaternion(u[t].qx, u[t].qy, u[t].qz, u[t].qw), o = new e.Vector3(0, 0, -1).applyQuaternion(a), s = i.clone().addScaledVector(o, .001);
|
|
1089
1121
|
r.push(i), r.push(s);
|
|
1090
1122
|
let c = new e.BufferGeometry();
|
|
1091
|
-
c.setFromPoints(r),
|
|
1123
|
+
c.setFromPoints(r), C.push(c), S.push(n);
|
|
1092
1124
|
}
|
|
1093
|
-
if (
|
|
1094
|
-
let t = c(
|
|
1095
|
-
t.rotateX(-Math.PI / 2), i.add(new e.Mesh(t,
|
|
1125
|
+
if (S.length > 0) {
|
|
1126
|
+
let t = c(S, !1);
|
|
1127
|
+
t.rotateX(-Math.PI / 2), i.add(new e.Mesh(t, T)), S.forEach((e) => e.dispose());
|
|
1096
1128
|
}
|
|
1097
|
-
if (
|
|
1098
|
-
let t = c(
|
|
1099
|
-
t.rotateX(-Math.PI / 2), i.add(new e.Line(t,
|
|
1129
|
+
if (C.length > 0) {
|
|
1130
|
+
let t = c(C, !1);
|
|
1131
|
+
t.rotateX(-Math.PI / 2), i.add(new e.Line(t, w)), C.forEach((e) => e.dispose());
|
|
1100
1132
|
}
|
|
1101
|
-
let
|
|
1102
|
-
|
|
1133
|
+
let E = _(p);
|
|
1134
|
+
E.rotateX(-Math.PI / 2), i.add(E), console.log("jsonpcdData", p), o.forEach((t) => {
|
|
1103
1135
|
if (t.doorAndBeamData && t.doorAndBeamData.length > 0) {
|
|
1104
1136
|
let n = new e.Group();
|
|
1105
1137
|
t.doorAndBeamData.forEach((t) => {
|
|
@@ -1321,7 +1353,7 @@ var T = (t, n, r, i = .05, a = .02) => {
|
|
|
1321
1353
|
let V = 0;
|
|
1322
1354
|
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
1355
|
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) continue;
|
|
1356
|
+
if (g > a.start.z - .15 && a.length < 1.5 || a.length, E < .5 || h < v.facePoints[2].z || k > w / 3 * 2) continue;
|
|
1325
1357
|
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
1358
|
if (Math.abs(v.facePoints[0].z - a.start.z) < .25) {
|
|
1327
1359
|
if (T > .2 && a.length > .5) {
|
|
@@ -1340,10 +1372,10 @@ var T = (t, n, r, i = .05, a = .02) => {
|
|
|
1340
1372
|
if (p && p.angle > 85) {
|
|
1341
1373
|
let t = y(r, o, new e.Line3(a, c).closestPointToPoint(r, !0, l), r);
|
|
1342
1374
|
if (!t || t.direction != "same") continue;
|
|
1343
|
-
let i =
|
|
1344
|
-
if (
|
|
1345
|
-
let
|
|
1346
|
-
|
|
1375
|
+
let i = r.clone().addScaledVector(s, 10), p = o.clone().addScaledVector(s.clone().negate(), 10), m = new e.Vector3().subVectors(c, a).normalize(), h = d(i, p, a.clone().addScaledVector(m.clone().negate(), .15), c.clone().addScaledVector(m, .15));
|
|
1376
|
+
if (h == null) continue;
|
|
1377
|
+
let g = r.distanceTo(h.point);
|
|
1378
|
+
g < u && (u = g, f = n);
|
|
1347
1379
|
} else if (p && p.angle < 5) {
|
|
1348
1380
|
let e = A(r, a), t = A(r, c), n = A(o, a), i = A(o, c);
|
|
1349
1381
|
if (e < .1 || t < .1 || n < .1 || i < .1) {
|
|
@@ -1368,10 +1400,10 @@ var T = (t, n, r, i = .05, a = .02) => {
|
|
|
1368
1400
|
if (p && p.angle > 85) {
|
|
1369
1401
|
let t = y(r, n, new e.Line3(o, c).closestPointToPoint(r, !0, l), r);
|
|
1370
1402
|
if (!t || t.direction != "same") continue;
|
|
1371
|
-
let i =
|
|
1372
|
-
if (
|
|
1373
|
-
let
|
|
1374
|
-
|
|
1403
|
+
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));
|
|
1404
|
+
if (h == null) continue;
|
|
1405
|
+
let g = r.distanceTo(h.point);
|
|
1406
|
+
g < u && (u = g, f = a);
|
|
1375
1407
|
}
|
|
1376
1408
|
if (p && p.angle < 5) {
|
|
1377
1409
|
let e = A(n, o), t = A(n, c), i = A(r, o), a = A(r, c);
|
|
@@ -1406,8 +1438,8 @@ var T = (t, n, r, i = .05, a = .02) => {
|
|
|
1406
1438
|
}
|
|
1407
1439
|
}
|
|
1408
1440
|
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 (
|
|
1441
|
+
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;
|
|
1442
|
+
if (Math.abs(v.facePoints[2].z - s[t].maxZ), l) {
|
|
1411
1443
|
let l = S(n.clone(), i.clone(), o.clone(), c.clone(), { distanceThreshold: 2 });
|
|
1412
1444
|
if (l && l.type == "collinear_gap") {
|
|
1413
1445
|
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);
|
|
@@ -1427,7 +1459,7 @@ var T = (t, n, r, i = .05, a = .02) => {
|
|
|
1427
1459
|
}
|
|
1428
1460
|
}
|
|
1429
1461
|
}
|
|
1430
|
-
(M || t || E > w / 2 && T > .8 && Math.abs(g - a.start.z) < .4) && (j = !0);
|
|
1462
|
+
console.log("isDoor", M), (M || t || E > w / 2 && T > .8 && Math.abs(g - a.start.z) < .4) && (j = !0);
|
|
1431
1463
|
}
|
|
1432
1464
|
}
|
|
1433
1465
|
if (j) {
|
|
@@ -1440,10 +1472,11 @@ var T = (t, n, r, i = .05, a = .02) => {
|
|
|
1440
1472
|
}
|
|
1441
1473
|
}
|
|
1442
1474
|
if (n) continue;
|
|
1443
|
-
|
|
1475
|
+
f.push({
|
|
1444
1476
|
start: L,
|
|
1445
1477
|
end: R
|
|
1446
|
-
})
|
|
1478
|
+
});
|
|
1479
|
+
{
|
|
1447
1480
|
let t = new e.Vector3().subVectors(L, R), n = new e.Vector3().subVectors(R, L), r = Infinity, i = Infinity, o = L.clone(), s = R.clone(), c = [];
|
|
1448
1481
|
if (a.originalPoints.forEach((e) => {
|
|
1449
1482
|
e.z - L.z < E / 2 + .1 && e.z - L.z > E / 2 - .1 && c.push(e);
|
|
@@ -1571,6 +1604,9 @@ var T = (t, n, r, i = .05, a = .02) => {
|
|
|
1571
1604
|
let M = t[m].rooftopPz - t[m].start.z, N = t[l].rooftopPz - t[l].start.z, P = M < N ? M : N;
|
|
1572
1605
|
t[m].doorAndBeamData || (t[m].doorAndBeamData = []), t[m].doorAndBeamData.push({
|
|
1573
1606
|
id: c,
|
|
1607
|
+
beamStart: D,
|
|
1608
|
+
beamEnd: D,
|
|
1609
|
+
beamHeight: P,
|
|
1574
1610
|
doorStart: D,
|
|
1575
1611
|
doorEnd: O,
|
|
1576
1612
|
doorHeight: P,
|