@soonspacejs/plugin-heat-cloud 2.13.5 → 2.13.7

package/dist/index.esm.js

@@ -1 +1,3360 @@
-import{TextureLoader as t,Texture as e,Matrix3 as n,Vector3 as i,Vector4 as a,Vector2 as o,AlphaFormat as r,RedFormat as s,RedIntegerFormat as l,RGFormat as c,RGIntegerFormat as u,RGBAFormat as d,RGBAIntegerFormat as m,LuminanceFormat as p,LuminanceAlphaFormat as f,DepthFormat as h,DepthStencilFormat as v,ShaderMaterial as g,GLSL3 as x,Matrix4 as _,DoubleSide as y,Mesh as w,BoxGeometry as M,Data3DTexture as z,FloatType as D,LinearFilter as S,UnsignedByteType as b,Color as C,FrontSide as N,UniformsUtils as U,ShaderLib as A,Points as P,PlaneGeometry as O,Box3 as I,Line3 as T,Quaternion as F,BackSide as R,RawShaderMaterial as V,BufferGeometry as E,Float32BufferAttribute as G,Uint8ClampedBufferAttribute as j}from"three";import L from"soonspacejs";var B,k=(t=>(t[t.x=0]="x",t[t.y=1]="y",t[t.z=2]="z",t))(k||{});(B=k||(k={})).toKey=function(t){return B[t]},B.toIndex=function(t){return B[t]},B.getCrossAxiss=function(t){return[(t+1)%3,(t+2)%3]};var H=(t=>(t[t.x=0]="x",t[t.y=1]="y",t[t.z=2]="z",t[t.w=3]="w",t))(H||{});(t=>{t.toKey=function(e){return t[e]},t.toIndex=function(e){return t[e]},t.getCrossAxiss=function(t){return[(t+1)%4,(t+2)%4]}})(H||(H={}));var q=Object.defineProperty,K=(t,e,n)=>(((t,e,n)=>{e in t?q(t,e,{enumerable:!0,configurable:!0,writable:!0,value:n}):t[e]=n})(t,"symbol"!=typeof e?e+"":e,n),n);function W(t,e,n,i="srgb"){const a=new OffscreenCanvas(e,n).getContext("2d"),o=a.createLinearGradient(0,0,256,1);for(const[e,n]of t)o.addColorStop(Number(e),n);return a.fillStyle=o,a.fillRect(0,0,e,n),a.getImageData(0,0,e,n,{colorSpace:i})}function Y(t,e){let{area:n,width:i,height:a,colorSpace:o}=e||{};o=o??"srgb",n=n||{};const r=n.x??0,s=n.y??0,l=X(t),c=n.width??l.width,u=n.height??l.height;i=i??c,a=a??u;const d=new OffscreenCanvas(i,a).getContext("2d");return d.drawImage(t,r,s,c,u,0,0,i,a),d.getImageData(0,0,i,a,{colorSpace:o})}function X(t){let e=0,n=0;return t instanceof HTMLImageElement?(e=t.naturalWidth,n=t.naturalHeight):t instanceof HTMLVideoElement?(e=t.videoWidth,n=t.videoHeight):(e=t.width,n=t.height),{width:e,height:n}}function Z(t){const{data:e,width:n,height:i,colorSpace:a}=t,o=e.length,r=new Uint8ClampedArray(o);for(let t=0;t<i;t++){const i=t*n,a=-(i+n),o=a+n,s=e.slice(4*a,4*o);r.set(s,4*i)}return new ImageData(r,n,i,{colorSpace:a})}function J(t,e){const{data:n,width:i}=t,a=e*i,o=a+i;return n.slice(4*a,4*o)}function Q(t,e,n){const{data:i,width:a}=t,o=4*(e*a+n);return i.slice(o,o+4)}function $(t){const e=document.createElement("canvas"),n=e.getContext("2d");return e.width=t.width,e.height=t.height,n.putImageData(t,0,0),e}function tt(t){return $(t).toDataURL()}function et(t){return t&&null!=t.data&&null!=t.width&&null!=t.height}function nt(t){return et(t)&&null!=t.depth}function it(t,e){for(const[n,i]of Object.entries(t))if(0<=i&&i<=e[n])return!0;return!1}function at(t,e){const{x:n,y:i}=e;return t.z*i*n+t.y*n+t.x}function ot(t,e){const{x:n,y:i}=e,a=Math.trunc(t/n),o=t-a*n,r=Math.trunc(a/i);return{x:o,y:a-r*i,z:r}}function rt(t,e,n=1){let{x:i,y:a,z:o}=e;const r=t.size;return i=Math.max(0,Math.min(r.x-1,i)),a=Math.max(0,Math.min(r.y-1,a)),o=Math.max(0,Math.min(r.z-1,o)),st(t,{x:i,y:a,z:o},n)}function st(t,e,n=1){const{data:i,size:a}=t,o=at(e,a);let r=[];if(n>0){const t=o*n;for(let e=0;e<n;e++)r.push(i[t+e])}return{index:o,value:r}}function lt(t,e){let{x:n,y:i,z:a}=e;const o=t.size;return n=Math.max(0,Math.min(o.x-1,n)),i=Math.max(0,Math.min(o.y-1,i)),a=Math.max(0,Math.min(o.z-1,a)),ct(t,{x:n,y:i,z:a})}function ct(t,e){const{data:n,size:i}=t,a=at(e,i);return{index:a,value:n[a]}}function ut(t,e){const{voidValue:n,valuesAccumulate:i,verifyVoid:a=(()=>!1)}=e;let{x:o,y:r,z:s}=t[0].size;const l=t.length;for(let e=1;e<l;e++){const n=t[e].size;o=Math.max(o,n.x),r=Math.max(r,n.y),s=Math.max(s,n.z)}const c=new Array(o*r*s),u=o*r;for(let e=0;e<s;e++){const s=e*u;for(let l=0;l<r;l++){const r=s+l*o;for(let s=0;s<o;s++){const o=r+s,u={x:s,y:l,z:e},d=[];for(const n of t){const t=n.size;if(s<t.x&&l<t.y&&e<t.z)continue;const{index:i,value:o}=ct(n,u);a(o,n,i)||d.push({value:o,index:i,data3D:n,coord:u})}const m=d.length;let p=n;m>1?p=i(d):1===m&&(p=d[0].value),c[o]=p}}}return{data:c,size:{x:o,y:r,z:s}}}function dt(t,e,n,i=1){const{size:a}=t,o=k.toKey(e),[r,s]=k.getCrossAxiss(e),l=k.toKey(r),c=k.toKey(s),u=a[l],d=a[c],m=[],p=a[o];n=Math.trunc(n),n=Math.max(0,Math.min(p-1,n));let f={x:0,y:0,z:0};f[o]=n;for(let e=0;e<d;e++){f[c]=e;for(let e=0;e<u;e++){f[l]=e;const{value:n}=st(t,f,i);m.push(...n)}}return{data:m,size:{x:u,y:d}}}function mt(t,e,n){e=Math.trunc(e);const{data:i,width:a,height:o,colorSpace:r}=n?Z(t):t,s=i.length,l=new Uint8ClampedArray(s*e);for(let t=0;t<e;t++)l.set(i,t*s);return{data:l,width:a,height:o,depth:e,colorSpace:r}}function pt(t,e,n){const i=n||{},a=i.reverseY,o=i.axis??"z",r=et(t)?t:Y(t,i);return"z"===o?mt(r,e,a):ft(r,o,e,a)}function ft(t,e,n,i){n=Math.trunc(n);const{data:a,width:o,height:r,colorSpace:s}=i?Z(t):t,l=k.toKey(e),c=(e+1)%3,u=(e+2)%3,d=k.toKey(c),m=k.toKey(u),p={[d]:o,[m]:r,[l]:n},f=new Uint8ClampedArray(a.length*n),h={x:0,y:0,z:0};for(let t=0;t<n;t++){h[l]=t;for(let t=0;t<r;t++){h[m]=t;for(let e=0;e<o;e++){h[d]=e;const n=4*(t*o+e),i=n+4,r=a.subarray(n,i),s=4*at(h,p);f.set(r,s)}}}return{data:f,width:p.x,height:p.y,depth:p.z,colorSpace:s}}function ht(t,e){return t.y*e.x+t.x}function vt(t,e){const{x:n}=e,i=Math.trunc(t/n);return{x:t-i*n,y:i}}function gt(t,e,n=1){let{x:i,y:a}=e;const o=t.size;return i=Math.max(0,Math.min(o.x-1,i)),a=Math.max(0,Math.min(o.y-1,a)),xt(t,{x:i,y:a},n)}function xt(t,e,n=1){const{data:i,size:a}=t,o=ht(e,a);let r=[];if(n>0){const t=o*n;for(let e=0;e<n;e++)r.push(i[t+e])}return{index:o,value:r}}function _t(t,e){let{x:n,y:i}=e;const a=t.size;return n=Math.max(0,Math.min(a.x-1,n)),i=Math.max(0,Math.min(a.y-1,i)),yt(t,{x:n,y:i})}function yt(t,e){const{data:n,size:i}=t,a=ht(e,i);return{index:a,value:n[a]}}function wt(t,e,n=1){const{data:i,size:a}=t,o=a.x,r=e*o,s=r+o;return("function"==typeof i.slice?i:Array.from(i)).slice(r*n,s*n)}class Mt{constructor(t){K(this,"isImageData3DSlice",!0),K(this,"_image3D"),K(this,"_canvas"),K(this,"_context"),K(this,"_axis",k.z),K(this,"_depth",0),K(this,"_sliceSize",null),K(this,"_slice"),this.image3D=t}get image3D(){return this._image3D}set image3D(t){this._image3D=t,this.updateSize()}get canvas(){let t=this._canvas;return t||(t=this._canvas=document.createElement("canvas"),t.width=this.sliceSize.x,t.height=this.sliceSize.y),t}get context(){return this._context??(this._context=this.canvas.getContext("2d"))}get axis(){return this._axis}set axis(t){this._axis=t,this.updateSize()}get depth(){return this._depth}set depth(t){this._depth=t,this.updateSlice()}get sliceSize(){let t=this._sliceSize;if(!t){const{axis:e,image3D:n}=this,{width:i,height:a,depth:o}=n,[r,s]=k.getCrossAxiss(e),l=[i,a,o];this._sliceSize=t={x:l[r],y:l[s]}}return t}get slice(){return this._slice??(this._slice=this.getSlice(this.axis,this.depth))}getSlice(t,e){const{data:n,width:i,height:a,depth:o,colorSpace:r}=this.image3D,{data:s,size:l}=dt({data:n,size:{x:i,y:a,z:o}},t,e,4),c=s instanceof Uint8ClampedArray?s:Uint8ClampedArray.from(s);return new ImageData(c,l.x,l.y,{colorSpace:r})}getDataURL(){return this.canvas.toDataURL()}getColor(t){const{x:e,y:n}=t;return Q(this.slice,n,e)}updateSize(){this._sliceSize=null;const t=this.canvas;t.width=this.sliceSize.x,t.height=this.sliceSize.y,this.updateSlice()}updateSlice(){this._slice=null,this.context.putImageData(this.slice,0,0)}}var zt=Object.defineProperty,Dt=(t,e,n)=>(((t,e,n)=>{e in t?zt(t,e,{enumerable:!0,configurable:!0,writable:!0,value:n}):t[e]=n})(t,"symbol"!=typeof e?e+"":e,n),n);function St(t){const e=t.length,n=[];if(0===e)return n;const[r,s]=function(t){return void 0!==t.w?[a,4]:void 0!==t.z?[i,3]:[o,2]}(t[0]);for(let i=0;i<e;i++){const a=t[i],o=(new r).copy(a),s=a.radius,l=n[i]??(n[i]=[]);for(let a=i+1;a<e;a++){const e=t[a];o.distanceToSquared(e)<=(e.radius+s)**2&&(l.push(a),(n[a]=[]).push(i))}}return n}let bt;function Ct(n){if("string"==typeof n)return(bt||(bt=new t)).load(n);if(n.isTexture)return n;const i=Array.isArray(n)?W(n,256,1):n,a=new e(i);return"srgb"===i.colorSpace&&(a.colorSpace="srgb"),a.needsUpdate=!0,a}function Nt(t,e,n){let i=function(t,e){const n=Math.sqrt(t.lengthSq()*e.lengthSq());if(0===n)return 0;let i=t.dot(e)/n;return i=Math.max(-1,Math.min(1,i)),Math.acos(i)}(t,e);return 0===i?i:t.clone().cross(e).dot(n)<0?-i:i}const Ut=180/Math.PI,At={yaw:[{name:"前",range:[-15,15]},{name:"左",range:[15,165]},{name:"右",range:[-165,-15]},{name:"后",range:[-180.1,-165]},{name:"后",range:[165,180.1]}],pitch:[{name:"前",range:[-15,15]},{name:"上",range:[15,165]},{name:"下",range:[-165,-15]},{name:"后",range:[-180.1,-165]},{name:"后",range:[165,180.1]}],roll:[{name:"上",range:[-15,15]},{name:"左",range:[15,165]},{name:"右",range:[-165,-15]},{name:"下",range:[-180.1,-165]},{name:"下",range:[165,180.1]}]};const Pt={degrees:!0,map:At,front:{x:0,y:0,z:1},up:{x:0,y:1,z:0}};let Ot=class{constructor(t){Dt(this,"_options"),Dt(this,"_listMap"),Dt(this,"_front"),Dt(this,"_up"),t&&(this.options=t)}static get options(){return this._options??(this.options=Pt)}static set options(t){this._options=Object.assign({},structuredClone(At),t)}get defaultOptions(){return this.constructor.options}get options(){return this._options??(this.options=this.defaultOptions)}set options(t){this._options=Object.assign({},structuredClone(this.defaultOptions),t),this._listMap=null,this._front=null,this._up=null}get degrees(){return this.options.degrees??(this.options.degrees=this.defaultOptions.degrees??!0)}set degrees(t){this.options.degrees=t}get map(){return this.options.map||(this.map=this.defaultOptions.map),this.options.map}set map(t){const e=structuredClone(this.defaultOptions.map),n=structuredClone(At);t=t?{yaw:t.yaw??e.yaw??n.yaw,pitch:t.pitch??e.pitch??n.pitch,roll:t.roll??e.roll??n.roll}:e,this.options.map=t,this._listMap=null}get listMap(){return this._listMap??(this._listMap=function(t){const e={};for(const[n,i]of Object.entries(t))e[n]=Array.isArray(i)?i:Object.entries(i).map((([t,e])=>({name:t,range:e})));return e}(this.map??{}))}get front(){return this._front||(this.front=this.options.front??this.defaultOptions.front??Pt.front),this._front}set front(t){this._front=(new i).copy(t)}get up(){return this._up||(this.up=this.options.up??this.defaultOptions.up??Pt.up),this._up}set up(t){this._up=(new i).copy(t)}computeAzimuth(t,e,n){const a=(new i).copy(t),o=(new i).copy(e??this.front),r=(new i).copy(n??this.up);let{yaw:s,pitch:l,roll:c}=function(t,e,n){const i=e.clone().negate(),a=n.clone().projectOnPlane(e),o=Nt(e,t.clone().projectOnPlane(a),a),r=e.clone().cross(a);return{yaw:o,pitch:Nt(e,t.clone().projectOnPlane(r),r),roll:Nt(a,t.clone().projectOnPlane(i),i)}}(a,o,r);return this.degrees&&(s*=Ut,l*=Ut,c*=Ut),{yaw:{angle:s,name:this.findAzimuthNames("yaw",s)},pitch:{angle:l,name:this.findAzimuthNames("pitch",l)},roll:{angle:c,name:this.findAzimuthNames("roll",c)}}}findAzimuthNames(t,e){const n=[],i=this.listMap[t];if(!i)return n;for(const{name:t,range:[a,o]}of i)(e>=a&&e<o||e<=a&&e>o)&&n.push(t);return n}};Dt(Ot,"_options"),new Ot;const It={[r]:1,[s]:1,[l]:1,[c]:3,[u]:2,[d]:4,[m]:4,[p]:1,[f]:2,[h]:1,[v]:2};function Tt(t){return It[t]}new n;var Ft=Object.defineProperty,Rt=(t,e,n)=>(((t,e,n)=>{e in t?Ft(t,e,{enumerable:!0,configurable:!0,writable:!0,value:n}):t[e]=n})(t,"symbol"!=typeof e?e+"":e,n),n);var Vt=(t=>(t[t.Fill=0]="Fill",t[t.Align=1]="Align",t[t.Raw=2]="Raw",t))(Vt||{});class Et extends g{constructor(t){const{map:e,opacity:n,accFactor:a,steps:r,alphaRange:s,fit:l,atomize:c,side:u,containerMin:d,containerMax:m,uniforms:p,...f}=t??{},h=n??1,v=a??1,g=r??100,_=(new o).copy(s??{x:0,y:.95}),y=l??0,w=c??!0,M=u??N,z=(new i).copy(d??{x:0,y:0,z:0}),D=new i;if(m)D.copy(m);else if(e){const{width:t,height:n,depth:a}=e.image,o=new i(t-1,n-1,a-1);D.addVectors(z,o)}const S={map:{value:e},containerMin:{value:z},containerMax:{value:D},fit:{value:y},opacity:{value:h},alphaRange:{value:_},accFactor:{value:v},steps:{value:g},isDoubleSide:{value:!1},atomize:{value:w},...p};super({glslVersion:x,uniforms:S,transparent:!0,vertexShader:"precision highp sampler3D;\n\nuniform sampler3D map;\nuniform vec3 containerMin;\nuniform vec3 containerMax;\n\n// 0: 充满；1: 对齐；2: 原始\nuniform int fit;\n\n\nout vec3 cameraPos;\nout vec3 lookDir;\nout vec3 displayMin;\nout vec3 displayMax;\n\n\n\nvoid main() {\n    cameraPos = vec3( inverse( modelMatrix ) * vec4( cameraPosition, 1.0 ) ).xyz;\n    lookDir = position - cameraPos;\n\n    vec3 containerSize = containerMax - containerMin;\n    vec3 origin = containerMin;\n    vec3 textSize =  vec3(textureSize(map,0));\n    displayMin = containerMin;\n    displayMax = containerMax;\n\n    switch (fit) {\n        case 1:\n            displayMax = displayMin + textSize;\n            break;\n        case 2:\n            origin = vec3(0.0);\n            break;\n        default:\n            textSize = containerSize;\n    }\n\n    vec3 scale = 1.0/textSize;\n    cameraPos = (cameraPos - origin)*scale;\n    lookDir *= scale;\n\n    displayMin = max((displayMin - origin)*scale,vec3(0.0));\n    displayMax = min((displayMax - origin)*scale,vec3(1.0));\n    \n    gl_Position = projectionMatrix * modelViewMatrix * vec4( position, 1.0 );\n}",...f}),Rt(this,"isVolumeMaterial",!0),Rt(this,"_side",N),this.opacity=h,this.side=M}get map(){return this.uniforms.map.value}set map(t){this.uniforms.map.value=t,this.uniformsNeedUpdate=!0}get containerMin(){return this.uniforms.containerMin.value}set containerMin(t){this.uniforms.containerMin.value.copy(t),this.uniformsNeedUpdate=!0}get containerMax(){return this.uniforms.containerMax.value}set containerMax(t){this.uniforms.containerMax.value.copy(t),this.uniformsNeedUpdate=!0}get containerSize(){return this.containerMax.clone().sub(this.containerMin)}get fit(){return this.uniforms.fit.value}set fit(t){this.uniforms.fit.value=t,this.uniformsNeedUpdate=!0}getFitTranslate(){return 2===this.fit?new i:this.containerMin.clone()}getFitScale(){var t;const e=new i(1,1,1),n=null==(t=this.map)?void 0:t.image,a=this.fit;if(!n||2===a||1===a)return e;const o=new i(n.width,n.height,n.depth);return e.copy(this.containerSize).divide(o),e}getFitMatrix(){var t;const e=new _,n=null==(t=this.map)?void 0:t.image,a=this.fit;if(!n||2===a)return e;if(0===a){const t=new i(n.width,n.height,n.depth),a=this.containerSize.divide(t);e.makeScale(a.x,a.y,a.z)}return e.setPosition(this.containerMin),e}getFitMatrixInvert(){return this.getFitMatrix().invert()}get side(){return this._side}set side(t){this._side=t,this.uniforms&&(this.uniforms.isDoubleSide.value=t===y,this.uniformsNeedUpdate=!0)}get opacity(){return this.uniforms.opacity.value}set opacity(t){this.uniforms&&(this.uniforms.opacity.value=t),this.uniformsNeedUpdate=!0}get alphaRange(){return this.uniforms.alphaRange.value}set alphaRange(t){this.uniforms.alphaRange.value.copy(t),this.uniformsNeedUpdate=!0}get atomize(){return this.uniforms.atomize.value}set atomize(t){this.uniforms.atomize.value=t,this.uniformsNeedUpdate=!0}get steps(){return this.uniforms.steps.value}set steps(t){this.uniforms.steps.value=t,this.uniformsNeedUpdate=!0}get accFactor(){return this.uniforms.accFactor.value}set accFactor(t){this.uniforms&&(this.uniforms.accFactor.value=t),this.uniformsNeedUpdate=!0}toMapPosition(t){const e=new i(t.x,t.y,t.z).applyMatrix4(this.getFitMatrixInvert());return e.set(Math.trunc(e.x),Math.trunc(e.y),Math.trunc(e.z)),e}toMapDepth(t,e){const a=k.toKey(t),o=new i;o[a]=1;const r=(new n).setFromMatrix4(this.getFitMatrix());return o.applyMatrix3(r),e/=o.length(),Math.trunc(e)}getData3DSlice(t,e){const n=this.map;return null!=n&&n.image?function(t,e,n){const{data:i,width:a,height:o,depth:r}=t.image;return dt({data:i,size:{x:a,y:o,z:r}},e,n,Tt(t.format))}(n,t,e=this.toMapDepth(t,e)):null}getItem(t){const e=this.map;return null!=e&&e.image?function(t,e){const{data:n,width:i,height:a,depth:o}=t.image;return rt({data:n,size:{x:i,y:a,z:o}},e,Tt(t.format))}(e,this.toMapPosition(t)).value:null}}function Gt(t){t.boundingBox||t.computeBoundingBox();const{x:e,y:n,z:i}=t.boundingBox.min;return t.translate(-e,-n,-i),t}class jt extends w{constructor(t){var e;const{width:n,height:i,depth:a}=(null==(e=t.map)?void 0:e.image)||{},o=new M(n,i,a);Gt(o),super(o,t),Rt(this,"isVolumeMesh",!0),Rt(this,"_geometry"),Rt(this,"autoUpdateMaterial",!0),Rt(this,"_material"),Rt(this,"autoUpdateGeometry",!0),Rt(this,"autoNormalize",!1),Object.defineProperties(this,{geometry:{get:()=>this._geometry,set:t=>{this._geometry=t,this.autoNormalize?this.normalize():this.autoUpdateMaterial&&this.updateMaterial()}},material:{get:()=>this._material,set:t=>{this._material=t,this.autoUpdateGeometry&&this.updateGeometry()}}}),this.geometry=o,this.material=t}updateMaterial(){const{geometry:t,material:e}=this;if(!e||!t)return!1;t.boundingBox||t.computeBoundingBox();const{min:n,max:i}=t.boundingBox;e.containerMin=n,e.containerMax=i}updateGeometry(){var t;const{geometry:e,material:n,autoUpdateMaterial:a}=this,o=null==(t=n.map)?void 0:t.image;if(!o)return!1;e.boundingBox||e.computeBoundingBox();const{width:r,height:s,depth:l}=o,c=e.boundingBox.getSize(new i);e.scale(r/c.x,s/c.y,l/c.z),a&&this.updateMaterial()}normalize(){const{geometry:t,autoUpdateMaterial:e}=this;t.boundingBox||t.computeBoundingBox();const{x:n,y:i,z:a}=t.boundingBox.min;t.translate(-n,-i,-a),e&&this.updateMaterial()}toMapPosition(t){const e=this.worldToLocal(new i(t.x,t.y,t.z));return this.material.toMapPosition(e)}toMaterialDepth(t,e){const a=k.toKey(t),o=new i;o[a]=1;const r=(new n).setFromMatrix4(this.matrixWorld);return o.applyMatrix3(r),e/=o.length()}toMapDepth(t,e){return e=this.toMaterialDepth(t,e),this.material.toMapDepth(t,e)}getData3DSlice(t,e){return e=this.toMaterialDepth(t,e),this.material.getData3DSlice(t,e)}getItem(t){const e=this.worldToLocal(new i(t.x,t.y,t.z));return this.material.getItem(e)}}class Lt extends Et{constructor(t){const{gradient:e,range:n,discardOut:i,voidRange:a,...r}=t??{},s=(new o).copy(n??{x:0,y:100}),l=i??!0,c=(new o).copy(a??{x:-100,y:-1});super({...r,fragmentShader:"precision highp float;\nprecision highp sampler3D;\n\nuniform sampler3D map;\nuniform sampler2D gradient;\n\nuniform bool atomize;\nuniform float steps;\n// 颜色累积系数\nuniform float accFactor;\nuniform vec2 range;\n// 是否丢弃超出范围的像素\nuniform bool discardOut;\n// 空值范围\nuniform vec2 voidRange;\nuniform float opacity;\nuniform vec2 alphaRange;\nuniform bool isDoubleSide;\n\nin vec3 cameraPos;\nin vec3 lookDir;\nin vec3 displayMin;\nin vec3 displayMax;\nout vec4 fragColor;\n\nvec4 colorBlend(vec4 near,vec4 far){\n    float nA = near.a;\n    if  (nA >= 1.0){\n        return near;\n    }\n    if  (nA <= 0.0){\n        return far;\n    }\n    \n    float fA = far.a;\n    float a = fA + nA - fA*nA;\n    vec3 color = (far.rgb * fA * (1.0 - nA) + near.rgb*nA)/a;\n    return vec4(color,a);\n}\n\n\n\n\nvec2 intersectBox( vec3 orig, vec3 dir ) {\n    vec3 inv_dir = 1.0 / dir;\n    vec3 tmin_tmp = (displayMin - orig ) * inv_dir;\n    vec3 tmax_tmp = (displayMax - orig ) * inv_dir;\n    vec3 tmin = min( tmin_tmp, tmax_tmp );\n    vec3 tmax = max( tmin_tmp, tmax_tmp );\n    float t0 = max( tmin.x, max( tmin.y, tmin.z ) );\n    float t1 = min( tmax.x, min( tmax.y, tmax.z ) );\n    return vec2( t0, t1 );\n}\n\nfloat getValue( vec3 point ) {\n    return texture( map, point ).r;\n}\n\n// 获取梯度颜色\nvec4 getGradientColor(float val) {\n    val = clamp( val, 0.0,1.0);\n    return texture2D(gradient, vec2(val, 0.5));\n}\n\n\n\n\nvoid main(){\n    vec3 rayDir = normalize( lookDir );\n\n    vec2 times = intersectBox( cameraPos, rayDir );\n    float tMin = times.x;\n    float tMax = times.y;\n    if ( tMin > tMax || tMax < 0.0 ) discard;\n    tMin = max( tMin, 0.0 );\n\n    vec3 point = cameraPos + tMin * rayDir;\n    \n    float step = 1.0/steps;\n    float opacityFactor = atomize ? step * accFactor : 1.0;\n\n\n    float invRangeLen = 1.0/(range.y - range.x);\n    vec3 stepDir = rayDir * step;\n\n    vec4 finalColor = vec4( 0,0,0,0);\n    float alphaMax = alphaRange.y;\n\n    for ( float t = tMin; t <= tMax; t += step,point += stepDir ) {\n        float val = getValue( point);\n        if (voidRange.x <= val && val <= voidRange.y ) continue;\n        val = (val - range.x) * invRangeLen;\n        if (discardOut && (val < 0.0 || val > 1.0)) continue;\n\n        vec4 gradientColor = getGradientColor(val);\n        gradientColor.a *= opacityFactor;\n        finalColor = colorBlend(finalColor,gradientColor);\n        if ( finalColor.a >= alphaMax ) break;\n    }\n\n    finalColor.a *= opacity;\n    if ( finalColor.a <= alphaRange.x ) discard;\n    if (isDoubleSide){\n        finalColor.a = 1.0 - sqrt(1.0 - finalColor.a);\n    }\n    fragColor = finalColor;\n\n}",uniforms:{gradient:{value:null==e?e:Ct(e)},range:{value:s},discardOut:{value:l},voidRange:{value:c}}}),Rt(this,"isGradientVolumeMaterial",!0)}get gradient(){return this.uniforms.gradient.value}set gradient(t){this.uniforms.gradient.value=t,this.uniformsNeedUpdate=!0}get range(){return this.uniforms.range.value}set range(t){this.uniforms.range.value.copy(t),this.uniformsNeedUpdate=!0}get discardOut(){return this.uniforms.discardOut.value}set discardOut(t){this.uniforms.discardOut.value=t,this.uniformsNeedUpdate=!0}get voidRange(){return this.uniforms.voidRange.value}set voidRange(t){this.uniforms.voidRange.value.copy(t),this.uniformsNeedUpdate=!0}}class Bt extends z{constructor(t,e,n,i,a){super(t,e,n,i),Rt(this,"isGradientData3DTexture",!0),this.format=s,this.type=a??D,this.minFilter=this.magFilter=S,this.unpackAlignment=1,this.needsUpdate=!0}getData3DSlice(t,e){const{data:n,width:i,height:a,depth:o}=this.image;return dt({data:n,size:{x:i,y:a,z:o}},t,e,1)}getItem(t){const{data:e,width:n,height:i,depth:a}=this.image;return rt({data:e,size:{x:n,y:i,z:a}},t).value}getValue(t){const{data:e,width:n,height:i,depth:a}=this.image;return lt({data:e,size:{x:n,y:i,z:a}},t).value}}function kt(t,e){const{data:n,size:i}=t,{x:a,y:o,z:r}=i,s=a*o*r,{voidValue:l=0,uint8:c}=e??{};let u;if(c){u=new Uint8Array(s);for(let t=0;t<s;t++){const e=n[t]??l;u[t]=Math.trunc(e)}}else{u=new Float32Array(s);for(let t=0;t<s;t++){const e=n[t]??l;u[t]=e}}return new Bt(u,a,o,r,c?b:D)}class Ht extends Et{constructor(t){super({...t,fragmentShader:"precision highp float;\nprecision highp sampler3D;\n\nuniform sampler3D map;\n\n\n\n// 颜色累积系数\nuniform float accFactor;\nuniform bool isDoubleSide;\nuniform bool atomize;\n\n\nuniform float opacity;\nuniform vec2 alphaRange;\nuniform float steps;\n\nin vec3 cameraPos;\nin vec3 lookDir;\nin vec3 displayMin;\nin vec3 displayMax;\nout vec4 fragColor;\n\n// 混合颜色\nvec4 colorBlend(vec4 near,vec4 far){\n    float nA = near.a;\n    if  (nA >= 1.0){\n        return near;\n    }\n    if  (nA <= 0.0){\n        return far;\n    }\n    \n    float fA = far.a;\n    float a = fA + nA - fA*nA;\n    vec3 color = (far.rgb * fA * (1.0 - nA) + near.rgb*nA)/a;\n    return vec4(color,a);\n}\n\n\n\n\n// 包围盒求交\nvec2 intersectBox( vec3 orig, vec3 dir ) {\n    vec3 inv_dir = 1.0 / dir;\n    vec3 tmin_tmp = (displayMin - orig ) * inv_dir;\n    vec3 tmax_tmp = (displayMax - orig ) * inv_dir;\n    vec3 tmin = min( tmin_tmp, tmax_tmp );\n    vec3 tmax = max( tmin_tmp, tmax_tmp );\n    float t0 = max( tmin.x, max( tmin.y, tmin.z ) );\n    float t1 = min( tmax.x, min( tmax.y, tmax.z ) );\n    return vec2( t0, t1 );\n}\n\n\n\nvoid main(){\n    vec3 rayDir = normalize( lookDir );\n\n    vec2 times = intersectBox( cameraPos, rayDir );\n    float tMin = times.x;\n    float tMax = times.y;\n    if ( tMin > tMax || tMax < 0.0 ) discard;\n    tMin = max( tMin, 0.0 );\n\n    vec3 point = cameraPos + tMin * rayDir;\n    \n    float step = 1.0/steps;\n    float opacityFactor = atomize ? step * accFactor : 1.0;\n\n    vec3 stepDir = rayDir * step;\n\n    vec4 finalColor = vec4( 0,0,0,0);\n    float alphaMax = alphaRange.y;\n\n    for ( float t = tMin; t <= tMax; t += step,point += stepDir ) {\n        vec4 textureColor = texture( map, point );\n        textureColor.a *= opacityFactor;\n        finalColor = colorBlend(finalColor,textureColor);\n        if ( finalColor.a >= alphaMax ) break;\n    }\n\n    finalColor.a *= opacity;\n    if ( finalColor.a <= alphaRange.x ) discard;\n    if (isDoubleSide){\n        finalColor.a = 1.0 - sqrt(1.0 - finalColor.a);\n    }\n    fragColor = finalColor;\n\n}"}),Rt(this,"isImageVolumeMaterial",!0)}}function qt(t,e,n){const i=n??{},{x:a,y:o}=i.range??{x:0,y:100},r=i.discardOut??!0,{x:s,y:l}=i.voidRange??{x:-100,y:-1},c=i.voidColor??[0,0,0,0],{width:u,height:d}=e,m=Math.trunc(d/2),p=u-1;function f(t){const n=Math.trunc(p*t);return Q(e,m,n)}const h=1/(o-a),{data:v,size:g}=t,x=v.length,_=[];for(let t=0;t<x;t++){let e=c,n=v[t];(n<=s||l<=n)&&(n=(n-a)*h,r&&(n<0||n>1)||(n=Math.min(Math.max(n,0),1),e=f(n))),_.push(...e)}return{data:_,size:g}}class Kt extends z{constructor(t,e,n,i){super(t,e,n,i),Rt(this,"isImageData3DTexture",!0),this.format=d,this.type=b,this.minFilter=this.magFilter=S,this.unpackAlignment=1,this.needsUpdate=!0}getData3DSlice(t,e){const{data:n,width:i,height:a,depth:o}=this.image;return dt({data:n,size:{x:i,y:a,z:o}},t,e,4)}getItem(t){return this.getColor(t)}getColor(t){const{data:e,width:n,height:i,depth:a}=this.image;return rt({data:e,size:{x:n,y:i,z:a}},t,4).value}}function Wt(t,e,n){const i=(null==n?void 0:n.reverseY)??!0,{data:a,width:o,height:r,depth:s,colorSpace:l}=pt(t,e,{...n,reverseY:i}),c=new Kt(a,o,r,s);return"srgb"===l&&(c.colorSpace="srgb"),c}function Yt(t,e){const n=e.gradient,i=Array.isArray(n)?W(n,256,1):et(n)?n:Y(n);let a=t;nt(t)&&(a={data:t.data,size:{x:t.width,y:t.height,z:t.depth}});const{data:o,size:r}=qt(a,i,e),s=Uint8ClampedArray.from(o);return new Kt(s,r.x,r.y,r.z)}class Xt extends g{constructor(t){const e=t??{},n=e.color??new C(1,1,1),i=e.solid??0,a=e.exp??1,o=e.opacity??1,r=e.side??N;super({glslVersion:x,uniforms:{color:{value:n},solid:{value:i},exp:{value:a},opacity:{value:o},isDoubleSide:{value:!1}},vertexShader:"out vec3 lookDir;\nout vec3 normalDir;\n\n#if defined( USE_COLOR_ALPHA ) || defined( USE_COLOR )\n    out vec4 verColor;\n#endif\n\nvoid main() {\n    vec3 cameraPos = vec3( inverse( modelMatrix ) * vec4( cameraPosition, 1.0 ) ).xyz;\n    lookDir = position - cameraPos;\n    normalDir = normal;\n\n    #if defined( USE_COLOR_ALPHA )\n        verColor = color;\n    #elif defined( USE_INSTANCING_COLOR )\n        overColor = vec3(instanceColor,1.0);\n    #elif defined( USE_COLOR )\n        overColor = vec3(color,1.0);\n    #endif\n    gl_Position = projectionMatrix * modelViewMatrix * vec4( position, 1.0 );\n}",fragmentShader:"uniform vec3 color;\nuniform float solid;\nuniform float exp;\n\nuniform float opacity;\nuniform bool isDoubleSide;\n\nin vec3 lookDir;\nin vec3 normalDir;\n\n#if defined( USE_COLOR_ALPHA ) || defined( USE_COLOR )\n    in  vec4 verColor;\n#endif\n\nout vec4 fragColor;\n\n\nvoid main(){\n    #if defined( USE_COLOR_ALPHA ) || defined( USE_COLOR )\n        fragColor = verColor;\n    #else\n        fragColor = vec4(color,1.0);\n    #endif\n\n    vec3 dir = normalize(lookDir);\n    vec3 nor = normalize(normalDir);\n    float depth = abs(dot(nor,dir));\n    float dist = sqrt(1.0 - pow(depth,2.0));\n\n    float finalAlpha = fragColor.a * depth;\n    if (dist > solid) {\n        float distAlpha = 1.0 - (dist - solid)/(1.0 - solid);\n        finalAlpha *= pow(distAlpha,exp);\n    }\n    \n    finalAlpha *= opacity;\n    if ( finalAlpha == 0.0 ) discard;\n    if (isDoubleSide){\n        finalAlpha = 1.0 - sqrt(1.0 - finalAlpha);\n    }\n    fragColor.a = finalAlpha;\n}",side:N,transparent:!0}),Rt(this,"isSphereFogMaterial",!0),Rt(this,"_side",N),this.opacity=o,this.side=r}get side(){return this._side}set side(t){this._side=t,this.uniforms&&(this.uniforms.isDoubleSide.value=t===y,this.uniformsNeedUpdate=!0)}get color(){return this.uniforms.color.value}set color(t){this.uniforms.color.value.copy(t),this.uniformsNeedUpdate=!0}get opacity(){return this.uniforms.opacity.value}set opacity(t){this.uniforms&&(this.uniforms.opacity.value=t),this.uniformsNeedUpdate=!0}get solid(){return this.uniforms.solid.value}set solid(t){this.uniforms&&(this.uniforms.solid.value=t),this.uniformsNeedUpdate=!0}get exp(){return this.uniforms.exp.value}set exp(t){this.uniforms&&(this.uniforms.exp.value=t),this.uniformsNeedUpdate=!0}}class Zt extends g{constructor(t,e){const{solid:n=0,exp:i=1,useRadius:a=!0,sizeAttenuation:o=!0,star:r=!1,...s}=t??{},{uniforms:l,...c}=e||{},u={star:{value:r},solid:{value:n},exp:{value:i}};super({uniforms:U.merge([A.points.uniforms,u,l]),transparent:!0,depthTest:!1,defines:{useRadius:a},...c}),Rt(this,"isFogPointsMaterial",!0),Rt(this,"sizeAttenuation",!0),Object.assign(this,{...s,solid:n,exp:i,sizeAttenuation:o})}get map(){return this.uniforms.map.value}set map(t){this.uniforms.map.value=t,this.uniformsNeedUpdate=!0}get star(){return this.uniforms.star.value}set star(t){this.uniforms.star.value=t,this.uniformsNeedUpdate=!0}get size(){return this.uniforms.size.value}set size(t){this.uniforms.size.value=t,this.uniformsNeedUpdate=!0}get opacity(){return this.uniforms.opacity.value}set opacity(t){this.uniforms&&(this.uniforms.opacity.value=t),this.uniformsNeedUpdate=!0}get solid(){return this.uniforms.solid.value}set solid(t){this.uniforms&&(this.uniforms.solid.value=t),this.uniformsNeedUpdate=!0}get exp(){return this.uniforms.exp.value}set exp(t){this.uniforms&&(this.uniforms.exp.value=t),this.uniformsNeedUpdate=!0}get scale(){return this.uniforms.scale.value}set scale(t){this.uniforms&&(this.uniforms.scale.value=t),this.uniformsNeedUpdate=!0}get fogDensity(){return this.uniforms.fogDensity.value}set fogDensity(t){this.uniforms&&(this.uniforms.fogDensity.value=t),this.uniformsNeedUpdate=!0}get fogNear(){return this.uniforms.fogNear.value}set fogNear(t){this.uniforms&&(this.uniforms.fogNear.value=t),this.uniformsNeedUpdate=!0}get fogFar(){return this.uniforms.fogFar.value}set fogFar(t){this.uniforms&&(this.uniforms.fogFar.value=t),this.uniformsNeedUpdate=!0}get fogColor(){return this.uniforms.fogColor.value}set fogColor(t){this.uniforms.fogColor.value.copy(t),this.uniformsNeedUpdate=!0}}class Jt extends Zt{constructor(t){super(t,{vertexShader:"#ifdef useRadius\n\tattribute float radius;\n#endif\n\nuniform float size;\nuniform float scale;\n\n#include <common>\n#include <color_pars_vertex>\n#include <fog_pars_vertex>\n#include <morphtarget_pars_vertex>\n#include <logdepthbuf_pars_vertex>\n#include <clipping_planes_pars_vertex>\n\n#ifdef USE_POINTS_UV\n\n\tvarying vec2 vUv;\n\tuniform mat3 uvTransform;\n\n#endif\n\nvoid main() {\n\n\t#ifdef USE_POINTS_UV\n\n\t\tvUv = ( uvTransform * vec3( uv, 1 ) ).xy;\n\n\t#endif\n\n\t#include <color_vertex>\n\t#include <morphcolor_vertex>\n\t#include <begin_vertex>\n\t#include <morphtarget_vertex>\n\t#include <project_vertex>\n\n\t#ifdef useRadius\n\t\tgl_PointSize = radius;\n\t#else\n\t\tgl_PointSize = size;\n\t#endif\n\n\t#ifdef USE_SIZEATTENUATION\n\n\t\tbool isPerspective = isPerspectiveMatrix( projectionMatrix );\n\n\t\tif ( isPerspective ) gl_PointSize *= ( scale / - mvPosition.z );\n\n\t#endif\n\n\t#include <logdepthbuf_vertex>\n\t#include <clipping_planes_vertex>\n\t#include <worldpos_vertex>\n\t#include <fog_vertex>\n\n}",fragmentShader:"uniform float solid;\nuniform float exp;\nuniform bool star;\n\nuniform vec3 diffuse;\nuniform float opacity;\n\n#include <common>\n#include <color_pars_fragment>\n#include <map_particle_pars_fragment>\n#include <alphatest_pars_fragment>\n#include <fog_pars_fragment>\n#include <logdepthbuf_pars_fragment>\n#include <clipping_planes_pars_fragment>\n\n\nvec4 getColor(vec4 originColor){\n\tfloat finalAlpha = originColor.a;\n\tfloat dist =  0.0;\n\tif (star){\n\t\tvec2 coord = abs(gl_PointCoord - vec2(0.5));\n\t\tdist =  max(coord.x,coord.y);\n\t}else{\n\t\tdist = distance(gl_PointCoord,vec2(0.5));\n\t\tif (dist>0.5) discard;\n\t}\n\tdist /= 0.5;\n\tfloat depth = sqrt(1.0 - pow(dist,2.0));\n\tfinalAlpha *= depth;\n    if (dist > solid) {\n        float distAlpha = 1.0 - (dist - solid)/(1.0 - solid);\n        finalAlpha *= pow(distAlpha,exp);\n    }\n\toriginColor.a = finalAlpha;\n\treturn originColor;\n}\n\nvoid main() {\n\n\t#include <clipping_planes_fragment>\n\n\tvec3 outgoingLight = vec3( 0.0 );\n\tvec4 diffuseColor = vec4( diffuse, opacity );\n\n\t#include <logdepthbuf_fragment>\n\t#include <map_particle_fragment>\n\t#include <color_fragment>\n\t#include <alphatest_fragment>\n\n\tdiffuseColor = getColor(diffuseColor);\n\toutgoingLight = diffuseColor.rgb;\n\n\t#include <output_fragment>\n\t#include <tonemapping_fragment>\n\t#include <encodings_fragment>\n\t#include <fog_fragment>\n\t#include <premultiplied_alpha_fragment>\n\n}"}),Rt(this,"isColorFogPointsMaterial",!0)}get color(){return this.uniforms.diffuse.value}set color(t){this.uniforms.diffuse.value.copy(t),this.uniformsNeedUpdate=!0}get alphaMap(){return this.uniforms.alphaMap.value}set alphaMap(t){this.uniforms&&(this.uniforms.alphaMap.value=t),this.uniformsNeedUpdate=!0}get uvTransform(){return this.uniforms.uvTransform.value}set uvTransform(t){this.uniforms&&this.uniforms.uvTransform.value.copy(t),this.uniformsNeedUpdate=!0}}function Qt(t){const e=new Jt(t),n=function(t){const{points:e,radius:n=10,color:i=[0,0,0,0]}=t,a=[],o=[],r=[];for(const{x:t,y:s,z:l,radius:c,color:u}of e){a.push(t,s,l),o.push(c??n);const e=u??i,d=e[3];e[3]=Math.trunc(255*d),r.push(...e)}const s=new E;return s.setAttribute("position",new G(a,3)),s.setAttribute("radius",new G(o,1)),s.setAttribute("color",new j(r,4,!0)),s}(t);return new P(n,e)}class $t extends Zt{constructor(t){const{range:e,...n}=t||{};super(n,{vertexShader:"attribute float value;\nuniform vec2 range;\nvarying float vValue;\n#ifdef useRadius\n\tattribute float radius;\n#endif\n\n\nuniform float size;\nuniform float scale;\n\n#include <common>\n// #include <color_pars_vertex>\n#include <fog_pars_vertex>\n#include <morphtarget_pars_vertex>\n#include <logdepthbuf_pars_vertex>\n#include <clipping_planes_pars_vertex>\n\n#ifdef USE_POINTS_UV\n\n\tvarying vec2 vUv;\n\tuniform mat3 uvTransform;\n\n#endif\n\nvoid main() {\n\tvValue = (value - range.x)/(range.y - range.x);\n\n\t#ifdef USE_POINTS_UV\n\n\t\tvUv = ( uvTransform * vec3( uv, 1 ) ).xy;\n\n\t#endif\n\n\t// #include <color_vertex>\n\t// #include <morphcolor_vertex>\n\t#include <begin_vertex>\n\t#include <morphtarget_vertex>\n\t#include <project_vertex>\n\n\t#ifdef useRadius\n\t\tgl_PointSize = radius;\n\t#else\n\t\tgl_PointSize = size;\n\t#endif\n\n\t#ifdef USE_SIZEATTENUATION\n\n\t\tbool isPerspective = isPerspectiveMatrix( projectionMatrix );\n\n\t\tif ( isPerspective ) gl_PointSize *= ( scale / - mvPosition.z );\n\n\t#endif\n\n\t#include <logdepthbuf_vertex>\n\t#include <clipping_planes_vertex>\n\t#include <worldpos_vertex>\n\t#include <fog_vertex>\n\n}",fragmentShader:"varying float vValue;\nuniform sampler2D map;\n\n\nuniform float solid;\nuniform float exp;\nuniform bool star;\n\n// uniform vec3 diffuse;\nuniform float opacity;\n\n#include <common>\n// #include <color_pars_fragment>\n// #include <map_particle_pars_fragment>\n#include <alphatest_pars_fragment>\n#include <fog_pars_fragment>\n#include <logdepthbuf_pars_fragment>\n#include <clipping_planes_pars_fragment>\n\n\n// 获取梯度颜色\nvec4 getGradientColor(float val) {\n\tval = clamp( val, 0.0,1.0);\n    return texture2D(map, vec2(val, 0.5));\n}\n\n\nvec4 getColor(){\n\tfloat val = vValue;\n\tfloat dist =  0.0;\n\tif (star){\n\t\tvec2 coord = abs(gl_PointCoord - vec2(0.5));\n\t\tdist =  max(coord.x,coord.y);\n\t}else{\n\t\tdist = distance(gl_PointCoord,vec2(0.5));\n\t\tif (dist>0.5) discard;\n\t}\n\tdist /= 0.5;\n\tfloat depth = sqrt(1.0 - pow(dist,2.0));\n    if (dist > solid) {\n        float solidDist = 1.0 - (dist - solid)/(1.0 - solid);\n        val *= pow(solidDist,exp);\n    }\n\tvec4 valColor = getGradientColor(val);\n\tvalColor.a *= depth * opacity;\n\treturn valColor;\n}\n\nvoid main() {\n\n\t#include <clipping_planes_fragment>\n\n\tvec3 outgoingLight = vec3( 0.0 );\n\tvec4 diffuseColor = getColor();\n\n\t#include <logdepthbuf_fragment>\n\t// #include <map_particle_fragment>\n\t// #include <color_fragment>\n\t#include <alphatest_fragment>\n\n\toutgoingLight = diffuseColor.rgb;\n\n\t#include <output_fragment>\n\t#include <tonemapping_fragment>\n\t#include <encodings_fragment>\n\t#include <fog_fragment>\n\t#include <premultiplied_alpha_fragment>\n\n}",uniforms:{range:{value:(new o).copy(e??{x:0,y:100})}}}),Rt(this,"isColorFogPointsMaterial",!0)}get range(){return this.uniforms.range.value}set range(t){this.uniforms.range.value.copy(t),this.uniformsNeedUpdate=!0}}function te(t){const e=t.gradient,n=null==e?e:Ct(e),i=new $t({...t,map:n}),a=function(t){const{points:e,radius:n=10,value:i=100}=t,a=[],o=[],r=[];for(const{x:t,y:s,z:l,radius:c,value:u}of e)a.push(t,s,l),o.push(c??n),r.push(u??i);const s=new E;return s.setAttribute("position",new G(a,3)),s.setAttribute("radius",new G(o,1)),s.setAttribute("value",new G(r,1)),s}(t);return new P(a,i)}class ee extends g{constructor(t){const e=t??{};let{map:n,opacity:i,axis:a,depth:o,side:r,...s}=e;i=i??1,a=a??k.z,o=o??0,r=r??y,super({...s,glslVersion:x,uniforms:{map:{value:n},opacity:{value:i},axis:{value:a},depth:{value:o}},transparent:!0,vertexShader:"\tout vec2 vUv;\n    \n\tvoid main() {\n\t\tvUv = uv;\n\t\tgl_Position = projectionMatrix * modelViewMatrix * vec4( position, 1.0 );\n\t}",fragmentShader:"precision highp float;\nprecision highp sampler3D;\n\nuniform sampler3D map;\n\nuniform float opacity;\nuniform int axis;\nuniform float depth;\n\nin vec2 vUv;\nout vec4 fragColor;\n\nvoid main(){\n    vec3 uv3 = vec3(0.0);\n    int depthSize = textureSize(map,0)[axis];\n    uv3[axis] = depth/float(depthSize);\n    int xIndex = (axis + 1)%3;\n    int yIndex = (axis + 2)%3;\n    uv3[xIndex] = vUv.x;\n    uv3[yIndex] = vUv.y;\n\n    fragColor = texture(map, uv3);\n    fragColor.a *= opacity;\n}",side:r}),Rt(this,"isSliceMaterial",!0),Rt(this,"_sliceSize",null),this.opacity=i}get map(){return this.uniforms.map.value}set map(t){this.uniforms.map.value=t,this._sliceSize=null,this.uniformsNeedUpdate=!0}get opacity(){return this.uniforms.opacity.value}set opacity(t){this.uniforms&&(this.uniforms.opacity.value=t),this.uniformsNeedUpdate=!0}get axis(){return this.uniforms.axis.value}set axis(t){this.uniforms.axis.value=t,this._sliceSize=null,this.uniformsNeedUpdate=!0}get depth(){return this.uniforms.depth.value}set depth(t){this.uniforms.depth.value=t,this.uniformsNeedUpdate=!0}get sliceSize(){let t=this._sliceSize;if(!t){this._sliceSize=t=new o(1,1);const{axis:e,map:n}=this;if(n){const{width:i,height:a,depth:o}=n.image,[r,s]=k.getCrossAxiss(e),l=[i,a,o];t.set(l[r],l[s])}}return t}}class ne extends w{constructor(t){const e=new ee(t),n=new O;super(n,e),Rt(this,"isSliceMesh",!0),Rt(this,"_geometry"),Rt(this,"_material"),Rt(this,"autoUpdateGeometry",!0),Object.defineProperties(this,{geometry:{get:()=>this._geometry,set:t=>{this._geometry=t,this.autoUpdateGeometry&&this.updateGeometry()}},material:{get:()=>this._material,set:t=>{this._material=t,this.autoUpdateGeometry&&this.updateGeometry()}}}),this.geometry=n,this.material=e}updateGeometry(){var t;const{geometry:e,material:n}=this;if(!n||!(null==(t=n.map)?void 0:t.image))return!1;e.boundingBox||e.computeBoundingBox();const a=n.sliceSize,o=e.boundingBox.getSize(new i);e.scale(a.x/o.x,a.y/o.y,1)}get map(){return this.material.map}set map(t){this.material.map=t,this.material=this.material}get axis(){return this.material.axis}set axis(t){this.material.axis=t,this.material=this.material}get depth(){return this.material.depth}set depth(t){this.material.depth=t}get sliceSize(){return this.material.sliceSize}}class ie extends Mt{constructor(t){super(nt(t)?t:t.image),Rt(this,"isImageData3DTextureSlice",!0)}set texture(t){this.image3D=nt(t)?t:t.image}}function ae({ratio:t,value:e}){return e-e*t}function oe(t,e){const{x:n,y:i}=e;let a=0;for(const{value:e,clim:{x:n,y:i}}of t)a+=(e-n)/(i-n);return a*(i-n)+n}function re(t,e,n){const{size:a,data3D:o,tags:r}=e,{clim:s,radius:l,value:c,hollow:u,valueGradient:d,valuesAccumulate:m}=t,p=(new i).copy(t),f=[],h=[],v={particles:f,values:h},g=l*u,x=p.clone().subScalar(l),_=p.clone().addScalar(l);x.max(new i(0,0,0));const y=(new i).copy(a).subScalar(1);_.min(y);let{x:w,y:M,z:z}=x;w=Math.trunc(w),M=Math.trunc(M),z=Math.trunc(z);let{x:D,y:S,z:b}=_;if(D=Math.trunc(D),S=Math.trunc(S),b=Math.trunc(b),w>D||M>S||z>b)return v;const{x:C,y:N}=a,U=C*N;for(let e=z;e<=b;e++){const a=e*U;for(let v=M;v<=S;v++){const x=a+v*C;t:for(let a=w;a<=D;a++){const _=x+a;if(r[_])continue;const y=new i(a,v,e),w=y.distanceTo(p);if(w>l||w<g)continue;r[_]=!0;const M=[],z={...{radius:l,hollow:u,hollowRadius:g,distance:w,clim:s,value:c},ratio:w/l,distance:w,point:y},D=d(z);M.push({...t,...z,value:D});for(const t of n){const e=y.distanceTo(t),{hollow:n,radius:i}=t;if(e<=i){const a=i*n;if(e<a)continue t;const o={...t,ratio:e/i,distance:e,hollowRadius:a,point:y},r=d(o);M.push({...o,value:r})}}const S=m(M,s);o[_]=S,f.push(y),h.push(S)}}}return v}function se(t){const{points:e,clim:n={x:0,y:100},radius:i=10,value:a=100,hollow:o=0,valueGradient:r=ae,valuesAccumulate:s=oe,size:l}=t,c={clim:n,radius:i,value:a,hollow:o,valueGradient:r,valuesAccumulate:s},u=e.map((t=>({...c,...t})));let d=l;l||(d=(new I).setFromPoints(u).max.addScalar(1));const m=St(u);let{x:p,y:f,z:h}=d;p=d.x=Math.trunc(p),f=d.y=Math.trunc(f),h=d.z=Math.trunc(h);const v=new Array(p*f*h),g={size:d,data3D:v,tags:new Array(p*f*h)},x=u.length;for(let t=0;t<x;t++){re(u[t],g,m[t].map((t=>u[t])))}return{data:v,size:d}}function le({startRadius:t,addedRadius:e,ratio:n}){return e*n+t}function ce({startHollow:t,addedHollow:e,ratio:n}){return e*n+t}function ue({startValue:t,addedValue:e,ratio:n}){return e*n+t}function de(t,e){const{x:n,y:i}=e;return t.reduce(((t,e)=>t+e-n),0)}function me(t){const{points:e,size:n,clim:a={x:0,y:100},radius:o=10,value:r=100,hollow:s=0,valueGradient:l=ae,radiusGradient:c=le,hollowGradient:u=ce,lineValueGradient:d=ue,valuesAccumulate:m=de}=t,p={radius:o,value:r,hollow:s},f=e.map((function(t){return Object.assign(new i,p,t)})),h=n?(new i).copy(n).subScalar(1):(new I).setFromPoints(f).max;h.x=Math.trunc(h.x),h.y=Math.trunc(h.y),h.z=Math.trunc(h.z);const v=h.clone().addScalar(1),{x:g,y:x,z:y}=v,w=g*x,M=new Array(w*y),z=[],D=[],S={particles:z,values:D,size:v,data:M},b=new i(0,0,0),C=e.length,N=Math.trunc(C/2);for(let t=0;t<N;t++){const e=2*t,n=f[e],p=f[e+1],v=new T(n,p),x=v.distance(),y=v.delta(new i).normalize(),C=new F;C.setFromUnitVectors(new i(0,1,0),y);const N=(new _).makeRotationFromQuaternion(C);N.setPosition(n);const{radius:U,value:A,hollow:P}=n,{radius:O,value:R,hollow:V}=p,E={length:x,startRadius:U,endRadius:O,addedRadius:O-U,defaultRadius:o,radius:o,startValue:A,endValue:R,addedValue:R-A,defaultValue:r,startHollow:P,endHollow:V,addedHollow:V-P,defaultHollow:s,clim:a},G=Math.max(U,O),j=new i(-G,0,-G),L=new i(G,x,G),B=new I(j,L);B.applyMatrix4(N);const{min:k,max:H}=B;k.max(b),H.min(h);let{x:q,y:K,z:W}=k;q=Math.trunc(q),K=Math.trunc(K),W=Math.trunc(W);let{x:Y,y:X,z:Z}=H;if(Y=Math.trunc(Y),X=Math.trunc(X),Z=Math.trunc(Z),q>Y||K>X||W>Z)return S;N.invert();for(let t=W;t<=Z;t++){const e=t*w;for(let n=K;n<=X;n++){const o=e+n*g;for(let e=q;e<=Y;e++){const r=o+e,s=new i(e,n,t),p=s.clone();p.applyMatrix4(N);const{x:f,y:h,z:v}=p;if(h<0||h>x)continue;const g=Math.hypot(f,v),_={...E,ratio:h/x,point:s},y=c(_);_.radius=y;const w=u(_),S=y*w;if(g>y||g<S)continue;const b=d(_);let C=l({...E,value:b,ratio:g/y,radius:y,hollow:w,hollowRadius:S,distance:g,point:s});const U=D[r];void 0!==U&&(C=m([U,C],a)),M[r]=C,z.push(s),D.push(C)}}}}return S}function pe(t,e){const a=(new i).setFromMatrixScale(e),o=Math.sqrt(a.lengthSq()/3),r=(new n).setFromMatrix4(e);let{points:s,size:l,radius:c,...u}=t;if(l&&(l=(new i).copy(l).applyMatrix3(r)),c&&(c*=o),s){const t=[];for(let{x:n,y:a,z:r,radius:l,...c}of s){l&&(l*=o);const s=new i(n,a,r).applyMatrix4(e);Object.assign(s,c,{radius:l}),t.push(s)}s=t}return{points:s,size:l,radius:c,...u}}function fe(t,e,n=1){let{points:a,size:o,radius:r,...s}=t;if(r&&(r*=n),a){const t=[];for(let{x:o,y:r,z:s,radius:l,...c}of a){l&&(l*=n);const a=new i(o,r,s).add(e).multiplyScalar(n);Object.assign(a,c,{radius:l}),t.push(a)}a=t}return o&&(o=(new i).copy(o).add(e).multiplyScalar(n)),{points:a,size:o,radius:r,...s}}function he(t,e){const{points:n,radius:a}=t;if(!n)return{options:t,position:new i,scale:1};const o=new I;for(const t of n){const e=(new i).copy(t),n=t.radius??a;n?(o.expandByPoint(e.clone().subScalar(n)),o.expandByPoint(e.addScalar(n))):o.expandByPoint(e)}const r=(new i).max(o.min),s=o.max.clone().sub(r).addScalar(1);t.size&&s.min(t.size);let{scale:l,maxSize:c}=e||{};l||(c=c??100,l=c/Math.max(s.x,s.y,s.z)),s.multiplyScalar(l);const u=fe(t,r.clone().negate(),l);return u.size=s,{options:u,position:r,scale:1/l}}const ve={uniforms:{u_size:{value:new i(1,1,1)},u_renderstyle:{value:0},u_renderthreshold:{value:.5},u_clim:{value:new o(1,1)},u_data:{value:null},u_cmdata:{value:null}},vertexShader:'\n\n\t\tvarying vec4 v_nearpos;\n\t\tvarying vec4 v_farpos;\n\t\tvarying vec3 v_position;\n\n\t\tvoid main() {\n\t\t\t\t// Prepare transforms to map to "camera view". See also:\n\t\t\t\t// https://threejs.org/docs/#api/renderers/webgl/WebGLProgram\n\t\t\t\tmat4 viewtransformf = modelViewMatrix;\n\t\t\t\tmat4 viewtransformi = inverse(modelViewMatrix);\n\n\t\t\t\t// Project local vertex coordinate to camera position. Then do a step\n\t\t\t\t// backward (in cam coords) to the near clipping plane, and project back. Do\n\t\t\t\t// the same for the far clipping plane. This gives us all the information we\n\t\t\t\t// need to calculate the ray and truncate it to the viewing cone.\n\t\t\t\tvec4 position4 = vec4(position, 1.0);\n\t\t\t\tvec4 pos_in_cam = viewtransformf * position4;\n\n\t\t\t\t// Intersection of ray and near clipping plane (z = -1 in clip coords)\n\t\t\t\tpos_in_cam.z = -pos_in_cam.w;\n\t\t\t\tv_nearpos = viewtransformi * pos_in_cam;\n\n\t\t\t\t// Intersection of ray and far clipping plane (z = +1 in clip coords)\n\t\t\t\tpos_in_cam.z = pos_in_cam.w;\n\t\t\t\tv_farpos = viewtransformi * pos_in_cam;\n\n\t\t\t\t// Set varyings and output pos\n\t\t\t\tv_position = position;\n\t\t\t\tgl_Position = projectionMatrix * viewMatrix * modelMatrix * position4;\n\t\t}',fragmentShader:"\n\n\t\t\t\tprecision highp float;\n\t\t\t\tprecision mediump sampler3D;\n\n\t\t\t\tuniform vec3 u_size;\n\t\t\t\tuniform int u_renderstyle;\n\t\t\t\tuniform float u_renderthreshold;\n\t\t\t\tuniform vec2 u_clim;\n\n\t\t\t\tuniform sampler3D u_data;\n\t\t\t\tuniform sampler2D u_cmdata;\n\n\t\t\t\tvarying vec3 v_position;\n\t\t\t\tvarying vec4 v_nearpos;\n\t\t\t\tvarying vec4 v_farpos;\n\n\t\t\t\t// The maximum distance through our rendering volume is sqrt(3).\n\t\t\t\tconst int MAX_STEPS = 887;\t// 887 for 512^3, 1774 for 1024^3\n\t\t\t\tconst int REFINEMENT_STEPS = 4;\n\t\t\t\tconst float relative_step_size = 1.0;\n\t\t\t\tconst vec4 ambient_color = vec4(0.2, 0.4, 0.2, 1.0);\n\t\t\t\tconst vec4 diffuse_color = vec4(0.8, 0.2, 0.2, 1.0);\n\t\t\t\tconst vec4 specular_color = vec4(1.0, 1.0, 1.0, 1.0);\n\t\t\t\tconst float shininess = 40.0;\n\n\t\t\t\tvoid cast_mip(vec3 start_loc, vec3 step, int nsteps, vec3 view_ray);\n\t\t\t\tvoid cast_iso(vec3 start_loc, vec3 step, int nsteps, vec3 view_ray);\n\n\t\t\t\tfloat sample1(vec3 texcoords);\n\t\t\t\tvec4 apply_colormap(float val);\n\t\t\t\tvec4 add_lighting(float val, vec3 loc, vec3 step, vec3 view_ray);\n\n\n\t\t\t\tvoid main() {\n\t\t\t\t\t\t// Normalize clipping plane info\n\t\t\t\t\t\tvec3 farpos = v_farpos.xyz / v_farpos.w;\n\t\t\t\t\t\tvec3 nearpos = v_nearpos.xyz / v_nearpos.w;\n\n\t\t\t\t\t\t// Calculate unit vector pointing in the view direction through this fragment.\n\t\t\t\t\t\tvec3 view_ray = normalize(nearpos.xyz - farpos.xyz);\n\n\t\t\t\t\t\t// Compute the (negative) distance to the front surface or near clipping plane.\n\t\t\t\t\t\t// v_position is the back face of the cuboid, so the initial distance calculated in the dot\n\t\t\t\t\t\t// product below is the distance from near clip plane to the back of the cuboid\n\t\t\t\t\t\tfloat distance = dot(nearpos - v_position, view_ray);\n\t\t\t\t\t\tdistance = max(distance, min((-0.5 - v_position.x) / view_ray.x,\n\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t(u_size.x - 0.5 - v_position.x) / view_ray.x));\n\t\t\t\t\t\tdistance = max(distance, min((-0.5 - v_position.y) / view_ray.y,\n\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t(u_size.y - 0.5 - v_position.y) / view_ray.y));\n\t\t\t\t\t\tdistance = max(distance, min((-0.5 - v_position.z) / view_ray.z,\n\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t(u_size.z - 0.5 - v_position.z) / view_ray.z));\n\n\t\t\t\t\t\t// Now we have the starting position on the front surface\n\t\t\t\t\t\tvec3 front = v_position + view_ray * distance;\n\n\t\t\t\t\t\t// Decide how many steps to take\n\t\t\t\t\t\tint nsteps = int(-distance / relative_step_size + 0.5);\n\t\t\t\t\t\tif ( nsteps < 1 )\n\t\t\t\t\t\t\t\tdiscard;\n\n\t\t\t\t\t\t// Get starting location and step vector in texture coordinates\n\t\t\t\t\t\tvec3 step = ((v_position - front) / u_size) / float(nsteps);\n\t\t\t\t\t\tvec3 start_loc = front / u_size;\n\n\t\t\t\t\t\t// For testing: show the number of steps. This helps to establish\n\t\t\t\t\t\t// whether the rays are correctly oriented\n\t\t\t\t\t\t//'gl_FragColor = vec4(0.0, float(nsteps) / 1.0 / u_size.x, 1.0, 1.0);\n\t\t\t\t\t\t//'return;\n\n\t\t\t\t\t\tif (u_renderstyle == 0)\n\t\t\t\t\t\t\t\tcast_mip(start_loc, step, nsteps, view_ray);\n\t\t\t\t\t\telse if (u_renderstyle == 1)\n\t\t\t\t\t\t\t\tcast_iso(start_loc, step, nsteps, view_ray);\n\n\t\t\t\t\t\tif (gl_FragColor.a < 0.05)\n\t\t\t\t\t\t\t\tdiscard;\n\t\t\t\t}\n\n\n\t\t\t\tfloat sample1(vec3 texcoords) {\n\t\t\t\t\t\t/* Sample float value from a 3D texture. Assumes intensity data. */\n\t\t\t\t\t\treturn texture(u_data, texcoords.xyz).r;\n\t\t\t\t}\n\n\n\t\t\t\tvec4 apply_colormap(float val) {\n\t\t\t\t\t\tval = (val - u_clim[0]) / (u_clim[1] - u_clim[0]);\n\t\t\t\t\t\treturn texture2D(u_cmdata, vec2(val, 0.5));\n\t\t\t\t}\n\n\n\t\t\t\tvoid cast_mip(vec3 start_loc, vec3 step, int nsteps, vec3 view_ray) {\n\n\t\t\t\t\t\tfloat max_val = -1e6;\n\t\t\t\t\t\tint max_i = 100;\n\t\t\t\t\t\tvec3 loc = start_loc;\n\n\t\t\t\t\t\t// Enter the raycasting loop. In WebGL 1 the loop index cannot be compared with\n\t\t\t\t\t\t// non-constant expression. So we use a hard-coded max, and an additional condition\n\t\t\t\t\t\t// inside the loop.\n\t\t\t\t\t\tfor (int iter=0; iter<MAX_STEPS; iter++) {\n\t\t\t\t\t\t\t\tif (iter >= nsteps)\n\t\t\t\t\t\t\t\t\t\tbreak;\n\t\t\t\t\t\t\t\t// Sample from the 3D texture\n\t\t\t\t\t\t\t\tfloat val = sample1(loc);\n\t\t\t\t\t\t\t\t// Apply MIP operation\n\t\t\t\t\t\t\t\tif (val > max_val) {\n\t\t\t\t\t\t\t\t\t\tmax_val = val;\n\t\t\t\t\t\t\t\t\t\tmax_i = iter;\n\t\t\t\t\t\t\t\t}\n\t\t\t\t\t\t\t\t// Advance location deeper into the volume\n\t\t\t\t\t\t\t\tloc += step;\n\t\t\t\t\t\t}\n\n\t\t\t\t\t\t// Refine location, gives crispier images\n\t\t\t\t\t\tvec3 iloc = start_loc + step * (float(max_i) - 0.5);\n\t\t\t\t\t\tvec3 istep = step / float(REFINEMENT_STEPS);\n\t\t\t\t\t\tfor (int i=0; i<REFINEMENT_STEPS; i++) {\n\t\t\t\t\t\t\t\tmax_val = max(max_val, sample1(iloc));\n\t\t\t\t\t\t\t\tiloc += istep;\n\t\t\t\t\t\t}\n\n\t\t\t\t\t\t// Resolve final color\n\t\t\t\t\t\tgl_FragColor = apply_colormap(max_val);\n\t\t\t\t}\n\n\n\t\t\t\tvoid cast_iso(vec3 start_loc, vec3 step, int nsteps, vec3 view_ray) {\n\n\t\t\t\t\t\tgl_FragColor = vec4(0.0);\t// init transparent\n\t\t\t\t\t\tvec4 color3 = vec4(0.0);\t// final color\n\t\t\t\t\t\tvec3 dstep = 1.5 / u_size;\t// step to sample derivative\n\t\t\t\t\t\tvec3 loc = start_loc;\n\n\t\t\t\t\t\tfloat low_threshold = u_renderthreshold - 0.02 * (u_clim[1] - u_clim[0]);\n\n\t\t\t\t\t\t// Enter the raycasting loop. In WebGL 1 the loop index cannot be compared with\n\t\t\t\t\t\t// non-constant expression. So we use a hard-coded max, and an additional condition\n\t\t\t\t\t\t// inside the loop.\n\t\t\t\t\t\tfor (int iter=0; iter<MAX_STEPS; iter++) {\n\t\t\t\t\t\t\t\tif (iter >= nsteps)\n\t\t\t\t\t\t\t\t\t\tbreak;\n\n\t\t\t\t\t\t\t\t// Sample from the 3D texture\n\t\t\t\t\t\t\t\tfloat val = sample1(loc);\n\n\t\t\t\t\t\t\t\tif (val > low_threshold) {\n\t\t\t\t\t\t\t\t\t\t// Take the last interval in smaller steps\n\t\t\t\t\t\t\t\t\t\tvec3 iloc = loc - 0.5 * step;\n\t\t\t\t\t\t\t\t\t\tvec3 istep = step / float(REFINEMENT_STEPS);\n\t\t\t\t\t\t\t\t\t\tfor (int i=0; i<REFINEMENT_STEPS; i++) {\n\t\t\t\t\t\t\t\t\t\t\t\tval = sample1(iloc);\n\t\t\t\t\t\t\t\t\t\t\t\tif (val > u_renderthreshold) {\n\t\t\t\t\t\t\t\t\t\t\t\t\t\tgl_FragColor = add_lighting(val, iloc, dstep, view_ray);\n\t\t\t\t\t\t\t\t\t\t\t\t\t\treturn;\n\t\t\t\t\t\t\t\t\t\t\t\t}\n\t\t\t\t\t\t\t\t\t\t\t\tiloc += istep;\n\t\t\t\t\t\t\t\t\t\t}\n\t\t\t\t\t\t\t\t}\n\n\t\t\t\t\t\t\t\t// Advance location deeper into the volume\n\t\t\t\t\t\t\t\tloc += step;\n\t\t\t\t\t\t}\n\t\t\t\t}\n\n\n\t\t\t\tvec4 add_lighting(float val, vec3 loc, vec3 step, vec3 view_ray)\n\t\t\t\t{\n\t\t\t\t\t// Calculate color by incorporating lighting\n\n\t\t\t\t\t\t// View direction\n\t\t\t\t\t\tvec3 V = normalize(view_ray);\n\n\t\t\t\t\t\t// calculate normal vector from gradient\n\t\t\t\t\t\tvec3 N;\n\t\t\t\t\t\tfloat val1, val2;\n\t\t\t\t\t\tval1 = sample1(loc + vec3(-step[0], 0.0, 0.0));\n\t\t\t\t\t\tval2 = sample1(loc + vec3(+step[0], 0.0, 0.0));\n\t\t\t\t\t\tN[0] = val1 - val2;\n\t\t\t\t\t\tval = max(max(val1, val2), val);\n\t\t\t\t\t\tval1 = sample1(loc + vec3(0.0, -step[1], 0.0));\n\t\t\t\t\t\tval2 = sample1(loc + vec3(0.0, +step[1], 0.0));\n\t\t\t\t\t\tN[1] = val1 - val2;\n\t\t\t\t\t\tval = max(max(val1, val2), val);\n\t\t\t\t\t\tval1 = sample1(loc + vec3(0.0, 0.0, -step[2]));\n\t\t\t\t\t\tval2 = sample1(loc + vec3(0.0, 0.0, +step[2]));\n\t\t\t\t\t\tN[2] = val1 - val2;\n\t\t\t\t\t\tval = max(max(val1, val2), val);\n\n\t\t\t\t\t\tfloat gm = length(N); // gradient magnitude\n\t\t\t\t\t\tN = normalize(N);\n\n\t\t\t\t\t\t// Flip normal so it points towards viewer\n\t\t\t\t\t\tfloat Nselect = float(dot(N, V) > 0.0);\n\t\t\t\t\t\tN = (2.0 * Nselect - 1.0) * N;\t// ==\tNselect * N - (1.0-Nselect)*N;\n\n\t\t\t\t\t\t// Init colors\n\t\t\t\t\t\tvec4 ambient_color = vec4(0.0, 0.0, 0.0, 0.0);\n\t\t\t\t\t\tvec4 diffuse_color = vec4(0.0, 0.0, 0.0, 0.0);\n\t\t\t\t\t\tvec4 specular_color = vec4(0.0, 0.0, 0.0, 0.0);\n\n\t\t\t\t\t\t// note: could allow multiple lights\n\t\t\t\t\t\tfor (int i=0; i<1; i++)\n\t\t\t\t\t\t{\n\t\t\t\t\t\t\t\t // Get light direction (make sure to prevent zero devision)\n\t\t\t\t\t\t\t\tvec3 L = normalize(view_ray);\t//lightDirs[i];\n\t\t\t\t\t\t\t\tfloat lightEnabled = float( length(L) > 0.0 );\n\t\t\t\t\t\t\t\tL = normalize(L + (1.0 - lightEnabled));\n\n\t\t\t\t\t\t\t\t// Calculate lighting properties\n\t\t\t\t\t\t\t\tfloat lambertTerm = clamp(dot(N, L), 0.0, 1.0);\n\t\t\t\t\t\t\t\tvec3 H = normalize(L+V); // Halfway vector\n\t\t\t\t\t\t\t\tfloat specularTerm = pow(max(dot(H, N), 0.0), shininess);\n\n\t\t\t\t\t\t\t\t// Calculate mask\n\t\t\t\t\t\t\t\tfloat mask1 = lightEnabled;\n\n\t\t\t\t\t\t\t\t// Calculate colors\n\t\t\t\t\t\t\t\tambient_color +=\tmask1 * ambient_color;\t// * gl_LightSource[i].ambient;\n\t\t\t\t\t\t\t\tdiffuse_color +=\tmask1 * lambertTerm;\n\t\t\t\t\t\t\t\tspecular_color += mask1 * specularTerm * specular_color;\n\t\t\t\t\t\t}\n\n\t\t\t\t\t\t// Calculate final color by componing different components\n\t\t\t\t\t\tvec4 final_color;\n\t\t\t\t\t\tvec4 color = apply_colormap(val);\n\t\t\t\t\t\tfinal_color = color * (ambient_color + diffuse_color) + specular_color;\n\t\t\t\t\t\tfinal_color.a = color.a;\n\t\t\t\t\t\treturn final_color;\n\t\t\t\t}"};var ge=(t=>(t[t.MIP=0]="MIP",t[t.ISO=1]="ISO",t))(ge||{});class xe extends g{constructor(t){const{style:e,threshold:n,clim:a,map:r,gradient:s}=t??{},l=new i(1,1,1);if(r){const{width:t,height:e,depth:n}=r.image;l.set(t,e,n)}const c=null==s?s:Ct(s),u=U.clone(ve.uniforms);u.u_size.value.copy(l),u.u_clim.value.copy(a??new o(1,1)),u.u_renderstyle.value=e??1,u.u_renderthreshold.value=n??.5,u.u_data.value=r??null,u.u_cmdata.value=c??null,super({...ve,uniforms:u,side:R}),Rt(this,"isExampleVolumeMaterial_1",!0),Rt(this,"_gradient")}get map(){return this.uniforms.u_data.value}set map(t){this.uniforms.u_data.value=t,this.uniformsNeedUpdate=!0}get gradientMap(){return this.uniforms.u_cmdata.value}get gradient(){return this._gradient}set gradient(t){this._gradient=t;const e=null==t?t:Ct(t);this.uniforms.u_cmdata.value=e??null,this.uniformsNeedUpdate=!0}get size(){return this.uniforms.u_size.value}set size(t){this.uniforms.u_size.value.copy(t),this.uniformsNeedUpdate=!0}get style(){return this.uniforms.u_renderstyle.value}set style(t){this.uniforms.u_renderstyle.value=t,this.uniformsNeedUpdate=!0}get clim(){return this.uniforms.u_clim.value}set clim(t){this.uniforms.u_clim.value=t,this.uniformsNeedUpdate=!0}get threshold(){return this.uniforms.u_renderthreshold.value}set threshold(t){this.uniforms.u_renderthreshold.value=t,this.uniformsNeedUpdate=!0}}const _e="in vec3 position;\n\nuniform mat4 modelMatrix;\nuniform mat4 modelViewMatrix;\nuniform mat4 projectionMatrix;\nuniform vec3 cameraPos;\n\nout vec3 vOrigin;\nout vec3 vDirection;\n\nvoid main() {\n    vec4 mvPosition = modelViewMatrix * vec4( position, 1.0 );\n    vOrigin = vec3( inverse( modelMatrix ) * vec4( cameraPos, 1.0 ) ).xyz;\n    vDirection = position - vOrigin;\n    gl_Position = projectionMatrix * mvPosition;\n}";class ye extends V{constructor(t){let{color:e,threshold:n,map:a,opacity:o,range:r,steps:s,frame:l,cameraPos:c}=t??{};e=e??new C(1,1,1),n=n??.25,o=o??.25,r=r??.1,s=s??100,l=l??100,c=c??new i,super({glslVersion:x,uniforms:{base:{value:e},map:{value:a},cameraPos:{value:c},threshold:{value:n},opacity:{value:o},range:{value:r},steps:{value:s},frame:{value:l}},vertexShader:_e,fragmentShader:"precision highp float;\nprecision highp sampler3D;\n\nuniform mat4 modelViewMatrix;\nuniform mat4 projectionMatrix;\n\nin vec3 vOrigin;\nin vec3 vDirection;\n\nout vec4 color;\n\nuniform vec3 base;\nuniform sampler3D map;\n\nuniform float threshold;\nuniform float range;\nuniform float opacity;\nuniform float steps;\nuniform float frame;\n\nuint wang_hash(uint seed)\n{\n        seed = (seed ^ 61u) ^ (seed >> 16u);\n        seed *= 9u;\n        seed = seed ^ (seed >> 4u);\n        seed *= 0x27d4eb2du;\n        seed = seed ^ (seed >> 15u);\n        return seed;\n}\n\nfloat randomFloat(inout uint seed)\n{\n        return float(wang_hash(seed)) / 4294967296.;\n}\n\nvec2 hitBox( vec3 orig, vec3 dir ) {\n    const vec3 box_min = vec3( - 0.5 );\n    const vec3 box_max = vec3( 0.5 );\n    vec3 inv_dir = 1.0 / dir;\n    vec3 tmin_tmp = ( box_min - orig ) * inv_dir;\n    vec3 tmax_tmp = ( box_max - orig ) * inv_dir;\n    vec3 tmin = min( tmin_tmp, tmax_tmp );\n    vec3 tmax = max( tmin_tmp, tmax_tmp );\n    float t0 = max( tmin.x, max( tmin.y, tmin.z ) );\n    float t1 = min( tmax.x, min( tmax.y, tmax.z ) );\n    return vec2( t0, t1 );\n}\n\nfloat sample1( vec3 p ) {\n    return texture( map, p ).r;\n}\n\nfloat shading( vec3 coord ) {\n    float step = 0.01;\n    return sample1( coord + vec3( - step ) ) - sample1( coord + vec3( step ) );\n}\n\nvoid main(){\n    vec3 rayDir = normalize( vDirection );\n    vec2 bounds = hitBox( vOrigin, rayDir );\n\n    if ( bounds.x > bounds.y ) discard;\n\n    bounds.x = max( bounds.x, 0.0 );\n\n    vec3 p = vOrigin + bounds.x * rayDir;\n    vec3 inc = 1.0 / abs( rayDir );\n    float delta = min( inc.x, min( inc.y, inc.z ) );\n    delta /= steps;\n\n    // Jitter\n\n    // Nice little seed from\n    // https://blog.demofox.org/2020/05/25/casual-shadertoy-path-tracing-1-basic-camera-diffuse-emissive/\n    uint seed = uint( gl_FragCoord.x ) * uint( 1973 ) + uint( gl_FragCoord.y ) * uint( 9277 ) + uint( frame ) * uint( 26699 );\n    vec3 size = vec3( textureSize( map, 0 ) );\n    float randNum = randomFloat( seed ) * 2.0 - 1.0;\n    p += rayDir * randNum * ( 1.0 / size );\n\n    //\n\n    vec4 ac = vec4( base, 0.0 );\n\n    for ( float t = bounds.x; t < bounds.y; t += delta ) {\n\n        float d = sample1( p + 0.5 );\n\n        d = smoothstep( threshold - range, threshold + range, d ) * opacity;\n\n        float col = shading( p + 0.5 ) * 3.0 + ( ( p.x + p.y ) * 0.25 ) + 0.2;\n\n        ac.rgb += ( 1.0 - ac.a ) * d * col;\n\n        ac.a += ( 1.0 - ac.a ) * d;\n\n        if ( ac.a >= 0.95 ) break;\n\n        p += rayDir * delta;\n\n    }\n\n    color = ac;\n\n    if ( color.a == 0.0 ) discard;\n\n}",side:R,transparent:!0}),Rt(this,"isExampleVolumeMaterial_Cloud",!0),this.opacity=o}get map(){return this.uniforms.map.value}set map(t){this.uniforms.map.value=t,this.uniformsNeedUpdate=!0}get color(){return this.uniforms.base.value}set color(t){this.uniforms.base.value=t,this.uniformsNeedUpdate=!0}get opacity(){return this.uniforms.opacity.value}set opacity(t){this.uniforms&&(this.uniforms.opacity.value=t),this.uniformsNeedUpdate=!0}get threshold(){return this.uniforms.threshold.value}set threshold(t){this.uniforms.threshold.value=t,this.uniformsNeedUpdate=!0}get range(){return this.uniforms.range.value}set range(t){this.uniforms.range.value=t,this.uniformsNeedUpdate=!0}get steps(){return this.uniforms.steps.value}set steps(t){this.uniforms.steps.value=t,this.uniformsNeedUpdate=!0}get frame(){return this.uniforms.frame.value}set frame(t){this.uniforms.frame.value=t,this.uniformsNeedUpdate=!0}get cameraPos(){return this.uniforms.cameraPos.value}set cameraPos(t){this.uniforms.cameraPos.value.copy(t),this.uniformsNeedUpdate=!0}}class we extends V{constructor(t){let{threshold:e,map:n,steps:a,cameraPos:o}=t??{};e=e??.25,a=a??100,o=o??new i,super({glslVersion:x,uniforms:{map:{value:n},cameraPos:{value:o},threshold:{value:e},steps:{value:a}},vertexShader:_e,fragmentShader:"precision highp float;\nprecision highp sampler3D;\n\nuniform mat4 modelViewMatrix;\nuniform mat4 projectionMatrix;\n\nin vec3 vOrigin;\nin vec3 vDirection;\n\nout vec4 color;\n\nuniform sampler3D map;\n\nuniform float threshold;\nuniform float steps;\n\nvec2 hitBox( vec3 orig, vec3 dir ) {\n    const vec3 box_min = vec3( - 0.5 );\n    const vec3 box_max = vec3( 0.5 );\n    vec3 inv_dir = 1.0 / dir;\n    vec3 tmin_tmp = ( box_min - orig ) * inv_dir;\n    vec3 tmax_tmp = ( box_max - orig ) * inv_dir;\n    vec3 tmin = min( tmin_tmp, tmax_tmp );\n    vec3 tmax = max( tmin_tmp, tmax_tmp );\n    float t0 = max( tmin.x, max( tmin.y, tmin.z ) );\n    float t1 = min( tmax.x, min( tmax.y, tmax.z ) );\n    return vec2( t0, t1 );\n}\n\nfloat sample1( vec3 p ) {\n    return texture( map, p ).r;\n}\n\n#define epsilon .0001\n\nvec3 normal( vec3 coord ) {\n    if ( coord.x < epsilon ) return vec3( 1.0, 0.0, 0.0 );\n    if ( coord.y < epsilon ) return vec3( 0.0, 1.0, 0.0 );\n    if ( coord.z < epsilon ) return vec3( 0.0, 0.0, 1.0 );\n    if ( coord.x > 1.0 - epsilon ) return vec3( - 1.0, 0.0, 0.0 );\n    if ( coord.y > 1.0 - epsilon ) return vec3( 0.0, - 1.0, 0.0 );\n    if ( coord.z > 1.0 - epsilon ) return vec3( 0.0, 0.0, - 1.0 );\n\n    float step = 0.01;\n    float x = sample1( coord + vec3( - step, 0.0, 0.0 ) ) - sample1( coord + vec3( step, 0.0, 0.0 ) );\n    float y = sample1( coord + vec3( 0.0, - step, 0.0 ) ) - sample1( coord + vec3( 0.0, step, 0.0 ) );\n    float z = sample1( coord + vec3( 0.0, 0.0, - step ) ) - sample1( coord + vec3( 0.0, 0.0, step ) );\n\n    return normalize( vec3( x, y, z ) );\n}\n\nvoid main(){\n\n    vec3 rayDir = normalize( vDirection );\n    vec2 bounds = hitBox( vOrigin, rayDir );\n\n    if ( bounds.x > bounds.y ) discard;\n\n    bounds.x = max( bounds.x, 0.0 );\n\n    vec3 p = vOrigin + bounds.x * rayDir;\n    vec3 inc = 1.0 / abs( rayDir );\n    float delta = min( inc.x, min( inc.y, inc.z ) );\n    delta /= steps;\n\n    for ( float t = bounds.x; t < bounds.y; t += delta ) {\n\n        float d = sample1( p + 0.5 );\n\n        if ( d > threshold ) {\n\n            color.rgb = normal( p + 0.5 ) * 0.5 + ( p * 1.5 + 0.25 );\n            color.a = 1.;\n            break;\n\n        }\n\n        p += rayDir * delta;\n\n    }\n\n    if ( color.a == 0.0 ) discard;\n\n}",side:R,transparent:!0}),Rt(this,"isExampleVolumeMaterial_Perlin",!0)}get map(){return this.uniforms.map.value}set map(t){this.uniforms.map.value=t,this.uniformsNeedUpdate=!0}get threshold(){return this.uniforms.threshold.value}set threshold(t){this.uniforms.threshold.value=t,this.uniformsNeedUpdate=!0}get steps(){return this.uniforms.steps.value}set steps(t){this.uniforms.steps.value=t,this.uniformsNeedUpdate=!0}get cameraPos(){return this.uniforms.cameraPos.value}set cameraPos(t){this.uniforms.cameraPos.value.copy(t),this.uniformsNeedUpdate=!0}}class Me{constructor(t){this.ssp=t,this.defaultColorGradient=[[0,"rgba(0,255,0,0)"],[.5,"rgba(64,255,255,0.5)"],[1,"rgba(255,64,255,1)"]],this.defaultGradientVolumeMaterialOptions={fit:Vt.Raw,accFactor:2,depthTest:!1,side:R,discardOut:!1,gradient:Ct(this.defaultColorGradient)}}_createTexture(t,e,n=!1){const i=he(Object.assign({points:t},e),Object.assign({maxSize:60},e)),{options:a,position:o,scale:r}=i,s=kt(n?me(a):se(a),Object.assign({},e)),l=new Lt(Object.assign(Object.assign(Object.assign({},this.defaultGradientVolumeMaterialOptions),{map:s}),e)),c=new jt(l);return c.position.copy(o),c.scale.set(r,r,r),this.ssp.addObject(c),c}createHeatCloud(t,e){return this._createTexture(t,e)}createLineHeat(t,e){return this._createTexture(t,e,!0)}async createImageExtrusion(t,e){var n;const i=new Ht(Object.assign({side:R,depthTest:!1},e)),a=Wt(await L.utils.imageLoader.loadAsync(t),null!==(n=null==e?void 0:e.depth)&&void 0!==n?n:200,e);i.map=a;const o=new jt(i);return this.ssp.addObject(o),o}createSliceMesh(t,e){let n=t.map;if(!n)return;if(t instanceof Lt){const{gradient:i,range:a,discardOut:o,voidRange:r}=t;n=Yt(n.image,Object.assign({gradient:null==i?void 0:i.image,range:a,discardOut:o,voidRange:r},e))}const i=new ne(Object.assign({map:n},e));return this.ssp.addObject(i),i}createImageSlice(t,e){const n=t.map;if(!n)return;let i=null;if(t instanceof Lt){const{gradient:a,range:o,discardOut:r,voidRange:s}=t;i=Yt(n.image,Object.assign({gradient:null==a?void 0:a.image,range:o,discardOut:r,voidRange:s},e))}if(t instanceof Ht&&(i=n),!i)return;return new ie(i)}}export{Jt as ColorFogPointsMaterial,xe as ExampleVolumeMaterial_1,ye as ExampleVolumeMaterial_Cloud,we as ExampleVolumeMaterial_Perlin,Zt as FogPointsMaterial,Bt as GradientData3DTexture,$t as GradientFogPointsMaterial,Lt as GradientVolumeMaterial,Me as HeatMapPlugin,Mt as ImageData3DSlice,Kt as ImageData3DTexture,ie as ImageData3DTextureSlice,Ht as ImageVolumeMaterial,ee as SliceMaterial,ne as SliceMesh,Xt as SphereFogMaterial,Vt as VolumeFit,Et as VolumeMaterial,jt as VolumeMesh,ge as VolumeRenderStyle,re as computeHeatData3DForPoint,Qt as createColorFogPoints,kt as createGradientData3DTexture,te as createGradientFogPoints,se as createHeatData3D,Wt as createImageData3DTextureByExtrudeImage,Yt as createImageData3DTextureFromGradient,me as createLineData3D,W as createLinearGradientImageData,Ct as createLinearGradientTexture,ht as data2DCoordToIndex,vt as data2DIndexToCoord,at as data3DCoordToIndex,ot as data3DIndexToCoord,Me as default,pt as extrudeImage,mt as extrudeImageData,ft as extrudeImageDataOnAxis,xt as getData2DItem,gt as getData2DItemSafe,wt as getData2DRow,yt as getData2DValue,_t as getData2DValueSafe,st as getData3DItem,rt as getData3DItemSafe,dt as getData3DSlice,ct as getData3DValue,lt as getData3DValueSafe,Q as getImageDateColor,Y as getImageDateOfImage,J as getImageDateRow,X as getNaturalSizeOfImageSource,qt as gradientData3DToImageData3D,ce as hollowGradient_Default,$ as imageDataToCanvas,tt as imageDataToUrl,et as isIImageData2,nt as isIImageData3,it as isOutOfSize,ue as lineValueGradient_Default,Gt as makeOriginOnBoundingBoxMinOfGeometry,ut as mergeData3Ds,de as numberValuesAccumulate_Default,he as optimizeGradientOptions,le as radiusGradient_Default,Z as reverseImageDateY,pe as transformGradientOptions,fe as translateScaleGradientOptions,ae as valueGradient_Default,oe as valuesAccumulate_Default};
+import { TextureLoader, Texture, Matrix3, Vector3, Vector4, Vector2, AlphaFormat, RedFormat, RedIntegerFormat, RGFormat, RGIntegerFormat, RGBAFormat, RGBAIntegerFormat, LuminanceFormat, LuminanceAlphaFormat, DepthFormat, DepthStencilFormat, ShaderMaterial, GLSL3, Matrix4, DoubleSide, Mesh, BoxGeometry, Data3DTexture, FloatType, LinearFilter, UnsignedByteType, Color, FrontSide, UniformsUtils, ShaderLib, Points, PlaneGeometry, Box3, Line3, Quaternion, BackSide, RawShaderMaterial, BufferGeometry, Float32BufferAttribute, Uint8ClampedBufferAttribute } from 'three';
+import SoonSpace from 'soonspacejs';
+
+var f$1 = /* @__PURE__ */ ((n) => (n[n.x = 0] = "x", n[n.y = 1] = "y", n[n.z = 2] = "z", n))(f$1 || {});
+((n) => {
+  function t(r) {
+    return n[r];
+  }
+  n.toKey = t;
+  function e(r) {
+    return n[r];
+  }
+  n.toIndex = e;
+  function o(r) {
+    const c = (r + 1) % 3, i = (r + 2) % 3;
+    return [c, i];
+  }
+  n.getCrossAxiss = o;
+})(f$1 || (f$1 = {}));
+var s$1 = /* @__PURE__ */ ((n) => (n[n.x = 0] = "x", n[n.y = 1] = "y", n[n.z = 2] = "z", n[n.w = 3] = "w", n))(s$1 || {});
+((n) => {
+  function t(r) {
+    return n[r];
+  }
+  n.toKey = t;
+  function e(r) {
+    return n[r];
+  }
+  n.toIndex = e;
+  function o(r) {
+    const c = (r + 1) % 4, i = (r + 2) % 4;
+    return [c, i];
+  }
+  n.getCrossAxiss = o;
+})(s$1 || (s$1 = {}));
+
+var f = /* @__PURE__ */ ((n) => (n[n.x = 0] = "x", n[n.y = 1] = "y", n[n.z = 2] = "z", n))(f || {});
+((n) => {
+  function t(r) {
+    return n[r];
+  }
+  n.toKey = t;
+  function e(r) {
+    return n[r];
+  }
+  n.toIndex = e;
+  function o(r) {
+    const c = (r + 1) % 3, i = (r + 2) % 3;
+    return [c, i];
+  }
+  n.getCrossAxiss = o;
+})(f || (f = {}));
+var s = /* @__PURE__ */ ((n) => (n[n.x = 0] = "x", n[n.y = 1] = "y", n[n.z = 2] = "z", n[n.w = 3] = "w", n))(s || {});
+((n) => {
+  function t(r) {
+    return n[r];
+  }
+  n.toKey = t;
+  function e(r) {
+    return n[r];
+  }
+  n.toIndex = e;
+  function o(r) {
+    const c = (r + 1) % 4, i = (r + 2) % 4;
+    return [c, i];
+  }
+  n.getCrossAxiss = o;
+})(s || (s = {}));
+
+function B$1(e, t, n, a = "srgb") {
+  const i = new OffscreenCanvas(t, n).getContext("2d"), o = i.createLinearGradient(0, 0, 256, 1);
+  for (const [c, r] of e)
+    o.addColorStop(Number(c), r);
+  return i.fillStyle = o, i.fillRect(0, 0, t, n), i.getImageData(0, 0, t, n, { colorSpace: a });
+}
+function _$1(e, t) {
+  const { x: n, y: a } = t;
+  return e.z * a * n + e.y * n + e.x;
+}
+function $$1(e, t, n = 1) {
+  let { x: a, y: s, z: i } = t;
+  const o = e.size;
+  return a = Math.max(0, Math.min(o.x - 1, a)), s = Math.max(0, Math.min(o.y - 1, s)), i = Math.max(0, Math.min(o.z - 1, i)), U$1(e, { x: a, y: s, z: i }, n);
+}
+function U$1(e, t, n = 1) {
+  const { data: a, size: s } = e, i = _$1(t, s);
+  let o = [];
+  if (n > 0) {
+    const c = i * n;
+    for (let r = 0; r < n; r++)
+      o.push(a[c + r]);
+  }
+  return { index: i, value: o };
+}
+function E$1(e, t, n, a = 1) {
+  const { size: s } = e, i = f.toKey(t), [o, c] = f.getCrossAxiss(t), r = f.toKey(o), u = f.toKey(c), h = s[r], g = s[u], l = [], f$1 = s[i];
+  n = Math.trunc(n), n = Math.max(0, Math.min(f$1 - 1, n));
+  let x = { x: 0, y: 0, z: 0 };
+  x[i] = n;
+  for (let m = 0; m < g; m++) {
+    x[u] = m;
+    for (let d = 0; d < h; d++) {
+      x[r] = d;
+      const { value: y } = U$1(e, x, a);
+      l.push(...y);
+    }
+  }
+  return { data: l, size: { x: h, y: g } };
+}
+
+var zt = Object.defineProperty;
+var Ut$1 = (t, n, e) => n in t ? zt(t, n, { enumerable: !0, configurable: !0, writable: !0, value: e }) : t[n] = e;
+var j$1 = (t, n, e) => (Ut$1(t, typeof n != "symbol" ? n + "" : n, e), e);
+function Pt$1(t) {
+  return t.w !== void 0 ? [Vector4, 4] : t.z !== void 0 ? [Vector3, 3] : [Vector2, 2];
+}
+function En(t) {
+  const n = t.length, e = [];
+  if (n === 0)
+    return e;
+  const [o, s] = Pt$1(t[0]);
+  for (let c = 0; c < n; c++) {
+    const r = t[c], i = new o().copy(r), l = r.radius, u = e[c] ?? (e[c] = []);
+    for (let a = c + 1; a < n; a++) {
+      const f = t[a];
+      i.distanceToSquared(f) <= (f.radius + l) ** 2 && (u.push(a), (e[a] = []).push(c));
+    }
+  }
+  return e;
+}
+let Tt$1;
+function Gn(t) {
+  if (typeof t == "string")
+    return (Tt$1 || (Tt$1 = new TextureLoader())).load(t);
+  if (t.isTexture)
+    return t;
+  const n = Array.isArray(t) ? B$1(t, 256, 1) : t, e = new Texture(n);
+  return n.colorSpace === "srgb" && (e.colorSpace = "srgb"), e.needsUpdate = !0, e;
+}
+function an(t, n) {
+  const e = Math.sqrt(t.lengthSq() * n.lengthSq());
+  if (e === 0)
+    return 0;
+  let o = t.dot(n) / e;
+  return o = Math.max(-1, Math.min(1, o)), Math.acos(o);
+}
+function J$1(t, n, e) {
+  let o = an(t, n);
+  return o === 0 ? o : t.clone().cross(n).dot(e) < 0 ? -o : o;
+}
+function fn$1(t, n, e) {
+  const o = n.clone().negate(), s = e.clone().projectOnPlane(n), c = t.clone().projectOnPlane(s), r = J$1(n, c, s), i = n.clone().cross(s), l = t.clone().projectOnPlane(i), u = J$1(n, l, i), a = t.clone().projectOnPlane(o), f = J$1(s, a, o);
+  return { yaw: r, pitch: u, roll: f };
+}
+const it$1 = 180 / Math.PI, ut = {
+  yaw: [
+    { name: "前", range: [-15, 15] },
+    { name: "左", range: [15, 165] },
+    { name: "右", range: [-165, -15] },
+    { name: "后", range: [-180.1, -165] },
+    { name: "后", range: [165, 180.1] }
+  ],
+  pitch: [
+    { name: "前", range: [-15, 15] },
+    { name: "上", range: [15, 165] },
+    { name: "下", range: [-165, -15] },
+    { name: "后", range: [-180.1, -165] },
+    { name: "后", range: [165, 180.1] }
+  ],
+  roll: [
+    { name: "上", range: [-15, 15] },
+    { name: "左", range: [15, 165] },
+    { name: "右", range: [-165, -15] },
+    { name: "下", range: [-180.1, -165] },
+    { name: "下", range: [165, 180.1] }
+  ]
+};
+function pn$1(t) {
+  const n = {};
+  for (const [e, o] of Object.entries(t))
+    n[e] = Array.isArray(o) ? o : Object.entries(o).map(([s, c]) => ({ name: s, range: c }));
+  return n;
+}
+const lt = {
+  degrees: !0,
+  map: ut,
+  front: { x: 0, y: 0, z: 1 },
+  up: { x: 0, y: 1, z: 0 }
+};
+let It$1 = class It {
+  constructor(n) {
+    j$1(this, "_options");
+    j$1(this, "_listMap");
+    j$1(this, "_front");
+    j$1(this, "_up");
+    n && (this.options = n);
+  }
+  static get options() {
+    return this._options ?? (this.options = lt);
+  }
+  static set options(n) {
+    this._options = Object.assign({}, structuredClone(ut), n);
+  }
+  /**
+   * 默认选项
+   */
+  get defaultOptions() {
+    return this.constructor.options;
+  }
+  get options() {
+    return this._options ?? (this.options = this.defaultOptions);
+  }
+  set options(n) {
+    this._options = Object.assign({}, structuredClone(this.defaultOptions), n), this._listMap = null, this._front = null, this._up = null;
+  }
+  /**
+   * 是否使用角度，而非弧度
+   */
+  get degrees() {
+    return this.options.degrees ?? (this.options.degrees = this.defaultOptions.degrees ?? !0);
+  }
+  set degrees(n) {
+    this.options.degrees = n;
+  }
+  get map() {
+    return this.options.map || (this.map = this.defaultOptions.map), this.options.map;
+  }
+  set map(n) {
+    const e = structuredClone(this.defaultOptions.map), o = structuredClone(ut);
+    n = n ? {
+      yaw: n.yaw ?? e.yaw ?? o.yaw,
+      pitch: n.pitch ?? e.pitch ?? o.pitch,
+      roll: n.roll ?? e.roll ?? o.roll
+    } : e, this.options.map = n, this._listMap = null;
+  }
+  get listMap() {
+    return this._listMap ?? (this._listMap = pn$1(this.map ?? {}));
+  }
+  get front() {
+    return this._front || (this.front = this.options.front ?? this.defaultOptions.front ?? lt.front), this._front;
+  }
+  set front(n) {
+    this._front = new Vector3().copy(n);
+  }
+  get up() {
+    return this._up || (this.up = this.options.up ?? this.defaultOptions.up ?? lt.up), this._up;
+  }
+  set up(n) {
+    this._up = new Vector3().copy(n);
+  }
+  /**
+   * 计算方位
+   * @param target 
+   * @param front 
+   * @param up 
+   */
+  computeAzimuth(n, e, o) {
+    const s = new Vector3().copy(n), c = new Vector3().copy(e ?? this.front), r = new Vector3().copy(o ?? this.up);
+    let { yaw: i, pitch: l, roll: u } = fn$1(s, c, r);
+    return this.degrees && (i *= it$1, l *= it$1, u *= it$1), {
+      yaw: {
+        angle: i,
+        name: this.findAzimuthNames("yaw", i)
+      },
+      pitch: {
+        angle: l,
+        name: this.findAzimuthNames("pitch", l)
+      },
+      roll: {
+        angle: u,
+        name: this.findAzimuthNames("roll", u)
+      }
+    };
+  }
+  /**
+   * 查找匹配的方位名字
+   * @param type - 类型
+   * @param angle - 角度
+   * @returns 匹配的名字列表
+   */
+  findAzimuthNames(n, e) {
+    const o = [], s = this.listMap[n];
+    if (!s)
+      return o;
+    for (const { name: c, range: [r, i] } of s)
+      (e >= r && e < i || e <= r && e > i) && o.push(c);
+    return o;
+  }
+};
+j$1(It$1, "_options");
+new It$1(); const mn$1 = {
+  [AlphaFormat]: 1,
+  [RedFormat]: 1,
+  [RedIntegerFormat]: 1,
+  [RGFormat]: 3,
+  [RGIntegerFormat]: 2,
+  [RGBAFormat]: 4,
+  [RGBAIntegerFormat]: 4,
+  [LuminanceFormat]: 1,
+  [LuminanceAlphaFormat]: 2,
+  [DepthFormat]: 1,
+  [DepthStencilFormat]: 2
+};
+function Mt(t) {
+  return mn$1[t];
+}
+function Yn(t, n, e) {
+  const { data: o, width: s, height: c, depth: r } = t.image, i = Mt(t.format);
+  return E$1({ data: o, size: { x: s, y: c, z: r } }, n, e, i);
+}
+function Zn(t, n) {
+  const { data: e, width: o, height: s, depth: c } = t.image, r = Mt(t.format);
+  return $$1({ data: e, size: { x: o, y: s, z: c } }, n, r);
+}
+new Matrix3();
+
+var T = Object.defineProperty;
+var R = (e, t, n) => t in e ? T(e, t, { enumerable: !0, configurable: !0, writable: !0, value: n }) : e[t] = n;
+var D = (e, t, n) => (R(e, typeof t != "symbol" ? t + "" : t, n), n);
+function B(e, t, n, a = "srgb") {
+  const i = new OffscreenCanvas(t, n).getContext("2d"), o = i.createLinearGradient(0, 0, 256, 1);
+  for (const [c, r] of e)
+    o.addColorStop(Number(c), r);
+  return i.fillStyle = o, i.fillRect(0, 0, t, n), i.getImageData(0, 0, t, n, { colorSpace: a });
+}
+function K(e, t) {
+  let { area: n, width: a, height: s, colorSpace: i } = t || {};
+  i = i ?? "srgb", n = n || {};
+  const o = n.x ?? 0, c = n.y ?? 0, r = b(e), u = n.width ?? r.width, h = n.height ?? r.height;
+  a = a ?? u, s = s ?? h;
+  const l = new OffscreenCanvas(a, s).getContext("2d");
+  return l.drawImage(e, o, c, u, h, 0, 0, a, s), l.getImageData(0, 0, a, s, { colorSpace: i });
+}
+function b(e) {
+  let t = 0, n = 0;
+  return e instanceof HTMLImageElement ? (t = e.naturalWidth, n = e.naturalHeight) : e instanceof HTMLVideoElement ? (t = e.videoWidth, n = e.videoHeight) : (t = e.width, n = e.height), { width: t, height: n };
+}
+function V(e) {
+  const { data: t, width: n, height: a, colorSpace: s } = e, i = t.length, o = new Uint8ClampedArray(i);
+  for (let c = 0; c < a; c++) {
+    const r = c * n, u = -(r + n), h = u + n, g = t.slice(u * 4, h * 4);
+    o.set(g, r * 4);
+  }
+  return new ImageData(o, n, a, { colorSpace: s });
+}
+function F(e, t) {
+  const { data: n, width: a } = e, s = t * a, i = s + a;
+  return n.slice(s * 4, i * 4);
+}
+function H(e, t, n) {
+  const { data: a, width: s } = e, o = (t * s + n) * 4;
+  return a.slice(o, o + 4);
+}
+function Y(e) {
+  const t = document.createElement("canvas"), n = t.getContext("2d");
+  return t.width = e.width, t.height = e.height, n.putImageData(e, 0, 0), t;
+}
+function J(e) {
+  return Y(e).toDataURL();
+}
+function L(e) {
+  return e && e.data != null && e.width != null && e.height != null;
+}
+function P(e) {
+  return L(e) && e.depth != null;
+}
+function Q(e, t) {
+  for (const [n, a] of Object.entries(e))
+    if (0 <= a && a <= t[n])
+      return !0;
+  return !1;
+}
+function _(e, t) {
+  const { x: n, y: a } = t;
+  return e.z * a * n + e.y * n + e.x;
+}
+function Z(e, t) {
+  const { x: n, y: a } = t, s = Math.trunc(e / n), i = e - s * n, o = Math.trunc(s / a), c = s - o * a;
+  return { x: i, y: c, z: o };
+}
+function $(e, t, n = 1) {
+  let { x: a, y: s, z: i } = t;
+  const o = e.size;
+  return a = Math.max(0, Math.min(o.x - 1, a)), s = Math.max(0, Math.min(o.y - 1, s)), i = Math.max(0, Math.min(o.z - 1, i)), U(e, { x: a, y: s, z: i }, n);
+}
+function U(e, t, n = 1) {
+  const { data: a, size: s } = e, i = _(t, s);
+  let o = [];
+  if (n > 0) {
+    const c = i * n;
+    for (let r = 0; r < n; r++)
+      o.push(a[c + r]);
+  }
+  return { index: i, value: o };
+}
+function tt(e, t) {
+  let { x: n, y: a, z: s } = t;
+  const i = e.size;
+  return n = Math.max(0, Math.min(i.x - 1, n)), a = Math.max(0, Math.min(i.y - 1, a)), s = Math.max(0, Math.min(i.z - 1, s)), N(e, { x: n, y: a, z: s });
+}
+function N(e, t) {
+  const { data: n, size: a } = e, s = _(t, a);
+  return { index: s, value: n[s] };
+}
+function et(e, t) {
+  const { voidValue: n, valuesAccumulate: a, verifyVoid: s = () => !1 } = t;
+  let { x: i, y: o, z: c } = e[0].size;
+  const r = e.length;
+  for (let l = 1; l < r; l++) {
+    const f = e[l].size;
+    i = Math.max(i, f.x), o = Math.max(o, f.y), c = Math.max(c, f.z);
+  }
+  const u = i * o * c, h = new Array(u), g = i * o;
+  for (let l = 0; l < c; l++) {
+    const f = l * g;
+    for (let x = 0; x < o; x++) {
+      const m = f + x * i;
+      for (let d = 0; d < i; d++) {
+        const y = m + d, v = { x: d, y: x, z: l }, M = [];
+        for (const I of e) {
+          const p = I.size;
+          if (d < p.x && x < p.y && l < p.z)
+            continue;
+          const { index: C, value: A } = N(I, v);
+          s(A, I, C) || M.push({ value: A, index: C, data3D: I, coord: v });
+        }
+        const S = M.length;
+        let w = n;
+        S > 1 ? w = a(M) : S === 1 && (w = M[0].value), h[y] = w;
+      }
+    }
+  }
+  return {
+    data: h,
+    size: { x: i, y: o, z: c }
+  };
+}
+function E(e, t, n, a = 1) {
+  const { size: s } = e, i = f$1.toKey(t), [o, c] = f$1.getCrossAxiss(t), r = f$1.toKey(o), u = f$1.toKey(c), h = s[r], g = s[u], l = [], f = s[i];
+  n = Math.trunc(n), n = Math.max(0, Math.min(f - 1, n));
+  let x = { x: 0, y: 0, z: 0 };
+  x[i] = n;
+  for (let m = 0; m < g; m++) {
+    x[u] = m;
+    for (let d = 0; d < h; d++) {
+      x[r] = d;
+      const { value: y } = U(e, x, a);
+      l.push(...y);
+    }
+  }
+  return { data: l, size: { x: h, y: g } };
+}
+function W(e, t, n) {
+  t = Math.trunc(t);
+  const { data: a, width: s, height: i, colorSpace: o } = n ? V(e) : e, c = a.length, r = c * t, u = new Uint8ClampedArray(r);
+  for (let h = 0; h < t; h++)
+    u.set(a, h * c);
+  return { data: u, width: s, height: i, depth: t, colorSpace: o };
+}
+function nt(e, t, n) {
+  const a = n || {}, s = a.reverseY, i = a.axis ?? "z", o = L(e) ? e : K(e, a);
+  return i === "z" ? W(o, t, s) : X(o, i, t, s);
+}
+function X(e, t, n, a) {
+  n = Math.trunc(n);
+  const { data: s, width: i, height: o, colorSpace: c } = a ? V(e) : e, r = f$1.toKey(t), u = (t + 1) % 3, h = (t + 2) % 3, g = f$1.toKey(u), l = f$1.toKey(h), f = {
+    [g]: i,
+    [l]: o,
+    [r]: n
+  }, x = new Uint8ClampedArray(s.length * n), m = { x: 0, y: 0, z: 0 };
+  for (let d = 0; d < n; d++) {
+    m[r] = d;
+    for (let y = 0; y < o; y++) {
+      m[l] = y;
+      for (let v = 0; v < i; v++) {
+        m[g] = v;
+        const S = (y * i + v) * 4, w = S + 4, I = s.subarray(S, w), C = _(m, f) * 4;
+        x.set(I, C);
+      }
+    }
+  }
+  return { data: x, width: f.x, height: f.y, depth: f.z, colorSpace: c };
+}
+function O(e, t) {
+  return e.y * t.x + e.x;
+}
+function at$1(e, t) {
+  const { x: n } = t, a = Math.trunc(e / n);
+  return { x: e - a * n, y: a };
+}
+function st(e, t, n = 1) {
+  let { x: a, y: s } = t;
+  const i = e.size;
+  return a = Math.max(0, Math.min(i.x - 1, a)), s = Math.max(0, Math.min(i.y - 1, s)), G(e, { x: a, y: s }, n);
+}
+function G(e, t, n = 1) {
+  const { data: a, size: s } = e, i = O(t, s);
+  let o = [];
+  if (n > 0) {
+    const c = i * n;
+    for (let r = 0; r < n; r++)
+      o.push(a[c + r]);
+  }
+  return { index: i, value: o };
+}
+function it(e, t) {
+  let { x: n, y: a } = t;
+  const s = e.size;
+  return n = Math.max(0, Math.min(s.x - 1, n)), a = Math.max(0, Math.min(s.y - 1, a)), j(e, { x: n, y: a });
+}
+function j(e, t) {
+  const { data: n, size: a } = e, s = O(t, a);
+  return { index: s, value: n[s] };
+}
+function ot$1(e, t, n = 1) {
+  const { data: a, size: s } = e, i = s.x, o = t * i, c = o + i;
+  return (typeof a.slice == "function" ? a : Array.from(a)).slice(o * n, c * n);
+}
+class ct {
+  constructor(t) {
+    D(this, "isImageData3DSlice", !0);
+    D(this, "_image3D");
+    D(this, "_canvas");
+    D(this, "_context");
+    D(this, "_axis", f$1.z);
+    D(this, "_depth", 0);
+    D(this, "_sliceSize", null);
+    D(this, "_slice");
+    this.image3D = t;
+  }
+  /**
+   * 3D图像数据
+   */
+  get image3D() {
+    return this._image3D;
+  }
+  set image3D(t) {
+    this._image3D = t, this.updateSize();
+  }
+  /**
+   * 呈现切片的 canvas 元素
+   * @remarks
+   * 当 {@link ImageData3DSlice.axis} 或 {@link ImageData3DSlice.depth} 变更时，它会自动更新切片
+   */
+  get canvas() {
+    let t = this._canvas;
+    return t || (t = this._canvas = document.createElement("canvas"), t.width = this.sliceSize.x, t.height = this.sliceSize.y), t;
+  }
+  /**
+   * canvas 的 2d 上下文对象
+   */
+  get context() {
+    return this._context ?? (this._context = this.canvas.getContext("2d"));
+  }
+  /**
+   * 轴
+   * @remarks
+   * 会沿该轴的垂直截面生成切片
+   * 
+   * @defaultValue Axis.z
+   */
+  get axis() {
+    return this._axis;
+  }
+  set axis(t) {
+    this._axis = t, this.updateSize();
+  }
+  /**
+   * 轴上的坐标值
+   * @remarks
+   * 会在轴的该位置处获取切片
+   * @defaultValue 0
+   */
+  get depth() {
+    return this._depth;
+  }
+  set depth(t) {
+    this._depth = t, this.updateSlice();
+  }
+  /**
+   * 切片的尺寸
+   * @remarks
+   * 不同轴上的切片的尺寸一般不一样
+   */
+  get sliceSize() {
+    let t = this._sliceSize;
+    if (!t) {
+      const { axis: n, image3D: a } = this, { width: s, height: i, depth: o } = a, [c, r] = f$1.getCrossAxiss(n), u = [s, i, o];
+      this._sliceSize = t = {
+        x: u[c],
+        y: u[r]
+      };
+    }
+    return t;
+  }
+  /**
+   * 获取切片数据
+   * @remarks
+   * 总是获取当前轴和深度处的切片
+   */
+  get slice() {
+    return this._slice ?? (this._slice = this.getSlice(this.axis, this.depth));
+  }
+  /**
+   * 获取切片数据
+   * @remarks
+   * 可以指定轴和深度；该方法不会更新  {@link ImageData3DSlice.axis} 和 {@link ImageData3DSlice.depth} 属性 
+   * @param axis - 轴
+   * @param depth - 深度
+   * @returns 
+   */
+  getSlice(t, n) {
+    const { data: a, width: s, height: i, depth: o, colorSpace: c } = this.image3D, { data: r, size: u } = E({ data: a, size: { x: s, y: i, z: o } }, t, n, 4), h = r instanceof Uint8ClampedArray ? r : Uint8ClampedArray.from(r);
+    return new ImageData(h, u.x, u.y, { colorSpace: c });
+  }
+  /**
+   * 获取切片图像的 url
+   * @returns 
+   */
+  getDataURL() {
+    return this.canvas.toDataURL();
+  }
+  /**
+   * 获取切片上指定坐标处的颜色
+   * @param coord - 坐标
+   * @returns 
+   */
+  getColor(t) {
+    const { x: n, y: a } = t;
+    return H(this.slice, a, n);
+  }
+  /**
+   * 更新切片尺寸相关的信息
+   */
+  updateSize() {
+    this._sliceSize = null;
+    const t = this.canvas;
+    t.width = this.sliceSize.x, t.height = this.sliceSize.y, this.updateSlice();
+  }
+  /**
+   * 更新切片相关的信息
+   */
+  updateSlice() {
+    this._slice = null, this.context.putImageData(this.slice, 0, 0);
+  }
+}
+
+var at = Object.defineProperty;
+var ot = (o, s, e) => s in o ? at(o, s, { enumerable: !0, configurable: !0, writable: !0, value: e }) : o[s] = e;
+var v = (o, s, e) => (ot(o, typeof s != "symbol" ? s + "" : s, e), e);
+const St = `precision highp sampler3D;
+
+uniform sampler3D map;
+uniform vec3 containerMin;
+uniform vec3 containerMax;
+
+// 0: 充满；1: 对齐；2: 原始
+uniform int fit;
+
+
+out vec3 cameraPos;
+out vec3 lookDir;
+out vec3 displayMin;
+out vec3 displayMax;
+
+
+
+void main() {
+    cameraPos = vec3( inverse( modelMatrix ) * vec4( cameraPosition, 1.0 ) ).xyz;
+    lookDir = position - cameraPos;
+
+    vec3 containerSize = containerMax - containerMin;
+    vec3 origin = containerMin;
+    vec3 textSize =  vec3(textureSize(map,0));
+    displayMin = containerMin;
+    displayMax = containerMax;
+
+    switch (fit) {
+        case 1:
+            displayMax = displayMin + textSize;
+            break;
+        case 2:
+            origin = vec3(0.0);
+            break;
+        default:
+            textSize = containerSize;
+    }
+
+    vec3 scale = 1.0/textSize;
+    cameraPos = (cameraPos - origin)*scale;
+    lookDir *= scale;
+
+    displayMin = max((displayMin - origin)*scale,vec3(0.0));
+    displayMax = min((displayMax - origin)*scale,vec3(1.0));
+    
+    gl_Position = projectionMatrix * modelViewMatrix * vec4( position, 1.0 );
+}`;
+var Dt = /* @__PURE__ */ ((o) => (o[o.Fill = 0] = "Fill", o[o.Align = 1] = "Align", o[o.Raw = 2] = "Raw", o))(Dt || {});
+class Xe extends ShaderMaterial {
+  constructor(e) {
+    const { map: t, opacity: n, accFactor: i, steps: a, alphaRange: r, fit: l, atomize: u, side: c, containerMin: p, containerMax: m, uniforms: f, ...h } = e ?? {}, _ = n ?? 1, M = i ?? 1, S = a ?? 100, D = new Vector2().copy(r ?? { x: 0, y: 0.95 }), w = l ?? 0, x = u ?? !0, y = c ?? FrontSide, g = new Vector3().copy(p ?? { x: 0, y: 0, z: 0 }), z = new Vector3();
+    if (m)
+      z.copy(m);
+    else if (t) {
+      const { width: U, height: T, depth: C } = t.image, V = new Vector3(U - 1, T - 1, C - 1);
+      z.addVectors(g, V);
+    }
+    const N = {
+      map: { value: t },
+      containerMin: { value: g },
+      containerMax: { value: z },
+      fit: { value: w },
+      opacity: { value: _ },
+      alphaRange: { value: D },
+      accFactor: { value: M },
+      steps: { value: S },
+      isDoubleSide: { value: !1 },
+      atomize: { value: x },
+      ...f
+    };
+    super({
+      glslVersion: GLSL3,
+      // @ts-ignore
+      uniforms: N,
+      transparent: !0,
+      vertexShader: St,
+      ...h
+    });
+    v(this, "isVolumeMaterial", !0);
+    v(this, "_side", FrontSide);
+    this.opacity = _, this.side = y;
+  }
+  /**
+   * 三维的纹理
+   */
+  get map() {
+    return this.uniforms.map.value;
+  }
+  set map(e) {
+    this.uniforms.map.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 渲染体积材质的容器的最小点
+   * @remarks
+   * 容器的最小点和最大点一般是被设置成 geometry 的AABB包围盒；但也可以不是
+   * 总之，体积材质的渲染依托于 geometry 的形状，它不会超出 geometry 的尺寸；
+   * 但可以控制容器的最小点和最大点来控制体积材质的渲染尺寸；
+   */
+  get containerMin() {
+    return this.uniforms.containerMin.value;
+  }
+  set containerMin(e) {
+    this.uniforms.containerMin.value.copy(e), this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 渲染体积材质的容器的最大点
+   * @remarks
+   * 容器的最小点和最大点一般是被设置成 geometry 的AABB包围盒；但也可以不是
+   * 总之，体积材质的渲染依托于 geometry 的形状，它不会超出 geometry 的尺寸；
+   * 但可以控制容器的最小点和最大点来控制体积材质的渲染尺寸；
+   */
+  get containerMax() {
+    return this.uniforms.containerMax.value;
+  }
+  set containerMax(e) {
+    this.uniforms.containerMax.value.copy(e), this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 渲染体积材质容器的尺寸
+   * @remarks
+   * 每次调用都会重新计算
+   */
+  get containerSize() {
+    return this.containerMax.clone().sub(this.containerMin);
+  }
+  /**
+   * 体积材质在容器内的填充模式
+   * @defaultValue VolumeFit.Fill
+   */
+  get fit() {
+    return this.uniforms.fit.value;
+  }
+  set fit(e) {
+    this.uniforms.fit.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 获取当前填模式下对 map 的应用的平移
+   * @returns
+   */
+  getFitTranslate() {
+    return this.fit === 2 ? new Vector3() : this.containerMin.clone();
+  }
+  /**
+   * 获取当前填模式下对 map 的应用的缩放
+   * @returns 
+   */
+  getFitScale() {
+    var a;
+    const e = new Vector3(1, 1, 1), t = (a = this.map) == null ? void 0 : a.image, n = this.fit;
+    if (!t || n === 2 || n === 1)
+      return e;
+    const i = new Vector3(t.width, t.height, t.depth);
+    return e.copy(this.containerSize).divide(i), e;
+  }
+  /**
+   * 获取当前填模式下对 map 的应用的变换矩阵
+   * @returns 
+   */
+  getFitMatrix() {
+    var i;
+    const e = new Matrix4(), t = (i = this.map) == null ? void 0 : i.image, n = this.fit;
+    if (!t || n === 2)
+      return e;
+    if (n === 0) {
+      const a = new Vector3(t.width, t.height, t.depth), l = this.containerSize.divide(a);
+      e.makeScale(l.x, l.y, l.z);
+    }
+    return e.setPosition(this.containerMin), e;
+  }
+  /**
+   * 获取当前填模式下对 map 的应用的变换矩阵的逆矩阵
+   * @returns 
+   */
+  getFitMatrixInvert() {
+    return this.getFitMatrix().invert();
+  }
+  /**
+   * Mesh三角形的渲染面
+   * @remarks
+   * 可设置为前面、后面、或者两面都渲染
+   * 如果只渲染前面，进入体积材质内部，体积材质就会消失
+   * 如果只渲染后面，体积初遮挡，就会消失
+   * 如果两面都渲染，则会隐藏看到容器的面
+   */
+  // @ts-ignore
+  get side() {
+    return this._side;
+  }
+  set side(e) {
+    this._side = e, this.uniforms && (this.uniforms.isDoubleSide.value = e === DoubleSide, this.uniformsNeedUpdate = !0);
+  }
+  /**
+   * 透明度
+   * @remarks
+   * 决定整体的透明度
+   */
+  // @ts-ignore
+  get opacity() {
+    return this.uniforms.opacity.value;
+  }
+  set opacity(e) {
+    this.uniforms && (this.uniforms.opacity.value = e), this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 有效透明度的取值范围
+   * @remarks
+   * 如果最终颜色的透明度小于或等于 alphaRange.x ，则被认为是完全透明
+   * 如果最终颜色的透明度在于或等于 alphaRange.y ，则被认为是完全不透明
+   * 取值范围应当在 0 - 1
+   * @defaultValue {x:0,y:0.95}
+   */
+  get alphaRange() {
+    return this.uniforms.alphaRange.value;
+  }
+  set alphaRange(e) {
+    this.uniforms.alphaRange.value.copy(e), this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 是否开启雾化的效果
+   * @remarks
+   * 开启后的渲染效果会更像雾
+   * @defaultValue true
+   */
+  get atomize() {
+    return this.uniforms.atomize.value;
+  }
+  set atomize(e) {
+    this.uniforms.atomize.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 体积材质渲染的采样数
+   * @remarks
+   * 采样数越高，占用的GPU资源就越大；所以该数据设置的适宜最好；
+   * @defaultValue 100
+   */
+  get steps() {
+    return this.uniforms.steps.value;
+  }
+  set steps(e) {
+    this.uniforms.steps.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 颜色累积系数
+   * @remarks
+   * 它是对体积材质进行积分时使用的系数；只在开启雾化效果{@link VolumeMaterial.atomize}时才有效；
+   * 值越小，最终呈现出的效果就越雾化；
+   * @defaultValue 1
+   */
+  get accFactor() {
+    return this.uniforms.accFactor.value;
+  }
+  set accFactor(e) {
+    this.uniforms && (this.uniforms.accFactor.value = e), this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 将材质空间下的坐标转为 map 空间下的坐标
+   * @param coord - 材质空间下的坐标
+   * @returns 
+   */
+  toMapPosition(e) {
+    const t = new Vector3(e.x, e.y, e.z).applyMatrix4(this.getFitMatrixInvert());
+    return t.set(Math.trunc(t.x), Math.trunc(t.y), Math.trunc(t.z)), t;
+  }
+  /**
+   * 将材质空间下的深度转为 map 空间下的深度
+   * @param axis - map坐标系下的坐标轴
+   * @param depth - 材质空间下的深度
+   * @returns 
+   */
+  toMapDepth(e, t) {
+    const n = f$1.toKey(e), i = new Vector3();
+    i[n] = 1;
+    const a = new Matrix3().setFromMatrix4(this.getFitMatrix());
+    return i.applyMatrix3(a), t /= i.length(), Math.trunc(t);
+  }
+  /**
+   * 获取3D数据的切片
+   * @remarks
+   * 切片就是指定轴的指定位置的垂直截面
+   * @param axis - map坐标系下的坐标轴
+   * @param depth - 材质空间下，在轴方向的坐标值
+   * @returns 
+   */
+  getData3DSlice(e, t) {
+    const n = this.map;
+    return n != null && n.image ? (t = this.toMapDepth(e, t), Yn(n, e, t)) : null;
+  }
+  /**
+   * 获取3D数据中指定坐标处的数据项目
+   * @param coord 
+   * @returns 
+   */
+  getItem(e) {
+    const t = this.map;
+    return t != null && t.image ? Zn(t, this.toMapPosition(e)).value : null;
+  }
+}
+function bt(o) {
+  o.boundingBox || o.computeBoundingBox();
+  const { x: s, y: e, z: t } = o.boundingBox.min;
+  return o.translate(-s, -e, -t), o;
+}
+class rn extends Mesh {
+  constructor(e) {
+    var r;
+    const { width: t, height: n, depth: i } = ((r = e.map) == null ? void 0 : r.image) || {}, a = new BoxGeometry(t, n, i);
+    bt(a);
+    super(a, e);
+    v(this, "isVolumeMesh", !0);
+    // @ts-ignore
+    // override get geometry(){
+    //     return this._geometry;
+    // }
+    // override set geometry(value){
+    //     this._geometry = value;
+    //     this.configMaterial();
+    // }
+    v(this, "_geometry");
+    /**
+     * 是否自动更新材质
+     * @remarks
+     * 当启动该选项后，在 VolumeMesh 监测到 geometry 变更时，会自动更新体积材质的相关选项，以适应新的 geometry
+     * @defaultValue true
+     */
+    v(this, "autoUpdateMaterial", !0);
+    // @ts-ignore
+    // override get material(){
+    //     return this._material;
+    // }
+    // override set material(value){
+    //     this._material = value;
+    //     this.configGeometry();
+    // }
+    v(this, "_material");
+    /**
+     * 是否自动更新几何体
+     * @remarks
+     * 当启动该选项后，在 VolumeMesh 监测到 material 变更时，会自动更新 geometry 的相关选项，以适应新的 体积材质
+     * @defaultValue true
+     */
+    v(this, "autoUpdateGeometry", !0);
+    /**
+     * 自动规范化
+     * @remarks
+     * 开启该选项，设置新的 geometry 时，对 geometry 自动执行规范化操作 
+     */
+    v(this, "autoNormalize", !1);
+    Object.defineProperties(this, {
+      geometry: {
+        get: () => this._geometry,
+        set: (l) => {
+          this._geometry = l, this.autoNormalize ? this.normalize() : this.autoUpdateMaterial && this.updateMaterial();
+        }
+      },
+      material: {
+        get: () => this._material,
+        set: (l) => {
+          this._material = l, this.autoUpdateGeometry && this.updateGeometry();
+        }
+      }
+    }), this.geometry = a, this.material = e;
+  }
+  /**
+   * 更新材质
+   * @remarks
+   * 更新体积材质的相关选项以使其匹配 geometry
+   */
+  updateMaterial() {
+    const { geometry: e, material: t } = this;
+    if (!(t && e))
+      return !1;
+    e.boundingBox || e.computeBoundingBox();
+    const { min: n, max: i } = e.boundingBox;
+    t.containerMin = n, t.containerMax = i;
+  }
+  /**
+   * 更新几何体
+   * @remarks
+   * 更新 geometry 的相关选项以使其匹配体积材质
+   */
+  updateGeometry() {
+    var c;
+    const { geometry: e, material: t, autoUpdateMaterial: n } = this, i = (c = t.map) == null ? void 0 : c.image;
+    if (!i)
+      return !1;
+    e.boundingBox || e.computeBoundingBox();
+    const { width: a, height: r, depth: l } = i, u = e.boundingBox.getSize(new Vector3());
+    e.scale(a / u.x, r / u.y, l / u.z), n && this.updateMaterial();
+  }
+  /**
+   * 规范化
+   * @remarks
+   * 将 geometry 的原点移动到包围盒的最小点处
+   */
+  normalize() {
+    const { geometry: e, autoUpdateMaterial: t } = this;
+    e.boundingBox || e.computeBoundingBox();
+    const { x: n, y: i, z: a } = e.boundingBox.min;
+    e.translate(-n, -i, -a), t && this.updateMaterial();
+  }
+  /**
+   * 将材质空间下的坐标转为 map 空间下的坐标
+   * @param coord - 材质空间下的坐标
+   * @returns 
+   */
+  toMapPosition(e) {
+    const t = this.worldToLocal(new Vector3(e.x, e.y, e.z));
+    return this.material.toMapPosition(t);
+  }
+  /**
+   * 将世界坐标系下的深度转为 material 空间下的深度
+   * @param axis - 材质空间下的坐标轴
+   * @param depth - 深度
+   * @returns 
+   */
+  toMaterialDepth(e, t) {
+    const n = f$1.toKey(e), i = new Vector3();
+    i[n] = 1;
+    const a = new Matrix3().setFromMatrix4(this.matrixWorld);
+    return i.applyMatrix3(a), t /= i.length(), t;
+  }
+  /**
+   * 将世界坐标系下的深度转为 map 空间下的深度
+   * @param axis - map空间下的坐标轴
+   * @param depth - 深度
+   * @returns 
+   */
+  toMapDepth(e, t) {
+    return t = this.toMaterialDepth(e, t), this.material.toMapDepth(e, t);
+  }
+  /**
+   * 获取3D数据的切片
+   * @remarks
+   * 切片就是指定轴的指定位置的垂直截面；
+   * 轴和深度都是局部坐标系下的
+   * @param axis - map空间下的坐标轴
+   * @param depth - 世界空间下在轴方向上的坐标值
+   * @returns 
+   */
+  getData3DSlice(e, t) {
+    return t = this.toMaterialDepth(e, t), this.material.getData3DSlice(e, t);
+  }
+  /**
+   * 获取3D数据中指定坐标处的数据项目
+   * @param coord - 世界坐标系下的坐标
+   * @returns 
+   */
+  getItem(e) {
+    const t = this.worldToLocal(new Vector3(e.x, e.y, e.z));
+    return this.material.getItem(t);
+  }
+}
+const Ct = `precision highp float;
+precision highp sampler3D;
+
+uniform sampler3D map;
+uniform sampler2D gradient;
+
+uniform bool atomize;
+uniform float steps;
+// 颜色累积系数
+uniform float accFactor;
+uniform vec2 range;
+// 是否丢弃超出范围的像素
+uniform bool discardOut;
+// 空值范围
+uniform vec2 voidRange;
+uniform float opacity;
+uniform vec2 alphaRange;
+uniform bool isDoubleSide;
+
+in vec3 cameraPos;
+in vec3 lookDir;
+in vec3 displayMin;
+in vec3 displayMax;
+out vec4 fragColor;
+
+vec4 colorBlend(vec4 near,vec4 far){
+    float nA = near.a;
+    if  (nA >= 1.0){
+        return near;
+    }
+    if  (nA <= 0.0){
+        return far;
+    }
+    
+    float fA = far.a;
+    float a = fA + nA - fA*nA;
+    vec3 color = (far.rgb * fA * (1.0 - nA) + near.rgb*nA)/a;
+    return vec4(color,a);
+}
+
+
+
+
+vec2 intersectBox( vec3 orig, vec3 dir ) {
+    vec3 inv_dir = 1.0 / dir;
+    vec3 tmin_tmp = (displayMin - orig ) * inv_dir;
+    vec3 tmax_tmp = (displayMax - orig ) * inv_dir;
+    vec3 tmin = min( tmin_tmp, tmax_tmp );
+    vec3 tmax = max( tmin_tmp, tmax_tmp );
+    float t0 = max( tmin.x, max( tmin.y, tmin.z ) );
+    float t1 = min( tmax.x, min( tmax.y, tmax.z ) );
+    return vec2( t0, t1 );
+}
+
+float getValue( vec3 point ) {
+    return texture( map, point ).r;
+}
+
+// 获取梯度颜色
+vec4 getGradientColor(float val) {
+    val = clamp( val, 0.0,1.0);
+    return texture2D(gradient, vec2(val, 0.5));
+}
+
+
+
+
+void main(){
+    vec3 rayDir = normalize( lookDir );
+
+    vec2 times = intersectBox( cameraPos, rayDir );
+    float tMin = times.x;
+    float tMax = times.y;
+    if ( tMin > tMax || tMax < 0.0 ) discard;
+    tMin = max( tMin, 0.0 );
+
+    vec3 point = cameraPos + tMin * rayDir;
+    
+    float step = 1.0/steps;
+    float opacityFactor = atomize ? step * accFactor : 1.0;
+
+
+    float invRangeLen = 1.0/(range.y - range.x);
+    vec3 stepDir = rayDir * step;
+
+    vec4 finalColor = vec4( 0,0,0,0);
+    float alphaMax = alphaRange.y;
+
+    for ( float t = tMin; t <= tMax; t += step,point += stepDir ) {
+        float val = getValue( point);
+        if (voidRange.x <= val && val <= voidRange.y ) continue;
+        val = (val - range.x) * invRangeLen;
+        if (discardOut && (val < 0.0 || val > 1.0)) continue;
+
+        vec4 gradientColor = getGradientColor(val);
+        gradientColor.a *= opacityFactor;
+        finalColor = colorBlend(finalColor,gradientColor);
+        if ( finalColor.a >= alphaMax ) break;
+    }
+
+    finalColor.a *= opacity;
+    if ( finalColor.a <= alphaRange.x ) discard;
+    if (isDoubleSide){
+        finalColor.a = 1.0 - sqrt(1.0 - finalColor.a);
+    }
+    fragColor = finalColor;
+
+}`;
+class sn extends Xe {
+  constructor(e) {
+    const { gradient: t, range: n, discardOut: i, voidRange: a, ...r } = e ?? {}, l = new Vector2().copy(n ?? { x: 0, y: 100 }), u = i ?? !0, c = new Vector2().copy(a ?? { x: -100, y: -1 }), m = {
+      gradient: { value: t == null ? t : Gn(t) },
+      range: { value: l },
+      discardOut: { value: u },
+      voidRange: { value: c }
+    };
+    super({ ...r, fragmentShader: Ct, uniforms: m });
+    v(this, "isGradientVolumeMaterial", !0);
+  }
+  /**
+   * 数值的梯度映射纹理
+   */
+  get gradient() {
+    return this.uniforms.gradient.value;
+  }
+  set gradient(e) {
+    this.uniforms.gradient.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 数值的映射区间
+   * @remarks
+   * x 为最小值，y 为最大值
+   * 
+   * @defaultValue {x:0,y:100}
+   */
+  get range() {
+    return this.uniforms.range.value;
+  }
+  set range(e) {
+    this.uniforms.range.value.copy(e), this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 是否丢弃超出映射范围的数值
+   * @remarks
+   * 当数值超出映射区间 {@link GradientVolumeMaterial.range}  时，是否丢弃；
+   * 
+   * 包含左右边界值
+   * @defaultValue true
+   */
+  get discardOut() {
+    return this.uniforms.discardOut.value;
+  }
+  set discardOut(e) {
+    this.uniforms.discardOut.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 空值的范围
+   * @remarks
+   * 当数值在此范围中时，会被认为是空的值，渲染时会被丢弃
+   * 
+   * 包含左右边界值
+   * 
+   * @defaultValue {x:-100,y:-1}
+   */
+  get voidRange() {
+    return this.uniforms.voidRange.value;
+  }
+  set voidRange(e) {
+    this.uniforms.voidRange.value.copy(e), this.uniformsNeedUpdate = !0;
+  }
+}
+class Pt extends Data3DTexture {
+  constructor(e, t, n, i, a) {
+    super(e, t, n, i);
+    v(this, "isGradientData3DTexture", !0);
+    this.format = RedFormat, this.type = a ?? FloatType, this.minFilter = this.magFilter = LinearFilter, this.unpackAlignment = 1, this.needsUpdate = !0;
+  }
+  /**
+   * 获取3D数据的切片
+   * @remarks
+   * 切片就是指定轴的指定位置的垂直截面
+   * @param axis - 轴
+   * @param depth - 在轴上的坐标值
+   * @returns 
+   */
+  getData3DSlice(e, t) {
+    const { data: n, width: i, height: a, depth: r } = this.image;
+    return E({ data: n, size: { x: i, y: a, z: r } }, e, t, 1);
+  }
+  /**
+   * 获取3D数据中指定坐标处的数据项目
+   * @remarks
+   * 在 ImageData3DTexture 中，与 {@link ImageData3DTexture.getColor} 返回的值一样；
+   * @param coord 
+   * @returns 
+   */
+  getItem(e) {
+    const { data: t, width: n, height: i, depth: a } = this.image;
+    return $({ data: t, size: { x: n, y: i, z: a } }, e).value;
+  }
+  /**
+   * 获取3D数据中指定坐标处的值
+   * @param coord - 坐标
+   * @returns 
+   */
+  getValue(e) {
+    const { data: t, width: n, height: i, depth: a } = this.image;
+    return tt({ data: t, size: { x: n, y: i, z: a } }, e).value;
+  }
+}
+function ln(o, s) {
+  const { data: e, size: t } = o, { x: n, y: i, z: a } = t, r = n * i * a, { voidValue: l = 0, uint8: u } = s ?? {};
+  let c;
+  if (u) {
+    c = new Uint8Array(r);
+    for (let m = 0; m < r; m++) {
+      const f = e[m] ?? l;
+      c[m] = Math.trunc(f);
+    }
+  } else {
+    c = new Float32Array(r);
+    for (let m = 0; m < r; m++) {
+      const f = e[m] ?? l;
+      c[m] = f;
+    }
+  }
+  const p = u ? UnsignedByteType : FloatType;
+  return new Pt(c, n, i, a, p);
+}
+const Nt = `precision highp float;
+precision highp sampler3D;
+
+uniform sampler3D map;
+
+
+
+// 颜色累积系数
+uniform float accFactor;
+uniform bool isDoubleSide;
+uniform bool atomize;
+
+
+uniform float opacity;
+uniform vec2 alphaRange;
+uniform float steps;
+
+in vec3 cameraPos;
+in vec3 lookDir;
+in vec3 displayMin;
+in vec3 displayMax;
+out vec4 fragColor;
+
+// 混合颜色
+vec4 colorBlend(vec4 near,vec4 far){
+    float nA = near.a;
+    if  (nA >= 1.0){
+        return near;
+    }
+    if  (nA <= 0.0){
+        return far;
+    }
+    
+    float fA = far.a;
+    float a = fA + nA - fA*nA;
+    vec3 color = (far.rgb * fA * (1.0 - nA) + near.rgb*nA)/a;
+    return vec4(color,a);
+}
+
+
+
+
+// 包围盒求交
+vec2 intersectBox( vec3 orig, vec3 dir ) {
+    vec3 inv_dir = 1.0 / dir;
+    vec3 tmin_tmp = (displayMin - orig ) * inv_dir;
+    vec3 tmax_tmp = (displayMax - orig ) * inv_dir;
+    vec3 tmin = min( tmin_tmp, tmax_tmp );
+    vec3 tmax = max( tmin_tmp, tmax_tmp );
+    float t0 = max( tmin.x, max( tmin.y, tmin.z ) );
+    float t1 = min( tmax.x, min( tmax.y, tmax.z ) );
+    return vec2( t0, t1 );
+}
+
+
+
+void main(){
+    vec3 rayDir = normalize( lookDir );
+
+    vec2 times = intersectBox( cameraPos, rayDir );
+    float tMin = times.x;
+    float tMax = times.y;
+    if ( tMin > tMax || tMax < 0.0 ) discard;
+    tMin = max( tMin, 0.0 );
+
+    vec3 point = cameraPos + tMin * rayDir;
+    
+    float step = 1.0/steps;
+    float opacityFactor = atomize ? step * accFactor : 1.0;
+
+    vec3 stepDir = rayDir * step;
+
+    vec4 finalColor = vec4( 0,0,0,0);
+    float alphaMax = alphaRange.y;
+
+    for ( float t = tMin; t <= tMax; t += step,point += stepDir ) {
+        vec4 textureColor = texture( map, point );
+        textureColor.a *= opacityFactor;
+        finalColor = colorBlend(finalColor,textureColor);
+        if ( finalColor.a >= alphaMax ) break;
+    }
+
+    finalColor.a *= opacity;
+    if ( finalColor.a <= alphaRange.x ) discard;
+    if (isDoubleSide){
+        finalColor.a = 1.0 - sqrt(1.0 - finalColor.a);
+    }
+    fragColor = finalColor;
+
+}`;
+class cn extends Xe {
+  constructor(e) {
+    super({ ...e, fragmentShader: Nt });
+    v(this, "isImageVolumeMaterial", !0);
+  }
+}
+function Ut(o, s, e) {
+  const t = e ?? {}, { x: n, y: i } = t.range ?? { x: 0, y: 100 }, a = t.discardOut ?? !0, { x: r, y: l } = t.voidRange ?? { x: -100, y: -1 }, u = t.voidColor ?? [0, 0, 0, 0], { width: c, height: p } = s, m = Math.trunc(p / 2), f = c - 1;
+  function h(x) {
+    const y = Math.trunc(f * x);
+    return H(s, m, y);
+  }
+  const _ = 1 / (i - n), { data: M, size: S } = o, D = M.length, w = [];
+  for (let x = 0; x < D; x++) {
+    let y = u, g = M[x];
+    (g <= r || l <= g) && (g = (g - n) * _, a && (g < 0 || g > 1) || (g = Math.min(Math.max(g, 0), 1), y = h(g))), w.push(...y);
+  }
+  return { data: w, size: S };
+}
+class $e extends Data3DTexture {
+  constructor(e, t, n, i) {
+    super(e, t, n, i);
+    v(this, "isImageData3DTexture", !0);
+    this.format = RGBAFormat, this.type = UnsignedByteType, this.minFilter = this.magFilter = LinearFilter, this.unpackAlignment = 1, this.needsUpdate = !0;
+  }
+  /**
+   * 获取3D数据的切片
+   * @remarks
+   * 切片就是指定轴的指定位置的垂直截面
+   * @param axis - 轴
+   * @param depth - 在轴上的坐标值
+   * @returns 
+   */
+  getData3DSlice(e, t) {
+    const { data: n, width: i, height: a, depth: r } = this.image;
+    return E({ data: n, size: { x: i, y: a, z: r } }, e, t, 4);
+  }
+  /**
+   * 获取3D数据中指定坐标处的数据项目
+   * @remarks
+   * 在 ImageData3DTexture 中，与 {@link ImageData3DTexture.getColor} 返回的值一样；
+   * @param coord 
+   * @returns 
+   */
+  getItem(e) {
+    return this.getColor(e);
+  }
+  /**
+   * 获取3D图像中指定坐标处的颜色
+   * @param coord - 坐标
+   * @returns 
+   */
+  getColor(e) {
+    const { data: t, width: n, height: i, depth: a } = this.image;
+    return $({ data: t, size: { x: n, y: i, z: a } }, e, 4).value;
+  }
+}
+function un(o, s, e) {
+  const t = (e == null ? void 0 : e.reverseY) ?? !0, { data: n, width: i, height: a, depth: r, colorSpace: l } = nt(o, s, { ...e, reverseY: t }), u = new $e(n, i, a, r);
+  return l === "srgb" && (u.colorSpace = "srgb"), u;
+}
+function mn(o, s) {
+  const e = s.gradient, t = Array.isArray(e) ? B(e, 256, 1) : L(e) ? e : K(e);
+  let n = o;
+  P(o) && (n = {
+    // @ts-ignore
+    data: o.data,
+    size: { x: o.width, y: o.height, z: o.depth }
+  });
+  const { data: i, size: a } = Ut(n, t, s), r = Uint8ClampedArray.from(i);
+  return new $e(r, a.x, a.y, a.z);
+}
+const At = `out vec3 lookDir;
+out vec3 normalDir;
+
+#if defined( USE_COLOR_ALPHA ) || defined( USE_COLOR )
+    out vec4 verColor;
+#endif
+
+void main() {
+    vec3 cameraPos = vec3( inverse( modelMatrix ) * vec4( cameraPosition, 1.0 ) ).xyz;
+    lookDir = position - cameraPos;
+    normalDir = normal;
+
+    #if defined( USE_COLOR_ALPHA )
+        verColor = color;
+    #elif defined( USE_INSTANCING_COLOR )
+        overColor = vec3(instanceColor,1.0);
+    #elif defined( USE_COLOR )
+        overColor = vec3(color,1.0);
+    #endif
+    gl_Position = projectionMatrix * modelViewMatrix * vec4( position, 1.0 );
+}`, Ft = `uniform vec3 color;
+uniform float solid;
+uniform float exp;
+
+uniform float opacity;
+uniform bool isDoubleSide;
+
+in vec3 lookDir;
+in vec3 normalDir;
+
+#if defined( USE_COLOR_ALPHA ) || defined( USE_COLOR )
+    in  vec4 verColor;
+#endif
+
+out vec4 fragColor;
+
+
+void main(){
+    #if defined( USE_COLOR_ALPHA ) || defined( USE_COLOR )
+        fragColor = verColor;
+    #else
+        fragColor = vec4(color,1.0);
+    #endif
+
+    vec3 dir = normalize(lookDir);
+    vec3 nor = normalize(normalDir);
+    float depth = abs(dot(nor,dir));
+    float dist = sqrt(1.0 - pow(depth,2.0));
+
+    float finalAlpha = fragColor.a * depth;
+    if (dist > solid) {
+        float distAlpha = 1.0 - (dist - solid)/(1.0 - solid);
+        finalAlpha *= pow(distAlpha,exp);
+    }
+    
+    finalAlpha *= opacity;
+    if ( finalAlpha == 0.0 ) discard;
+    if (isDoubleSide){
+        finalAlpha = 1.0 - sqrt(1.0 - finalAlpha);
+    }
+    fragColor.a = finalAlpha;
+}`;
+class dn extends ShaderMaterial {
+  constructor(e) {
+    const t = e ?? {}, n = t.color ?? new Color(1, 1, 1), i = t.solid ?? 0, a = t.exp ?? 1, r = t.opacity ?? 1, l = t.side ?? FrontSide;
+    super({
+      glslVersion: GLSL3,
+      // @ts-ignore
+      uniforms: {
+        color: { value: n },
+        solid: { value: i },
+        exp: { value: a },
+        opacity: { value: r },
+        isDoubleSide: { value: !1 }
+      },
+      vertexShader: At,
+      fragmentShader: Ft,
+      side: FrontSide,
+      transparent: !0
+    });
+    v(this, "isSphereFogMaterial", !0);
+    v(this, "_side", FrontSide);
+    this.opacity = r, this.side = l;
+  }
+  /**
+   * Mesh三角形的渲染面
+   * @remarks
+   * 可设置为前面、后面、或者两面都渲染
+   * 如果只渲染前面，进入体积材质内部，体积材质就会消失
+   * 如果只渲染后面，体积初遮挡，就会消失
+   * 如果两面都渲染，则会隐藏看到容器的面
+   */
+  // @ts-ignore
+  get side() {
+    return this._side;
+  }
+  set side(e) {
+    this._side = e, this.uniforms && (this.uniforms.isDoubleSide.value = e === DoubleSide, this.uniformsNeedUpdate = !0);
+  }
+  /**
+   * 球形雾的颜色
+   * @remarks
+   * 当开启顶点颜色 {@link SphereFogMaterial.vertexColors} 时，此设置将无论；
+   */
+  get color() {
+    return this.uniforms.color.value;
+  }
+  set color(e) {
+    this.uniforms.color.value.copy(e), this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 透明度
+   * @remarks
+   * 决定整体的透明度
+   * @defaultValue 1
+   */
+  // @ts-ignore
+  get opacity() {
+    return this.uniforms.opacity.value;
+  }
+  set opacity(e) {
+    this.uniforms && (this.uniforms.opacity.value = e), this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 实心因子
+   * @remarks
+   * 在实心范围内，球形雾不会衰减
+   * @defaultValue 0
+   */
+  get solid() {
+    return this.uniforms.solid.value;
+  }
+  set solid(e) {
+    this.uniforms && (this.uniforms.solid.value = e), this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 雾的衰减指数
+   * @defaultValue 1
+   */
+  get exp() {
+    return this.uniforms.exp.value;
+  }
+  set exp(e) {
+    this.uniforms && (this.uniforms.exp.value = e), this.uniformsNeedUpdate = !0;
+  }
+}
+class Ye extends ShaderMaterial {
+  constructor(e, t) {
+    const { solid: n = 0, exp: i = 1, useRadius: a = !0, sizeAttenuation: r = !0, star: l = !1, ...u } = e ?? {}, { uniforms: c, ...p } = t || {}, m = {
+      star: { value: l },
+      solid: { value: n },
+      exp: { value: i }
+    };
+    super({
+      // glslVersion: GLSL3,
+      // @ts-ignore
+      uniforms: UniformsUtils.merge([ShaderLib.points.uniforms, m, c]),
+      transparent: !0,
+      depthTest: !1,
+      // depthWrite:false,
+      defines: {
+        useRadius: a
+      },
+      ...p
+    });
+    v(this, "isFogPointsMaterial", !0);
+    /**
+     * 点的大小是否因相机深度而衰减
+     * @remarks
+     * 仅限透视摄像头
+     * 
+     * @defaultValue true
+     */
+    v(this, "sizeAttenuation", !0);
+    Object.assign(this, { ...u, solid: n, exp: i, sizeAttenuation: r });
+  }
+  /**
+   * 使用来自Texture的数据设置点的颜色
+   * @remarks
+   * 可以选择包括一个alpha通道，通常与 .transparent或.alphaTest。
+   */
+  get map() {
+    return this.uniforms.map.value;
+  }
+  set map(e) {
+    this.uniforms.map.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 是否渲染成星星的样式
+   * @remarks
+   * 否则，将被渲染成球形
+   */
+  get star() {
+    return this.uniforms.star.value;
+  }
+  set star(e) {
+    this.uniforms.star.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 点的大小
+   * @defaultValue 1.0
+   */
+  get size() {
+    return this.uniforms.size.value;
+  }
+  set size(e) {
+    this.uniforms.size.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 透明度
+   * @remarks
+   * 决定整体的透明度
+   * @defaultValue 1
+   */
+  // @ts-ignore
+  get opacity() {
+    return this.uniforms.opacity.value;
+  }
+  set opacity(e) {
+    this.uniforms && (this.uniforms.opacity.value = e), this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 实心因子
+   * @remarks
+   * 在实心范围内，雾点不会衰减
+   * @defaultValue 0
+   */
+  get solid() {
+    return this.uniforms.solid.value;
+  }
+  set solid(e) {
+    this.uniforms && (this.uniforms.solid.value = e), this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 雾的衰减指数
+   * @defaultValue 1
+   */
+  get exp() {
+    return this.uniforms.exp.value;
+  }
+  set exp(e) {
+    this.uniforms && (this.uniforms.exp.value = e), this.uniformsNeedUpdate = !0;
+  }
+  get scale() {
+    return this.uniforms.scale.value;
+  }
+  set scale(e) {
+    this.uniforms && (this.uniforms.scale.value = e), this.uniformsNeedUpdate = !0;
+  }
+  get fogDensity() {
+    return this.uniforms.fogDensity.value;
+  }
+  set fogDensity(e) {
+    this.uniforms && (this.uniforms.fogDensity.value = e), this.uniformsNeedUpdate = !0;
+  }
+  get fogNear() {
+    return this.uniforms.fogNear.value;
+  }
+  set fogNear(e) {
+    this.uniforms && (this.uniforms.fogNear.value = e), this.uniformsNeedUpdate = !0;
+  }
+  get fogFar() {
+    return this.uniforms.fogFar.value;
+  }
+  set fogFar(e) {
+    this.uniforms && (this.uniforms.fogFar.value = e), this.uniformsNeedUpdate = !0;
+  }
+  get fogColor() {
+    return this.uniforms.fogColor.value;
+  }
+  set fogColor(e) {
+    this.uniforms.fogColor.value.copy(e), this.uniformsNeedUpdate = !0;
+  }
+}
+const Rt = `#ifdef useRadius
+	attribute float radius;
+#endif
+
+uniform float size;
+uniform float scale;
+
+#include <common>
+#include <color_pars_vertex>
+#include <fog_pars_vertex>
+#include <morphtarget_pars_vertex>
+#include <logdepthbuf_pars_vertex>
+#include <clipping_planes_pars_vertex>
+
+#ifdef USE_POINTS_UV
+
+	varying vec2 vUv;
+	uniform mat3 uvTransform;
+
+#endif
+
+void main() {
+
+	#ifdef USE_POINTS_UV
+
+		vUv = ( uvTransform * vec3( uv, 1 ) ).xy;
+
+	#endif
+
+	#include <color_vertex>
+	#include <morphcolor_vertex>
+	#include <begin_vertex>
+	#include <morphtarget_vertex>
+	#include <project_vertex>
+
+	#ifdef useRadius
+		gl_PointSize = radius;
+	#else
+		gl_PointSize = size;
+	#endif
+
+	#ifdef USE_SIZEATTENUATION
+
+		bool isPerspective = isPerspectiveMatrix( projectionMatrix );
+
+		if ( isPerspective ) gl_PointSize *= ( scale / - mvPosition.z );
+
+	#endif
+
+	#include <logdepthbuf_vertex>
+	#include <clipping_planes_vertex>
+	#include <worldpos_vertex>
+	#include <fog_vertex>
+
+}`, It = `uniform float solid;
+uniform float exp;
+uniform bool star;
+
+uniform vec3 diffuse;
+uniform float opacity;
+
+#include <common>
+#include <color_pars_fragment>
+#include <map_particle_pars_fragment>
+#include <alphatest_pars_fragment>
+#include <fog_pars_fragment>
+#include <logdepthbuf_pars_fragment>
+#include <clipping_planes_pars_fragment>
+
+
+vec4 getColor(vec4 originColor){
+	float finalAlpha = originColor.a;
+	float dist =  0.0;
+	if (star){
+		vec2 coord = abs(gl_PointCoord - vec2(0.5));
+		dist =  max(coord.x,coord.y);
+	}else{
+		dist = distance(gl_PointCoord,vec2(0.5));
+		if (dist>0.5) discard;
+	}
+	dist /= 0.5;
+	float depth = sqrt(1.0 - pow(dist,2.0));
+	finalAlpha *= depth;
+    if (dist > solid) {
+        float distAlpha = 1.0 - (dist - solid)/(1.0 - solid);
+        finalAlpha *= pow(distAlpha,exp);
+    }
+	originColor.a = finalAlpha;
+	return originColor;
+}
+
+void main() {
+
+	#include <clipping_planes_fragment>
+
+	vec3 outgoingLight = vec3( 0.0 );
+	vec4 diffuseColor = vec4( diffuse, opacity );
+
+	#include <logdepthbuf_fragment>
+	#include <map_particle_fragment>
+	#include <color_fragment>
+	#include <alphatest_fragment>
+
+	diffuseColor = getColor(diffuseColor);
+	outgoingLight = diffuseColor.rgb;
+
+	#include <output_fragment>
+	#include <tonemapping_fragment>
+	#include <encodings_fragment>
+	#include <fog_fragment>
+	#include <premultiplied_alpha_fragment>
+
+}`;
+class Ot extends Ye {
+  constructor(e) {
+    super(e, { vertexShader: Rt, fragmentShader: It });
+    v(this, "isColorFogPointsMaterial", !0);
+  }
+  /**
+   * 雾点的颜色
+   * @remarks
+   * 当开启顶点颜色 {@link SphereFogMaterial.vertexColors} 时，此设置将无论；
+   */
+  get color() {
+    return this.uniforms.diffuse.value;
+  }
+  set color(e) {
+    this.uniforms.diffuse.value.copy(e), this.uniformsNeedUpdate = !0;
+  }
+  get alphaMap() {
+    return this.uniforms.alphaMap.value;
+  }
+  set alphaMap(e) {
+    this.uniforms && (this.uniforms.alphaMap.value = e), this.uniformsNeedUpdate = !0;
+  }
+  get uvTransform() {
+    return this.uniforms.uvTransform.value;
+  }
+  set uvTransform(e) {
+    this.uniforms && this.uniforms.uvTransform.value.copy(e), this.uniformsNeedUpdate = !0;
+  }
+}
+function Tt(o) {
+  const { points: s, radius: e = 10, color: t = [0, 0, 0, 0] } = o, n = [], i = [], a = [];
+  for (const { x: l, y: u, z: c, radius: p, color: m } of s) {
+    n.push(l, u, c), i.push(p ?? e);
+    const f = m ?? t, h = f[3];
+    f[3] = Math.trunc(h * 255), a.push(...f);
+  }
+  const r = new BufferGeometry();
+  return r.setAttribute("position", new Float32BufferAttribute(n, 3)), r.setAttribute("radius", new Float32BufferAttribute(i, 1)), r.setAttribute("color", new Uint8ClampedBufferAttribute(a, 4, !0)), r;
+}
+function pn(o) {
+  const s = new Ot(o), e = Tt(o);
+  return new Points(e, s);
+}
+const Vt = `attribute float value;
+uniform vec2 range;
+varying float vValue;
+#ifdef useRadius
+	attribute float radius;
+#endif
+
+
+uniform float size;
+uniform float scale;
+
+#include <common>
+// #include <color_pars_vertex>
+#include <fog_pars_vertex>
+#include <morphtarget_pars_vertex>
+#include <logdepthbuf_pars_vertex>
+#include <clipping_planes_pars_vertex>
+
+#ifdef USE_POINTS_UV
+
+	varying vec2 vUv;
+	uniform mat3 uvTransform;
+
+#endif
+
+void main() {
+	vValue = (value - range.x)/(range.y - range.x);
+
+	#ifdef USE_POINTS_UV
+
+		vUv = ( uvTransform * vec3( uv, 1 ) ).xy;
+
+	#endif
+
+	// #include <color_vertex>
+	// #include <morphcolor_vertex>
+	#include <begin_vertex>
+	#include <morphtarget_vertex>
+	#include <project_vertex>
+
+	#ifdef useRadius
+		gl_PointSize = radius;
+	#else
+		gl_PointSize = size;
+	#endif
+
+	#ifdef USE_SIZEATTENUATION
+
+		bool isPerspective = isPerspectiveMatrix( projectionMatrix );
+
+		if ( isPerspective ) gl_PointSize *= ( scale / - mvPosition.z );
+
+	#endif
+
+	#include <logdepthbuf_vertex>
+	#include <clipping_planes_vertex>
+	#include <worldpos_vertex>
+	#include <fog_vertex>
+
+}`, Et = `varying float vValue;
+uniform sampler2D map;
+
+
+uniform float solid;
+uniform float exp;
+uniform bool star;
+
+// uniform vec3 diffuse;
+uniform float opacity;
+
+#include <common>
+// #include <color_pars_fragment>
+// #include <map_particle_pars_fragment>
+#include <alphatest_pars_fragment>
+#include <fog_pars_fragment>
+#include <logdepthbuf_pars_fragment>
+#include <clipping_planes_pars_fragment>
+
+
+// 获取梯度颜色
+vec4 getGradientColor(float val) {
+	val = clamp( val, 0.0,1.0);
+    return texture2D(map, vec2(val, 0.5));
+}
+
+
+vec4 getColor(){
+	float val = vValue;
+	float dist =  0.0;
+	if (star){
+		vec2 coord = abs(gl_PointCoord - vec2(0.5));
+		dist =  max(coord.x,coord.y);
+	}else{
+		dist = distance(gl_PointCoord,vec2(0.5));
+		if (dist>0.5) discard;
+	}
+	dist /= 0.5;
+	float depth = sqrt(1.0 - pow(dist,2.0));
+    if (dist > solid) {
+        float solidDist = 1.0 - (dist - solid)/(1.0 - solid);
+        val *= pow(solidDist,exp);
+    }
+	vec4 valColor = getGradientColor(val);
+	valColor.a *= depth * opacity;
+	return valColor;
+}
+
+void main() {
+
+	#include <clipping_planes_fragment>
+
+	vec3 outgoingLight = vec3( 0.0 );
+	vec4 diffuseColor = getColor();
+
+	#include <logdepthbuf_fragment>
+	// #include <map_particle_fragment>
+	// #include <color_fragment>
+	#include <alphatest_fragment>
+
+	outgoingLight = diffuseColor.rgb;
+
+	#include <output_fragment>
+	#include <tonemapping_fragment>
+	#include <encodings_fragment>
+	#include <fog_fragment>
+	#include <premultiplied_alpha_fragment>
+
+}`;
+class Gt extends Ye {
+  constructor(e) {
+    const { range: t, ...n } = e || {}, a = {
+      range: { value: new Vector2().copy(t ?? { x: 0, y: 100 }) }
+    };
+    super(n, { vertexShader: Vt, fragmentShader: Et, uniforms: a });
+    v(this, "isColorFogPointsMaterial", !0);
+  }
+  /**
+   * 数值的映射区间
+   * @remarks
+   * x 为最小值，y 为最大值
+   * 
+   * @defaultValue {x:0,y:100}
+   */
+  get range() {
+    return this.uniforms.range.value;
+  }
+  set range(e) {
+    this.uniforms.range.value.copy(e), this.uniformsNeedUpdate = !0;
+  }
+}
+function Bt(o) {
+  const { points: s, radius: e = 10, value: t = 100 } = o, n = [], i = [], a = [];
+  for (const { x: l, y: u, z: c, radius: p, value: m } of s)
+    n.push(l, u, c), i.push(p ?? e), a.push(m ?? t);
+  const r = new BufferGeometry();
+  return r.setAttribute("position", new Float32BufferAttribute(n, 3)), r.setAttribute("radius", new Float32BufferAttribute(i, 1)), r.setAttribute("value", new Float32BufferAttribute(a, 1)), r;
+}
+function fn(o) {
+  const s = o.gradient, e = s == null ? s : Gn(s), t = new Gt({ ...o, map: e }), n = Bt(o);
+  return new Points(n, t);
+}
+const Lt = `	out vec2 vUv;
+    
+	void main() {
+		vUv = uv;
+		gl_Position = projectionMatrix * modelViewMatrix * vec4( position, 1.0 );
+	}`, kt = `precision highp float;
+precision highp sampler3D;
+
+uniform sampler3D map;
+
+uniform float opacity;
+uniform int axis;
+uniform float depth;
+
+in vec2 vUv;
+out vec4 fragColor;
+
+void main(){
+    vec3 uv3 = vec3(0.0);
+    int depthSize = textureSize(map,0)[axis];
+    uv3[axis] = depth/float(depthSize);
+    int xIndex = (axis + 1)%3;
+    int yIndex = (axis + 2)%3;
+    uv3[xIndex] = vUv.x;
+    uv3[yIndex] = vUv.y;
+
+    fragColor = texture(map, uv3);
+    fragColor.a *= opacity;
+}`;
+class jt extends ShaderMaterial {
+  constructor(e) {
+    const t = e ?? {};
+    let { map: n, opacity: i, axis: a, depth: r, side: l, ...u } = t;
+    i = i ?? 1, a = a ?? f$1.z, r = r ?? 0, l = l ?? DoubleSide;
+    super({
+      ...u,
+      glslVersion: GLSL3,
+      // @ts-ignore
+      uniforms: {
+        map: { value: n },
+        opacity: { value: i },
+        axis: { value: a },
+        depth: { value: r }
+      },
+      transparent: !0,
+      vertexShader: Lt,
+      fragmentShader: kt,
+      side: l
+    });
+    v(this, "isSliceMaterial", !0);
+    v(this, "_sliceSize", null);
+    this.opacity = i;
+  }
+  /**
+   * 三维的纹理
+   */
+  get map() {
+    return this.uniforms.map.value;
+  }
+  set map(e) {
+    this.uniforms.map.value = e, this._sliceSize = null, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 透明度
+   * @remarks
+   * 决定整体的透明度
+   * @defaultValue 1
+   */
+  // @ts-ignore
+  get opacity() {
+    return this.uniforms.opacity.value;
+  }
+  set opacity(e) {
+    this.uniforms && (this.uniforms.opacity.value = e), this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 轴
+   * @remarks
+   * 切片就是垂直于该轴的截面
+   * 
+   * @defaultValue Axis.z
+   */
+  get axis() {
+    return this.uniforms.axis.value;
+  }
+  set axis(e) {
+    this.uniforms.axis.value = e, this._sliceSize = null, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 轴上的坐标值
+   * @remarks
+   * 会在轴的该位置处获取切片
+   * @defaultValue 0
+   */
+  get depth() {
+    return this.uniforms.depth.value;
+  }
+  set depth(e) {
+    this.uniforms.depth.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 切片的尺寸
+   * @remarks
+   * 不同轴上的切片的尺寸一般不一样
+   */
+  get sliceSize() {
+    let e = this._sliceSize;
+    if (!e) {
+      this._sliceSize = e = new Vector2(1, 1);
+      const { axis: t, map: n } = this;
+      if (n) {
+        const { width: i, height: a, depth: r } = n.image, [l, u] = f$1.getCrossAxiss(t), c = [i, a, r];
+        e.set(c[l], c[u]);
+      }
+    }
+    return e;
+  }
+}
+class vn extends Mesh {
+  constructor(e) {
+    const t = new jt(e), n = new PlaneGeometry();
+    super(n, t);
+    v(this, "isSliceMesh", !0);
+    // @ts-ignore
+    // override get geometry(){
+    //     return this._geometry;
+    // }
+    // override set geometry(value){
+    //     this._geometry = value;
+    //     this.configMaterial();
+    // }
+    v(this, "_geometry");
+    // @ts-ignore
+    // override get material(){
+    //     return this._material;
+    // }
+    // override set material(value){
+    //     this._material = value;
+    //     this.configGeometry();
+    // }
+    v(this, "_material");
+    /**
+     * 是否自动更新几何体
+     * @remarks
+     * 当启动该选项后，在 SliceMesh 监测到切片尺寸相关的变更时，会自动更新 geometry 的尺寸，以适应切片的尺寸
+     * @defaultValue true
+     */
+    v(this, "autoUpdateGeometry", !0);
+    Object.defineProperties(this, {
+      geometry: {
+        get: () => this._geometry,
+        set: (i) => {
+          this._geometry = i, this.autoUpdateGeometry && this.updateGeometry();
+        }
+      },
+      material: {
+        get: () => this._material,
+        set: (i) => {
+          this._material = i, this.autoUpdateGeometry && this.updateGeometry();
+        }
+      }
+    }), this.geometry = n, this.material = t;
+  }
+  /**
+   * 更新几何体
+   */
+  updateGeometry() {
+    var r;
+    const { geometry: e, material: t } = this;
+    if (!t || !((r = t.map) == null ? void 0 : r.image))
+      return !1;
+    e.boundingBox || e.computeBoundingBox();
+    const i = t.sliceSize, a = e.boundingBox.getSize(new Vector3());
+    e.scale(i.x / a.x, i.y / a.y, 1);
+  }
+  /**
+   * 三维的纹理
+   */
+  get map() {
+    return this.material.map;
+  }
+  set map(e) {
+    this.material.map = e, this.material = this.material;
+  }
+  /**
+   * 轴
+   * @remarks
+   * 切片就是垂直于该轴的截面
+   * 
+   * @defaultValue Axis.z
+   */
+  get axis() {
+    return this.material.axis;
+  }
+  set axis(e) {
+    this.material.axis = e, this.material = this.material;
+  }
+  /**
+   * 轴上的坐标值
+   * @remarks
+   * 会在轴的该位置处获取切片
+   * @defaultValue 0
+   */
+  get depth() {
+    return this.material.depth;
+  }
+  set depth(e) {
+    this.material.depth = e;
+  }
+  /**
+   * 切片的尺寸
+   * @remarks
+   * 不同轴上的切片的尺寸一般不一样
+   */
+  get sliceSize() {
+    return this.material.sliceSize;
+  }
+}
+class hn extends ct {
+  constructor(e) {
+    const t = P(e) ? e : e.image;
+    super(t);
+    v(this, "isImageData3DTextureSlice", !0);
+  }
+  /**
+   * 设置材质或 IImageData3 类数据
+   */
+  set texture(e) {
+    this.image3D = P(e) ? e : e.image;
+  }
+}
+function Ze({ ratio: o, value: s }) {
+  return s - s * o;
+}
+function Ht(o, s) {
+  const { x: e, y: t } = s;
+  let n = 0;
+  for (const { value: i, clim: { x: a, y: r } } of o)
+    n += (i - a) / (r - a);
+  return n * (t - e) + e;
+}
+function qt(o, s, e) {
+  const { size: t, data3D: n, tags: i } = s, { clim: a, radius: r, value: l, hollow: u, valueGradient: c, valuesAccumulate: p } = o, m = new Vector3().copy(o), f = [], h = [], _ = {
+    particles: f,
+    values: h
+  }, M = r * u, S = m.clone().subScalar(r), D = m.clone().addScalar(r);
+  S.max(new Vector3(0, 0, 0));
+  const w = new Vector3().copy(t).subScalar(1);
+  D.min(w);
+  let { x, y, z: g } = S;
+  x = Math.trunc(x), y = Math.trunc(y), g = Math.trunc(g);
+  let { x: z, y: N, z: U } = D;
+  if (z = Math.trunc(z), N = Math.trunc(N), U = Math.trunc(U), x > z || y > N || g > U)
+    return _;
+  const { x: T, y: C } = t, V = T * C;
+  for (let A = g; A <= U; A++) {
+    const he = A * V;
+    for (let R = y; R <= N; R++) {
+      const J = he + R * T;
+      e:
+        for (let P = x; P <= z; P++) {
+          const L = J + P;
+          if (i[L])
+            continue;
+          const F = new Vector3(P, R, A), b = F.distanceTo(m);
+          if (b > r || b < M)
+            continue;
+          i[L] = !0;
+          const k = {
+            radius: r,
+            hollow: u,
+            hollowRadius: M,
+            distance: b,
+            clim: a,
+            value: l
+          }, E = [], j = {
+            ...k,
+            ratio: b / r,
+            distance: b,
+            point: F
+          }, H = c(j);
+          E.push({
+            ...o,
+            ...j,
+            value: H
+          });
+          for (const G of e) {
+            const B = F.distanceTo(G), { hollow: ge, radius: X } = G;
+            if (B <= X) {
+              const $ = X * ge;
+              if (B < $)
+                continue e;
+              const I = {
+                ...G,
+                ratio: B / X,
+                distance: B,
+                hollowRadius: $,
+                point: F
+              }, xe = c(I);
+              E.push({ ...I, value: xe });
+            }
+          }
+          const q = p(E, a);
+          n[L] = q, f.push(F), h.push(q);
+        }
+    }
+  }
+  return _;
+}
+function gn(o) {
+  const { points: s, clim: e = { x: 0, y: 100 }, radius: t = 10, value: n = 100, hollow: i = 0, valueGradient: a = Ze, valuesAccumulate: r = Ht, size: l } = o, u = { clim: e, radius: t, value: n, hollow: i, valueGradient: a, valuesAccumulate: r }, c = s.map((x) => ({ ...u, ...x }));
+  let p = l;
+  l || (p = new Box3().setFromPoints(c).max.addScalar(1));
+  const m = En(c);
+  let { x: f, y: h, z: _ } = p;
+  f = p.x = Math.trunc(f), h = p.y = Math.trunc(h), _ = p.z = Math.trunc(_);
+  const M = new Array(f * h * _), S = new Array(f * h * _), D = { size: p, data3D: M, tags: S }, w = c.length;
+  for (let x = 0; x < w; x++) {
+    const y = c[x], g = m[x].map((z) => c[z]);
+    qt(y, D, g);
+  }
+  return {
+    data: M,
+    size: p
+  };
+}
+function Xt({ startRadius: o, addedRadius: s, ratio: e }) {
+  return s * e + o;
+}
+function $t({ startHollow: o, addedHollow: s, ratio: e }) {
+  return s * e + o;
+}
+function Yt({ startValue: o, addedValue: s, ratio: e }) {
+  return s * e + o;
+}
+function Zt(o, s) {
+  const { x: e, y: t } = s;
+  return o.reduce((n, i) => n + i - e, 0);
+}
+function xn(o) {
+  const { points: s, size: e, clim: t = { x: 0, y: 100 }, radius: n = 10, value: i = 100, hollow: a = 0, valueGradient: r = Ze, radiusGradient: l = Xt, hollowGradient: u = $t, lineValueGradient: c = Yt, valuesAccumulate: p = Zt } = o, m = { radius: n, value: i, hollow: a }, f = s.map(function(C) {
+    return Object.assign(new Vector3(), m, C);
+  }), h = e ? new Vector3().copy(e).subScalar(1) : new Box3().setFromPoints(f).max;
+  h.x = Math.trunc(h.x), h.y = Math.trunc(h.y), h.z = Math.trunc(h.z);
+  const _ = h.clone().addScalar(1), { x: M, y: S, z: D } = _, w = M * S, x = new Array(w * D), y = [], g = [], z = {
+    particles: y,
+    values: g,
+    size: _,
+    data: x
+  }, N = new Vector3(0, 0, 0), U = s.length, T = Math.trunc(U / 2);
+  for (let C = 0; C < T; C++) {
+    const V = C * 2, A = f[V], he = V + 1, R = f[he], J = new Line3(A, R), P = J.distance(), L = J.delta(new Vector3()).normalize(), F = new Quaternion();
+    F.setFromUnitVectors(new Vector3(0, 1, 0), L);
+    const b = new Matrix4().makeRotationFromQuaternion(F);
+    b.setPosition(A);
+    const { radius: k, value: E, hollow: j } = A, { radius: H, value: q, hollow: G } = R, B = H - k, ge = q - E, X = G - j, $ = { length: P, startRadius: k, endRadius: H, addedRadius: B, defaultRadius: n, radius: n, startValue: E, endValue: q, addedValue: ge, defaultValue: i, startHollow: j, endHollow: G, addedHollow: X, defaultHollow: a, clim: t }, I = Math.max(k, H), xe = new Vector3(-I, 0, -I), Ke = new Vector3(I, P, I), Se = new Box3(xe, Ke);
+    Se.applyMatrix4(b);
+    const { min: De, max: be } = Se;
+    De.max(N), be.min(h);
+    let { x: ee, y: te, z: ne } = De;
+    ee = Math.trunc(ee), te = Math.trunc(te), ne = Math.trunc(ne);
+    let { x: ie, y: ae, z: oe } = be;
+    if (ie = Math.trunc(ie), ae = Math.trunc(ae), oe = Math.trunc(oe), ee > ie || te > ae || ne > oe)
+      return z;
+    b.invert();
+    for (let re = ne; re <= oe; re++) {
+      const Qe = re * w;
+      for (let se = te; se <= ae; se++) {
+        const Je = Qe + se * M;
+        for (let le = ee; le <= ie; le++) {
+          const Ce = Je + le, ce = new Vector3(le, se, re), Pe = ce.clone();
+          Pe.applyMatrix4(b);
+          const { x: et, y: _e, z: tt } = Pe;
+          if (_e < 0 || _e > P)
+            continue;
+          const ue = Math.hypot(et, tt), me = {
+            ...$,
+            ratio: _e / P,
+            point: ce
+          }, Y = l(me);
+          me.radius = Y;
+          const Ne = u(me), Ue = Y * Ne;
+          if (ue > Y || ue < Ue)
+            continue;
+          const nt = c(me), it = {
+            ...$,
+            value: nt,
+            ratio: ue / Y,
+            radius: Y,
+            hollow: Ne,
+            hollowRadius: Ue,
+            distance: ue,
+            point: ce
+          };
+          let de = r(it);
+          const Ae = g[Ce];
+          Ae !== void 0 && (de = p([Ae, de], t)), x[Ce] = de, y.push(ce), g.push(de);
+        }
+      }
+    }
+  }
+  return z;
+}
+function _n(o, s) {
+  const e = new Vector3().setFromMatrixScale(s), t = Math.sqrt(e.lengthSq() / 3), n = new Matrix3().setFromMatrix4(s);
+  let { points: i, size: a, radius: r, ...l } = o;
+  if (a && (a = new Vector3().copy(a).applyMatrix3(n)), r && (r *= t), i) {
+    const u = [];
+    for (let { x: c, y: p, z: m, radius: f, ...h } of i) {
+      f && (f *= t);
+      const _ = new Vector3(c, p, m).applyMatrix4(s);
+      Object.assign(_, h, { radius: f }), u.push(_);
+    }
+    i = u;
+  }
+  return { points: i, size: a, radius: r, ...l };
+}
+function Wt(o, s, e = 1) {
+  let { points: t, size: n, radius: i, ...a } = o;
+  if (i && (i *= e), t) {
+    const r = [];
+    for (let { x: l, y: u, z: c, radius: p, ...m } of t) {
+      p && (p *= e);
+      const f = new Vector3(l, u, c).add(s).multiplyScalar(e);
+      Object.assign(f, m, { radius: p }), r.push(f);
+    }
+    t = r;
+  }
+  return n && (n = new Vector3().copy(n).add(s).multiplyScalar(e)), { points: t, size: n, radius: i, ...a };
+}
+function yn(o, s) {
+  const { points: e, radius: t } = o;
+  if (!e)
+    return { options: o, position: new Vector3(), scale: 1 };
+  const n = new Box3();
+  for (const p of e) {
+    const m = new Vector3().copy(p), f = p.radius ?? t;
+    f ? (n.expandByPoint(m.clone().subScalar(f)), n.expandByPoint(m.addScalar(f))) : n.expandByPoint(m);
+  }
+  const i = new Vector3().max(n.min), a = n.max.clone().sub(i).addScalar(1);
+  o.size && a.min(o.size);
+  let { scale: r, maxSize: l } = s || {};
+  r || (l = l ?? 100, r = l / Math.max(a.x, a.y, a.z)), a.multiplyScalar(r);
+  const u = i.clone().negate(), c = Wt(o, u, r);
+  return c.size = a, { options: c, position: i, scale: 1 / r };
+}
+const Fe = {
+  uniforms: {
+    u_size: { value: new Vector3(1, 1, 1) },
+    u_renderstyle: { value: 0 },
+    u_renderthreshold: { value: 0.5 },
+    u_clim: { value: new Vector2(1, 1) },
+    u_data: { value: null },
+    u_cmdata: { value: null }
+  },
+  vertexShader: (
+    /* glsl */
+    `
+
+		varying vec4 v_nearpos;
+		varying vec4 v_farpos;
+		varying vec3 v_position;
+
+		void main() {
+				// Prepare transforms to map to "camera view". See also:
+				// https://threejs.org/docs/#api/renderers/webgl/WebGLProgram
+				mat4 viewtransformf = modelViewMatrix;
+				mat4 viewtransformi = inverse(modelViewMatrix);
+
+				// Project local vertex coordinate to camera position. Then do a step
+				// backward (in cam coords) to the near clipping plane, and project back. Do
+				// the same for the far clipping plane. This gives us all the information we
+				// need to calculate the ray and truncate it to the viewing cone.
+				vec4 position4 = vec4(position, 1.0);
+				vec4 pos_in_cam = viewtransformf * position4;
+
+				// Intersection of ray and near clipping plane (z = -1 in clip coords)
+				pos_in_cam.z = -pos_in_cam.w;
+				v_nearpos = viewtransformi * pos_in_cam;
+
+				// Intersection of ray and far clipping plane (z = +1 in clip coords)
+				pos_in_cam.z = pos_in_cam.w;
+				v_farpos = viewtransformi * pos_in_cam;
+
+				// Set varyings and output pos
+				v_position = position;
+				gl_Position = projectionMatrix * viewMatrix * modelMatrix * position4;
+		}`
+  ),
+  fragmentShader: (
+    /* glsl */
+    `
+
+				precision highp float;
+				precision mediump sampler3D;
+
+				uniform vec3 u_size;
+				uniform int u_renderstyle;
+				uniform float u_renderthreshold;
+				uniform vec2 u_clim;
+
+				uniform sampler3D u_data;
+				uniform sampler2D u_cmdata;
+
+				varying vec3 v_position;
+				varying vec4 v_nearpos;
+				varying vec4 v_farpos;
+
+				// The maximum distance through our rendering volume is sqrt(3).
+				const int MAX_STEPS = 887;	// 887 for 512^3, 1774 for 1024^3
+				const int REFINEMENT_STEPS = 4;
+				const float relative_step_size = 1.0;
+				const vec4 ambient_color = vec4(0.2, 0.4, 0.2, 1.0);
+				const vec4 diffuse_color = vec4(0.8, 0.2, 0.2, 1.0);
+				const vec4 specular_color = vec4(1.0, 1.0, 1.0, 1.0);
+				const float shininess = 40.0;
+
+				void cast_mip(vec3 start_loc, vec3 step, int nsteps, vec3 view_ray);
+				void cast_iso(vec3 start_loc, vec3 step, int nsteps, vec3 view_ray);
+
+				float sample1(vec3 texcoords);
+				vec4 apply_colormap(float val);
+				vec4 add_lighting(float val, vec3 loc, vec3 step, vec3 view_ray);
+
+
+				void main() {
+						// Normalize clipping plane info
+						vec3 farpos = v_farpos.xyz / v_farpos.w;
+						vec3 nearpos = v_nearpos.xyz / v_nearpos.w;
+
+						// Calculate unit vector pointing in the view direction through this fragment.
+						vec3 view_ray = normalize(nearpos.xyz - farpos.xyz);
+
+						// Compute the (negative) distance to the front surface or near clipping plane.
+						// v_position is the back face of the cuboid, so the initial distance calculated in the dot
+						// product below is the distance from near clip plane to the back of the cuboid
+						float distance = dot(nearpos - v_position, view_ray);
+						distance = max(distance, min((-0.5 - v_position.x) / view_ray.x,
+																				(u_size.x - 0.5 - v_position.x) / view_ray.x));
+						distance = max(distance, min((-0.5 - v_position.y) / view_ray.y,
+																				(u_size.y - 0.5 - v_position.y) / view_ray.y));
+						distance = max(distance, min((-0.5 - v_position.z) / view_ray.z,
+																				(u_size.z - 0.5 - v_position.z) / view_ray.z));
+
+						// Now we have the starting position on the front surface
+						vec3 front = v_position + view_ray * distance;
+
+						// Decide how many steps to take
+						int nsteps = int(-distance / relative_step_size + 0.5);
+						if ( nsteps < 1 )
+								discard;
+
+						// Get starting location and step vector in texture coordinates
+						vec3 step = ((v_position - front) / u_size) / float(nsteps);
+						vec3 start_loc = front / u_size;
+
+						// For testing: show the number of steps. This helps to establish
+						// whether the rays are correctly oriented
+						//'gl_FragColor = vec4(0.0, float(nsteps) / 1.0 / u_size.x, 1.0, 1.0);
+						//'return;
+
+						if (u_renderstyle == 0)
+								cast_mip(start_loc, step, nsteps, view_ray);
+						else if (u_renderstyle == 1)
+								cast_iso(start_loc, step, nsteps, view_ray);
+
+						if (gl_FragColor.a < 0.05)
+								discard;
+				}
+
+
+				float sample1(vec3 texcoords) {
+						/* Sample float value from a 3D texture. Assumes intensity data. */
+						return texture(u_data, texcoords.xyz).r;
+				}
+
+
+				vec4 apply_colormap(float val) {
+						val = (val - u_clim[0]) / (u_clim[1] - u_clim[0]);
+						return texture2D(u_cmdata, vec2(val, 0.5));
+				}
+
+
+				void cast_mip(vec3 start_loc, vec3 step, int nsteps, vec3 view_ray) {
+
+						float max_val = -1e6;
+						int max_i = 100;
+						vec3 loc = start_loc;
+
+						// Enter the raycasting loop. In WebGL 1 the loop index cannot be compared with
+						// non-constant expression. So we use a hard-coded max, and an additional condition
+						// inside the loop.
+						for (int iter=0; iter<MAX_STEPS; iter++) {
+								if (iter >= nsteps)
+										break;
+								// Sample from the 3D texture
+								float val = sample1(loc);
+								// Apply MIP operation
+								if (val > max_val) {
+										max_val = val;
+										max_i = iter;
+								}
+								// Advance location deeper into the volume
+								loc += step;
+						}
+
+						// Refine location, gives crispier images
+						vec3 iloc = start_loc + step * (float(max_i) - 0.5);
+						vec3 istep = step / float(REFINEMENT_STEPS);
+						for (int i=0; i<REFINEMENT_STEPS; i++) {
+								max_val = max(max_val, sample1(iloc));
+								iloc += istep;
+						}
+
+						// Resolve final color
+						gl_FragColor = apply_colormap(max_val);
+				}
+
+
+				void cast_iso(vec3 start_loc, vec3 step, int nsteps, vec3 view_ray) {
+
+						gl_FragColor = vec4(0.0);	// init transparent
+						vec4 color3 = vec4(0.0);	// final color
+						vec3 dstep = 1.5 / u_size;	// step to sample derivative
+						vec3 loc = start_loc;
+
+						float low_threshold = u_renderthreshold - 0.02 * (u_clim[1] - u_clim[0]);
+
+						// Enter the raycasting loop. In WebGL 1 the loop index cannot be compared with
+						// non-constant expression. So we use a hard-coded max, and an additional condition
+						// inside the loop.
+						for (int iter=0; iter<MAX_STEPS; iter++) {
+								if (iter >= nsteps)
+										break;
+
+								// Sample from the 3D texture
+								float val = sample1(loc);
+
+								if (val > low_threshold) {
+										// Take the last interval in smaller steps
+										vec3 iloc = loc - 0.5 * step;
+										vec3 istep = step / float(REFINEMENT_STEPS);
+										for (int i=0; i<REFINEMENT_STEPS; i++) {
+												val = sample1(iloc);
+												if (val > u_renderthreshold) {
+														gl_FragColor = add_lighting(val, iloc, dstep, view_ray);
+														return;
+												}
+												iloc += istep;
+										}
+								}
+
+								// Advance location deeper into the volume
+								loc += step;
+						}
+				}
+
+
+				vec4 add_lighting(float val, vec3 loc, vec3 step, vec3 view_ray)
+				{
+					// Calculate color by incorporating lighting
+
+						// View direction
+						vec3 V = normalize(view_ray);
+
+						// calculate normal vector from gradient
+						vec3 N;
+						float val1, val2;
+						val1 = sample1(loc + vec3(-step[0], 0.0, 0.0));
+						val2 = sample1(loc + vec3(+step[0], 0.0, 0.0));
+						N[0] = val1 - val2;
+						val = max(max(val1, val2), val);
+						val1 = sample1(loc + vec3(0.0, -step[1], 0.0));
+						val2 = sample1(loc + vec3(0.0, +step[1], 0.0));
+						N[1] = val1 - val2;
+						val = max(max(val1, val2), val);
+						val1 = sample1(loc + vec3(0.0, 0.0, -step[2]));
+						val2 = sample1(loc + vec3(0.0, 0.0, +step[2]));
+						N[2] = val1 - val2;
+						val = max(max(val1, val2), val);
+
+						float gm = length(N); // gradient magnitude
+						N = normalize(N);
+
+						// Flip normal so it points towards viewer
+						float Nselect = float(dot(N, V) > 0.0);
+						N = (2.0 * Nselect - 1.0) * N;	// ==	Nselect * N - (1.0-Nselect)*N;
+
+						// Init colors
+						vec4 ambient_color = vec4(0.0, 0.0, 0.0, 0.0);
+						vec4 diffuse_color = vec4(0.0, 0.0, 0.0, 0.0);
+						vec4 specular_color = vec4(0.0, 0.0, 0.0, 0.0);
+
+						// note: could allow multiple lights
+						for (int i=0; i<1; i++)
+						{
+								 // Get light direction (make sure to prevent zero devision)
+								vec3 L = normalize(view_ray);	//lightDirs[i];
+								float lightEnabled = float( length(L) > 0.0 );
+								L = normalize(L + (1.0 - lightEnabled));
+
+								// Calculate lighting properties
+								float lambertTerm = clamp(dot(N, L), 0.0, 1.0);
+								vec3 H = normalize(L+V); // Halfway vector
+								float specularTerm = pow(max(dot(H, N), 0.0), shininess);
+
+								// Calculate mask
+								float mask1 = lightEnabled;
+
+								// Calculate colors
+								ambient_color +=	mask1 * ambient_color;	// * gl_LightSource[i].ambient;
+								diffuse_color +=	mask1 * lambertTerm;
+								specular_color += mask1 * specularTerm * specular_color;
+						}
+
+						// Calculate final color by componing different components
+						vec4 final_color;
+						vec4 color = apply_colormap(val);
+						final_color = color * (ambient_color + diffuse_color) + specular_color;
+						final_color.a = color.a;
+						return final_color;
+				}`
+  )
+};
+var Kt = /* @__PURE__ */ ((o) => (o[o.MIP = 0] = "MIP", o[o.ISO = 1] = "ISO", o))(Kt || {});
+class Mn extends ShaderMaterial {
+  constructor(e) {
+    const { style: t, threshold: n, clim: i, map: a, gradient: r } = e ?? {}, l = new Vector3(1, 1, 1);
+    if (a) {
+      const { width: p, height: m, depth: f } = a.image;
+      l.set(p, m, f);
+    }
+    const u = r == null ? r : Gn(r), c = UniformsUtils.clone(Fe.uniforms);
+    c.u_size.value.copy(l), c.u_clim.value.copy(i ?? new Vector2(1, 1)), c.u_renderstyle.value = t ?? 1, c.u_renderthreshold.value = n ?? 0.5, c.u_data.value = a ?? null, c.u_cmdata.value = u ?? null;
+    super({
+      ...Fe,
+      // @ts-ignore
+      uniforms: c,
+      side: BackSide
+    });
+    v(this, "isExampleVolumeMaterial_1", !0);
+    v(this, "_gradient");
+  }
+  /**
+   * 三维纹理
+   */
+  get map() {
+    return this.uniforms.u_data.value;
+  }
+  set map(e) {
+    this.uniforms.u_data.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 颜色映射纹理
+   * @remarks
+   * 会将垂直居中的那条水平像素从左到右映射到 clim 表示的限制区间，
+   * 即 clim.x 对应最左边的像素，clim.y 对应最右边的像素
+   */
+  get gradientMap() {
+    return this.uniforms.u_cmdata.value;
+  }
+  /**
+   * 颜色映射选项
+   */
+  get gradient() {
+    return this._gradient;
+  }
+  set gradient(e) {
+    this._gradient = e;
+    const t = e == null ? e : Gn(e);
+    this.uniforms.u_cmdata.value = t ?? null, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 三维纹理的尺寸
+   * @remarks
+   * 三维纹理的尺寸 要 和 Geometry 的尺寸是一样的，不但比列一样，大小也是一样的
+   * 但 mesh 本身可以进行几何变换；
+   * Geometry 左下角的坐标要从 (0,0,0) 开始；
+   */
+  get size() {
+    return this.uniforms.u_size.value;
+  }
+  set size(e) {
+    this.uniforms.u_size.value.copy(e), this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 渲染样式
+   */
+  get style() {
+    return this.uniforms.u_renderstyle.value;
+  }
+  set style(e) {
+    this.uniforms.u_renderstyle.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 映射区间
+   * @remarks
+   * x 为最小值，y 为最大值
+   */
+  get clim() {
+    return this.uniforms.u_clim.value;
+  }
+  set clim(e) {
+    this.uniforms.u_clim.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 渲染最低阈值
+   * @remarks
+   * 只对 ISO 样式有效；
+   * 
+   * 当三维纹理中的值低于该值时，不会被渲染，即视为空
+   */
+  get threshold() {
+    return this.uniforms.u_renderthreshold.value;
+  }
+  set threshold(e) {
+    this.uniforms.u_renderthreshold.value = e, this.uniformsNeedUpdate = !0;
+  }
+}
+const We = `in vec3 position;
+
+uniform mat4 modelMatrix;
+uniform mat4 modelViewMatrix;
+uniform mat4 projectionMatrix;
+uniform vec3 cameraPos;
+
+out vec3 vOrigin;
+out vec3 vDirection;
+
+void main() {
+    vec4 mvPosition = modelViewMatrix * vec4( position, 1.0 );
+    vOrigin = vec3( inverse( modelMatrix ) * vec4( cameraPos, 1.0 ) ).xyz;
+    vDirection = position - vOrigin;
+    gl_Position = projectionMatrix * mvPosition;
+}`, Qt = `precision highp float;
+precision highp sampler3D;
+
+uniform mat4 modelViewMatrix;
+uniform mat4 projectionMatrix;
+
+in vec3 vOrigin;
+in vec3 vDirection;
+
+out vec4 color;
+
+uniform vec3 base;
+uniform sampler3D map;
+
+uniform float threshold;
+uniform float range;
+uniform float opacity;
+uniform float steps;
+uniform float frame;
+
+uint wang_hash(uint seed)
+{
+        seed = (seed ^ 61u) ^ (seed >> 16u);
+        seed *= 9u;
+        seed = seed ^ (seed >> 4u);
+        seed *= 0x27d4eb2du;
+        seed = seed ^ (seed >> 15u);
+        return seed;
+}
+
+float randomFloat(inout uint seed)
+{
+        return float(wang_hash(seed)) / 4294967296.;
+}
+
+vec2 hitBox( vec3 orig, vec3 dir ) {
+    const vec3 box_min = vec3( - 0.5 );
+    const vec3 box_max = vec3( 0.5 );
+    vec3 inv_dir = 1.0 / dir;
+    vec3 tmin_tmp = ( box_min - orig ) * inv_dir;
+    vec3 tmax_tmp = ( box_max - orig ) * inv_dir;
+    vec3 tmin = min( tmin_tmp, tmax_tmp );
+    vec3 tmax = max( tmin_tmp, tmax_tmp );
+    float t0 = max( tmin.x, max( tmin.y, tmin.z ) );
+    float t1 = min( tmax.x, min( tmax.y, tmax.z ) );
+    return vec2( t0, t1 );
+}
+
+float sample1( vec3 p ) {
+    return texture( map, p ).r;
+}
+
+float shading( vec3 coord ) {
+    float step = 0.01;
+    return sample1( coord + vec3( - step ) ) - sample1( coord + vec3( step ) );
+}
+
+void main(){
+    vec3 rayDir = normalize( vDirection );
+    vec2 bounds = hitBox( vOrigin, rayDir );
+
+    if ( bounds.x > bounds.y ) discard;
+
+    bounds.x = max( bounds.x, 0.0 );
+
+    vec3 p = vOrigin + bounds.x * rayDir;
+    vec3 inc = 1.0 / abs( rayDir );
+    float delta = min( inc.x, min( inc.y, inc.z ) );
+    delta /= steps;
+
+    // Jitter
+
+    // Nice little seed from
+    // https://blog.demofox.org/2020/05/25/casual-shadertoy-path-tracing-1-basic-camera-diffuse-emissive/
+    uint seed = uint( gl_FragCoord.x ) * uint( 1973 ) + uint( gl_FragCoord.y ) * uint( 9277 ) + uint( frame ) * uint( 26699 );
+    vec3 size = vec3( textureSize( map, 0 ) );
+    float randNum = randomFloat( seed ) * 2.0 - 1.0;
+    p += rayDir * randNum * ( 1.0 / size );
+
+    //
+
+    vec4 ac = vec4( base, 0.0 );
+
+    for ( float t = bounds.x; t < bounds.y; t += delta ) {
+
+        float d = sample1( p + 0.5 );
+
+        d = smoothstep( threshold - range, threshold + range, d ) * opacity;
+
+        float col = shading( p + 0.5 ) * 3.0 + ( ( p.x + p.y ) * 0.25 ) + 0.2;
+
+        ac.rgb += ( 1.0 - ac.a ) * d * col;
+
+        ac.a += ( 1.0 - ac.a ) * d;
+
+        if ( ac.a >= 0.95 ) break;
+
+        p += rayDir * delta;
+
+    }
+
+    color = ac;
+
+    if ( color.a == 0.0 ) discard;
+
+}`;
+class zn extends RawShaderMaterial {
+  constructor(e) {
+    let { color: t, threshold: n, map: i, opacity: a, range: r, steps: l, frame: u, cameraPos: c } = e ?? {};
+    t = t ?? new Color(1, 1, 1), n = n ?? 0.25, a = a ?? 0.25, r = r ?? 0.1, l = l ?? 100, u = u ?? 100, c = c ?? new Vector3();
+    super({
+      glslVersion: GLSL3,
+      // @ts-ignore
+      uniforms: {
+        base: { value: t },
+        map: { value: i },
+        cameraPos: { value: c },
+        threshold: { value: n },
+        opacity: { value: a },
+        range: { value: r },
+        steps: { value: l },
+        frame: { value: u }
+      },
+      vertexShader: We,
+      fragmentShader: Qt,
+      side: BackSide,
+      transparent: !0
+    });
+    v(this, "isExampleVolumeMaterial_Cloud", !0);
+    this.opacity = a;
+  }
+  /**
+   * 三维纹理
+   */
+  get map() {
+    return this.uniforms.map.value;
+  }
+  set map(e) {
+    this.uniforms.map.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 颜色
+   */
+  get color() {
+    return this.uniforms.base.value;
+  }
+  set color(e) {
+    this.uniforms.base.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 透明度
+   */
+  // @ts-ignore
+  get opacity() {
+    return this.uniforms.opacity.value;
+  }
+  set opacity(e) {
+    this.uniforms && (this.uniforms.opacity.value = e), this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 渲染最低阈值
+   */
+  get threshold() {
+    return this.uniforms.threshold.value;
+  }
+  set threshold(e) {
+    this.uniforms.threshold.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 范围
+   */
+  get range() {
+    return this.uniforms.range.value;
+  }
+  set range(e) {
+    this.uniforms.range.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 体积材质渲染的采样数
+   */
+  get steps() {
+    return this.uniforms.steps.value;
+  }
+  set steps(e) {
+    this.uniforms.steps.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 渲染帧次
+   */
+  get frame() {
+    return this.uniforms.frame.value;
+  }
+  set frame(e) {
+    this.uniforms.frame.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 相机的位置
+   */
+  get cameraPos() {
+    return this.uniforms.cameraPos.value;
+  }
+  set cameraPos(e) {
+    this.uniforms.cameraPos.value.copy(e), this.uniformsNeedUpdate = !0;
+  }
+}
+const Jt = `precision highp float;
+precision highp sampler3D;
+
+uniform mat4 modelViewMatrix;
+uniform mat4 projectionMatrix;
+
+in vec3 vOrigin;
+in vec3 vDirection;
+
+out vec4 color;
+
+uniform sampler3D map;
+
+uniform float threshold;
+uniform float steps;
+
+vec2 hitBox( vec3 orig, vec3 dir ) {
+    const vec3 box_min = vec3( - 0.5 );
+    const vec3 box_max = vec3( 0.5 );
+    vec3 inv_dir = 1.0 / dir;
+    vec3 tmin_tmp = ( box_min - orig ) * inv_dir;
+    vec3 tmax_tmp = ( box_max - orig ) * inv_dir;
+    vec3 tmin = min( tmin_tmp, tmax_tmp );
+    vec3 tmax = max( tmin_tmp, tmax_tmp );
+    float t0 = max( tmin.x, max( tmin.y, tmin.z ) );
+    float t1 = min( tmax.x, min( tmax.y, tmax.z ) );
+    return vec2( t0, t1 );
+}
+
+float sample1( vec3 p ) {
+    return texture( map, p ).r;
+}
+
+#define epsilon .0001
+
+vec3 normal( vec3 coord ) {
+    if ( coord.x < epsilon ) return vec3( 1.0, 0.0, 0.0 );
+    if ( coord.y < epsilon ) return vec3( 0.0, 1.0, 0.0 );
+    if ( coord.z < epsilon ) return vec3( 0.0, 0.0, 1.0 );
+    if ( coord.x > 1.0 - epsilon ) return vec3( - 1.0, 0.0, 0.0 );
+    if ( coord.y > 1.0 - epsilon ) return vec3( 0.0, - 1.0, 0.0 );
+    if ( coord.z > 1.0 - epsilon ) return vec3( 0.0, 0.0, - 1.0 );
+
+    float step = 0.01;
+    float x = sample1( coord + vec3( - step, 0.0, 0.0 ) ) - sample1( coord + vec3( step, 0.0, 0.0 ) );
+    float y = sample1( coord + vec3( 0.0, - step, 0.0 ) ) - sample1( coord + vec3( 0.0, step, 0.0 ) );
+    float z = sample1( coord + vec3( 0.0, 0.0, - step ) ) - sample1( coord + vec3( 0.0, 0.0, step ) );
+
+    return normalize( vec3( x, y, z ) );
+}
+
+void main(){
+
+    vec3 rayDir = normalize( vDirection );
+    vec2 bounds = hitBox( vOrigin, rayDir );
+
+    if ( bounds.x > bounds.y ) discard;
+
+    bounds.x = max( bounds.x, 0.0 );
+
+    vec3 p = vOrigin + bounds.x * rayDir;
+    vec3 inc = 1.0 / abs( rayDir );
+    float delta = min( inc.x, min( inc.y, inc.z ) );
+    delta /= steps;
+
+    for ( float t = bounds.x; t < bounds.y; t += delta ) {
+
+        float d = sample1( p + 0.5 );
+
+        if ( d > threshold ) {
+
+            color.rgb = normal( p + 0.5 ) * 0.5 + ( p * 1.5 + 0.25 );
+            color.a = 1.;
+            break;
+
+        }
+
+        p += rayDir * delta;
+
+    }
+
+    if ( color.a == 0.0 ) discard;
+
+}`;
+class wn extends RawShaderMaterial {
+  constructor(e) {
+    let { threshold: t, map: n, steps: i, cameraPos: a } = e ?? {};
+    t = t ?? 0.25, i = i ?? 100, a = a ?? new Vector3();
+    super({
+      glslVersion: GLSL3,
+      // @ts-ignore
+      uniforms: {
+        map: { value: n },
+        cameraPos: { value: a },
+        threshold: { value: t },
+        steps: { value: i }
+      },
+      vertexShader: We,
+      fragmentShader: Jt,
+      side: BackSide,
+      transparent: !0
+    });
+    v(this, "isExampleVolumeMaterial_Perlin", !0);
+  }
+  /**
+   * 三维纹理
+   */
+  get map() {
+    return this.uniforms.map.value;
+  }
+  set map(e) {
+    this.uniforms.map.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 渲染最低阈值
+   */
+  get threshold() {
+    return this.uniforms.threshold.value;
+  }
+  set threshold(e) {
+    this.uniforms.threshold.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 体积材质渲染的采样数
+   */
+  get steps() {
+    return this.uniforms.steps.value;
+  }
+  set steps(e) {
+    this.uniforms.steps.value = e, this.uniformsNeedUpdate = !0;
+  }
+  /**
+   * 相机的位置
+   */
+  get cameraPos() {
+    return this.uniforms.cameraPos.value;
+  }
+  set cameraPos(e) {
+    this.uniforms.cameraPos.value.copy(e), this.uniformsNeedUpdate = !0;
+  }
+}
+
+class HeatMapPlugin {
+    constructor(ssp) {
+        this.ssp = ssp;
+        this.defaultColorGradient = [
+            [0, 'rgba(0,255,0,0)'],
+            [0.5, 'rgba(64,255,255,0.5)'],
+            [1, 'rgba(255,64,255,1)']
+        ];
+        this.defaultGradientVolumeMaterialOptions = {
+            fit: Dt.Raw,
+            accFactor: 2,
+            depthTest: false,
+            side: BackSide,
+            discardOut: false,
+            // 数值的梯度映射纹理
+            gradient: Gn(this.defaultColorGradient),
+        };
+    }
+    _createTexture(points, options, isLine = false) {
+        //优化
+        const optimized = yn(Object.assign({ points }, options), Object.assign({ maxSize: 60 }, options));
+        const { options: optimizedOptions, position, scale, } = optimized;
+        //创建热力3D数据
+        const heatData3D = isLine ? xn(optimizedOptions) : gn(optimizedOptions);
+        //创建热力3D纹理
+        const heatTexture = ln(heatData3D, Object.assign({}, options));
+        //创建热力材质
+        const heatMaterial = new sn(Object.assign(Object.assign(Object.assign({}, this.defaultGradientVolumeMaterialOptions), { map: heatTexture }), options));
+        const volumMesh = new rn(heatMaterial);
+        volumMesh.position.copy(position);
+        volumMesh.scale.set(scale, scale, scale);
+        this.ssp.addObject(volumMesh);
+        return volumMesh;
+    }
+    /**
+     * 创建热力云
+     * return material
+     */
+    createHeatCloud(points, options) {
+        return this._createTexture(points, options);
+    }
+    //创建线状热力
+    createLineHeat(points, options) {
+        return this._createTexture(points, options, true);
+    }
+    //创建图片挤压
+    async createImageExtrusion(imageUrl, options) {
+        var _a;
+        const material = new cn(Object.assign({ side: BackSide, depthTest: false }, options));
+        const image = await SoonSpace.utils.imageLoader.loadAsync(imageUrl);
+        const texture = un(image, (_a = options === null || options === void 0 ? void 0 : options.depth) !== null && _a !== void 0 ? _a : 200, options);
+        material.map = texture;
+        const volumMesh = new rn(material);
+        this.ssp.addObject(volumMesh);
+        return volumMesh;
+    }
+    //创建切片
+    createSliceMesh(material, options) {
+        let texture = material.map;
+        if (!texture)
+            return;
+        if (material instanceof sn) {
+            const { gradient, range, discardOut, voidRange, } = material;
+            texture = mn(texture.image, Object.assign({ gradient: gradient === null || gradient === void 0 ? void 0 : gradient.image, range,
+                discardOut,
+                voidRange }, options));
+        }
+        const sliceMesh = new vn(Object.assign({ map: texture }, options));
+        this.ssp.addObject(sliceMesh);
+        return sliceMesh;
+    }
+    //创建切片dom
+    createImageSlice(material, options) {
+        const texture = material.map;
+        if (!texture)
+            return;
+        let imageTexture = null;
+        if (material instanceof sn) {
+            const { gradient, range, discardOut, voidRange, } = material;
+            imageTexture = mn(texture.image, Object.assign({ gradient: gradient === null || gradient === void 0 ? void 0 : gradient.image, range,
+                discardOut,
+                voidRange }, options));
+        }
+        if (material instanceof cn) {
+            imageTexture = texture;
+        }
+        if (!imageTexture)
+            return;
+        const imageSlice = new hn(imageTexture);
+        return imageSlice;
+    }
+}
+
+export { Ot as ColorFogPointsMaterial, Mn as ExampleVolumeMaterial_1, zn as ExampleVolumeMaterial_Cloud, wn as ExampleVolumeMaterial_Perlin, Ye as FogPointsMaterial, Pt as GradientData3DTexture, Gt as GradientFogPointsMaterial, sn as GradientVolumeMaterial, HeatMapPlugin, ct as ImageData3DSlice, $e as ImageData3DTexture, hn as ImageData3DTextureSlice, cn as ImageVolumeMaterial, jt as SliceMaterial, vn as SliceMesh, dn as SphereFogMaterial, Dt as VolumeFit, Xe as VolumeMaterial, rn as VolumeMesh, Kt as VolumeRenderStyle, qt as computeHeatData3DForPoint, pn as createColorFogPoints, ln as createGradientData3DTexture, fn as createGradientFogPoints, gn as createHeatData3D, un as createImageData3DTextureByExtrudeImage, mn as createImageData3DTextureFromGradient, xn as createLineData3D, B as createLinearGradientImageData, Gn as createLinearGradientTexture, O as data2DCoordToIndex, at$1 as data2DIndexToCoord, _ as data3DCoordToIndex, Z as data3DIndexToCoord, HeatMapPlugin as default, nt as extrudeImage, W as extrudeImageData, X as extrudeImageDataOnAxis, G as getData2DItem, st as getData2DItemSafe, ot$1 as getData2DRow, j as getData2DValue, it as getData2DValueSafe, U as getData3DItem, $ as getData3DItemSafe, E as getData3DSlice, N as getData3DValue, tt as getData3DValueSafe, H as getImageDateColor, K as getImageDateOfImage, F as getImageDateRow, b as getNaturalSizeOfImageSource, Ut as gradientData3DToImageData3D, $t as hollowGradient_Default, Y as imageDataToCanvas, J as imageDataToUrl, L as isIImageData2, P as isIImageData3, Q as isOutOfSize, Yt as lineValueGradient_Default, bt as makeOriginOnBoundingBoxMinOfGeometry, et as mergeData3Ds, Zt as numberValuesAccumulate_Default, yn as optimizeGradientOptions, Xt as radiusGradient_Default, V as reverseImageDateY, _n as transformGradientOptions, Wt as translateScaleGradientOptions, Ze as valueGradient_Default, Ht as valuesAccumulate_Default };
+//# sourceMappingURL=index.esm.js.map