diff --git a/assets/index-DgVXdzB4.js b/assets/index-jh3uDHfm.js
similarity index 99%
rename from assets/index-DgVXdzB4.js
rename to assets/index-jh3uDHfm.js
index 915027a..aba57c6 100644
--- a/assets/index-DgVXdzB4.js
+++ b/assets/index-jh3uDHfm.js
@@ -1700,4 +1700,4 @@ void main(void) {
             <\/script>
           </body>
         </html>`}saveHTML(u="untitled.html",m="gl1",A){const S=this.generateHTML(m,A);NVUtilities.download(S,u,"application/html")}json(){return this.document.opts=this.opts,this.document.scene=this.scene,this.document.volumes=this.volumes,this.document.meshes=this.meshes,this.drawScene(),this.document.previewImageDataURL=this.canvas.toDataURL(),this.document.json()}async saveDocument(u="untitled.nvd"){this.document.title=u,log.debug("saveDocument",this.volumes[0]),this.drawScene(),this.document.previewImageDataURL=this.canvas.toDataURL(),this.document.volumes=this.volumes,this.document.meshes=this.meshes,this.document.download(u)}async loadVolumes(u){if(this.loadingText="loading...",this.drawScene(),this.thumbnailVisible)return this.deferredVolumes=u,this;this.volumes=[],this.gl.clearColor(0,0,0,1),this.gl.clear(this.gl.COLOR_BUFFER_BIT);for(let m=0;m<u.length;m++){u[m].colorMap!==void 0&&(u[m].colormap=u[m].colorMap);const A={url:u[m].url,headers:u[m].headers,name:u[m].name,colormap:u[m].colormap,colormapNegative:u[m].colormapNegative,opacity:u[m].opacity,urlImgData:u[m].urlImgData,cal_min:u[m].cal_min,cal_max:u[m].cal_max,trustCalMinMax:this.opts.trustCalMinMax,isManifest:u[m].isManifest,frame4D:u[m].frame4D,limitFrames4D:u[m].limitFrames4D||this.opts.limitFrames4D,colorbarVisible:u[m].colorbarVisible};await this.addVolumeFromUrl(A)}return this}async addMeshFromUrl(u){const m=this.getFileExt(u.url);if(m==="JCON"||m==="JSON"){const v=await(await fetch(u.url,{})).json(),w=this.loadConnectomeAsMesh(v);return this.mediaUrlMap.set(w,u.url),this.onMeshAddedFromUrl(u,w),this.addMesh(w),w}const A=await NVMesh2.loadFromUrl({...u,gl:this.gl});return this.mediaUrlMap.set(A,u.url),this.onMeshAddedFromUrl(u,A),this.addMesh(A),A}async loadMeshes(u){if(this.loadingText="loading...",this.drawScene(),this.thumbnailVisible)return this.deferredMeshes=u,this;this.initialized,this.meshes=[],this.gl.clearColor(0,0,0,1),this.gl.clear(this.gl.COLOR_BUFFER_BIT);for(let m=0;m<u.length;m++)await this.addMeshFromUrl(u[m]);return this.updateGLVolume(),this.drawScene(),this}async loadConnectomeFromUrl(u,m={}){const S=await(await fetch(u,{headers:m})).json();return this.loadConnectome(S)}async loadFreeSurferConnectomeFromUrl(u,m={}){const S=await(await fetch(u,{headers:m})).json();return this.loadFreeSurferConnectome(S)}async loadFreeSurferConnectome(u){const m=NVConnectome.convertFreeSurferConnectome(u);return this.loadConnectome(m)}handleNodeAdded(u){const m=u.detail.node,A=[1,1,1,1];this.addLabel(m.name,{textColor:A,bulletScale:1,bulletColor:A,lineWidth:0,lineColor:A,lineTerminator:"none",textScale:1},[m.x,m.y,m.z]),this.drawScene()}loadConnectomeAsMesh(u){let m=u;if("data_type"in u&&u.data_type==="fs_pointset")m=NVConnectome.convertFreeSurferConnectome(u),log.warn("converted FreeSurfer connectome",m);else if("nodes"in u){const A=u.nodes;"names"in A&&"X"in A&&"Y"in A&&"Z"in A&&"Color"in A&&"Size"in A&&(m=NVConnectome.convertLegacyConnectome(u),log.warn("converted legacy connectome",m))}else throw new Error("not a known connectome format");return new NVConnectome(this.gl,m)}loadConnectome(u){this.loadingText="loading...",this.drawScene(),this.meshes=[],this.gl.clearColor(0,0,0,1),this.gl.clear(this.gl.COLOR_BUFFER_BIT);const m=this.loadConnectomeAsMesh(u);return this.addMesh(m),this.drawScene(),this}createEmptyDrawing(){if(this.back===null||!this.back.dims||Math.min(Math.min(this.back.dims[1],this.back.dims[2]),this.back.dims[3])<1)return;const m=this.back.dims[1]*this.back.dims[2]*this.back.dims[3];this.drawBitmap=new Uint8Array(m),this.drawClearAllUndoBitmaps(),this.drawAddUndoBitmap(),this.drawTexture=this.r8Tex(this.drawTexture,TEXTURE7_DRAW,this.back.dims,!0),this.refreshDrawing(!1)}r16Tex(u,m,A,S){u&&this.gl.deleteTexture(u),u=this.gl.createTexture(),this.gl.activeTexture(m),this.gl.bindTexture(this.gl.TEXTURE_3D,u),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_MIN_FILTER,this.gl.NEAREST),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_MAG_FILTER,this.gl.NEAREST),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_WRAP_R,this.gl.CLAMP_TO_EDGE),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_WRAP_S,this.gl.CLAMP_TO_EDGE),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_WRAP_T,this.gl.CLAMP_TO_EDGE),this.gl.pixelStorei(this.gl.UNPACK_ALIGNMENT,1),this.gl.texStorage3D(this.gl.TEXTURE_3D,1,this.gl.R16I,A[1],A[2],A[3]);const v=A[1]*A[2]*A[3];return S.length!==v&&(S=new Int16Array(v)),this.gl.texSubImage3D(this.gl.TEXTURE_3D,0,0,0,0,A[1],A[2],A[3],this.gl.RED_INTEGER,this.gl.SHORT,S),u}drawGrowCut(){if(!this.back||!this.back.dims)throw new Error("back not defined");const u=this.back.hdr,m=this.gl,A=u.dims[1]*u.dims[2]*u.dims[3];if(!this.drawBitmap||this.drawBitmap.length!==A){log.debug("bitmap dims are wrong");return}const S=m.createFramebuffer();m.bindFramebuffer(m.FRAMEBUFFER,S),m.disable(m.CULL_FACE),m.viewport(0,0,this.back.dims[1],this.back.dims[2]),m.disable(m.BLEND);let v=img2ras16(this.back);const w=this.r16Tex(null,TEXTURE11_GC_BACK,this.back.dims,v);for(let Y=1;Y<A;Y++)v[Y]=this.drawBitmap[Y];const D=this.r16Tex(null,TEXTURE14_GC_LABEL0,this.back.dims,v),R=this.r16Tex(null,TEXTURE15_GC_LABEL1,this.back.dims,v),P=1e4;for(let Y=1;Y<A;Y++)v[Y]>0&&(v[Y]=P);const L=this.r16Tex(null,TEXTURE12_GC_STRENGTH0,this.back.dims,v),d=this.r16Tex(null,TEXTURE13_GC_STRENGTH1,this.back.dims,v);m.bindVertexArray(this.genericVAO);const E=this.growCutShader;E.use(m);const e=128;m.uniform1i(E.uniforms.finalPass,0),m.uniform1i(E.uniforms.backTex,11);for(let Y=0;Y<e;Y++){m.uniform1i(E.uniforms.labelTex,14),m.uniform1i(E.uniforms.strengthTex,12);for(let G=0;G<this.back.dims[3];G++){const J=1/this.back.dims[3]*(G+.5);m.uniform1f(E.uniforms.coordZ,J),m.framebufferTextureLayer(m.FRAMEBUFFER,m.COLOR_ATTACHMENT0,R,0,G),m.framebufferTextureLayer(m.FRAMEBUFFER,m.COLOR_ATTACHMENT1,d,0,G),m.drawBuffers([m.COLOR_ATTACHMENT0,m.COLOR_ATTACHMENT1]),m.checkFramebufferStatus(m.FRAMEBUFFER)!==m.FRAMEBUFFER_COMPLETE&&log.error("Incomplete framebuffer"),m.drawArrays(m.TRIANGLE_STRIP,0,4)}Y===e-1&&m.uniform1i(E.uniforms.finalPass,1),m.uniform1i(E.uniforms.labelTex,15),m.uniform1i(E.uniforms.strengthTex,13);for(let G=0;G<this.back.dims[3];G++){const J=1/this.back.dims[3]*(G+.5);m.uniform1f(E.uniforms.coordZ,J),m.framebufferTextureLayer(m.FRAMEBUFFER,m.COLOR_ATTACHMENT0,D,0,G),m.framebufferTextureLayer(m.FRAMEBUFFER,m.COLOR_ATTACHMENT1,L,0,G),m.drawBuffers([m.COLOR_ATTACHMENT0,m.COLOR_ATTACHMENT1]),m.checkFramebufferStatus(m.FRAMEBUFFER)!==m.FRAMEBUFFER_COMPLETE&&log.error("Incomplete framebuffer"),m.drawArrays(m.TRIANGLE_STRIP,0,4)}}m.drawBuffers([m.COLOR_ATTACHMENT0]);const N=m.COLOR_ATTACHMENT1,q=D;m.readBuffer(N);const z=m.getParameter(m.IMPLEMENTATION_COLOR_READ_FORMAT),e0=m.getParameter(m.IMPLEMENTATION_COLOR_READ_TYPE);(z!==m.RED_INTEGER||e0!==m.SHORT)&&log.debug("readPixels will fail."),v=new Int16Array;const H=this.back.dims[1]*this.back.dims[2],j=new Int16Array(H);for(let Y=0;Y<this.back.dims[3];Y++)m.framebufferTextureLayer(m.FRAMEBUFFER,N,q,0,Y),m.readPixels(0,0,this.back.dims[1],this.back.dims[2],z,e0,j),v=Int16Array.from([...v,...j]);let y=v[0];for(let Y=0;Y<v.length;Y++)y=Math.max(y,v[Y]);for(let Y=1;Y<A;Y++)this.drawBitmap[Y]=v[Y];m.deleteTexture(w),m.deleteTexture(L),m.deleteTexture(d),m.deleteTexture(D),m.deleteTexture(R),m.bindVertexArray(this.unusedVAO),m.viewport(0,0,m.canvas.width,m.canvas.height),m.bindFramebuffer(m.FRAMEBUFFER,null),m.deleteFramebuffer(S),this.drawAddUndoBitmap(),this.refreshDrawing(!0)}drawPt(u,m,A,S){var R;if(!((R=this.back)!=null&&R.dims))throw new Error("back.dims not set");const v=this.back.dims[1],w=this.back.dims[2],D=this.back.dims[3];u=Math.min(Math.max(u,0),v-1),m=Math.min(Math.max(m,0),w-1),A=Math.min(Math.max(A,0),D-1),this.drawBitmap[u+m*v+A*v*w]=S}drawPenLine(u,m,A){const S=Math.abs(u[0]-m[0]),v=Math.abs(u[1]-m[1]),w=Math.abs(u[2]-m[2]);let D=-1,R=-1,P=-1;m[0]>u[0]&&(D=1),m[1]>u[1]&&(R=1),m[2]>u[2]&&(P=1);let L=u[0],d=u[1],E=u[2];const e=m[0],N=m[1],q=m[2];if(S>=v&&S>=w){let z=2*v-S,e0=2*w-S;for(;L!==e;)L+=D,z>=0&&(d+=R,z-=2*S),e0>=0&&(E+=P,e0-=2*S),z+=2*v,e0+=2*w,this.drawPt(L,d,E,A)}else if(v>=S&&v>=w){let z=2*S-v,e0=2*w-v;for(;d!==N;)d+=R,z>=0&&(L+=D,z-=2*v),e0>=0&&(E+=P,e0-=2*v),z+=2*S,e0+=2*w,this.drawPt(L,d,E,A)}else{let z=2*v-w,e0=2*S-w;for(;E!==q;)E+=P,z>=0&&(d+=R,z-=2*w),e0>=0&&(L+=D,e0-=2*w),z+=2*v,e0+=2*S,this.drawPt(L,d,E,A)}}drawFloodFillCore(u,m,A=6){var L;if(!((L=this.back)!=null&&L.dims))throw new Error("back.dims undefined");const S=[this.back.dims[1],this.back.dims[2],this.back.dims[3]],v=S[0],w=v*S[1];function D(d){return d[0]+d[1]*v+d[2]*w}function R(d){const E=Math.floor(d/w),e=Math.floor((d-E*w)/v);return[Math.floor(d%v),e,E]}const P=[];for(P.push(m),u[m]=2;P.length>0;){let d=function(N){const q=e.slice();if(q[0]+=N[0],q[1]+=N[1],q[2]+=N[2],q[0]<0||q[1]<0||q[2]<0||q[0]>=S[0]||q[1]>=S[1]||q[2]>=S[2])return;const z=D(q);u[z]===1&&(u[z]=2,P.push(z))};const E=P[0];P.shift();const e=R(E);d([0,0,-1]),d([0,0,1]),d([0,-1,0]),d([0,1,0]),d([-1,0,0]),d([1,0,0]),!(A<=6)&&(d([-1,-1,0]),d([1,1,0]),d([-1,1,0]),d([1,1,0]),d([0,-1,-1]),d([0,1,-1]),d([-1,0,-1]),d([1,0,-1]),d([0,-1,1]),d([0,1,1]),d([-1,0,1]),d([1,0,1]),!(A<=18)&&(d([-1,-1,-1]),d([1,-1,-1]),d([-1,1,-1]),d([1,1,-1]),d([-1,-1,1]),d([1,-1,1]),d([-1,1,1]),d([1,1,1])))}}drawFloodFill(u,m=0,A=0,S=NaN,v=NaN,w=6){var q;if(!this.drawBitmap)throw new Error("drawBitmap undefined");if(!((q=this.back)!=null&&q.dims))throw new Error("back.dims undefined");m=Math.abs(m);const D=[this.back.dims[1],this.back.dims[2],this.back.dims[3]];if(u[0]<0||u[1]<0||u[2]<0||u[0]>=D[0]||u[1]>=D[1]||u[2]>=D[2])return;const R=D[0],P=R*D[1],L=P*D[2],d=this.drawBitmap.slice();if(d.length!==P*D[2])return;function E(z){return z[0]+z[1]*R+z[2]*P}const e=E(u),N=d[e];if(N===m){A!==0?log.debug("drawFloodFill selected voxel is not part of a drawing"):log.debug("drawFloodFill selected voxel is already desired color");return}for(let z=1;z<L;z++)d[z]=0,this.drawBitmap[z]===N&&(d[z]=1);if(this.drawFloodFillCore(d,e,w),A!==0){const z=this.volumes[0].img2RAS();let e0=z[e],H=e0;if(isFinite(v)&&isFinite(S))e0=v,H=S;else{for(let j=1;j<L;j++)d[j]===2&&(e0=Math.max(e0,z[j]),H=Math.min(H,z[j]));A===Number.POSITIVE_INFINITY&&(e0=A),A===Number.NEGATIVE_INFINITY&&(H=A)}log.debug("Intensity range of selected cluster :",H,e0);for(let j=1;j<L;j++)d[j]=0,z[j]>=H&&z[j]<=e0&&(d[j]=1);this.drawFloodFillCore(d,e,w),m=N}for(let z=1;z<L;z++)d[z]===2&&(this.drawBitmap[z]=m);this.drawAddUndoBitmap(),this.refreshDrawing(!1)}drawPenFilled(){var N;const u=this.drawPenFillPts.length;if(u<2){this.drawPenFillPts=[];return}const m=this.drawPenAxCorSag;let A=0,S=1;if(m===1&&(S=2),m===2&&(A=1,S=2),!((N=this.back)!=null&&N.dims))throw new Error("back.dims undefined");const v=[this.back.dims[A+1],this.back.dims[S+1]],w=new Uint8Array(v[0]*v[1]);let D=1;function R(q,z){const e0=Math.abs(q[0]-z[0]),H=Math.abs(q[1]-z[1]);w[q[0]+q[1]*v[0]]=D,w[z[0]+z[1]*v[0]]=D;let j=-1,y=-1;z[0]>q[0]&&(j=1),z[1]>q[1]&&(y=1);let Y=q[0],G=q[1];const J=z[0],i0=z[1];if(e0>=H){let Q=2*H-e0;for(;Y!==J;)Y+=j,Q>=0&&(G+=y,Q-=2*e0),Q+=2*H,w[Y+G*v[0]]=D}else{let Q=2*e0-H;for(;G!==i0;)G+=y,Q>=0&&(Y+=j,Q-=2*H),Q+=2*e0,w[Y+G*v[0]]=D}}const P=[this.drawPenFillPts[0][A],this.drawPenFillPts[0][S]];let L=P;for(let q=1;q<u;q++){const z=[this.drawPenFillPts[q][A],this.drawPenFillPts[q][S]];R(L,z),L=z}R(P,L);const d=[];function E(q){if(q[0]<0||q[1]<0||q[0]>=v[0]||q[1]>=v[1])return;const z=q[0]+q[1]*v[0];w[z]===0&&(d.push(q),w[z]=2)}for(let q=0;q<v[0];q++)E([q,0]);for(let q=0;q<v[0];q++)E([q,v[1]-1]);for(let q=0;q<v[1];q++)E([0,q]);for(let q=0;q<v[1];q++)E([v[0]-1,q]);for(;d.length>0;){const q=d.shift();E([q[0]-1,q[1]]),E([q[0]+1,q[1]]),E([q[0],q[1]-1]),E([q[0],q[1]+1])}D=this.opts.penValue;const e=this.drawPenFillPts[0][3-(A+S)];if(!this.drawBitmap)throw new Error("drawBitmap undefined");if(m===0){const q=e*v[0]*v[1];for(let z=0;z<v[0]*v[1];z++)w[z]!==2&&(this.drawBitmap[z+q]=D)}else{let q=1;const z=this.back.dims[1]*this.back.dims[2];let e0=e*this.back.dims[1];m===2&&(q=this.back.dims[1],e0=e);let H=0;for(let j=0;j<v[1];j++)for(let y=0;y<v[0];y++)w[H]!==2&&(this.drawBitmap[y*q+j*z+e0]=D),H++}if(!this.drawFillOverwrites&&this.drawUndoBitmaps[this.currentDrawUndoBitmap].length>0){const q=this.drawBitmap.length,z=decodeRLE(this.drawUndoBitmaps[this.currentDrawUndoBitmap],q);for(let e0=0;e0<q;e0++)z[e0]!==0&&(this.drawBitmap[e0]=z[e0])}this.drawPenFillPts=[],this.drawAddUndoBitmap(),this.refreshDrawing(!1)}closeDrawing(){this.drawClearAllUndoBitmaps(),this.rgbaTex(this.drawTexture,TEXTURE7_DRAW,[2,2,2,2],!0),this.drawBitmap=null,this.drawScene()}refreshDrawing(u=!0){var S;if(!((S=this.back)!=null&&S.dims))throw new Error("back.dims undefined");if(!this.drawBitmap)throw new Error("drawBitmap undefined");const m=this.back.dims.slice(),A=this.back.dims[1]*this.back.dims[2]*this.back.dims[3];this.drawBitmap.length===8?(m[1]=2,m[2]=2,m[3]=2):A!==this.drawBitmap.length&&log.warn("Drawing bitmap must match the background image"),this.gl.activeTexture(TEXTURE7_DRAW),this.gl.bindTexture(this.gl.TEXTURE_3D,this.drawTexture),this.gl.texSubImage3D(this.gl.TEXTURE_3D,0,0,0,0,m[1],m[2],m[3],this.gl.RED,this.gl.UNSIGNED_BYTE,this.drawBitmap),u&&this.drawScene()}r8Tex(u,m,A,S=!1){if(u&&this.gl.deleteTexture(u),u=this.gl.createTexture(),this.gl.activeTexture(m),this.gl.bindTexture(this.gl.TEXTURE_3D,u),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_MIN_FILTER,this.gl.NEAREST),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_MAG_FILTER,this.gl.NEAREST),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_WRAP_R,this.gl.CLAMP_TO_EDGE),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_WRAP_S,this.gl.CLAMP_TO_EDGE),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_WRAP_T,this.gl.CLAMP_TO_EDGE),this.gl.pixelStorei(this.gl.UNPACK_ALIGNMENT,1),this.gl.texStorage3D(this.gl.TEXTURE_3D,1,this.gl.R8,A[1],A[2],A[3]),S){const v=new Uint8Array(A[1]*A[2]*A[3]);this.gl.texSubImage3D(this.gl.TEXTURE_3D,0,0,0,0,A[1],A[2],A[3],this.gl.RED,this.gl.UNSIGNED_BYTE,v)}return u}rgbaTex(u,m,A,S=!1){if(u&&this.gl.deleteTexture(u),u=this.gl.createTexture(),this.gl.activeTexture(m),this.gl.bindTexture(this.gl.TEXTURE_3D,u),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_MIN_FILTER,this.gl.LINEAR),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_MAG_FILTER,this.gl.LINEAR),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_WRAP_R,this.gl.CLAMP_TO_EDGE),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_WRAP_S,this.gl.CLAMP_TO_EDGE),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_WRAP_T,this.gl.CLAMP_TO_EDGE),this.gl.pixelStorei(this.gl.UNPACK_ALIGNMENT,1),this.gl.texStorage3D(this.gl.TEXTURE_3D,1,this.gl.RGBA8,A[1],A[2],A[3]),S){const v=new Uint8Array(A[1]*A[2]*A[3]*4);this.gl.texSubImage3D(this.gl.TEXTURE_3D,0,0,0,0,A[1],A[2],A[3],this.gl.RGBA,this.gl.UNSIGNED_BYTE,v)}return u}requestCORSIfNotSameOrigin(u,m){new URL(m,window.location.href).origin!==window.location.origin&&(u.crossOrigin="")}async loadPngAsTexture(u,m){return new Promise((A,S)=>{const v=new Image;v.onload=()=>{if(!this.bmpShader)return;let w;m===4?(this.bmpTexture!==null&&this.gl.deleteTexture(this.bmpTexture),this.bmpTexture=this.gl.createTexture(),w=this.bmpTexture,this.bmpTextureWH=v.width/v.height,this.gl.activeTexture(TEXTURE4_THUMBNAIL),this.bmpShader.use(this.gl),this.gl.uniform1i(this.bmpShader.uniforms.bmpTexture,4)):m===5?(this.gl.activeTexture(TEXTURE5_MATCAP),this.matCapTexture!==null&&this.gl.deleteTexture(this.matCapTexture),this.matCapTexture=this.gl.createTexture(),w=this.matCapTexture):(this.fontShader.use(this.gl),this.gl.activeTexture(TEXTURE3_FONT),this.gl.uniform1i(this.fontShader.uniforms.fontTexture,3),this.fontTexture!==null&&this.gl.deleteTexture(this.fontTexture),this.fontTexture=this.gl.createTexture(),w=this.fontTexture),this.gl.bindTexture(this.gl.TEXTURE_2D,w),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_WRAP_S,this.gl.CLAMP_TO_EDGE),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_WRAP_T,this.gl.CLAMP_TO_EDGE),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_MIN_FILTER,this.gl.LINEAR),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_MAG_FILTER,this.gl.LINEAR),this.gl.texImage2D(this.gl.TEXTURE_2D,0,this.gl.RGBA,this.gl.RGBA,this.gl.UNSIGNED_BYTE,v),A(w),m!==4&&this.drawScene()},v.onerror=S,this.requestCORSIfNotSameOrigin(v,u),v.src=u})}async loadFontTexture(u){return this.loadPngAsTexture(u,3)}async loadBmpTexture(u){return this.loadPngAsTexture(u,4)}async loadMatCapTexture(u){return this.loadPngAsTexture(u,5)}initFontMets(){if(!this.fontMetrics)throw new Error("fontMetrics undefined");this.fontMets={distanceRange:this.fontMetrics.atlas.distanceRange,size:this.fontMetrics.atlas.size,mets:{}};for(let A=0;A<256;A++)this.fontMets.mets[A]={xadv:0,uv_lbwh:[0,0,0,0],lbwh:[0,0,0,0]};const u=this.fontMetrics.atlas.width,m=this.fontMetrics.atlas.height;for(let A=0;A<this.fontMetrics.glyphs.length;A++){const S=this.fontMetrics.glyphs[A],v=S.unicode;if(this.fontMets.mets[v].xadv=S.advance,S.planeBounds===void 0)continue;let w=S.atlasBounds.left/u,D=(m-S.atlasBounds.top)/m,R=(S.atlasBounds.right-S.atlasBounds.left)/u,P=(S.atlasBounds.top-S.atlasBounds.bottom)/m;this.fontMets.mets[v].uv_lbwh=[w,D,R,P],w=S.planeBounds.left,D=S.planeBounds.bottom,R=S.planeBounds.right-S.planeBounds.left,P=S.planeBounds.top-S.planeBounds.bottom,this.fontMets.mets[v].lbwh=[w,D,R,P]}}async loadFont(u=Roboto_Regular_default,m=Roboto_Regular_default2){await this.loadFontTexture(u);const A=await fetch(m);if(!A.ok)throw Error(A.statusText);const S=await A.text();this.fontMetrics=JSON.parse(S),this.initFontMets(),this.fontShader.use(this.gl),this.drawScene()}async loadDefaultMatCap(){return this.loadMatCapTexture(Shiny_default)}async loadDefaultFont(){await this.loadFontTexture(this.DEFAULT_FONT_GLYPH_SHEET),this.fontMetrics=this.DEFAULT_FONT_METRICS,this.initFontMets()}async initText(){this.fontShader=new Shader(this.gl,vertFontShader,fragFontShader),this.fontShader.use(this.gl),await this.loadDefaultFont(),await this.loadDefaultMatCap(),this.drawLoadingText(this.loadingText)}meshShaderNameToNumber(u="Phong"){const m=u.toLowerCase();for(let A=0;A<this.meshShaders.length;A++)if(this.meshShaders[A].Name.toLowerCase()===m)return A}setMeshShader(u,m=2){let A=0;if(typeof m=="number"?A=m:A=this.meshShaderNameToNumber(m),A===void 0)throw new Error("shaderIndex undefined");A=Math.min(A,this.meshShaders.length-1),A=Math.max(A,0);const S=this.getMeshIndexByID(u);if(S>=this.meshes.length){log.debug("Unable to change shader until mesh is loaded (maybe you need async)");return}this.meshes[S].meshShaderIndex=A,this.updateGLVolume(),this.onMeshShaderChanged(S,A)}createCustomMeshShader(u,m="Custom"){if(!u)throw new Error("Need fragment shader");const A=this.meshShaderNameToNumber(m);A>=0&&(this.gl.deleteProgram(this.meshShaders[A].shader.program),this.meshShaders.splice(A,1));const S=new Shader(this.gl,vertMeshShader,u);return S.use(this.gl),{Name:m,Frag:u,shader:S}}setCustomMeshShader(u="",m="Custom"){const A=this.createCustomMeshShader(u,m);return this.meshShaders.push(A),this.onCustomMeshShaderAdded(u,m),this.meshShaders.length-1}meshShaderNames(u=!0){const m=[];for(let A=0;A<this.meshShaders.length;A++)m.push(this.meshShaders[A].Name);return u===!0?m.sort():m}initRenderShader(u,m=0){u.use(this.gl),this.gl.uniform1i(u.uniforms.volume,0),this.gl.uniform1i(u.uniforms.colormap,1),this.gl.uniform1i(u.uniforms.overlay,2),this.gl.uniform1i(u.uniforms.drawing,7),this.gl.uniform1fv(u.uniforms.renderDrawAmbientOcclusion,[this.renderDrawAmbientOcclusion,1]),this.gl.uniform1f(u.uniforms.gradientAmount,m)}async init(){const u=this.gl.getExtension("WEBGL_debug_renderer_info");if(u){const w=this.gl.getParameter(u.UNMASKED_VENDOR_WEBGL),D=this.gl.getParameter(u.UNMASKED_RENDERER_WEBGL);log.info("renderer vendor: ",w),log.info("renderer: ",D)}else log.info("debug_renderer_info unavailable");const m=this.gl.getParameter(this.gl.RENDERER);log.info("firefox renderer: ",m),this.gl.clearDepth(0),this.gl.enable(this.gl.CULL_FACE),this.gl.cullFace(this.gl.FRONT),this.gl.enable(this.gl.BLEND),this.gl.blendFunc(this.gl.SRC_ALPHA,this.gl.ONE_MINUS_SRC_ALPHA),this.volumeTexture=this.rgbaTex(this.volumeTexture,TEXTURE0_BACK_VOL,[2,2,2,2],!0),this.overlayTexture=this.rgbaTex(this.overlayTexture,TEXTURE2_OVERLAY_VOL,[2,2,2,2],!0),this.drawTexture=this.r8Tex(this.drawTexture,TEXTURE7_DRAW,[2,2,2,2],!0);const A=[1,1,0,1,0,0,0,1,0,0,0,0],S=this.gl;this.cuboidVertexBuffer=S.createBuffer(),S.bindBuffer(S.ARRAY_BUFFER,this.cuboidVertexBuffer),S.bufferData(S.ARRAY_BUFFER,new Float32Array(A),S.STATIC_DRAW),this.genericVAO=S.createVertexArray(),S.bindVertexArray(this.genericVAO),S.bindBuffer(S.ARRAY_BUFFER,this.cuboidVertexBuffer),S.enableVertexAttribArray(0),S.vertexAttribPointer(0,3,S.FLOAT,!1,0,0),S.bindVertexArray(this.unusedVAO),this.pickingMeshShader=new Shader(S,vertMeshShader,fragMeshDepthShader),this.pickingMeshShader.use(S),this.pickingImageShader=new Shader(S,vertRenderShader,fragVolumePickingShader),this.pickingImageShader.use(S),S.uniform1i(this.pickingImageShader.uniforms.volume,0),S.uniform1i(this.pickingImageShader.uniforms.colormap,1),S.uniform1i(this.pickingImageShader.uniforms.overlay,2),S.uniform1i(this.pickingImageShader.uniforms.drawing,7),this.sliceMMShader=new Shader(S,vertSliceMMShader,fragSliceMMShader),this.sliceMMShader.use(S),S.uniform1i(this.sliceMMShader.uniforms.volume,0),S.uniform1i(this.sliceMMShader.uniforms.colormap,1),S.uniform1i(this.sliceMMShader.uniforms.overlay,2),S.uniform1i(this.sliceMMShader.uniforms.drawing,7),S.uniform1f(this.sliceMMShader.uniforms.drawOpacity,this.drawOpacity),this.sliceV1Shader=new Shader(S,vertSliceMMShader,fragSliceV1Shader),this.sliceV1Shader.use(S),S.uniform1i(this.sliceV1Shader.uniforms.volume,0),S.uniform1i(this.sliceV1Shader.uniforms.colormap,1),S.uniform1i(this.sliceV1Shader.uniforms.overlay,2),S.uniform1i(this.sliceV1Shader.uniforms.drawing,7),S.uniform1f(this.sliceV1Shader.uniforms.drawOpacity,this.drawOpacity),this.orientCubeShader=new Shader(S,vertOrientCubeShader,fragOrientCubeShader),this.orientCubeShaderVAO=S.createVertexArray(),S.bindVertexArray(this.orientCubeShaderVAO);const v=S.createBuffer();S.enableVertexAttribArray(0),S.enableVertexAttribArray(1),S.bindBuffer(S.ARRAY_BUFFER,v),S.bufferData(S.ARRAY_BUFFER,orientCube,S.STATIC_DRAW),S.vertexAttribPointer(0,3,S.FLOAT,!1,24,0),S.enableVertexAttribArray(1),S.vertexAttribPointer(1,3,S.FLOAT,!1,24,12),S.bindVertexArray(this.unusedVAO),this.rectShader=new Shader(S,vertRectShader,fragRectShader),this.rectShader.use(S),this.lineShader=new Shader(S,vertLineShader,fragRectShader),this.lineShader.use(S),this.line3DShader=new Shader(S,vertLine3DShader,fragRectShader),this.line3DShader.use(S),this.circleShader=new Shader(S,vertCircleShader,fragCircleShader),this.circleShader.use(S),this.renderVolumeShader=new Shader(S,vertRenderShader,fragRenderShader),this.initRenderShader(this.renderVolumeShader),this.renderSliceShader=new Shader(S,vertRenderShader,fragRenderSliceShader),this.initRenderShader(this.renderSliceShader),this.renderGradientShader=new Shader(S,vertRenderShader,fragRenderGradientShader),this.initRenderShader(this.renderGradientShader,.3),S.uniform1i(this.renderGradientShader.uniforms.matCap,5),S.uniform1i(this.renderGradientShader.uniforms.gradient,6),this.renderShader=this.renderVolumeShader,this.colorbarShader=new Shader(S,vertColorbarShader,fragColorbarShader),this.colorbarShader.use(S),S.uniform1i(this.colorbarShader.uniforms.colormap,1),this.blurShader=new Shader(S,blurVertShader,blurFragShader),this.sobelShader=new Shader(S,blurVertShader,sobelFragShader),this.growCutShader=new Shader(S,vertGrowCutShader,fragGrowCutShader),this.passThroughShader=new Shader(S,vertPassThroughShader,fragPassThroughShader),this.orientShaderAtlasU=new Shader(S,vertOrientShader,fragOrientShaderU.concat(fragOrientShaderAtlas)),this.orientShaderAtlasI=new Shader(S,vertOrientShader,fragOrientShaderI.concat(fragOrientShaderAtlas)),this.orientShaderU=new Shader(S,vertOrientShader,fragOrientShaderU.concat(fragOrientShader)),this.orientShaderI=new Shader(S,vertOrientShader,fragOrientShaderI.concat(fragOrientShader)),this.orientShaderF=new Shader(S,vertOrientShader,fragOrientShaderF.concat(fragOrientShader)),this.orientShaderRGBU=new Shader(S,vertOrientShader,fragOrientShaderU.concat(fragRGBOrientShader)),this.surfaceShader=new Shader(S,vertSurfaceShader,fragSurfaceShader),this.surfaceShader.use(S),this.fiberShader=new Shader(S,vertFiberShader,fragFiberShader),this.pickingImageShader.use(S);for(let w=0;w<this.meshShaders.length;w++){const D=this.meshShaders[w];D.Name==="Flat"?D.shader=new Shader(S,vertFlatMeshShader,fragFlatMeshShader):D.shader=new Shader(S,vertMeshShader,D.Frag),D.shader.use(S),D.shader.isMatcap=D.Name==="Matcap",D.shader.isMatcap&&S.uniform1i(D.shader.uniforms.matCap,5)}return this.bmpShader=new Shader(S,vertBmpShader,fragBmpShader),await this.initText(),this.opts.thumbnail.length>0&&(await this.loadBmpTexture(this.opts.thumbnail),this.thumbnailVisible=!0),this.updateGLVolume(),this.initialized=!0,this.resizeListener(),this.drawScene(),this}gradientGL(u){const m=this.gl,A=[0,0,0,0,1,0,1,0,0,1,1,0],S=m.createVertexArray();m.bindVertexArray(S);const v=m.createBuffer();m.bindBuffer(m.ARRAY_BUFFER,v),m.bufferData(m.ARRAY_BUFFER,new Float32Array(A),m.STATIC_DRAW),m.enableVertexAttribArray(0),m.vertexAttribPointer(0,3,m.FLOAT,!1,0,0);const w=m.createFramebuffer();m.bindFramebuffer(m.FRAMEBUFFER,w),m.disable(m.CULL_FACE),m.viewport(0,0,u.dims[1],u.dims[2]),m.disable(m.BLEND);const D=this.rgbaTex(null,TEXTURE8_GRADIENT_TEMP,u.dims),R=this.blurShader;R.use(m),m.activeTexture(TEXTURE0_BACK_VOL),m.bindTexture(m.TEXTURE_3D,this.volumeTexture);const P=.7;m.uniform1i(R.uniforms.intensityVol,0),m.uniform1f(R.uniforms.dX,P/u.dims[1]),m.uniform1f(R.uniforms.dY,P/u.dims[2]),m.uniform1f(R.uniforms.dZ,P/u.dims[3]),m.bindVertexArray(S);for(let E=0;E<u.dims[3]-1;E++){const e=1/u.dims[3]*(E+.5);m.uniform1f(R.uniforms.coordZ,e),m.framebufferTextureLayer(m.FRAMEBUFFER,m.COLOR_ATTACHMENT0,D,0,E),m.clear(m.DEPTH_BUFFER_BIT),m.drawArrays(m.TRIANGLE_STRIP,0,A.length/3)}const L=this.sobelShader;L.use(m),m.activeTexture(TEXTURE8_GRADIENT_TEMP),m.bindTexture(m.TEXTURE_3D,D),m.uniform1i(L.uniforms.intensityVol,8);const d=.7;m.uniform1f(L.uniforms.dX,d/u.dims[1]),m.uniform1f(L.uniforms.dY,d/u.dims[2]),m.uniform1f(L.uniforms.dZ,d/u.dims[3]),m.uniform1f(L.uniforms.coordZ,.5),m.bindVertexArray(S),m.activeTexture(TEXTURE0_BACK_VOL),this.gradientTexture!==null&&m.deleteTexture(this.gradientTexture),this.gradientTexture=this.rgbaTex(this.gradientTexture,TEXTURE6_GRADIENT,u.dims);for(let E=0;E<u.dims[3]-1;E++){const e=1/u.dims[3]*(E+.5);m.uniform1f(L.uniforms.coordZ,e),m.framebufferTextureLayer(m.FRAMEBUFFER,m.COLOR_ATTACHMENT0,this.gradientTexture,0,E),m.clear(m.DEPTH_BUFFER_BIT),m.drawArrays(m.TRIANGLE_STRIP,0,A.length/3)}m.deleteFramebuffer(w),m.deleteTexture(D),m.deleteBuffer(v),m.bindFramebuffer(m.FRAMEBUFFER,null)}updateGLVolume(){let u=0;const m=this.volumes.length;this.refreshColormaps();for(let A=0;A<m;A++)this.volumes[A].toRAS&&(this.refreshLayers(this.volumes[A],u),u++);if(this.furthestVertexFromOrigin=0,m>0&&(this.furthestVertexFromOrigin=this.volumeObject3D.furthestVertexFromOrigin),this.meshes)for(let A=0;A<this.meshes.length;A++)this.furthestVertexFromOrigin=Math.max(this.furthestVertexFromOrigin,this.meshes[A].furthestVertexFromOrigin);this.onVolumeUpdated&&this.onVolumeUpdated(),this.drawScene()}getDescriptives(u=0,m=[],A=!1){const S=this.volumes[u].hdr;let v=S.scl_slope;isNaN(v)&&(v=1);let w=S.scl_inter;isNaN(w)&&(w=1);const D=this.volumes[u].img,R=D.length,P=new Float32Array(R);for(let J=0;J<R;J++)P[J]=D[J]*v+w;const L=new Uint8Array(R);for(let J=0;J<R;J++)L[J]=1;if(m.length>0)for(let J=0;J<m.length;J++){const i0=this.volumes[m[J]].img;if(i0.length!==R){log.debug("Mask resolution does not match image. Skipping masking layer "+m[J]);continue}for(let Q=0;Q<R;Q++)(i0[Q]===0||isNaN(i0[Q]))&&(L[Q]=0)}else if(m.length<1&&A)for(let J=0;J<R;J++)(this.drawBitmap[J]===0||isNaN(this.drawBitmap[J]))&&(L[J]=0);let d=0,E=0,e=0,N=Number.NEGATIVE_INFINITY,q=Number.POSITIVE_INFINITY,z=0,e0=0,H=0;for(let J=0;J<R;J++){if(L[J]<1)continue;const i0=P[J];d++;let Q=E+(i0-E)/d;e=e+(i0-E)*(i0-Q),E=Q,i0!==0&&(z++,Q=e0+(i0-e0)/z,H=H+(i0-e0)*(i0-Q),e0=Q,q=Math.min(i0,N),N=Math.max(i0,N))}const j=Math.sqrt(e/(d-1)),y=Math.sqrt(H/(z-1)),Y=q,G=N;return d!==z&&(q=Math.min(0,N),N=Math.max(0,N)),{mean:E,stdev:j,nvox:d,volumeMM3:d*S.pixDims[1]*S.pixDims[2]*S.pixDims[3],volumeML:d*S.pixDims[1]*S.pixDims[2]*S.pixDims[3]*.001,min:q,max:N,meanNot0:e0,stdevNot0:y,nvoxNot0:z,minNot0:Y,maxNot0:G,cal_min:this.volumes[u].cal_min,cal_max:this.volumes[u].cal_max,robust_min:this.volumes[u].robust_min,robust_max:this.volumes[u].robust_max}}refreshLayers(u,m){var Y;if(this.refreshColormaps(),this.volumes.length<1)return;const A=u.hdr;let S=u.img;u.frame4D>0&&u.frame4D<u.nFrame4D&&(S=u.img.slice(u.frame4D*u.nVox3D,(u.frame4D+1)*u.nVox3D));const v=u.opacity;if(m>1&&v===0)return;let w=null;if(!this.back)throw new Error("back undefined");this.gl.bindVertexArray(this.unusedVAO),this.crosshairs3D&&(this.crosshairs3D.mm[0]=NaN);let D=clone$2(u.toRAS);if(m===0){this.volumeObject3D=u.toNiivueObject3D(this.VOLUME_ID,this.gl),invert(D,D),this.back.matRAS=u.matRAS,this.back.dims=u.dimsRAS,this.back.pixDims=u.pixDimsRAS,w=this.rgbaTex(this.volumeTexture,TEXTURE0_BACK_VOL,u.dimsRAS);const{volScale:G,vox:J}=this.sliceScale(!0);if(this.volScale=G,this.vox=J,this.volumeObject3D.scale=G,!this.renderShader)throw new Error("renderShader undefined");this.renderShader.use(this.gl),this.gl.uniform3fv(this.renderShader.uniforms.texVox,J),this.gl.uniform3fv(this.renderShader.uniforms.volScale,G);const i0=this.pickingImageShader;i0.use(this.gl),this.gl.uniform1i(i0.uniforms.volume,0),this.gl.uniform1i(i0.uniforms.colormap,1),this.gl.uniform1i(i0.uniforms.overlay,2),this.gl.uniform3fv(i0.uniforms.volScale,G),log.debug(this.volumeObject3D)}else{((Y=this.back)==null?void 0:Y.dims)===void 0&&log.error("Fatal error: Unable to render overlay: background dimensions not defined!");const G=this.mm2frac(u.mm000,0,!0);let J=this.mm2frac(u.mm100,0,!0),i0=this.mm2frac(u.mm010,0,!0),Q=this.mm2frac(u.mm001,0,!0);J=subtract$1(J,J,G),i0=subtract$1(i0,i0,G),Q=subtract$1(Q,Q,G),D=fromValues$3(J[0],i0[0],Q[0],G[0],J[1],i0[1],Q[1],G[1],J[2],i0[2],Q[2],G[2],0,0,0,1),invert(D,D),m===1?(w=this.rgbaTex(this.overlayTexture,TEXTURE2_OVERLAY_VOL,this.back.dims),this.overlayTexture=w,this.overlayTextureID=w):w=this.overlayTextureID}const R=this.gl.createFramebuffer();this.gl.bindFramebuffer(this.gl.FRAMEBUFFER,R),this.gl.disable(this.gl.CULL_FACE),this.gl.viewport(0,0,this.back.dims[1],this.back.dims[2]),this.gl.disable(this.gl.BLEND);const P=this.gl.createTexture();this.gl.activeTexture(TEXTURE9_ORIENT),this.gl.bindTexture(this.gl.TEXTURE_3D,P),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_MIN_FILTER,this.gl.NEAREST),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_MAG_FILTER,this.gl.NEAREST),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_WRAP_R,this.gl.CLAMP_TO_EDGE),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_WRAP_S,this.gl.CLAMP_TO_EDGE),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_WRAP_T,this.gl.CLAMP_TO_EDGE),this.gl.pixelStorei(this.gl.UNPACK_ALIGNMENT,1);let L=this.orientShaderU;if(!A)throw new Error("hdr undefined");if(!S)throw new Error("img undefined");if(A.datatypeCode===2)A.intent_code===1002&&(L=this.orientShaderAtlasU),this.gl.texStorage3D(this.gl.TEXTURE_3D,1,this.gl.R8UI,A.dims[1],A.dims[2],A.dims[3]),this.gl.texSubImage3D(this.gl.TEXTURE_3D,0,0,0,0,A.dims[1],A.dims[2],A.dims[3],this.gl.RED_INTEGER,this.gl.UNSIGNED_BYTE,S);else if(A.datatypeCode===4)L=this.orientShaderI,A.intent_code===1002&&(L=this.orientShaderAtlasI),this.gl.texStorage3D(this.gl.TEXTURE_3D,1,this.gl.R16I,A.dims[1],A.dims[2],A.dims[3]),this.gl.texSubImage3D(this.gl.TEXTURE_3D,0,0,0,0,A.dims[1],A.dims[2],A.dims[3],this.gl.RED_INTEGER,this.gl.SHORT,S);else if(A.datatypeCode===16)this.gl.texStorage3D(this.gl.TEXTURE_3D,1,this.gl.R32F,A.dims[1],A.dims[2],A.dims[3]),this.gl.texSubImage3D(this.gl.TEXTURE_3D,0,0,0,0,A.dims[1],A.dims[2],A.dims[3],this.gl.RED,this.gl.FLOAT,S),L=this.orientShaderF;else if(A.datatypeCode===64){let G=new Float32Array;G=Float32Array.from(S),this.gl.texStorage3D(this.gl.TEXTURE_3D,1,this.gl.R32F,A.dims[1],A.dims[2],A.dims[3]),this.gl.texSubImage3D(this.gl.TEXTURE_3D,0,0,0,0,A.dims[1],A.dims[2],A.dims[3],this.gl.RED,this.gl.FLOAT,G),L=this.orientShaderF}else A.datatypeCode===128?(L=this.orientShaderRGBU,L.use(this.gl),this.gl.uniform1i(L.uniforms.hasAlpha,0),this.gl.texStorage3D(this.gl.TEXTURE_3D,1,this.gl.RGB8UI,A.dims[1],A.dims[2],A.dims[3]),this.gl.texSubImage3D(this.gl.TEXTURE_3D,0,0,0,0,A.dims[1],A.dims[2],A.dims[3],this.gl.RGB_INTEGER,this.gl.UNSIGNED_BYTE,S)):A.datatypeCode===512?(A.intent_code===1002&&(L=this.orientShaderAtlasU),this.gl.texStorage3D(this.gl.TEXTURE_3D,1,this.gl.R16UI,A.dims[1],A.dims[2],A.dims[3]),this.gl.texSubImage3D(this.gl.TEXTURE_3D,0,0,0,0,A.dims[1],A.dims[2],A.dims[3],this.gl.RED_INTEGER,this.gl.UNSIGNED_SHORT,S)):A.datatypeCode===2304&&(L=this.orientShaderRGBU,L.use(this.gl),this.gl.uniform1i(L.uniforms.hasAlpha,1),this.gl.texStorage3D(this.gl.TEXTURE_3D,1,this.gl.RGBA8UI,A.dims[1],A.dims[2],A.dims[3]),this.gl.texSubImage3D(this.gl.TEXTURE_3D,0,0,0,0,A.dims[1],A.dims[2],A.dims[3],this.gl.RGBA_INTEGER,this.gl.UNSIGNED_BYTE,S));u.global_min===void 0&&u.calMinMax();let d=null;if(this.gl.bindVertexArray(this.genericVAO),m>1){if(!this.back.dims)throw new Error("back.dims undefined");d=this.rgbaTex(d,TEXTURE10_BLEND,this.back.dims),this.gl.bindTexture(this.gl.TEXTURE_3D,d);const G=this.passThroughShader;G.use(this.gl),this.gl.uniform1i(G.uniforms.in3D,2);for(let J=0;J<this.back.dims[3];J++){const i0=1/this.back.dims[3]*(J+.5);this.gl.uniform1f(G.uniforms.coordZ,i0),this.gl.framebufferTextureLayer(this.gl.FRAMEBUFFER,this.gl.COLOR_ATTACHMENT0,d,0,J),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4)}}else d=this.rgbaTex(d,TEXTURE10_BLEND,[2,2,2,2]);L.use(this.gl),this.gl.activeTexture(TEXTURE1_COLORMAPS);let E=null;if(u.colormapLabel!==null&&u.colormapLabel.lut.length>7){const G=u.colormapLabel.max-u.colormapLabel.min+1;E=this.createColormapTexture(E,1,G),this.gl.texSubImage2D(this.gl.TEXTURE_2D,0,0,0,G,1,this.gl.RGBA,this.gl.UNSIGNED_BYTE,u.colormapLabel.lut),this.gl.uniform1f(L.uniforms.cal_min,u.colormapLabel.min-.5),this.gl.uniform1f(L.uniforms.cal_max,u.colormapLabel.max+.5),this.gl.bindTexture(this.gl.TEXTURE_2D,E)}else this.gl.bindTexture(this.gl.TEXTURE_2D,this.colormapTexture),this.gl.uniform1f(L.uniforms.cal_min,u.cal_min),this.gl.uniform1f(L.uniforms.cal_max,u.cal_max);this.gl.uniform1i(L.uniforms.isAlphaThreshold,u.alphaThreshold),this.gl.uniform1i(L.uniforms.isAdditiveBlend,this.opts.isAdditiveBlend?1:0);let e=Number.POSITIVE_INFINITY,N=Number.NEGATIVE_INFINITY;if(u.colormapNegative.length>0&&(e=Math.min(-u.cal_min,-u.cal_max),N=Math.max(-u.cal_min,-u.cal_max),isFinite(u.cal_minNeg)&&isFinite(u.cal_maxNeg)&&(e=Math.min(u.cal_minNeg,u.cal_maxNeg),N=Math.max(u.cal_minNeg,u.cal_maxNeg))),!L)throw new Error("orientShader undefined");this.gl.uniform1f(L.uniforms.layer??null,m),this.gl.uniform1f(L.uniforms.cal_minNeg??null,e),this.gl.uniform1f(L.uniforms.cal_maxNeg??null,N),this.gl.bindTexture(this.gl.TEXTURE_3D,P),this.gl.uniform1i(L.uniforms.intensityVol??null,9),this.gl.uniform1i(L.uniforms.blend3D??null,10),this.gl.uniform1i(L.uniforms.colormap??null,1),this.gl.uniform1f(L.uniforms.scl_inter??null,A.scl_inter),this.gl.uniform1f(L.uniforms.scl_slope??null,A.scl_slope),this.gl.uniform1f(L.uniforms.opacity??null,v),this.gl.uniform1i(L.uniforms.modulationVol??null,7);let q=null;if(u.modulationImage!==null&&u.modulationImage>=0&&u.modulationImage<this.volumes.length){log.debug("modulating",this.volumes);const G=this.volumes[u.modulationImage].hdr;if(G.dims[1]===A.dims[1]&&G.dims[2]===A.dims[2]&&G.dims[3]===A.dims[3]){u.modulateAlpha?(this.gl.uniform1i(L.uniforms.modulation,2),this.gl.uniform1f(L.uniforms.opacity,1)):this.gl.uniform1i(L.uniforms.modulation,1),q=this.r8Tex(q,this.gl.TEXTURE7,A.dims,!0),this.gl.activeTexture(this.gl.TEXTURE7),this.gl.bindTexture(this.gl.TEXTURE_3D,q);const J=A.dims[1]*A.dims[2]*A.dims[3],i0=new Uint8Array(J),Q=this.volumes[u.modulationImage].cal_min,f=1/(this.volumes[u.modulationImage].cal_max-Q),o0=this.volumes[u.modulationImage].img.buffer;let l0=new Uint8Array(o0);switch(G.datatypeCode){case 4:l0=new Int16Array(o0);break;case 16:l0=new Float32Array(o0);break;case 64:l0=new Float64Array(o0);break;case 128:l0=new Uint8Array(o0);break;case 512:l0=new Uint16Array(o0);break}log.debug(this.volumes[u.modulationImage]);const c0=this.volumes[u.modulationImage].colormapNegative.length>0;let d0=this.volumes[u.modulationImage].cal_min,S0=this.volumes[u.modulationImage].cal_max;isFinite(this.volumes[u.modulationImage].cal_minNeg)&&isFinite(this.volumes[u.modulationImage].cal_maxNeg)&&(d0=this.volumes[u.modulationImage].cal_minNeg,S0=this.volumes[u.modulationImage].cal_minNeg),d0=Math.abs(d0),S0=Math.abs(S0),d0>S0&&([d0,S0]=[S0,d0]);const p0=1/(S0-d0);let T0=Math.abs(u.modulateAlpha);T0=Math.max(T0,1);const x0=this.volumes[u.modulationImage].frame4D*J;for(let M0=0;M0<J;M0++){const ee=l0[M0+x0]*G.scl_slope+G.scl_inter;let ie=(ee-Q)*f;c0&&ee<0&&(ie=(Math.abs(ee)-d0)*p0),ie=Math.min(Math.max(ie,0),1),ie=Math.pow(ie,T0)*255,i0[M0]=ie}this.gl.texSubImage3D(this.gl.TEXTURE_3D,0,0,0,0,A.dims[1],A.dims[2],A.dims[3],this.gl.RED,this.gl.UNSIGNED_BYTE,i0)}else log.debug("Modulation image dimensions do not match target")}else this.gl.uniform1i(L.uniforms.modulation,0);if(this.gl.uniformMatrix4fv(L.uniforms.mtx,!1,D),!this.back.dims)throw new Error("back.dims undefined");let z=0;A.intent_code===1002&&(z=this.opts.atlasOutline),this.gl.uniform4fv(L.uniforms.xyzaFrac,[1/this.back.dims[1],1/this.back.dims[2],1/this.back.dims[3],z]),log.debug("back dims: ",this.back.dims);for(let G=0;G<this.back.dims[3];G++){const J=1/this.back.dims[3]*(G+.5);this.gl.uniform1f(L.uniforms.coordZ,J),this.gl.framebufferTextureLayer(this.gl.FRAMEBUFFER,this.gl.COLOR_ATTACHMENT0,w,0,G),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4)}if(this.gl.bindVertexArray(this.unusedVAO),this.gl.deleteTexture(P),this.gl.deleteTexture(q),this.gl.deleteTexture(d),this.gl.viewport(0,0,this.gl.canvas.width,this.gl.canvas.height),this.gl.bindFramebuffer(this.gl.FRAMEBUFFER,null),this.gl.deleteFramebuffer(R),m===0&&(this.volumeTexture=w,this.gradientTextureAmount>0?this.gradientGL(A):(this.gradientTexture!==null&&this.gl.deleteTexture(this.gradientTexture),this.gradientTexture=null)),!this.renderShader)throw new Error("renderShader undefined");this.renderShader.use(this.gl);const e0=this.sliceScale(!0),H=e0.vox,j=e0.volScale;if(this.gl.uniform1f(this.renderShader.uniforms.overlays,this.overlays),this.gl.uniform4fv(this.renderShader.uniforms.clipPlaneColor,this.opts.clipPlaneColor),this.gl.uniform1f(this.renderShader.uniforms.backOpacity,this.volumes[0].opacity),this.gl.uniform1f(this.renderShader.uniforms.renderOverlayBlend,this.opts.renderOverlayBlend),this.gl.uniform4fv(this.renderShader.uniforms.clipPlane,this.scene.clipPlane),this.gl.uniform3fv(this.renderShader.uniforms.texVox,H),this.gl.uniform3fv(this.renderShader.uniforms.volScale,j),!this.pickingImageShader)throw new Error("pickingImageShader undefined");this.pickingImageShader.use(this.gl),this.gl.uniform1f(this.pickingImageShader.uniforms.overlays,this.overlays.length),this.gl.uniform3fv(this.pickingImageShader.uniforms.texVox,H);let y=this.sliceMMShader;if(this.opts.isV1SliceShader&&(y=this.sliceV1Shader),!y)throw new Error("slice shader undefined");y.use(this.gl),this.gl.uniform1f(y.uniforms.overlays,this.overlays.length),this.gl.uniform1f(y.uniforms.drawOpacity,this.drawOpacity),E!==null&&(this.gl.deleteTexture(E),this.gl.activeTexture(TEXTURE1_COLORMAPS),this.gl.bindTexture(this.gl.TEXTURE_2D,this.colormapTexture)),this.gl.uniform1i(y.uniforms.drawing,7),this.gl.activeTexture(TEXTURE7_DRAW),this.gl.bindTexture(this.gl.TEXTURE_3D,this.drawTexture),this.updateInterpolation(m)}colormaps(){return cmapper.colormaps()}addColormap(u,m){cmapper.addColormap(u,m)}setColormap(u,m){const A=this.getVolumeIndexByID(u);this.volumes[A].colormap=m,this.updateGLVolume()}idx(u,m,A,S){return A*S[0]*S[1]+m*S[0]+u}check_previous_slice(u,m,A,S,v,w,D,R){const P=new Uint32Array(27);let L=0;if(!v)return 0;const d=u[this.idx(A,S,v,w)];if(D>=6){const E=this.idx(A,S,v-1,w);d===u[E]&&(P[L++]=m[E])}if(D>=18){if(A){const E=this.idx(A-1,S,v-1,w);d===u[E]&&(P[L++]=m[E])}if(S){const E=this.idx(A,S-1,v-1,w);d===u[E]&&(P[L++]=m[E])}if(A<w[0]-1){const E=this.idx(A+1,S,v-1,w);d===u[E]&&(P[L++]=m[E])}if(S<w[1]-1){const E=this.idx(A,S+1,v-1,w);d===u[E]&&(P[L++]=m[E])}}if(D===26){if(A&&S){const E=this.idx(A-1,S-1,v-1,w);d===u[E]&&(P[L++]=m[E])}if(A<w[0]-1&&S){const E=this.idx(A+1,S-1,v-1,w);d===u[E]&&(P[L++]=m[E])}if(A&&S<w[1]-1){const E=this.idx(A-1,S+1,v-1,w);d===u[E]&&(P[L++]=m[E])}if(A<w[0]-1&&S<w[1]-1){const E=this.idx(A+1,S+1,v-1,w);d===u[E]&&(P[L++]=m[E])}}return L?(this.fill_tratab(R,P,L),P[0]):0}do_initial_labelling(u,m,A){let S=1;const v=8192;let w=v,D=new Uint32Array(w).fill(0);const R=new Uint32Array(m[0]*m[1]*m[2]).fill(0),P=new Uint32Array(27);for(let L=0;L<m[2];L++)for(let d=0;d<m[1];d++)for(let E=0;E<m[0];E++){let e=0;const N=u[this.idx(E,d,L,m)];if(N!==0){if(P[0]=this.check_previous_slice(u,R,E,d,L,m,A,D),P[0]&&(e+=1),A>=6){if(E){const q=this.idx(E-1,d,L,m);N===u[q]&&(P[e++]=R[q])}if(d){const q=this.idx(E,d-1,L,m);N===u[q]&&(P[e++]=R[q])}}if(A>=18){if(d&&E){const q=this.idx(E-1,d-1,L,m);N===u[q]&&(P[e++]=R[q])}if(d&&E<m[0]-1){const q=this.idx(E+1,d-1,L,m);N===u[q]&&(P[e++]=R[q])}}if(e)R[this.idx(E,d,L,m)]=P[0],this.fill_tratab(D,P,e);else{if(R[this.idx(E,d,L,m)]=S,S>=w){w+=v;const q=new Uint32Array(w);q.set(D),D=q}D[S-1]=S,S++}}}for(let L=0;L<S-1;L++){let d=L;for(;D[d]!==d+1;)d=D[d]-1;D[L]=d+1}return[S-1,D,R]}fill_tratab(u,m,A){let S=0;const v=new Uint32Array(A+5).fill(0);let D=2147483647;for(let R=0;R<A;R++){let P=m[R];for(S=0;u[P-1]!==P;)if(P=u[P-1],S++,S>100){log.info(`
-Ooh no!!`);break}v[R]=P,D=Math.min(D,P)}for(let R=0;R<A;R++)u[v[R]-1]=D}translate_labels(u,m,A,S){const v=m[0]*m[1]*m[2];let w=0;const D=new Uint32Array(v).fill(0);for(let L=0;L<S;L++)w=Math.max(w,A[L]);const R=new Uint32Array(w).fill(0);let P=0;for(let L=0;L<v;L++)u[L]&&(R[A[u[L]-1]-1]||(P+=1,R[A[u[L]-1]-1]=P),D[L]=R[A[u[L]-1]-1]);return[P,D]}largest_original_cluster_labels(u,m,A){const S=u.length,v=new Uint32Array(m+1).fill(0),w=new Uint32Array(m+1).fill(0);for(let P=0;P<S;P++){const L=u[P],d=A[P];v[d]=L,w[d]++}let D=0;for(let P=0;P<m+1;P++){const L=v[P];D=Math.max(D,L);for(let d=0;d<m+1;d++)d!==P&&L===v[d]&&(w[P]<w[d]||w[P]===w[d]&&P<d)&&(v[P]=0)}const R=new Uint32Array(S).fill(0);for(let P=0;P<S;P++)R[P]=v[A[P]];return[D,R]}bwlabel(u,m,A=26,S=!1,v=!1){const w=Date.now(),D=m[0]*m[1]*m[2],R=new Uint32Array(D).fill(0);if(![6,18,26].includes(A))return log.info("bwlabel: conn must be 6, 18 or 26."),[0,R];if(m[0]<2||m[1]<2||m[2]<1)return log.info("bwlabel: img must be 2 or 3-dimensional"),[0,R];if(S)for(let N=0;N<D;N++)u[N]!==0&&(R[N]=1);else R.set(u);let[P,L,d]=this.do_initial_labelling(R,m,A);L===void 0&&(L=new Uint32Array);const[E,e]=this.translate_labels(d,m,L,P);if(log.info(A+" neighbor clustering into "+E+" regions in "+(Date.now()-w)+"ms"),v){const[N,q]=this.largest_original_cluster_labels(R,E,e);return[N,q]}return[E,e]}async createConnectedLabelImage(u,m=26,A=!1,S=!1){var d,E;const v=this.getVolumeIndexByID(u),w=Uint32Array.from(((d=this.volumes[v].dims)==null?void 0:d.slice(1,4))??[]),D=Uint32Array.from(((E=this.volumes[v].img)==null?void 0:E.slice())??[]),[R,P]=this.bwlabel(D,w,m,A,S),L=this.volumes[v].clone();L.opacity=.5,L.colormap="random";for(let e=0;e<L.img.length;e++)L.img[e]=P[e];return L.cal_min=0,L.cal_max=R,L}async scalecropUint8(u,m=0,A=255,S,v){const w=u.length,D=new Uint8Array(w);for(let R=0;R<w;R++){let P=u[R];P=m+v*(P-S),P=Math.max(P,m),P=Math.min(P,A),D[R]=P}return D}async scalecropFloat32(u,m=0,A=1,S,v){const w=u.length,D=new Float32Array(w);for(let R=0;R<w;R++){let P=u[R];P=m+v*(P-S),P=Math.max(P,m),P=Math.min(P,A),D[R]=P}return D}getScale(u,m=0,A=255,S=0,v=.999){let w=u.global_min,D=u.global_max;if(u.hdr.datatypeCode===2)return[w,1];if((!isFinite(S)||!isFinite(v))&&isFinite(u.cal_min)&&isFinite(u.cal_max)&&u.cal_max>u.cal_min){w=u.cal_min,D=u.cal_max;const j=(A-m)/(D-w);return log.info(" Robust Rescale:  min: "+w+"  max: "+D+" scale: "+j),console.log("Robust Rescale:  min: "+w+"  max: "+D+" scale: "+j),[w,j]}const R=u.img,P=u.hdr.dims[1]*u.hdr.dims[2]*u.hdr.dims[3];if(u.hdr.scl_slope!==1||u.hdr.scl_inter!==0){const j=u.img,y=new Float32Array(u.img.length);for(let Y=0;Y<P;Y++)y[Y]=j[Y]*u.hdr.scl_slope+u.hdr.scl_inter}if(w<0&&log.warn("WARNING: Input image has value(s) below 0.0 !"),log.info(" Input:    min: "+w+"  max: "+D),S===0&&v===1)return[w,1];let L=0;for(let j=0;j<P;j++)Math.abs(R[j])>=1e-15&&L++;const d=1e3,E=(D-w)/d,e=new Array(d).fill(0);for(let j=0;j<P;j++){const y=R[j];let Y=Math.floor((y-w)/E);Y=Math.min(Y,d-1),e[Y]++}const N=new Array(d).fill(0);N[0]=e[0];for(let j=1;j<d;j++)N[j]=N[j-1]+e[j];let q=Math.floor(S*P),z=0;for(;z<d&&!(N[z]>=q);)z++;const e0=w;for(w=z*E+e0,q=P-Math.floor((1-v)*L),z=0;z<d-1&&!(N[z+1]>=q);)z++;D=z*E+e0;let H=1;return w!==D&&(H=(A-m)/(D-w)),log.info(" Rescale:  min: "+w+"  max: "+D+" scale: "+H),[w,H]}conformVox2Vox(u,m,A=256,S=1,v=!1){const w=m.flat(),D=fromValues$3(w[0],w[1],w[2],w[3],w[4],w[5],w[6],w[7],w[8],w[9],w[10],w[11],w[12],w[13],w[14],w[15]),R=fromValues$1(u[1]/2,u[2]/2,u[3]/2,1),P=create$1(),L=create$3();transpose(L,D),transformMat4(P,R,L);const d=fromValues$2(P[0],P[1],P[2]),E=fromValues$2(S,S,S);let e=fromValues$3(-1,0,0,0,0,0,1,0,0,-1,0,0,0,0,0,1);v&&(e=fromValues$3(1,0,0,0,0,1,0,0,0,0,1,0,0,0,0,1)),transpose(e,e);const N=fromValues$1(A,A,A,1),q=create$3();scale$3(q,e,E);const z=fromValues$1(N[0],N[1],N[2],1);transformMat4(z,z,q),scale$1(z,z,.5);const e0=create$2();subtract$1(e0,d,fromValues$2(z[0],z[1],z[2]));const H=create$3();transpose(H,q),H[3]=e0[0],H[7]=e0[1],H[11]=e0[2];const j=create$3();invert(j,H);const y=create$3();mul(y,D,j);const Y=create$3();return invert(Y,y),[H,y,Y]}async createNiftiArray(u=[256,256,256],m=[1,1,1],A=[1,0,0,-128,0,1,0,-128,0,0,1,-128,0,0,0,1],S=2,v=new Uint8Array){return await NVImage.createNiftiArray(u,m,A,S,v)}async niftiArray2NVImage(u=new Uint8Array){return await NVImage.loadFromUrl({url:u})}async loadFromUrl(u){return await NVImage.loadFromUrl({url:u})}async conform(u,m=!1,A=!0,S=!1,v=!1){const R=this.conformVox2Vox(u.hdr.dims,u.hdr.affine.flat(),256,1,m),P=R[0],L=R[2],d=256*256*256,E=new Float32Array(d),e=new Float32Array(u.img),N=u.hdr.dims[1]*u.hdr.dims[2]*u.hdr.dims[3];if(u.hdr.scl_slope!==1||u.hdr.scl_inter!==0)for(let o0=0;o0<N;o0++)e[o0]=e[o0]*u.hdr.scl_slope+u.hdr.scl_inter;const q=u.hdr.dims[1],z=u.hdr.dims[2],e0=u.hdr.dims[3],H=q*z;let j=-1;function y(o0,l0,c0){return o0+l0*q+c0*H}const Y=L[0],G=L[4],J=L[8];if(A)for(let o0=0;o0<256;o0++)for(let l0=0;l0<256;l0++){const c0=l0*L[1]+o0*L[2]+L[3],d0=l0*L[5]+o0*L[6]+L[7],S0=l0*L[9]+o0*L[10]+L[11];for(let p0=0;p0<256;p0++){const T0=p0*Y+c0,x0=p0*G+d0,M0=p0*J+S0,ee=Math.floor(T0),ie=Math.floor(x0),Fe=Math.floor(M0);if(j++,ee<0||ie<0||Fe<0)continue;const Oe=Math.ceil(T0),Ee=Math.ceil(x0),I0=Math.ceil(M0);if(Oe>=q||Ee>=z||I0>=e0)continue;const Be=T0-ee,Ge=x0-ie,ti=M0-Fe,ai=1-Be,ri=1-Ge,qe=1-ti,be=y(ee,ie,Fe);let J0=0;J0+=e[be]*ai*ri*qe,J0+=e[be+H]*ai*ri*ti,J0+=e[be+q]*ai*Ge*qe,J0+=e[be+q+H]*ai*Ge*ti,J0+=e[be+1]*Be*ri*qe,J0+=e[be+1+H]*Be*ri*ti,J0+=e[be+1+q]*Be*Ge*qe,J0+=e[be+1+q+H]*Be*Ge*ti,E[j]=J0}}else for(let o0=0;o0<256;o0++)for(let l0=0;l0<256;l0++){const c0=l0*L[1]+o0*L[2]+L[3],d0=l0*L[5]+o0*L[6]+L[7],S0=l0*L[9]+o0*L[10]+L[11];for(let p0=0;p0<256;p0++){const T0=Math.round(p0*Y+c0),x0=Math.round(p0*G+d0),M0=Math.round(p0*J+S0);j++,!(T0<0||x0<0||M0<0)&&(T0>=q||x0>=z||M0>=e0||(E[j]=e[y(T0,x0,M0)]))}}let i0=0;v&&(i0=NaN);let Q=new Uint8Array;if(S){const o0=await this.getScale(u,0,1,i0),l0=await this.scalecropFloat32(E,0,1,o0[0],o0[1]);Q=await this.createNiftiArray([256,256,256],[1,1,1],Array.from(P),16,new Uint8Array(l0.buffer))}else{const o0=await this.getScale(u,0,255,i0),l0=await this.scalecropUint8(E,0,255,o0[0],o0[1]);Q=await this.createNiftiArray([256,256,256],[1,1,1],Array.from(P),2,l0)}return await this.niftiArray2NVImage(Q)}setRenderDrawAmbientOcclusion(u){if(!this.renderShader)throw new Error("renderShader undefined");this.renderDrawAmbientOcclusion=u,this.renderShader.use(this.gl),this.gl.uniform1fv(this.renderShader.uniforms.renderDrawAmbientOcclusion,[this.renderDrawAmbientOcclusion,1]),this.drawScene()}setColorMap(u,m){this.setColormap(u,m)}setColormapNegative(u,m){const A=this.getVolumeIndexByID(u);this.volumes[A].colormapNegative=m,this.updateGLVolume()}setModulationImage(u,m,A=0){const S=this.getVolumeIndexByID(u);let v=null;m.length>0&&(v=this.getVolumeIndexByID(m)),this.volumes[S].modulationImage=v,this.volumes[S].modulateAlpha=A,this.updateGLVolume()}setGamma(u=1){cmapper.gamma=u,this.updateGLVolume()}async loadDeferred4DVolumes(u){const m=this.getVolumeIndexByID(u),A=this.volumes[m];if(A.nTotalFrame4D<=A.nFrame4D)return;let S;A.fileObject?S=await NVImage.loadFromFile({file:A.fileObject}):S=await NVImage.loadFromUrl({url:A.url}),S&&(A.img=S.img.slice(),A.nTotalFrame4D=S.nTotalFrame4D,A.nFrame4D=S.nFrame4D,this.updateGLVolume())}setFrame4D(u,m){const A=this.getVolumeIndexByID(u),S=this.volumes[A];m>S.nFrame4D-1&&(m=S.nFrame4D-1),m<0&&(m=0),m!==S.frame4D&&(S.frame4D=m,this.updateGLVolume(),this.onFrameChange(S,m),this.createOnLocationChange())}getFrame4D(u){const m=this.getVolumeIndexByID(u);return this.volumes[m].frame4D}colormapFromKey(u){return cmapper.colormapFromKey(u)}colormap(u="",m=!1){return cmapper.colormap(u,m)}createColormapTexture(u=null,m=0,A=256){return u!==null&&this.gl.deleteTexture(u),m<1||A<1?null:(u=this.gl.createTexture(),this.gl.activeTexture(TEXTURE1_COLORMAPS),this.gl.bindTexture(this.gl.TEXTURE_2D,u),this.gl.texStorage2D(this.gl.TEXTURE_2D,1,this.gl.RGBA8,A,m),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_MIN_FILTER,this.gl.LINEAR),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_MAG_FILTER,this.gl.LINEAR),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_WRAP_R,this.gl.CLAMP_TO_EDGE),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_WRAP_S,this.gl.CLAMP_TO_EDGE),this.gl.pixelStorei(this.gl.UNPACK_ALIGNMENT,1),u)}addColormapList(u="",m=NaN,A=NaN,S=!1,v=!1,w=!0,D=!1){u.length<1&&(w=!1),this.colormapLists.push({name:u,min:m,max:A,alphaThreshold:S,negative:v,visible:w,invert:D})}refreshColormaps(){if(this.colormapLists=[],this.volumes.length<1&&this.meshes.length<1)return;const u=this.volumes.length;if(u>0)for(let w=0;w<u;w++){const D=this.volumes[w],R=negMinMax(D.cal_min,D.cal_max,D.cal_minNeg,D.cal_maxNeg);this.addColormapList(D.colormapNegative,R[0],R[1],D.alphaThreshold!==void 0,!0,D.colorbarVisible,D.colormapInvert),this.addColormapList(D.colormap,D.cal_min,D.cal_max,D.alphaThreshold!==void 0,!1,D.colorbarVisible,D.colormapInvert)}const m=this.meshes.length;if(m>0)for(let w=0;w<m;w++){const D=this.meshes[w];if(!D.colorbarVisible)continue;const R=D.layers.length;if("edgeColormap"in D&&"edges"in D&&D.edges!==void 0){const P=negMinMax(D.edgeMin,D.edgeMax,NaN,NaN);this.addColormapList(D.edgeColormapNegative,P[0],P[1],!1,!0,!0,D.colormapInvert),this.addColormapList(D.edgeColormap,D.edgeMin,D.edgeMax,!1,!1,!0,D.colormapInvert)}if(!(R<1))for(let P=0;P<R;P++){const L=this.meshes[w].layers[P];if(!L.colorbarVisible||L.colormap.length<1)continue;const d=negMinMax(L.cal_min,L.cal_max,L.cal_minNeg,L.cal_maxNeg);this.addColormapList(L.colormapNegative,d[0],d[1],L.alphaThreshold,!0,!0,L.colormapInvert),this.addColormapList(L.colormap,L.cal_min,L.cal_max,L.alphaThreshold,!1,!0,L.colormapInvert)}}const A=this.colormapLists.length;if(A<1)return;this.colormapTexture=this.createColormapTexture(this.colormapTexture,A+1);let S=new Uint8ClampedArray;function v(w){const D=new Uint8ClampedArray(S.length+w.length);D.set(S),D.set(w,S.length),S=D}for(let w=0;w<A;w++)v(Array.from(this.colormap(this.colormapLists[w].name,this.colormapLists[w].invert)));return v(Array.from(this.drawLut.lut)),this.gl.texSubImage2D(this.gl.TEXTURE_2D,0,0,0,256,A+1,this.gl.RGBA,this.gl.UNSIGNED_BYTE,S),this}sliceScale(u=!1){var w;let m=this.screenFieldOfViewMM(0);u&&(m=this.screenFieldOfViewVox(0));const A=Math.max(m[0],Math.max(m[1],m[2])),S=[m[0]/A,m[1]/A,m[2]/A];if(!((w=this.back)!=null&&w.dims))throw new Error("back.dims undefined");const v=[this.back.dims[1],this.back.dims[2],this.back.dims[3]];return{volScale:S,vox:v,longestAxis:A,dimsMM:m}}tileIndex(u,m){for(let A=0;A<this.screenSlices.length;A++){const S=this.screenSlices[A].leftTopWidthHeight;if(u>S[0]&&m>S[1]&&u<S[0]+S[2]&&m<S[1]+S[3])return A}return-1}inRenderTile(u,m){const A=this.tileIndex(u,m);return A>=0&&this.screenSlices[A].axCorSag===4?A:-1}sliceScroll3D(u=0){if(u!==0){if(this.volumes.length>0&&this.scene.clipPlaneDepthAziElev[0]<1.8){const m=this.scene.clipPlaneDepthAziElev.slice();return u>0&&(m[0]=Math.min(1.5,m[0]+.025)),u<0&&(m[0]=Math.max(-1.5,m[0]-.025)),m[0]!==this.scene.clipPlaneDepthAziElev[0]?(this.scene.clipPlaneDepthAziElev=m,this.setClipPlane(this.scene.clipPlaneDepthAziElev)):void 0}u>0&&(this.scene.volScaleMultiplier=Math.min(2,this.scene.volScaleMultiplier*1.1)),u<0&&(this.scene.volScaleMultiplier=Math.max(.5,this.scene.volScaleMultiplier*.9)),this.drawScene()}}deleteThumbnail(){this.bmpTexture&&(this.gl.deleteTexture(this.bmpTexture),this.bmpTexture=null,this.thumbnailVisible=!1)}inGraphTile(u,m){if(this.graph.opacity<=0||this.volumes.length<1||this.volumes[0].nFrame4D<1||!this.graph.plotLTWH||this.graph.plotLTWH[2]<1||this.graph.plotLTWH[3]<1)return!1;const A=[(u-this.graph.LTWH[0])/this.graph.LTWH[2],(m-this.graph.LTWH[1])/this.graph.LTWH[3]];return A[0]>0&&A[1]>0&&A[0]<=1&&A[1]<=1}mouseClick(u,m,A=0,S=!0){if(u*=this.uiData.dpr,m*=this.uiData.dpr,this.canvas.focus(),this.thumbnailVisible){this.thumbnailVisible=!1,Promise.all([this.loadVolumes(this.deferredVolumes),this.loadMeshes(this.deferredMeshes)]).catch(v=>{throw v});return}if(this.inGraphTile(u,m)){if(!this.graph.plotLTWH)throw new Error("plotLTWH undefined");const v=[(u-this.graph.plotLTWH[0])/this.graph.plotLTWH[2],(m-this.graph.plotLTWH[1])/this.graph.plotLTWH[3]];if(v[0]>0&&v[1]>0&&v[0]<=1&&v[1]<=1){const w=Math.round(v[0]*(this.volumes[0].nFrame4D-1));this.setFrame4D(this.volumes[0].id,w);return}v[0]>.5&&v[1]>1&&this.loadDeferred4DVolumes(this.volumes[0].id).catch(w=>{throw w});return}if(this.inRenderTile(u,m)>=0){this.sliceScroll3D(A),this.drawScene();return}if(!(this.screenSlices.length<1||this.gl.canvas.height<1||this.gl.canvas.width<1))for(let v=0;v<this.screenSlices.length;v++){const w=this.screenSlices[v].axCorSag;if(this.drawPenAxCorSag>=0&&this.drawPenAxCorSag!==w||w>2)continue;const D=this.screenXY2TextureFrac(u,m,v,!1);if(!(D[0]<0)){if(!S){this.scene.crosshairPos[2-w]=A,this.drawScene();return}if(A!==0){let R=1;A<0&&(R=-1);const P=[0,0,0];P[2-w]=R,this.moveCrosshairInVox(P[0],P[1],P[2]),this.drawScene(),this.createOnLocationChange(w);return}if(this.opts.isForceMouseClickToVoxelCenters?this.scene.crosshairPos=clone$1(this.vox2frac(this.frac2vox(D))):this.scene.crosshairPos=clone$1(D),this.opts.drawingEnabled){const R=this.frac2vox(this.scene.crosshairPos);if(!isFinite(this.opts.penValue)||this.opts.penValue<0||Object.is(this.opts.penValue,-0)){isFinite(this.opts.penValue)?this.drawFloodFill(R,Math.abs(this.opts.penValue),this.opts.floodFillNeighbors):this.drawFloodFill(R,0,this.opts.penValue,this.opts.floodFillNeighbors);return}if(isNaN(this.drawPenLocation[0]))this.drawPenAxCorSag=w,this.drawPenFillPts=[],this.drawPt(...R,this.opts.penValue);else{if(R[0]===this.drawPenLocation[0]&&R[1]===this.drawPenLocation[1]&&R[2]===this.drawPenLocation[2])return;this.drawPenLine(R,this.drawPenLocation,this.opts.penValue)}this.drawPenLocation=R,this.opts.isFilledPen&&this.drawPenFillPts.push(R),this.refreshDrawing(!1)}this.drawScene(),this.createOnLocationChange(w);return}}}drawRuler(){let u=[],m=[];for(let R=0;R<this.screenSlices.length;R++)if(this.screenSlices[R].axCorSag!==4&&this.screenSlices[R].fovMM.length>1){m=this.screenSlices[R].leftTopWidthHeight,u=this.screenSlices[R].fovMM;break}if(m.length<4)return;const S=100/u[0]*m[2],v=m[0]+.5*m[2]-.5*S,w=m[1]+m[3]-2*this.opts.rulerWidth,D=[v,w,v+S,w];this.drawRuler10cm(D)}drawRuler10cm(u){if(!this.lineShader)throw new Error("lineShader undefined");this.gl.bindVertexArray(this.genericVAO),this.lineShader.use(this.gl),this.gl.uniform4fv(this.lineShader.uniforms.lineColor,this.opts.rulerColor),this.gl.uniform2fv(this.lineShader.uniforms.canvasWidthHeight,[this.gl.canvas.width,this.gl.canvas.height]),this.gl.uniform1f(this.lineShader.uniforms.thickness,this.opts.rulerWidth),this.gl.uniform4fv(this.lineShader.uniforms.startXYendXY,u),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4);const m=-.1*(u[0]-u[2]),A=u[1],S=A-2*this.opts.rulerWidth,v=A-4*this.opts.rulerWidth;for(let w=0;w<11;w++){const D=u[0]+w*m,R=[D,A,D,S];w%5===0&&(R[3]=v),this.gl.uniform4fv(this.lineShader.uniforms.startXYendXY,R),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4)}this.gl.bindVertexArray(this.unusedVAO)}screenXY2mm(u,m,A=-1){let S;for(let v=0;v<this.screenSlices.length;v++){let w=v;if(A>=0&&(w=A),this.screenSlices[w].axCorSag>2)continue;const R=this.screenSlices[w].leftTopWidthHeight;if(u<R[0]||m<R[1]||u>R[0]+R[2]||m>R[1]+R[3]||(S=this.screenXY2TextureFrac(u,m,w,!1),S[0]<0))continue;const P=this.frac2mm(S);return fromValues$1(P[0],P[1],P[2],w)}return fromValues$1(NaN,NaN,NaN,NaN)}dragForPanZoom(u){const m=this.screenXY2mm(u[2],u[3]);if(isNaN(m[0]))return;const A=this.screenXY2mm(u[0],u[1],m[3]);if(isNaN(A[0])||isNaN(m[0])||isNaN(m[3]))return;const S=create$1(),v=this.uiData.pan2DxyzmmAtMouseDown[3];sub(S,m,A),this.scene.pan2Dxyzmm[0]=this.uiData.pan2DxyzmmAtMouseDown[0]+v*S[0],this.scene.pan2Dxyzmm[1]=this.uiData.pan2DxyzmmAtMouseDown[1]+v*S[1],this.scene.pan2Dxyzmm[2]=this.uiData.pan2DxyzmmAtMouseDown[2]+v*S[2]}dragForCenterButton(u){this.dragForPanZoom(u)}dragForSlicer3D(u){let m=this.uiData.pan2DxyzmmAtMouseDown[3];const A=u[3]-u[1];m+=A*.01,m=Math.max(m,.1),m=Math.min(m,10);const v=this.scene.pan2Dxyzmm[3]-m;this.opts.yoke3Dto2DZoom&&(this.scene.volScaleMultiplier=m),this.scene.pan2Dxyzmm[3]=m;const w=this.frac2mm(this.scene.crosshairPos);this.scene.pan2Dxyzmm[0]+=v*w[0],this.scene.pan2Dxyzmm[1]+=v*w[1],this.scene.pan2Dxyzmm[2]+=v*w[2]}drawMeasurementTool(u){const m=this.gl;if(m.bindVertexArray(this.genericVAO),m.depthFunc(m.ALWAYS),m.enable(m.BLEND),m.blendFunc(m.SRC_ALPHA,m.ONE_MINUS_SRC_ALPHA),!this.lineShader)throw new Error("lineShader undefined");this.lineShader.use(this.gl),m.uniform4fv(this.lineShader.uniforms.lineColor,this.opts.rulerColor),m.uniform2fv(this.lineShader.uniforms.canvasWidthHeight,[m.canvas.width,m.canvas.height]),m.uniform1f(this.lineShader.uniforms.thickness,this.opts.rulerWidth),m.uniform4fv(this.lineShader.uniforms.startXYendXY,u),m.drawArrays(m.TRIANGLE_STRIP,0,4);const A=this.opts.rulerColor;A[3]=1,m.uniform4fv(this.lineShader.uniforms.lineColor,A);const S=this.opts.rulerWidth;m.uniform1f(this.lineShader.uniforms.thickness,S*2);let v=[u[0],u[1]-S,u[0],u[1]+S];m.uniform4fv(this.lineShader.uniforms.startXYendXY,v),m.drawArrays(m.TRIANGLE_STRIP,0,4),v=[u[2],u[3]-S,u[2],u[3]+S],m.uniform4fv(this.lineShader.uniforms.startXYendXY,v),m.drawArrays(m.TRIANGLE_STRIP,0,4);let w=this.canvasPos2frac([u[0],u[1]]),D=this.canvasPos2frac([u[2],u[3]]);if(w[0]>=0&&D[0]>=0){const R=this.frac2mm(w);w=fromValues$2(R[0],R[1],R[2]);const P=this.frac2mm(D);D=fromValues$2(P[0],P[1],P[2]);const L=create$2();sub$1(L,w,D);const d=len(L);let E=2;d>9&&(E=1),d>99&&(E=0);const e=d.toFixed(E);this.drawTextBetween(u,e,1,A)}m.bindVertexArray(this.unusedVAO)}drawRect(u,m=[1,0,0,-1]){if(m[3]<0&&(m=this.opts.crosshairColor),!this.rectShader)throw new Error("rectShader undefined");this.rectShader.use(this.gl),this.gl.enable(this.gl.BLEND),this.gl.uniform4fv(this.rectShader.uniforms.lineColor,m),this.gl.uniform2fv(this.rectShader.uniforms.canvasWidthHeight,[this.gl.canvas.width,this.gl.canvas.height]),this.gl.uniform4f(this.rectShader.uniforms.leftTopWidthHeight,u[0],u[1],u[2],u[3]),this.gl.bindVertexArray(this.genericVAO),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4),this.gl.bindVertexArray(this.unusedVAO)}drawCircle(u,m=this.opts.fontColor,A=1){if(!this.circleShader)throw new Error("circleShader undefined");this.circleShader.use(this.gl),this.gl.enable(this.gl.BLEND),this.gl.uniform4fv(this.circleShader.uniforms.circleColor,m),this.gl.uniform2fv(this.circleShader.uniforms.canvasWidthHeight,[this.gl.canvas.width,this.gl.canvas.height]),this.gl.uniform4f(this.circleShader.uniforms.leftTopWidthHeight,u[0],u[1],u[2],u[3]),this.gl.uniform1f(this.circleShader.uniforms.fillPercent,A),this.gl.uniform4fv(this.circleShader.uniforms.circleColor,m),this.gl.bindVertexArray(this.genericVAO),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4),this.gl.bindVertexArray(this.unusedVAO)}drawSelectionBox(u){this.drawRect(u,this.opts.selectionBoxColor)}effectiveCanvasHeight(){return this.gl.canvas.height-this.colorbarHeight}effectiveCanvasWidth(){return this.gl.canvas.width-this.getLegendPanelWidth()}getAllLabels(){const S=this.meshes.filter(w=>w.type==="connectome").flatMap(w=>w.nodes).map(w=>w.label).filter(w=>w!==void 0);return[...this.document.labels,...S]}getBulletMarginWidth(){let u=0;const m=this.getAllLabels();if(m.length===0)return 0;const A=m.length===1?m[0].style.bulletScale:m.reduce((w,D)=>w.style.bulletScale>D.style.bulletScale?w:D).style.bulletScale,S=m.length===1?m[0]:m.reduce((w,D)=>{const R=this.opts.textHeight*this.gl.canvas.height*w.style.textScale,P=this.opts.textHeight*this.gl.canvas.height*D.style.textScale;return this.textHeight(R,w.text)>this.textHeight(P,D.text)?w:D}),v=this.opts.textHeight*this.gl.canvas.height*S.style.textScale;return u=this.textHeight(v,S.text)*A,u+=v,u}getLegendPanelWidth(){const u=this.getAllLabels();if(!this.opts.showLegend||u.length===0)return 0;const A=this.opts.textHeight*this.gl.canvas.height*1;let S=0;const v=u.reduce((P,L)=>{const d=this.opts.textHeight*this.gl.canvas.height*P.style.textScale,E=this.opts.textHeight*this.gl.canvas.height*L.style.textScale;return this.textWidth(d,P.text)>this.textWidth(E,L.text)?P:L}),w=this.opts.textHeight*this.gl.canvas.height*v.style.textScale,D=this.textWidth(w,v.text),R=this.getBulletMarginWidth();return D&&(S=R+D,S+=A*2),S}getLegendPanelHeight(){const u=this.getAllLabels();let m=0;const S=this.opts.textHeight*this.gl.canvas.height*1;for(const v of u){const w=this.opts.textHeight*this.gl.canvas.height*v.style.textScale,D=this.textHeight(w,v.text);m+=D}return m&&(m+=S/2*(u.length+1)),m}reserveColorbarPanel(){let u=Math.max(this.opts.textHeight,.01);u=u*Math.min(this.gl.canvas.height,this.gl.canvas.width);const m=3*u,A=[0,this.gl.canvas.height-m,this.gl.canvas.width,m];return this.colorbarHeight=A[3]+1,A}drawColorbarCore(u=0,m=[0,0,0,0],A=!1,S=0,v=1,w){if(m[2]<=0||m[3]<=0)return;let D=Math.max(this.opts.textHeight,.01);D=D*Math.min(this.gl.canvas.height,this.gl.canvas.width);let R=D;const P=3*D;let L=D;if(m[3]<P){if(m[3]<3)return;R=1,L=m[3]-2}this.gl.disable(this.gl.DEPTH_TEST),this.colorbarHeight=m[3]+1;const d=[m[0]+R,m[1],m[2]-2*R,L],E=[d[0]-1,d[1]-1,d[2]+2,d[3]+2];if(this.drawRect(E,this.opts.crosshairColor),!this.colorbarShader)throw new Error("colorbarShader undefined");this.colorbarShader.use(this.gl),this.gl.activeTexture(TEXTURE1_COLORMAPS),this.gl.bindTexture(this.gl.TEXTURE_2D,this.colormapTexture),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_MIN_FILTER,this.gl.NEAREST),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_MAG_FILTER,this.gl.NEAREST);const e=u;if(this.gl.uniform1f(this.colorbarShader.uniforms.layer,e),this.gl.uniform2fv(this.colorbarShader.uniforms.canvasWidthHeight,[this.gl.canvas.width,this.gl.canvas.height]),this.gl.disable(this.gl.CULL_FACE),A){const J=[d[0]+d[2],d[1],-d[2],d[3]];this.gl.uniform4fv(this.colorbarShader.uniforms.leftTopWidthHeight,J)}else this.gl.uniform4fv(this.colorbarShader.uniforms.leftTopWidthHeight,d);this.gl.bindVertexArray(this.genericVAO),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4),this.gl.bindVertexArray(this.unusedVAO),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_MIN_FILTER,this.gl.LINEAR),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_MAG_FILTER,this.gl.LINEAR);let N=0;if(w&&v<0&&A?(N=v,v=0):w&&S>0&&(N=S,S=0),S===v||D<1)return;const q=Math.abs(v-S);let[z,e0]=tickSpacing(S,v);e0<S&&(e0+=z);function H(J){return J.toFixed(6).replace(/\.?0*$/,"")}let j=e0;const y=[0,d[1]+d[3]-D*.5,2,D*.75],Y=y[1]+y[3],G=1;for(;j<=v;){y[0]=d[0]+(j-S)/q*d[2],this.drawRect(y);const J=H(G*j);this.drawTextBelow([y[0],Y],J),j+=z}if(N!==0){const J=[d[0]+(N-S)/q*d[2],d[1]-d[3]*.25,2,d[3]*1.5];this.drawRect(J)}}drawColorbar(){const u=this.colormapLists,m=u.length;if(m<1)return;let A=0;for(let R=0;R<m;R++)u[R].visible&&A++;if(A<1)return;let S=this.reserveColorbarPanel(),v=Math.max(this.opts.textHeight,.01);v=v*Math.min(this.gl.canvas.height,this.gl.canvas.width);const w=3*v;let D=S[2]/A;(S[2]<=0||S[3]<=0)&&(D=this.gl.canvas.width/A,S=[0,this.gl.canvas.height-w,D,w]),S[2]=D;for(let R=0;R<m;R++)u[R].visible&&(this.drawColorbarCore(R,S,u[R].negative,u[R].min,u[R].max,u[R].alphaThreshold),S[0]+=D)}textWidth(u,m){if(!m)return 0;let A=0;const S=new TextEncoder().encode(m);for(let v=0;v<m.length;v++)A+=u*this.fontMets.mets[S[v]].xadv;return A}textHeight(u,m){if(!m)return 0;const A=new Set(Array.from(m)),S=new TextEncoder().encode(Array.from(A).join("")),w=Object.values(this.fontMets.mets).filter((D,R)=>S.includes(R)).reduce((D,R)=>D.lbwh[3]>R.lbwh[3]?D:R).lbwh[3];return u*w}drawChar(u,m,A){if(!this.fontShader)throw new Error("fontShader undefined");const S=this.fontMets.mets[A],v=u[0]+m*S.lbwh[0],w=-(m*S.lbwh[1]),D=m*S.lbwh[2],R=m*S.lbwh[3],P=u[1]+(w-R)+m;return this.gl.uniform4f(this.fontShader.uniforms.leftTopWidthHeight,v,P,D,R),this.gl.uniform4fv(this.fontShader.uniforms.uvLeftTopWidthHeight,S.uv_lbwh),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4),m*S.xadv}drawLoadingText(u){if(!this.canvas)throw new Error("canvas undefined");this.gl.viewport(0,0,this.gl.canvas.width,this.gl.canvas.height),this.gl.enable(this.gl.CULL_FACE),this.gl.enable(this.gl.BLEND),this.drawTextBelow([this.canvas.width/2,this.canvas.height/2],u,3)}drawText(u,m,A=1,S=null){if(this.opts.textHeight<=0)return;if(!this.fontShader)throw new Error("fontShader undefined");this.fontShader.use(this.gl);const v=this.opts.textHeight*Math.min(this.gl.canvas.height,this.gl.canvas.width)*A;this.gl.enable(this.gl.BLEND),this.gl.uniform2f(this.fontShader.uniforms.canvasWidthHeight,this.gl.canvas.width,this.gl.canvas.height),S===null&&(S=this.opts.fontColor),this.gl.uniform4fv(this.fontShader.uniforms.fontColor,S);let w=v/this.fontMets.size*this.fontMets.distanceRange;w=Math.max(w,1),this.gl.uniform1f(this.fontShader.uniforms.screenPxRange,w);const D=new TextEncoder().encode(m);this.gl.bindVertexArray(this.genericVAO);for(let R=0;R<m.length;R++)u[0]+=this.drawChar(u,v,D[R]);this.gl.bindVertexArray(this.unusedVAO)}drawTextRight(u,m,A=1,S=null){this.opts.textHeight<=0||(u[1]-=.5*this.opts.textHeight*this.gl.canvas.height,this.drawText(u,m,A,S))}drawTextLeft(u,m,A=1,S=null){if(this.opts.textHeight<=0)return;const v=this.opts.textHeight*this.gl.canvas.height*A;u[0]-=this.textWidth(v,m),u[1]-=.5*v,this.drawText(u,m,A,S)}drawTextRightBelow(u,m,A=1,S=null){this.opts.textHeight<=0||this.drawText(u,m,A,S)}drawTextBetween(u,m,A=1,S=null){if(this.opts.textHeight<=0)return;const v=[(u[0]+u[2])*.5,(u[1]+u[3])*.5],w=this.opts.textHeight*this.gl.canvas.height*A,D=this.textWidth(w,m);v[0]-=.5*D,v[1]-=.5*w;const R=[v[0]-1,v[1]-1,D+2,w+2];let P=S;P===null&&(P=this.opts.crosshairColor),P&&P[0]+P[1]+P[2]>.8?P=[0,0,0,.5]:P=[1,1,1,.5],this.drawRect(R,P),this.drawText(v,m,A,S)}drawTextBelow(u,m,A=1,S=null){if(this.opts.textHeight<=0)return;if(!this.canvas)throw new Error("canvas undefined");let v=this.opts.textHeight*this.gl.canvas.height*A,w=this.textWidth(v,m);w>this.canvas.width&&(A*=(this.canvas.width-2)/w,v=this.opts.textHeight*this.gl.canvas.height*A,w=this.textWidth(v,m)),u[0]-=.5*this.textWidth(v,m),u[0]=Math.max(u[0],1),u[0]=Math.min(u[0],this.canvas.width-w-1),this.drawText(u,m,A,S)}updateInterpolation(u,m=!1){let A=this.gl.LINEAR;!m&&this.opts.isNearestInterpolation&&(A=this.gl.NEAREST),u===0?this.gl.activeTexture(TEXTURE0_BACK_VOL):this.gl.activeTexture(TEXTURE2_OVERLAY_VOL),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_MIN_FILTER,A),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_MAG_FILTER,A)}setAtlasOutline(u){this.opts.atlasOutline=u,this.updateGLVolume(),this.drawScene()}setInterpolation(u){this.opts.isNearestInterpolation=u;const m=this.volumes.length;if(!(m<1)){for(let A=0;A<m;A++)this.updateInterpolation(A);this.drawScene()}}calculateMvpMatrix2D(u,m,A,S=1/0,v=0,w=0,D=0,R){this.gl.viewport(u[0],this.gl.canvas.height-(u[1]+u[3]),u[2],u[3]);let L=m[0],d=A[0],E=[L,m[1]],e=[d-L,A[1]-m[1]];R&&(E=[A[0],m[1]],e=[m[0]-A[0],A[1]-m[1]],L=-A[0],d=-m[0]);const N=2*Math.max(Math.abs(m[2]),Math.abs(A[2])),q=create$3();let z=.01,e0=N*8;if(S!==1/0){let J=R;D===0&&(w===0||w===180)&&(J=!J);let i0=N*1.8-v;J||(i0=N*1.8+v),z=i0-S,e0=i0+S}ortho(q,L,d,m[1],A[1],z,e0);const H=create$3();H[0]=-1;const j=fromValues$2(0,0,-N*1.8);translate(H,H,j),rotateX(H,H,deg2rad(270-D)),rotateZ(H,H,deg2rad(w-180));const y=create$3();invert(y,H);const Y=create$3();transpose(Y,y);const G=create$3();return multiply(G,q,H),{modelViewProjectionMatrix:G,modelMatrix:H,normalMatrix:Y,leftTopMM:E,fovMM:e}}swizzleVec3MM(u,m){return m===1?u=swizzleVec3(u,[0,2,1]):m===2&&(u=swizzleVec3(u,[1,2,0])),u}screenFieldOfViewVox(u=0){const m=clone$1(this.volumeObject3D.fieldOfViewDeObliqueMM);return this.swizzleVec3MM(m,u)}screenFieldOfViewMM(u=0,m=!1){if(!m&&!this.opts.isSliceMM)return this.screenFieldOfViewVox(u);const A=this.volumeObject3D.extentsMin,S=this.volumeObject3D.extentsMax;let v=fromValues$2(A[0],A[1],A[2]),w=fromValues$2(S[0],S[1],S[2]);v=this.swizzleVec3MM(v,u),w=this.swizzleVec3MM(w,u);const D=create$2();return subtract$1(D,w,v),D}screenFieldOfViewExtendedVox(u=0){const m=this.volumes[0].extentsMinOrtho,A=this.volumes[0].extentsMaxOrtho;let S=fromValues$2(m[0],m[1],m[2]),v=fromValues$2(A[0],A[1],A[2]);const w=create$3();S=this.swizzleVec3MM(S,u),v=this.swizzleVec3MM(v,u);const D=create$2();return subtract$1(D,v,S),{mnMM:S,mxMM:v,rotation:w,fovMM:D}}screenFieldOfViewExtendedMM(u=0){if(!this.volumeObject3D)throw new Error("volumeObject3D undefined");const m=this.volumeObject3D.extentsMin,A=this.volumeObject3D.extentsMax;let S=fromValues$2(m[0],m[1],m[2]),v=fromValues$2(A[0],A[1],A[2]);const w=create$3();S=this.swizzleVec3MM(S,u),v=this.swizzleVec3MM(v,u);const D=create$2();return subtract$1(D,v,S),{mnMM:S,mxMM:v,rotation:w,fovMM:D}}drawSliceOrientationText(u,m){this.gl.viewport(0,0,this.gl.canvas.width,this.gl.canvas.height);let A="S";m===0&&(A="A");let S=this.opts.isRadiologicalConvention?"R":"L";if(m===2&&(S=this.opts.sagittalNoseLeft?"A":"P"),this.opts.isCornerOrientationText){this.drawTextRightBelow([u[0],u[1]],S+A);return}this.drawTextBelow([u[0]+u[2]*.5,u[1]],A),this.drawTextRight([u[0],u[1]+u[3]*.5],S)}xyMM2xyzMM(u,m){let A=2;u===1&&(A=1),u===2&&(A=0);let S=[0,0,0],v=[1,1,0],w=[1,0,1];S[A]=m,v[A]=m,w[A]=m,S=this.frac2mm(S),v=this.frac2mm(v),w=this.frac2mm(w),S=this.swizzleVec3MM(fromValues$2(S[0],S[1],S[2]),u),v=this.swizzleVec3MM(fromValues$2(v[0],v[1],v[2]),u),w=this.swizzleVec3MM(fromValues$2(w[0],w[1],w[2]),u);const D=(v[0]-S[0])*(w[1]-S[1])-(w[0]-S[0])*(v[1]-S[1]);let R=(v[0]-S[0])*(w[2]-S[2])-(w[0]-S[0])*(v[2]-S[2]);R/=D;let P=(v[1]-S[1])*(w[2]-S[2])-(w[1]-S[1])*(v[2]-S[2]);P/=D;const L=[0,0,0,0,0];return L[0]=S[0],L[1]=S[1],L[2]=S[2],L[3]=P,L[4]=R,L}draw2D(u,m,A=NaN){let S=this.volumes[0].frac2mm.slice(),v=this.screenFieldOfViewExtendedMM(m),w=create$3();this.opts.isSliceMM||(S=this.volumes[0].frac2mmOrtho.slice(),w=clone$2(this.volumes[0].mm2ortho),v=this.screenFieldOfViewExtendedVox(m));let D=this.opts.isRadiologicalConvention&&m<2;A===1/0||A===-1/0?(D=A!==1/0,m===1&&(D=!D)):this.opts.sagittalNoseLeft&&m===2&&(D=!D);let R=0,P=0;m===2?P=D?90:-90:m===1?P=D?180:0:(P=D?180:0,R=D?-90:90);const L=this.gl;if(u[2]===0||u[3]===0){const e0=L.canvas.width/v.fovMM[0],H=L.canvas.height/v.fovMM[1],j=Math.min(e0,H),y=e0/j,Y=H/j;v.fovMM[0]*=y,v.fovMM[1]*=Y;let G=(v.mnMM[0]+v.mxMM[0])*.5;v.mnMM[0]=G-v.fovMM[0]*.5,v.mxMM[0]=G+v.fovMM[0]*.5,G=(v.mnMM[1]+v.mxMM[1])*.5,v.mnMM[1]=G-v.fovMM[1]*.5,v.mxMM[1]=G+v.fovMM[1]*.5,u=[0,0,L.canvas.width,L.canvas.height]}if(isNaN(A)){const e0=this.scene.pan2Dxyzmm,H=this.swizzleVec3MM(fromValues$2(e0[0],e0[1],e0[2]),m),j=this.scene.pan2Dxyzmm[3];v.mnMM[0]-=H[0],v.mxMM[0]-=H[0],v.mnMM[1]-=H[1],v.mxMM[1]-=H[1],v.mnMM[0]/=j,v.mxMM[0]/=j,v.mnMM[1]/=j,v.mxMM[1]/=j}let d=2;m===1&&(d=1),m===2&&(d=0);let E=this.scene.crosshairPos[d],e=this.frac2mm(this.scene.crosshairPos);!isNaN(A)&&A!==1/0&&A!==-1/0&&(e=this.frac2mm([.5,.5,.5]),e[d]=A,E=this.mm2frac(e)[d]);const N=e[d];L.clear(L.DEPTH_BUFFER_BIT);let q=this.calculateMvpMatrix2D(u,v.mnMM,v.mxMM,1/0,0,P,R,D);if(A===1/0||A===-1/0){const e0=u.slice();this.draw3D(u,q.modelViewProjectionMatrix,q.modelMatrix,q.normalMatrix,P,R);const H=this.screenSlices[this.screenSlices.length-1];H.leftTopWidthHeight=e0,H.axCorSag=m,H.sliceFrac=1/0,H.AxyzMxy=this.xyMM2xyzMM(m,E),H.leftTopMM=q.leftTopMM,H.fovMM=q.fovMM;return}L.enable(L.DEPTH_TEST),L.blendFunc(L.SRC_ALPHA,L.ONE_MINUS_SRC_ALPHA),L.disable(L.BLEND),L.depthFunc(L.GREATER),L.disable(L.CULL_FACE);let z=this.sliceMMShader;if(this.opts.isV1SliceShader&&(z=this.sliceV1Shader),!z)throw new Error("slice Shader undefined");if(z.use(this.gl),L.uniform1f(z.uniforms.overlayOutlineWidth,this.overlayOutlineWidth),L.uniform1f(z.uniforms.overlayAlphaShader,this.overlayAlphaShader),L.uniform1i(z.uniforms.isAlphaClipDark,this.isAlphaClipDark?1:0),L.uniform1i(z.uniforms.backgroundMasksOverlays,this.backgroundMasksOverlays),L.uniform1f(z.uniforms.drawOpacity,this.drawOpacity),L.enable(L.BLEND),L.blendFunc(L.SRC_ALPHA,L.ONE_MINUS_SRC_ALPHA),L.uniform1f(z.uniforms.opacity,this.volumes[0].opacity),L.uniform1i(z.uniforms.axCorSag,m),L.uniform1f(z.uniforms.slice,E),L.uniformMatrix4fv(z.uniforms.frac2mm,!1,S),L.uniformMatrix4fv(z.uniforms.mvpMtx,!1,q.modelViewProjectionMatrix.slice()),L.bindVertexArray(this.genericVAO),L.drawArrays(L.TRIANGLE_STRIP,0,4),L.bindVertexArray(this.unusedVAO),this.screenSlices.push({leftTopWidthHeight:u,axCorSag:m,sliceFrac:E,AxyzMxy:this.xyMM2xyzMM(m,E),leftTopMM:q.leftTopMM,screen2frac:[],fovMM:q.fovMM}),isNaN(A)&&this.drawCrosshairs3D(!0,1,q.modelViewProjectionMatrix,!0,this.opts.isSliceMM),this.opts.meshThicknessOn2D>0){this.opts.meshThicknessOn2D!==1/0&&(q=this.calculateMvpMatrix2D(u,v.mnMM,v.mxMM,this.opts.meshThicknessOn2D,N,P,R,D));const e0=clone$2(q.modelViewProjectionMatrix);multiply(e0,e0,w),this.drawMesh3D(!0,1,e0,q.modelMatrix,q.normalMatrix)}isNaN(A)&&this.drawCrosshairs3D(!1,.15,q.modelViewProjectionMatrix,!0,this.opts.isSliceMM),this.drawSliceOrientationText(u,m),this.readyForSync=!0}calculateMvpMatrix(u,m=[0,0,0,0],A,S){(m[2]===0||m[3]===0)&&(m=[0,0,this.gl.canvas.width,this.gl.canvas.height]);const v=m[2]/m[3];let w=this.furthestFromPivot;const D=this.pivot3D,R=create$3();w=.8*w/this.scene.volScaleMultiplier,v<1?ortho(R,-w,w,-w/v,w/v,w*.01,w*8):ortho(R,-w*v,w*v,-w,w,w*.01,w*8);const P=create$3();P[0]=-1;const L=fromValues$2(0,0,-w*1.8);translate(P,P,L),this.position&&translate(P,P,this.position),rotateX(P,P,deg2rad(270-S)),rotateZ(P,P,deg2rad(A-180)),translate(P,P,[-D[0],-D[1],-D[2]]);const d=create$3();invert(d,P);const E=create$3();transpose(E,d);const e=create$3();return multiply(e,R,P),[e,P,E]}calculateModelMatrix(u,m){if(!this.back)throw new Error("back undefined");const A=create$3();if(A[0]=-1,rotateX(A,A,deg2rad(270-m)),rotateZ(A,A,deg2rad(u-180)),this.back.obliqueRAS){const S=clone$2(this.back.obliqueRAS);multiply(A,A,S)}return A}calculateRayDirection(u,m){const A=this.calculateModelMatrix(u,m),S=fromValues$3(1,0,0,0,0,-1,0,0,0,0,-1,0,0,0,0,1),v=create$3();multiply(v,S,A);const w=create$3();invert(w,v);const D=fromValues$1(0,0,-1,1);transformMat4(D,D,w);const R=fromValues$2(D[0],D[1],D[2]);normalize$1(R,R);const P=5e-5;return Math.abs(R[0])<P&&(R[0]=P),Math.abs(R[1])<P&&(R[1]=P),Math.abs(R[2])<P&&(R[2]=P),R}sceneExtentsMinMax(u=!0){let m=fromValues$2(0,0,0),A=fromValues$2(0,0,0);if(this.volumes.length>0){if(!this.volumeObject3D)throw new Error("volumeObject3D undefined");m=fromValues$2(this.volumeObject3D.extentsMin[0],this.volumeObject3D.extentsMin[1],this.volumeObject3D.extentsMin[2]),A=fromValues$2(this.volumeObject3D.extentsMax[0],this.volumeObject3D.extentsMax[1],this.volumeObject3D.extentsMax[2]),u||(m=fromValues$2(this.volumes[0].extentsMinOrtho[0],this.volumes[0].extentsMinOrtho[1],this.volumes[0].extentsMinOrtho[2]),A=fromValues$2(this.volumes[0].extentsMaxOrtho[0],this.volumes[0].extentsMaxOrtho[1],this.volumes[0].extentsMaxOrtho[2]))}if(this.meshes.length>0){if(this.volumes.length<1){const v=this.meshes[0].extentsMin,w=this.meshes[0].extentsMax;m=fromValues$2(v[0],v[1],v[2]),A=fromValues$2(w[0],w[1],w[2])}for(let v=0;v<this.meshes.length;v++){const w=this.meshes[v].extentsMin,D=this.meshes[v].extentsMax,R=fromValues$2(w[0],w[1],w[2]);min(m,m,R);const P=fromValues$2(D[0],D[1],D[2]);max$2(A,A,P)}}const S=create$2();return subtract$1(S,A,m),[m,A,S]}setPivot3D(){const[u,m]=this.sceneExtentsMinMax(),A=create$2();add$1(A,u,m),scale$2(A,A,.5),this.pivot3D=[A[0],A[1],A[2]],subtract$1(A,m,u),this.extentsMin=u,this.extentsMax=m,this.furthestFromPivot=length$1(A)*.5}getMaxVols(){if(this.volumes.length<1)return 0;let u=0;for(let m=0;m<this.volumes.length;m++)u=Math.max(u,this.volumes[m].nFrame4D);return u}detectPartialllyLoaded4D(){if(this.volumes.length<1)return!1;for(let u=0;u<this.volumes.length;u++)if(this.volumes[u].nFrame4D<this.volumes[u].hdr.dims[4])return!0;return!1}drawGraph(){if(this.getMaxVols()<2)return;const u=this.graph;let m=0;if(this.graph.autoSizeMultiplanar&&this.opts.sliceType===3)for(let l0=0;l0<this.screenSlices.length;l0++){const c0=this.screenSlices[l0].axCorSag;if(c0===0&&(m=this.screenSlices[l0].leftTopWidthHeight[1]),c0!==2)continue;const d0=this.screenSlices[l0].leftTopWidthHeight.slice();d0[1]===m?(u.LTWH[0]=d0[0]+d0[2],u.LTWH[1]=d0[1]):(u.LTWH[0]=d0[0],u.LTWH[1]=d0[1]+d0[3]),u.LTWH[2]=d0[2],u.LTWH[3]=d0[2]}if(u.opacity<=0||u.LTWH[2]<=5||u.LTWH[3]<=5||u.LTWH[0]+u.LTWH[2]>this.gl.canvas.width||u.LTWH[1]+u.LTWH[3]>this.gl.canvas.height)return;u.backColor=[.15,.15,.15,u.opacity],u.lineColor=[1,1,1,1],this.opts.backColor[0]+this.opts.backColor[1]+this.opts.backColor[2]>1.5&&(u.backColor=[.95,.95,.95,u.opacity],u.lineColor=[0,0,0,1]),u.textColor=u.lineColor.slice(),u.lineThickness=4,u.lineAlpha=1,u.lines=[];const A=[];if(u.vols.length<1)this.volumes[0]!=null&&A.push(0);else for(let l0=0;l0<u.vols.length;l0++){const c0=u.vols[l0];this.volumes[c0]==null||this.volumes[c0].nFrame4D<2||A.push(c0)}if(A.length<1)return;const S=this.volumes[A[0]].nFrame4D;if(this.graph.selectedColumn=this.volumes[A[0]].frame4D,S<2){log.debug("Unable to generate a graph: Selected volume is 3D not 4D");return}for(let l0=0;l0<A.length;l0++){u.lines[l0]=[];const c0=this.frac2vox(this.scene.crosshairPos),d0=this.volumes[A[l0]];let S0=d0.nFrame4D;S0=Math.min(S0,S);for(let p0=0;p0<S0;p0++){const T0=d0.getValue(c0[0],c0[1],c0[2],p0);u.lines[l0].push(T0)}}u.lineRGB=[[1,0,0],[0,.7,0],[0,0,1],[1,1,0],[1,0,1],[0,1,1],[1,1,1],[0,0,0]];let v=u.lines[0][0],w=u.lines[0][0];for(let l0=0;l0<u.lines.length;l0++)for(let c0=0;c0<u.lines[l0].length;c0++){const d0=u.lines[l0][c0];v=Math.min(d0,v),w=Math.max(d0,w)}if(u.normalizeValues&&w>v){const l0=w-v;for(let c0=0;c0<u.lines.length;c0++)for(let d0=0;d0<u.lines[c0].length;d0++){const S0=u.lines[c0][d0];u.lines[c0][d0]=(S0-v)/l0}v=0,w=1}v>=w&&(w=v+1),this.drawRect(u.LTWH,u.backColor);const[D,R,P]=tickSpacing(v,w),L=Math.max(0,-1*Math.floor(Math.log(D)/Math.log(10)));v=Math.min(R,v),w=Math.max(P,w);function d(l0){return l0.toFixed(6).replace(/\.?0*$/,"")}const e=.07*(Math.min(u.LTWH[2],u.LTWH[3])/(this.fontMets.size*this.uiData.dpr));let N=this.opts.textHeight*this.gl.canvas.height*e;N<16&&(N=0);let q=0,z=R;if(N>0)for(;z<=w;){const l0=z.toFixed(L),c0=this.textWidth(N,l0);q=Math.max(c0,q),z+=D}const e0=.05,H=Math.abs(u.LTWH[2]),j=Math.abs(u.LTWH[3]),y=[u.LTWH[0]+e0*H+q,u.LTWH[1]+e0*j,u.LTWH[2]-q-2*e0*H,u.LTWH[3]-N-2*e0*j];this.graph.LTWH=u.LTWH,this.graph.plotLTWH=y,this.drawRect(y,this.opts.backColor);const Y=w-v,G=y[3]/Y,J=y[2]/(u.lines[0].length-1),i0=y[1]+y[3];z=R+.5*D;const Q=u.lineColor.slice();for(Q[3]=.25*u.lineColor[3];z<=w;){const l0=i0-(z-v)*G;this.drawLine([y[0],l0,y[0]+y[2],l0],.5*u.lineThickness,Q),z+=D}z=R;const f=.5*u.lineThickness;for(;z<=w;){const l0=i0-(z-v)*G;this.drawLine([y[0]-f,l0,y[0]+y[2]+u.lineThickness,l0],u.lineThickness,u.lineColor);const c0=z.toFixed(L);N>0&&this.drawTextLeft([y[0]-6,l0],c0,e,u.textColor),z+=D}let o0=1;for(;u.lines[0].length/o0>20;)o0*=5;for(let l0=0;l0<u.lines[0].length;l0+=o0){const c0=l0*J+y[0];let d0=u.lineThickness;if(l0%2===1)d0*=.5,this.drawLine([c0,y[1],c0,y[1]+y[3]],d0,Q);else{const S0=d(l0);N>0&&this.drawTextBelow([c0,2+y[1]+y[3]],S0,e,u.textColor),this.drawLine([c0,y[1],c0,y[1]+y[3]],d0,u.lineColor)}}for(let l0=0;l0<u.lines.length;l0++){let c0=[1,0,0,u.lineAlpha];l0<u.lineRGB.length&&(c0=[u.lineRGB[l0][0],u.lineRGB[l0][1],u.lineRGB[l0][2],u.lineAlpha]);for(let d0=1;d0<u.lines[l0].length;d0++){const S0=(d0-1)*J,p0=d0*J,T0=(u.lines[l0][d0-1]-v)*G,x0=(u.lines[l0][d0]-v)*G,M0=[y[0]+S0,y[1]+y[3]-T0,y[0]+p0,y[1]+y[3]-x0];this.drawLine(M0,u.lineThickness,c0)}}if(u.selectedColumn>=0&&u.selectedColumn<u.lines[0].length){const l0=u.selectedColumn*J+y[0];this.drawLine([l0,y[1],l0,y[1]+y[3]],u.lineThickness,[u.lineRGB[3][0],u.lineRGB[3][1],u.lineRGB[3][2],1])}this.detectPartialllyLoaded4D()&&this.drawTextBelow([y[0]+y[2],y[1]+y[3]+N*.5],"...",e,u.textColor)}depthPicker(u,m){if(!this.uiData.mouseDepthPicker)return;this.uiData.mouseDepthPicker=!1;const A=this.gl,S=this.mousePos[0]*u[2]/u[2],v=A.canvas.height-this.mousePos[1]*u[3]/u[3]-1,w=new Uint8Array(4);if(A.readPixels(S,v,1,1,A.RGBA,A.UNSIGNED_BYTE,w),this.selectedObjectId=w[3],this.selectedObjectId===this.VOLUME_ID){this.scene.crosshairPos=new Float32Array(w.slice(0,3)).map(E=>E/255);return}const D=unpackFloatFromVec4i(w);if(D>1)return;const R=(this.mousePos[0]-u[0])/u[2],P=(A.canvas.height-this.mousePos[1]-u[1])/u[3],L=unProject(R,P,D,m),d=this.mm2frac(L,0,!0);d[0]<0||d[0]>1||d[1]<0||d[1]>1||d[2]<0||d[2]>1||(this.scene.crosshairPos=this.mm2frac(L,0,!0))}drawImage3D(u,m,A){if(this.volumes.length===0)return;const S=this.gl,v=this.calculateRayDirection(m,A),w=this.volumeObject3D;if(w){S.enable(S.BLEND),S.blendFunc(S.SRC_ALPHA,S.ONE_MINUS_SRC_ALPHA),S.enable(S.CULL_FACE),S.cullFace(S.FRONT);let D=this.renderShader;if(this.uiData.mouseDepthPicker&&(D=this.pickingImageShader),D.use(this.gl),S.uniform1i(D.uniforms.backgroundMasksOverlays,this.backgroundMasksOverlays),this.gradientTextureAmount>0){S.activeTexture(TEXTURE6_GRADIENT),S.bindTexture(S.TEXTURE_3D,this.gradientTexture);const R=this.calculateModelMatrix(m,A),P=create$3();invert(P,R);const L=create$3();transpose(L,P),S.uniformMatrix4fv(D.uniforms.normMtx,!1,L)}this.drawBitmap&&this.drawBitmap.length>8?S.uniform2f(D.uniforms.renderDrawAmbientOcclusionXY,this.renderDrawAmbientOcclusion,this.drawOpacity):S.uniform2f(D.uniforms.renderDrawAmbientOcclusionXY,this.renderDrawAmbientOcclusion,0),S.uniformMatrix4fv(D.uniforms.mvpMtx,!1,u),S.uniformMatrix4fv(D.uniforms.matRAS,!1,this.back.matRAS),S.uniform3fv(D.uniforms.rayDir,v),this.gradientTextureAmount<0?S.uniform4fv(D.uniforms.clipPlane,[this.scene.crosshairPos[0],this.scene.crosshairPos[1],this.scene.crosshairPos[2],30]):S.uniform4fv(D.uniforms.clipPlane,this.scene.clipPlane),S.uniform1f(D.uniforms.drawOpacity,1),S.bindVertexArray(w.vao),S.drawElements(w.mode,w.indexCount,S.UNSIGNED_SHORT,0),S.bindVertexArray(this.unusedVAO)}}drawOrientationCube(u,m=0,A=0){if(!this.opts.isOrientCube)return;const S=.05*Math.min(u[2],u[3]);if(S<5)return;const v=this.gl;v.enable(v.CULL_FACE),v.cullFace(v.BACK),this.orientCubeShader.use(this.gl),v.bindVertexArray(this.orientCubeShaderVAO);const w=create$3(),D=create$3();ortho(D,0,v.canvas.width,0,v.canvas.height,-10*S,10*S);let R=0;u[1]===0&&(R=v.canvas.height-this.effectiveCanvasHeight()),translate(w,w,[1.8*S+u[0],R+1.8*S+u[1],0]),scale$3(w,w,[S,S,S]),rotateX(w,w,deg2rad(270-A)),rotateZ(w,w,deg2rad(-m));const P=create$3();multiply(P,D,w),v.uniformMatrix4fv(this.orientCubeShader.uniforms.u_matrix,!1,P),v.drawArrays(v.TRIANGLE_STRIP,0,168),v.bindVertexArray(this.unusedVAO),this.gl.disable(this.gl.CULL_FACE)}createOnLocationChange(u=NaN){const[m,A,S]=this.sceneExtentsMinMax(!0),v=Math.max(Math.max(S[0],S[1]),S[2]);function w(E){return Math.max(0,-Math.ceil(Math.log10(Math.abs(E))))}let D=w(v*.001);const R=this.frac2mm(this.scene.crosshairPos,0,!0);function P(E,e=0){return parseFloat(E.toFixed(e))}let L=P(R[0],D)+"×"+P(R[1],D)+"×"+P(R[2],D);if(this.volumes.length>0&&this.volumes[0].nFrame4D>0&&(L+="×"+P(this.volumes[0].frame4D)),this.volumes.length>0){let E=" = ";for(let q=0;q<this.volumes.length;q++){const z=this.volumes[q].mm2vox(R);let e0=this.volumes[q].getValue(z[0],z[1],z[2],this.volumes[q].frame4D);if(D=3,this.volumes[q].colormapLabel!==null){const H=Math.round(e0);H>=0&&H<this.volumes[q].colormapLabel.labels.length?E+=this.volumes[q].colormapLabel.labels[H]:E+="undefined("+P(e0,D)+")"}else E+=P(e0,D);this.volumes[q].imaginary&&(e0=this.volumes[q].getValue(z[0],z[1],z[2],this.volumes[q].frame4D,!0),e0>=0&&(E+="+"),E+=P(e0,D)),E+="   "}L+=E;const e=this.back.dimsRAS,N=e[1]*e[2]*e[3];if(this.drawBitmap&&this.drawBitmap.length===N){const q=this.frac2vox(this.scene.crosshairPos),z=q[0]+q[1]*e[1]+q[2]*e[1]*e[2];L+=" "+this.drawLut.labels[this.drawBitmap[z]]}}const d={mm:this.frac2mm(this.scene.crosshairPos,0,!0),axCorSag:u,vox:this.frac2vox(this.scene.crosshairPos),frac:this.scene.crosshairPos,xy:[this.mousePos[0],this.mousePos[1]],values:this.volumes.map(E=>{const e=this.frac2mm(this.scene.crosshairPos,0,!0),N=E.mm2vox(e),q=E.getValue(N[0],N[1],N[2],E.frame4D);return{name:E.name,value:q,id:E.id,mm:e,vox:N}}),string:L};this.onLocationChange(d)}addLabel(u,m,A){const S={textColor:this.opts.legendTextColor,textScale:1,textAlignment:"left",lineWidth:0,lineColor:this.opts.legendTextColor,lineTerminator:"none",bulletScale:0,bulletColor:this.opts.legendTextColor},v=m?{...S,...m}:{...S},w=new NVLabel3D(u,v,A);return this.document.labels.push(w),w}calculateScreenPoint(u,m,A){const S=create$1();return transformMat4(S,[...u,1],m),S[3]!==0&&(S[0]=(S[0]/S[3]+1)*.5*A[2],S[1]=(1-S[1]/S[3])*.5*A[3],S[2]/=S[3],S[0]+=A[0],S[1]+=A[1]),S}getLabelAtPoint(u){log.debug("screenPoint",u);const m=this.getLegendPanelHeight(),A=this.getLegendPanelWidth(),S=this.gl.canvas.width-A;let v=(this.canvas.height-m)/2;if(log.debug("panelrect",S,v,S+A,v+m),u[0]<S||u[1]<v||u[0]>S+A||u[1]>v+m)return null;const D=this.opts.textHeight*Math.min(this.gl.canvas.height,this.gl.canvas.width)*1,R=this.getAllLabels();for(const P of R){const L=this.opts.textHeight*this.gl.canvas.height*P.style.textScale,d=this.textHeight(L,P.text);if(u[1]>=v&&u[1]<=v+d+D/2)return P;v+=d,v+=D/2}return null}drawLabelLine(u,m,A,S,v=!1){const w=Array.isArray(u.points)&&Array.isArray(u.points[0])?u.points:[u.points];for(const D of w){const R=this.calculateScreenPoint(D,A,S);v?this.drawDottedLine([...m,R[0],R[1]],u.style.lineWidth,u.style.lineColor):this.draw3DLine(m,[R[0],R[1],R[2]],u.style.lineWidth,u.style.lineColor)}}draw3DLabel(u,m,A,S,v=0,w,D=!1){const R=u.text,P=m[0],L=m[1],d=this.opts.textHeight*Math.min(this.gl.canvas.height,this.gl.canvas.width)*1,E=this.textHeight(u.style.textScale,R)*d;if(u.style.lineWidth>0&&Array.isArray(u.points)&&this.drawLabelLine(u,[P,L+E],A,S,D),u.style.bulletScale){const N=u.style.bulletScale*E,q=E-N,z=L+q/2+N/2,e0=P+(v-N)/2;this.drawCircle([e0,z,N,N],u.style.bulletColor)}let e=P;if(u.style.textAlignment!=="left"){const N=this.textWidth(u.style.textScale,u.text)*d;if(u.style.textAlignment==="right")e=P+w-d*1.5-N;else{const q=w-(v||d);e+=(q-N)/2}}else e+=v;this.drawText([e,L],R,u.style.textScale,u.style.textColor)}draw3DLabels(u,m,A=!1){const S=this.getAllLabels();if(!this.opts.showLegend||S.length===0)return;if(!this.canvas)throw new Error("canvas undefined");const v=this.gl;v.disable(v.CULL_FACE),v.viewport(0,0,this.canvas.width,this.canvas.height);const D=this.opts.textHeight*Math.min(this.gl.canvas.height,this.gl.canvas.width)*1,R=this.getBulletMarginWidth(),P=this.getLegendPanelHeight(),L=this.getLegendPanelWidth(),d=v.canvas.width-L;let E=(this.canvas.height-P)/2;this.drawRect([v.canvas.width-L,E,L-D,P],this.opts.legendBackgroundColor);const e=v.getParameter(v.BLEND),N=v.getParameter(v.DEPTH_FUNC);A||(v.disable(v.BLEND),v.depthFunc(v.GREATER));for(const q of S){this.draw3DLabel(q,[d,E],u,m,R,L,A);const z=this.opts.textHeight*this.gl.canvas.height*q.style.textScale,e0=this.textHeight(z,q.text);E+=e0,E+=D/2}A||(v.depthFunc(N),e&&v.enable(v.BLEND))}draw3D(u=[0,0,0,0],m=null,A=null,S=null,v=null,w=0){const D=v!==null;this.setPivot3D(),D||(v=this.scene.renderAzimuth,w=this.scene.renderElevation);const R=this.gl;m===null&&([m,A,S]=this.calculateMvpMatrix(null,u,v,w));let P=[...u];if(u[2]===0||u[3]===0?(u=[0,0,R.canvas.width,R.canvas.height],P=[...u],this.screenSlices.push({leftTopWidthHeight:u,axCorSag:4,sliceFrac:0,AxyzMxy:[],leftTopMM:[],fovMM:[isRadiological(A),0]})):(this.screenSlices.push({leftTopWidthHeight:u.slice(),axCorSag:4,sliceFrac:0,AxyzMxy:[],leftTopMM:[],fovMM:[isRadiological(A),0]}),u[1]=R.canvas.height-u[3]-u[1]),R.enable(R.DEPTH_TEST),R.depthFunc(R.ALWAYS),R.depthMask(!0),R.clearDepth(0),this.draw3DLabels(m,P,!1),R.viewport(u[0],u[1],u[2],u[3]),this.volumes.length>0&&(this.updateInterpolation(0,!0),this.updateInterpolation(1,!0),this.drawImage3D(m,v,w)),this.updateInterpolation(0),this.updateInterpolation(1),D||this.drawCrosshairs3D(!0,1,m),this.drawMesh3D(!0,1,m,A,S),this.uiData.mouseDepthPicker){this.depthPicker(u,m),this.createOnLocationChange(),this.draw3D(u,m,A,S,v,w);return}this.opts.meshXRay>0&&this.drawMesh3D(!1,this.opts.meshXRay,m,A,S),this.draw3DLabels(m,P,!1),R.viewport(u[0],u[1],u[2],u[3]),D||this.drawCrosshairs3D(!1,.15,m),R.viewport(0,0,R.canvas.width,R.canvas.height),this.drawOrientationCube(u,v,w);const L="azimuth: "+this.scene.renderAzimuth.toFixed(0)+" elevation: "+this.scene.renderElevation.toFixed(0);return this.readyForSync=!0,this.sync(),this.draw3DLabels(m,P,!0),L}drawMesh3D(u=!0,m=1,A,S,v){if(this.meshes.length<1)return;const w=this.gl;A||([A,S,v]=this.calculateMvpMatrix(this.volumeObject3D,void 0,this.scene.renderAzimuth,this.scene.renderElevation)),w.enable(w.DEPTH_TEST),w.blendFunc(w.SRC_ALPHA,w.ONE_MINUS_SRC_ALPHA),w.disable(w.BLEND),w.depthFunc(w.GREATER),w.disable(w.CULL_FACE),u?(w.disable(w.BLEND),w.depthFunc(w.GREATER)):(w.enable(w.BLEND),w.depthFunc(w.ALWAYS),w.enable(w.CULL_FACE)),w.cullFace(w.BACK);let D=this.meshShaders[0].shader,R=!1;for(let P=0;P<this.meshes.length;P++)if(this.meshes[P].visible!==!1&&(D=this.meshShaders[this.meshes[P].meshShaderIndex].shader,this.uiData.mouseDepthPicker&&(D=this.pickingMeshShader),D.use(this.gl),w.uniformMatrix4fv(D.uniforms.mvpMtx,!1,A),w.uniformMatrix4fv(D.uniforms.normMtx,!1,v),w.uniform1f(D.uniforms.opacity,m),!(this.meshes[P].indexCount<3))){if(this.meshes[P].offsetPt0&&(this.meshes[P].fiberSides<3||this.meshes[P].fiberRadius<=0)){R=!0;continue}D.isMatcap&&(w.activeTexture(TEXTURE5_MATCAP),w.bindTexture(w.TEXTURE_2D,this.matCapTexture)),w.bindVertexArray(this.meshes[P].vao),w.drawElements(w.TRIANGLES,this.meshes[P].indexCount,w.UNSIGNED_INT,0),w.bindVertexArray(this.unusedVAO)}if(!R){w.enable(w.BLEND),w.depthFunc(w.ALWAYS);return}D=this.fiberShader,D.use(this.gl),w.uniformMatrix4fv(D.uniforms.mvpMtx,!1,A),w.uniform1f(D.uniforms.opacity,m);for(let P=0;P<this.meshes.length;P++)this.meshes[P].indexCount<3||this.meshes[P].offsetPt0&&(this.meshes[P].fiberSides>=3&&this.meshes[P].fiberRadius>0||(w.bindVertexArray(this.meshes[P].vaoFiber),w.drawElements(w.LINE_STRIP,this.meshes[P].indexCount,w.UNSIGNED_INT,0),w.bindVertexArray(this.unusedVAO)));w.enable(w.BLEND),w.depthFunc(w.ALWAYS),this.readyForSync=!0}drawCrosshairs3D(u=!0,m=1,A=null,S=!1,v=!0){if(!this.opts.show3Dcrosshair&&!S||this.opts.crosshairWidth<=0&&S)return;const w=this.gl,D=this.frac2mm(this.scene.crosshairPos,0,v);if(this.crosshairs3D===null||this.crosshairs3D.mm[0]!==D[0]||this.crosshairs3D.mm[1]!==D[1]||this.crosshairs3D.mm[2]!==D[2]){this.crosshairs3D!==null&&(w.deleteBuffer(this.crosshairs3D.indexBuffer),w.deleteBuffer(this.crosshairs3D.vertexBuffer));const[L,d,E]=this.sceneExtentsMinMax(v);let e=1;if(this.volumes.length>0){if(!this.back)throw new Error("back undefined");e=.5*Math.min(Math.min(this.back.pixDims[1],this.back.pixDims[2]),this.back.pixDims[3])}else(E[0]<50||E[0]>1e3)&&(e=E[0]*.02);e*=this.opts.crosshairWidth,this.crosshairs3D=NiivueObject3D.generateCrosshairs(this.gl,1,D,L,d,e,20,this.opts.crosshairGap),this.crosshairs3D.mm=D}if(!this.surfaceShader)throw new Error("surfaceShader undefined");const R=this.surfaceShader;R.use(this.gl),A==null&&([A]=this.calculateMvpMatrix(this.crosshairs3D,void 0,this.scene.renderAzimuth,this.scene.renderElevation)),w.uniformMatrix4fv(R.uniforms.mvpMtx,!1,A),w.bindBuffer(w.ELEMENT_ARRAY_BUFFER,this.crosshairs3D.indexBuffer),w.enable(w.DEPTH_TEST);const P=[...this.opts.crosshairColor];u?(w.disable(w.BLEND),w.depthFunc(w.GREATER)):(w.enable(w.BLEND),w.blendFunc(w.SRC_ALPHA,w.ONE_MINUS_SRC_ALPHA),w.depthFunc(w.ALWAYS)),P[3]=m,w.uniform4fv(R.uniforms.surfaceColor,P),w.bindVertexArray(this.crosshairs3D.vao),w.drawElements(w.TRIANGLES,this.crosshairs3D.indexCount,w.UNSIGNED_INT,0),w.bindVertexArray(this.unusedVAO)}mm2frac(u,m=0,A=!1){if(this.volumes.length<1){const S=fromValues$2(.1,.5,.5),[v,w,D]=this.sceneExtentsMinMax();return S[0]=(u[0]-v[0])/D[0],S[1]=(u[1]-v[1])/D[1],S[2]=(u[2]-v[2])/D[2],isFinite(S)||(isFinite(S[0])||(S[0]=.5),isFinite(S[1])||(S[1]=.5),isFinite(S[2])||(S[2]=.5),this.meshes.length<1&&log.error("mm2frac() not finite: objects not (yet) loaded.")),S}return this.volumes[m].convertMM2Frac(u,A||this.opts.isSliceMM)}vox2frac(u,m=0){return this.volumes[m].convertVox2Frac(u)}frac2vox(u,m=0){return this.volumes.length<=m?[0,0,0]:this.volumes[m].convertFrac2Vox(u)}moveCrosshairInVox(u,m,A){const S=this.frac2vox(this.scene.crosshairPos);S[0]+=u,S[1]+=m,S[2]+=A,S[0]=clamp(S[0],0,this.volumes[0].dimsRAS[1]-1),S[1]=clamp(S[1],0,this.volumes[0].dimsRAS[2]-1),S[2]=clamp(S[2],0,this.volumes[0].dimsRAS[3]-1),this.scene.crosshairPos=this.vox2frac(S),this.createOnLocationChange(),this.drawScene()}frac2mm(u,m=0,A=!1){const S=fromValues$1(u[0],u[1],u[2],1);if(this.volumes.length>0)return this.volumes[m].convertFrac2MM(u,A||this.opts.isSliceMM);{const[v,w]=this.sceneExtentsMinMax(),D=(R,P,L)=>R*(1-L)+P*L;S[0]=D(v[0],w[0],u[0]),S[1]=D(v[1],w[1],u[1]),S[2]=D(v[2],w[2],u[2])}return S}screenXY2TextureFrac(u,m,A,S=!0){const v=fromValues$2(-1,-1,-1),w=this.screenSlices[A].axCorSag;if(w>2)return v;const D=this.screenSlices[A].leftTopWidthHeight.slice();let R=!1;D[2]<0&&(R=!0,D[0]+=D[2],D[2]=-D[2]);let P=(u-D[0])/D[2];R&&(P=1-P);const L=1-(m-D[1])/D[3];if(P<0||P>1||L<0||L>1||this.screenSlices[A].AxyzMxy.length<4)return v;let d=fromValues$2(0,0,0);d[0]=this.screenSlices[A].leftTopMM[0]+P*this.screenSlices[A].fovMM[0],d[1]=this.screenSlices[A].leftTopMM[1]+L*this.screenSlices[A].fovMM[1];const E=this.screenSlices[A].AxyzMxy;d[2]=E[2]+E[4]*(d[1]-E[1])-E[3]*(d[0]-E[0]),w===1&&(d=swizzleVec3(d,[0,2,1])),w===2&&(d=swizzleVec3(d,[2,0,1]));const e=this.mm2frac(d);return S&&(e[0]<0||e[0]>1||e[1]<0||e[1]>1||e[2]<0||e[2]>1)?v:e}canvasPos2frac(u){for(let m=0;m<this.screenSlices.length;m++){const A=this.screenXY2TextureFrac(u[0],u[1],m);if(A[0]>=0)return A}return[-1,-1,-1]}scaleSlice(u,m,A=0,S=0){const v=this.effectiveCanvasWidth()-A,w=this.effectiveCanvasHeight()-S;let D=v/u;m*D>w&&(D=w/m);const R=u*D,P=m*D;return[(v-R)*.5,(w-P)*.5,R,P,D]}drawThumbnail(){if(!this.bmpShader)throw new Error("bmpShader undefined");this.bmpShader.use(this.gl),this.gl.uniform2f(this.bmpShader.uniforms.canvasWidthHeight,this.gl.canvas.width,this.gl.canvas.height);let u=this.gl.canvas.height,m=this.gl.canvas.height*this.bmpTextureWH;m>this.gl.canvas.width&&(u=this.gl.canvas.width/this.bmpTextureWH,m=this.gl.canvas.width),this.gl.uniform4f(this.bmpShader.uniforms.leftTopWidthHeight,0,0,m,u),this.gl.bindVertexArray(this.genericVAO),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4),this.gl.bindVertexArray(this.unusedVAO)}drawLine(u,m=1,A=[1,0,0,-1]){if(this.gl.bindVertexArray(this.genericVAO),!this.lineShader)throw new Error("lineShader undefined");this.lineShader.use(this.gl),A[3]<0&&(A=this.opts.crosshairColor),this.gl.uniform4fv(this.lineShader.uniforms.lineColor,A),this.gl.uniform2fv(this.lineShader.uniforms.canvasWidthHeight,[this.gl.canvas.width,this.gl.canvas.height]),this.gl.uniform1f(this.lineShader.uniforms.thickness,m),this.gl.uniform4fv(this.lineShader.uniforms.startXYendXY,u),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4),this.gl.bindVertexArray(this.unusedVAO)}draw3DLine(u,m,A=1,S=[1,0,0,-1]){if(this.gl.bindVertexArray(this.genericVAO),!this.line3DShader)throw new Error("line3DShader undefined");this.line3DShader.use(this.gl),S[3]<0&&(S=this.opts.crosshairColor),this.gl.uniform4fv(this.line3DShader.uniforms.lineColor,S),this.gl.uniform2fv(this.line3DShader.uniforms.canvasWidthHeight,[this.gl.canvas.width,this.gl.canvas.height]),this.gl.uniform1f(this.line3DShader.uniforms.thickness,A),this.gl.uniform2fv(this.line3DShader.uniforms.startXY,u),this.gl.uniform3fv(this.line3DShader.uniforms.endXYZ,m),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4),this.gl.bindVertexArray(this.unusedVAO)}drawDottedLine(u,m=1,A=[1,0,0,-1]){if(this.gl.bindVertexArray(this.genericVAO),!this.lineShader)throw new Error("lineShader undefined");this.lineShader.use(this.gl);const S=A[3]<0?[...this.opts.crosshairColor]:[...A];S[3]=.3;const v=fromValues(u[2]-u[0],u[3]-u[1]),w=length(v);normalize(v,v);const R=this.opts.textHeight*Math.min(this.gl.canvas.height,this.gl.canvas.width)*1;scale(v,v,R/2);const P=length(v);let L=Math.floor(w/P);w%P&&L++;const d=[u[0],u[1]];this.gl.uniform4fv(this.lineShader.uniforms.lineColor,S),this.gl.uniform2fv(this.lineShader.uniforms.canvasWidthHeight,[this.gl.canvas.width,this.gl.canvas.height]),this.gl.uniform1f(this.lineShader.uniforms.thickness,m);for(let E=0;E<L-1;E++){if(E%2){d[0]+=v[0],d[1]+=v[1];continue}const e=[d[0],d[1],d[0]+v[0],d[1]+v[1]];this.gl.uniform4fv(this.lineShader.uniforms.startXYendXY,e),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4),d[0]+=v[0],d[1]+=v[1]}this.gl.bindVertexArray(this.unusedVAO)}drawGraphLine(u,m=[1,0,0,.5],A=2){this.drawLine(u,A,m)}drawCrossLinesMM(u,m,A,S,v){if(u<0||this.screenSlices.length<=u)return;const w=this.screenSlices[u];let D=w.sliceFrac;D===1/0&&log.warn("Rendering approximate cross lines in world view mode"),D===1/0&&(D=.5);let P=S.slice(),L=v.slice();const d=Math.max(1,this.opts.crosshairWidth);m===1&&(P=A.slice()),m===2&&(P=A.slice(),L=S.slice());function E(e){const N=fromValues(0,0);return N[0]=w.leftTopWidthHeight[0]+(e[0]-w.leftTopMM[0])/w.fovMM[0]*w.leftTopWidthHeight[2],N[1]=w.leftTopWidthHeight[1]+w.leftTopWidthHeight[3]-(e[1]-w.leftTopMM[1])/w.fovMM[1]*w.leftTopWidthHeight[3],N}if(P.length>0&&m===0){const e=D,N=1;for(let q=0;q<P.length;q++){const z=this.frac2mm([.5,.5,.5]);z[N]=P[q];let e0=this.mm2frac(z);e0=e0[N];let H=this.frac2mm([0,e0,e]);H=swizzleVec3(H,[0,1,2]);let j=this.frac2mm([1,e0,e]);j=swizzleVec3(j,[0,1,2]),H=E(H),j=E(j),this.drawLine([H[0],H[1],j[0],j[1]],d)}}if(P.length>0&&m===1){const e=D,N=2;for(let q=0;q<P.length;q++){const z=this.frac2mm([.5,.5,.5]);z[N]=P[q];let e0=this.mm2frac(z);e0=e0[N];let H=this.frac2mm([0,e,e0]);H=swizzleVec3(H,[0,2,1]);let j=this.frac2mm([1,e,e0]);j=swizzleVec3(j,[0,2,1]),H=E(H),j=E(j),this.drawLine([H[0],H[1],j[0],j[1]],d)}}if(P.length>0&&m===2){const e=D,N=2;for(let q=0;q<P.length;q++){const z=this.frac2mm([.5,.5,.5]);z[N]=P[q];let e0=this.mm2frac(z);e0=e0[N];let H=this.frac2mm([e,0,e0]);H=swizzleVec3(H,[1,2,0]);let j=this.frac2mm([e,1,e0]);j=swizzleVec3(j,[1,2,0]),H=E(H),j=E(j),this.drawLine([H[0],H[1],j[0],j[1]],d)}}if(L.length>0&&m===0){const e=D,N=0;for(let q=0;q<L.length;q++){const z=this.frac2mm([.5,.5,.5]);z[N]=L[q];let e0=this.mm2frac(z);e0=e0[N];let H=this.frac2mm([e0,0,e]);H=swizzleVec3(H,[0,1,2]);let j=this.frac2mm([e0,1,e]);j=swizzleVec3(j,[0,1,2]),H=E(H),j=E(j),this.drawLine([H[0],H[1],j[0],j[1]],d)}}if(L.length>0&&m===1){const e=D,N=0;for(let q=0;q<L.length;q++){const z=this.frac2mm([.5,.5,.5]);z[N]=L[q];let e0=this.mm2frac(z);e0=e0[N];let H=this.frac2mm([e0,e,0]);H=swizzleVec3(H,[0,2,1]);let j=this.frac2mm([e0,e,1]);j=swizzleVec3(j,[0,2,1]),H=E(H),j=E(j),this.drawLine([H[0],H[1],j[0],j[1]],d)}}if(L.length>0&&m===2){const e=D,N=1;for(let q=0;q<L.length;q++){const z=this.frac2mm([.5,.5,.5]);z[N]=L[q];let e0=this.mm2frac(z);e0=e0[N];let H=this.frac2mm([e,e0,0]);H=swizzleVec3(H,[1,2,0]);let j=this.frac2mm([e,e0,1]);j=swizzleVec3(j,[1,2,0]),H=E(H),j=E(j),this.drawLine([H[0],H[1],j[0],j[1]],d)}}}drawCrossLines(u,m,A,S,v){if(u<0||this.screenSlices.length<=u)return;if(this.opts.isSliceMM)return this.drawCrossLinesMM(u,m,A,S,v);if(this.screenSlices[u].sliceFrac===1/0)return this.drawCrossLinesMM(u,m,A,S,v);const w=this.screenSlices[u];let D=S.slice(),R=v.slice();if(m===1&&(D=A.slice()),m===2&&(D=A.slice(),R=S.slice()),D.length>0){const P=w.leftTopWidthHeight.slice();let L=2;m===0&&(L=1);const d=this.frac2mm([.5,.5,.5]);for(let E=0;E<D.length;E++){d[L]=D[E];const e=this.mm2frac(d);this.drawRect([P[0],P[1]+P[3]-e[L]*P[3],P[2],1])}}if(R.length>0){const P=w.leftTopWidthHeight.slice(),L=w.fovMM[0]<0;let d=0;m===2&&(d=1);const E=this.frac2mm([.5,.5,.5]);for(let e=0;e<R.length;e++){E[d]=R[e];const N=this.mm2frac(E);L?this.drawRect([P[0]+(P[2]-N[d]*P[2]),P[1],1,P[3]]):this.drawRect([P[0]+N[d]*P[2],P[1],1,P[3]])}}}drawMosaic(u){if(this.volumes.length===0){log.debug("Unable to draw mosaic until voxel-based image is loaded");return}this.screenSlices=[];const m=this.screenFieldOfViewMM(0,!0),A=this.screenFieldOfViewMM(0);u=u.replaceAll(";"," ;").trim();const S=[],v=[],w=[],D=u.split(/\s+/);let R=1;const P=this.opts.textHeight;for(let L=0;L<2;L++){let d=!1,E=!1;d=!1;let e=0,N=0,q=0,z=0,e0=!1,H=0;for(let Y=0;Y<D.length;Y++){const G=D[Y];if(G.includes("X")){E=!0;continue}if(G.includes("L")){e0=!G.includes("-");continue}if(G.includes("V")||G.includes("H")){Y++;continue}G.includes("A")&&(H=0),G.includes("C")&&(H=1),G.includes("S")&&(H=2),G.includes("R")&&(d=!0),G.includes(";")&&(q+=e,z=Math.max(z,N),e=0,N=0);const J=parseFloat(G);if(isNaN(J))continue;let i0=0,Q=0,f=A;if(d&&(f=m),H===2?i0=f[1]:i0=f[0],H===0?Q=f[1]:Q=f[2],L===0)d||(H===0&&S.push(J),H===1&&v.push(J),H===2&&w.push(J));else{const o0=[R*N,R*q,R*i0,R*Q];if(this.opts.textHeight=e0?P:0,d){let l0=J<0?-1/0:1/0;Object.is(J,-0)&&(l0=-1/0),this.draw2D(o0,H,l0)}else this.draw2D(o0,H,J);E&&this.drawCrossLines(this.screenSlices.length-1,H,S,v,w),d=!1,E=!1}N+=i0,e=Math.max(e,Q)}if(q+=e,z=Math.max(z,N),z<=0||q<=0)break;const j=this.gl.canvas.width/z,y=this.effectiveCanvasHeight()/q;R=Math.min(j,y)}this.opts.textHeight=P}drawSceneCore(){if(!this.initialized)return;if(this.colorbarHeight=0,this.gl.clearColor(this.opts.backColor[0],this.opts.backColor[1],this.opts.backColor[2],this.opts.backColor[3]),this.gl.clear(this.gl.COLOR_BUFFER_BIT|this.gl.DEPTH_BUFFER_BIT),this.bmpTexture&&this.thumbnailVisible){this.drawThumbnail();return}let u="";if(this.volumes.length===0||typeof this.volumes[0].dims>"u"){if(this.meshes.length>0){this.screenSlices=[],this.opts.sliceType=4,this.draw3D(),this.opts.isColorbar&&this.drawColorbar();return}this.drawLoadingText(this.loadingText);return}if(this.back===null)return;if(this.uiData.isDragging&&this.scene.clipPlaneDepthAziElev[0]<1.8&&this.inRenderTile(this.uiData.dragStart[0],this.uiData.dragStart[1])>=0){const v=this.uiData.dragStart[0]-this.uiData.dragEnd[0],w=this.uiData.dragStart[1]-this.uiData.dragEnd[1],D=this.uiData.dragClipPlaneStartDepthAziElev.slice();if(D[1]-=v,D[1]=D[1]%360,D[2]+=w,D[1]!==this.scene.clipPlaneDepthAziElev[1]||D[2]!==this.scene.clipPlaneDepthAziElev[2])return this.scene.clipPlaneDepthAziElev=D,this.setClipPlane(this.scene.clipPlaneDepthAziElev)}if(this.sliceMosaicString.length<1&&this.opts.sliceType===4){this.opts.isColorbar&&this.reserveColorbarPanel(),this.screenSlices=[],this.draw3D(),this.opts.isColorbar&&this.drawColorbar();return}this.opts.isColorbar&&this.reserveColorbarPanel();const m=this.getMaxVols(),A=this.opts.sliceType===3&&m>1&&this.graph.autoSizeMultiplanar&&this.graph.opacity>0;if(this.sliceMosaicString.length>0)this.drawMosaic(this.sliceMosaicString);else if(this.gl.viewport(0,0,this.gl.canvas.width,this.gl.canvas.height),this.screenSlices=[],this.opts.sliceType===0)this.draw2D([0,0,0,0],0);else if(this.opts.sliceType===1)this.draw2D([0,0,0,0],1);else if(this.opts.sliceType===2)this.draw2D([0,0,0,0],2);else{const v=isFinite(this.drawPenLocation[0])&&this.opts.drawingEnabled,{volScale:w}=this.sliceScale();typeof this.opts.multiplanarPadPixels!="number"&&log.debug("multiplanarPadPixels must be numeric");const D=parseFloat(`${this.opts.multiplanarPadPixels}`),R=this.scaleSlice(w[0]+w[1],w[1]+w[2],D*1,D*1),P=Math.max(Math.max(w[1],w[2]),w[0]),L=this.scaleSlice(w[0]+w[0]+w[1],Math.max(w[1],w[2]),D*2),d=this.scaleSlice(w[0]+w[0]+w[1]+P,Math.max(w[1],w[2]),D*3),E=this.scaleSlice(P,w[1]+w[2]+w[2],0,D*2),e=this.scaleSlice(P,w[1]+w[2]+w[2]+P,0,D*3);let N=!v&&(m<2||!A),q=!1,z=!1,e0=!1;if(this.opts.multiplanarLayout===1?q=!0:this.opts.multiplanarLayout===2?z=!0:this.opts.multiplanarLayout===3?e0=!0:E[4]>L[4]&&E[4]>R[4]?q=!0:L[4]>R[4]?e0=!0:z=!0,q){let H=E;this.opts.multiplanarForceRender||e[4]>=E[4]?H=e:N=!1;const j=w[0]*H[4],y=w[1]*H[4],Y=w[2]*H[4],G=P*H[4];this.draw2D([H[0],H[1],j,y],0),this.draw2D([H[0],H[1]+y+D,j,Y],1),this.draw2D([H[0],H[1]+y+D+Y+D,y,Y],2),N&&this.draw3D([H[0],H[1]+y+Y+Y+D*3,G,G])}else if(e0){let H=L;this.opts.multiplanarForceRender||d[4]>=L[4]?H=d:N=!1;const j=w[0]*H[4],y=w[1]*H[4],Y=w[2]*H[4];this.draw2D([H[0],H[1],j,y],0),this.draw2D([H[0]+j+D,H[1],j,Y],1),this.draw2D([H[0]+j+j+D*2,H[1],y,Y],2),N&&this.draw3D([H[0]+j+j+y+D*3,H[1],H[3],H[3]])}else if(z){const H=R,j=w[0]*H[4],y=w[1]*H[4],Y=w[2]*H[4];this.draw2D([H[0],H[1]+Y+D,j,y],0),this.draw2D([H[0],H[1],j,Y],1),this.draw2D([H[0]+j+D,H[1],y,Y],2),N&&this.draw3D([H[0]+j+D,H[1]+Y+D,y,y])}}if(this.opts.isRuler&&this.drawRuler(),this.opts.isColorbar&&this.drawColorbar(),A&&this.drawGraph(),this.uiData.isDragging){if(this.uiData.mouseButtonCenterDown){this.dragForCenterButton([this.uiData.dragStart[0],this.uiData.dragStart[1],this.uiData.dragEnd[0],this.uiData.dragEnd[1]]);return}if(this.opts.dragMode===4){this.dragForSlicer3D([this.uiData.dragStart[0],this.uiData.dragStart[1],this.uiData.dragEnd[0],this.uiData.dragEnd[1]]);return}if(this.opts.dragMode===3){this.dragForPanZoom([this.uiData.dragStart[0],this.uiData.dragStart[1],this.uiData.dragEnd[0],this.uiData.dragEnd[1]]);return}if(this.inRenderTile(this.uiData.dragStart[0],this.uiData.dragStart[1])>=0)return;if(this.opts.dragMode===2){this.drawMeasurementTool([this.uiData.dragStart[0],this.uiData.dragStart[1],this.uiData.dragEnd[0],this.uiData.dragEnd[1]]);return}const v=Math.abs(this.uiData.dragStart[0]-this.uiData.dragEnd[0]),w=Math.abs(this.uiData.dragStart[1]-this.uiData.dragEnd[1]);this.drawSelectionBox([Math.min(this.uiData.dragStart[0],this.uiData.dragEnd[0]),Math.min(this.uiData.dragStart[1],this.uiData.dragEnd[1]),v,w])}const S=this.frac2mm([this.scene.crosshairPos[0],this.scene.crosshairPos[1],this.scene.crosshairPos[2]]);return u=S[0].toFixed(2)+"×"+S[1].toFixed(2)+"×"+S[2].toFixed(2),this.readyForSync=!0,this.sync(),u}drawScene(){if(this.isBusy){this.needsRefresh=!0;return}this.isBusy=!1,this.needsRefresh=!1;let u=this.drawSceneCore();return this._gl!==null&&this.gl.finish(),this.needsRefresh&&(u=this.drawScene()),u}get gl(){if(!this._gl)throw new Error("unable to get WebGL context. Maybe the browser doesn't support WebGL2.");return this._gl}set gl(u){this._gl=u}};async function main(){async function T(S){await A.loadFromFile(S[0]),A.setColormap(A.volumes[2].id,"red"),lesionSlider.oninput()}openBtn.onclick=async function(){let S=document.createElement("input");S.style.display="none",S.type="file",document.body.appendChild(S),await A.removeVolume(A.volumes[2]),S.onchange=async function(v){await T(v.target.files)},S.click()},predictBtn.onclick=function(){window.alert("Outcome prediction: 0.70289")},aboutBtn.onclick=function(){window.alert("Drag and drop NIfTI images. Use pulldown menu to choose brainchop model")};function u(S){document.getElementById("intensity").innerHTML=S.string}const m={backColor:[.4,.4,.4,1],show3Dcrosshair:!0,onLocationChange:u};maskSlider.oninput=function(){A.setOpacity(1,this.value/255)},lesionSlider.oninput=function(){A.setOpacity(2,this.value/255)};const A=new Niivue(m);A.attachToCanvas(gl1),A.opts.dragMode=A.dragModes.pan,A.opts.multiplanarForceRender=!0,A.opts.yoke3Dto2DZoom=!0,A.opts.crosshairGap=11,A.setInterpolation(!0),await A.loadVolumes([{url:"./betsct1_unsmooth.nii.gz"},{url:"./mask_vox.nii.gz",colormap:"blue"},{url:"./M2095_lesion.nii.gz",colormap:"red"}]),maskSlider.oninput(),lesionSlider.oninput()}main();
+Ooh no!!`);break}v[R]=P,D=Math.min(D,P)}for(let R=0;R<A;R++)u[v[R]-1]=D}translate_labels(u,m,A,S){const v=m[0]*m[1]*m[2];let w=0;const D=new Uint32Array(v).fill(0);for(let L=0;L<S;L++)w=Math.max(w,A[L]);const R=new Uint32Array(w).fill(0);let P=0;for(let L=0;L<v;L++)u[L]&&(R[A[u[L]-1]-1]||(P+=1,R[A[u[L]-1]-1]=P),D[L]=R[A[u[L]-1]-1]);return[P,D]}largest_original_cluster_labels(u,m,A){const S=u.length,v=new Uint32Array(m+1).fill(0),w=new Uint32Array(m+1).fill(0);for(let P=0;P<S;P++){const L=u[P],d=A[P];v[d]=L,w[d]++}let D=0;for(let P=0;P<m+1;P++){const L=v[P];D=Math.max(D,L);for(let d=0;d<m+1;d++)d!==P&&L===v[d]&&(w[P]<w[d]||w[P]===w[d]&&P<d)&&(v[P]=0)}const R=new Uint32Array(S).fill(0);for(let P=0;P<S;P++)R[P]=v[A[P]];return[D,R]}bwlabel(u,m,A=26,S=!1,v=!1){const w=Date.now(),D=m[0]*m[1]*m[2],R=new Uint32Array(D).fill(0);if(![6,18,26].includes(A))return log.info("bwlabel: conn must be 6, 18 or 26."),[0,R];if(m[0]<2||m[1]<2||m[2]<1)return log.info("bwlabel: img must be 2 or 3-dimensional"),[0,R];if(S)for(let N=0;N<D;N++)u[N]!==0&&(R[N]=1);else R.set(u);let[P,L,d]=this.do_initial_labelling(R,m,A);L===void 0&&(L=new Uint32Array);const[E,e]=this.translate_labels(d,m,L,P);if(log.info(A+" neighbor clustering into "+E+" regions in "+(Date.now()-w)+"ms"),v){const[N,q]=this.largest_original_cluster_labels(R,E,e);return[N,q]}return[E,e]}async createConnectedLabelImage(u,m=26,A=!1,S=!1){var d,E;const v=this.getVolumeIndexByID(u),w=Uint32Array.from(((d=this.volumes[v].dims)==null?void 0:d.slice(1,4))??[]),D=Uint32Array.from(((E=this.volumes[v].img)==null?void 0:E.slice())??[]),[R,P]=this.bwlabel(D,w,m,A,S),L=this.volumes[v].clone();L.opacity=.5,L.colormap="random";for(let e=0;e<L.img.length;e++)L.img[e]=P[e];return L.cal_min=0,L.cal_max=R,L}async scalecropUint8(u,m=0,A=255,S,v){const w=u.length,D=new Uint8Array(w);for(let R=0;R<w;R++){let P=u[R];P=m+v*(P-S),P=Math.max(P,m),P=Math.min(P,A),D[R]=P}return D}async scalecropFloat32(u,m=0,A=1,S,v){const w=u.length,D=new Float32Array(w);for(let R=0;R<w;R++){let P=u[R];P=m+v*(P-S),P=Math.max(P,m),P=Math.min(P,A),D[R]=P}return D}getScale(u,m=0,A=255,S=0,v=.999){let w=u.global_min,D=u.global_max;if(u.hdr.datatypeCode===2)return[w,1];if((!isFinite(S)||!isFinite(v))&&isFinite(u.cal_min)&&isFinite(u.cal_max)&&u.cal_max>u.cal_min){w=u.cal_min,D=u.cal_max;const j=(A-m)/(D-w);return log.info(" Robust Rescale:  min: "+w+"  max: "+D+" scale: "+j),console.log("Robust Rescale:  min: "+w+"  max: "+D+" scale: "+j),[w,j]}const R=u.img,P=u.hdr.dims[1]*u.hdr.dims[2]*u.hdr.dims[3];if(u.hdr.scl_slope!==1||u.hdr.scl_inter!==0){const j=u.img,y=new Float32Array(u.img.length);for(let Y=0;Y<P;Y++)y[Y]=j[Y]*u.hdr.scl_slope+u.hdr.scl_inter}if(w<0&&log.warn("WARNING: Input image has value(s) below 0.0 !"),log.info(" Input:    min: "+w+"  max: "+D),S===0&&v===1)return[w,1];let L=0;for(let j=0;j<P;j++)Math.abs(R[j])>=1e-15&&L++;const d=1e3,E=(D-w)/d,e=new Array(d).fill(0);for(let j=0;j<P;j++){const y=R[j];let Y=Math.floor((y-w)/E);Y=Math.min(Y,d-1),e[Y]++}const N=new Array(d).fill(0);N[0]=e[0];for(let j=1;j<d;j++)N[j]=N[j-1]+e[j];let q=Math.floor(S*P),z=0;for(;z<d&&!(N[z]>=q);)z++;const e0=w;for(w=z*E+e0,q=P-Math.floor((1-v)*L),z=0;z<d-1&&!(N[z+1]>=q);)z++;D=z*E+e0;let H=1;return w!==D&&(H=(A-m)/(D-w)),log.info(" Rescale:  min: "+w+"  max: "+D+" scale: "+H),[w,H]}conformVox2Vox(u,m,A=256,S=1,v=!1){const w=m.flat(),D=fromValues$3(w[0],w[1],w[2],w[3],w[4],w[5],w[6],w[7],w[8],w[9],w[10],w[11],w[12],w[13],w[14],w[15]),R=fromValues$1(u[1]/2,u[2]/2,u[3]/2,1),P=create$1(),L=create$3();transpose(L,D),transformMat4(P,R,L);const d=fromValues$2(P[0],P[1],P[2]),E=fromValues$2(S,S,S);let e=fromValues$3(-1,0,0,0,0,0,1,0,0,-1,0,0,0,0,0,1);v&&(e=fromValues$3(1,0,0,0,0,1,0,0,0,0,1,0,0,0,0,1)),transpose(e,e);const N=fromValues$1(A,A,A,1),q=create$3();scale$3(q,e,E);const z=fromValues$1(N[0],N[1],N[2],1);transformMat4(z,z,q),scale$1(z,z,.5);const e0=create$2();subtract$1(e0,d,fromValues$2(z[0],z[1],z[2]));const H=create$3();transpose(H,q),H[3]=e0[0],H[7]=e0[1],H[11]=e0[2];const j=create$3();invert(j,H);const y=create$3();mul(y,D,j);const Y=create$3();return invert(Y,y),[H,y,Y]}async createNiftiArray(u=[256,256,256],m=[1,1,1],A=[1,0,0,-128,0,1,0,-128,0,0,1,-128,0,0,0,1],S=2,v=new Uint8Array){return await NVImage.createNiftiArray(u,m,A,S,v)}async niftiArray2NVImage(u=new Uint8Array){return await NVImage.loadFromUrl({url:u})}async loadFromUrl(u){return await NVImage.loadFromUrl({url:u})}async conform(u,m=!1,A=!0,S=!1,v=!1){const R=this.conformVox2Vox(u.hdr.dims,u.hdr.affine.flat(),256,1,m),P=R[0],L=R[2],d=256*256*256,E=new Float32Array(d),e=new Float32Array(u.img),N=u.hdr.dims[1]*u.hdr.dims[2]*u.hdr.dims[3];if(u.hdr.scl_slope!==1||u.hdr.scl_inter!==0)for(let o0=0;o0<N;o0++)e[o0]=e[o0]*u.hdr.scl_slope+u.hdr.scl_inter;const q=u.hdr.dims[1],z=u.hdr.dims[2],e0=u.hdr.dims[3],H=q*z;let j=-1;function y(o0,l0,c0){return o0+l0*q+c0*H}const Y=L[0],G=L[4],J=L[8];if(A)for(let o0=0;o0<256;o0++)for(let l0=0;l0<256;l0++){const c0=l0*L[1]+o0*L[2]+L[3],d0=l0*L[5]+o0*L[6]+L[7],S0=l0*L[9]+o0*L[10]+L[11];for(let p0=0;p0<256;p0++){const T0=p0*Y+c0,x0=p0*G+d0,M0=p0*J+S0,ee=Math.floor(T0),ie=Math.floor(x0),Fe=Math.floor(M0);if(j++,ee<0||ie<0||Fe<0)continue;const Oe=Math.ceil(T0),Ee=Math.ceil(x0),I0=Math.ceil(M0);if(Oe>=q||Ee>=z||I0>=e0)continue;const Be=T0-ee,Ge=x0-ie,ti=M0-Fe,ai=1-Be,ri=1-Ge,qe=1-ti,be=y(ee,ie,Fe);let J0=0;J0+=e[be]*ai*ri*qe,J0+=e[be+H]*ai*ri*ti,J0+=e[be+q]*ai*Ge*qe,J0+=e[be+q+H]*ai*Ge*ti,J0+=e[be+1]*Be*ri*qe,J0+=e[be+1+H]*Be*ri*ti,J0+=e[be+1+q]*Be*Ge*qe,J0+=e[be+1+q+H]*Be*Ge*ti,E[j]=J0}}else for(let o0=0;o0<256;o0++)for(let l0=0;l0<256;l0++){const c0=l0*L[1]+o0*L[2]+L[3],d0=l0*L[5]+o0*L[6]+L[7],S0=l0*L[9]+o0*L[10]+L[11];for(let p0=0;p0<256;p0++){const T0=Math.round(p0*Y+c0),x0=Math.round(p0*G+d0),M0=Math.round(p0*J+S0);j++,!(T0<0||x0<0||M0<0)&&(T0>=q||x0>=z||M0>=e0||(E[j]=e[y(T0,x0,M0)]))}}let i0=0;v&&(i0=NaN);let Q=new Uint8Array;if(S){const o0=await this.getScale(u,0,1,i0),l0=await this.scalecropFloat32(E,0,1,o0[0],o0[1]);Q=await this.createNiftiArray([256,256,256],[1,1,1],Array.from(P),16,new Uint8Array(l0.buffer))}else{const o0=await this.getScale(u,0,255,i0),l0=await this.scalecropUint8(E,0,255,o0[0],o0[1]);Q=await this.createNiftiArray([256,256,256],[1,1,1],Array.from(P),2,l0)}return await this.niftiArray2NVImage(Q)}setRenderDrawAmbientOcclusion(u){if(!this.renderShader)throw new Error("renderShader undefined");this.renderDrawAmbientOcclusion=u,this.renderShader.use(this.gl),this.gl.uniform1fv(this.renderShader.uniforms.renderDrawAmbientOcclusion,[this.renderDrawAmbientOcclusion,1]),this.drawScene()}setColorMap(u,m){this.setColormap(u,m)}setColormapNegative(u,m){const A=this.getVolumeIndexByID(u);this.volumes[A].colormapNegative=m,this.updateGLVolume()}setModulationImage(u,m,A=0){const S=this.getVolumeIndexByID(u);let v=null;m.length>0&&(v=this.getVolumeIndexByID(m)),this.volumes[S].modulationImage=v,this.volumes[S].modulateAlpha=A,this.updateGLVolume()}setGamma(u=1){cmapper.gamma=u,this.updateGLVolume()}async loadDeferred4DVolumes(u){const m=this.getVolumeIndexByID(u),A=this.volumes[m];if(A.nTotalFrame4D<=A.nFrame4D)return;let S;A.fileObject?S=await NVImage.loadFromFile({file:A.fileObject}):S=await NVImage.loadFromUrl({url:A.url}),S&&(A.img=S.img.slice(),A.nTotalFrame4D=S.nTotalFrame4D,A.nFrame4D=S.nFrame4D,this.updateGLVolume())}setFrame4D(u,m){const A=this.getVolumeIndexByID(u),S=this.volumes[A];m>S.nFrame4D-1&&(m=S.nFrame4D-1),m<0&&(m=0),m!==S.frame4D&&(S.frame4D=m,this.updateGLVolume(),this.onFrameChange(S,m),this.createOnLocationChange())}getFrame4D(u){const m=this.getVolumeIndexByID(u);return this.volumes[m].frame4D}colormapFromKey(u){return cmapper.colormapFromKey(u)}colormap(u="",m=!1){return cmapper.colormap(u,m)}createColormapTexture(u=null,m=0,A=256){return u!==null&&this.gl.deleteTexture(u),m<1||A<1?null:(u=this.gl.createTexture(),this.gl.activeTexture(TEXTURE1_COLORMAPS),this.gl.bindTexture(this.gl.TEXTURE_2D,u),this.gl.texStorage2D(this.gl.TEXTURE_2D,1,this.gl.RGBA8,A,m),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_MIN_FILTER,this.gl.LINEAR),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_MAG_FILTER,this.gl.LINEAR),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_WRAP_R,this.gl.CLAMP_TO_EDGE),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_WRAP_S,this.gl.CLAMP_TO_EDGE),this.gl.pixelStorei(this.gl.UNPACK_ALIGNMENT,1),u)}addColormapList(u="",m=NaN,A=NaN,S=!1,v=!1,w=!0,D=!1){u.length<1&&(w=!1),this.colormapLists.push({name:u,min:m,max:A,alphaThreshold:S,negative:v,visible:w,invert:D})}refreshColormaps(){if(this.colormapLists=[],this.volumes.length<1&&this.meshes.length<1)return;const u=this.volumes.length;if(u>0)for(let w=0;w<u;w++){const D=this.volumes[w],R=negMinMax(D.cal_min,D.cal_max,D.cal_minNeg,D.cal_maxNeg);this.addColormapList(D.colormapNegative,R[0],R[1],D.alphaThreshold!==void 0,!0,D.colorbarVisible,D.colormapInvert),this.addColormapList(D.colormap,D.cal_min,D.cal_max,D.alphaThreshold!==void 0,!1,D.colorbarVisible,D.colormapInvert)}const m=this.meshes.length;if(m>0)for(let w=0;w<m;w++){const D=this.meshes[w];if(!D.colorbarVisible)continue;const R=D.layers.length;if("edgeColormap"in D&&"edges"in D&&D.edges!==void 0){const P=negMinMax(D.edgeMin,D.edgeMax,NaN,NaN);this.addColormapList(D.edgeColormapNegative,P[0],P[1],!1,!0,!0,D.colormapInvert),this.addColormapList(D.edgeColormap,D.edgeMin,D.edgeMax,!1,!1,!0,D.colormapInvert)}if(!(R<1))for(let P=0;P<R;P++){const L=this.meshes[w].layers[P];if(!L.colorbarVisible||L.colormap.length<1)continue;const d=negMinMax(L.cal_min,L.cal_max,L.cal_minNeg,L.cal_maxNeg);this.addColormapList(L.colormapNegative,d[0],d[1],L.alphaThreshold,!0,!0,L.colormapInvert),this.addColormapList(L.colormap,L.cal_min,L.cal_max,L.alphaThreshold,!1,!0,L.colormapInvert)}}const A=this.colormapLists.length;if(A<1)return;this.colormapTexture=this.createColormapTexture(this.colormapTexture,A+1);let S=new Uint8ClampedArray;function v(w){const D=new Uint8ClampedArray(S.length+w.length);D.set(S),D.set(w,S.length),S=D}for(let w=0;w<A;w++)v(Array.from(this.colormap(this.colormapLists[w].name,this.colormapLists[w].invert)));return v(Array.from(this.drawLut.lut)),this.gl.texSubImage2D(this.gl.TEXTURE_2D,0,0,0,256,A+1,this.gl.RGBA,this.gl.UNSIGNED_BYTE,S),this}sliceScale(u=!1){var w;let m=this.screenFieldOfViewMM(0);u&&(m=this.screenFieldOfViewVox(0));const A=Math.max(m[0],Math.max(m[1],m[2])),S=[m[0]/A,m[1]/A,m[2]/A];if(!((w=this.back)!=null&&w.dims))throw new Error("back.dims undefined");const v=[this.back.dims[1],this.back.dims[2],this.back.dims[3]];return{volScale:S,vox:v,longestAxis:A,dimsMM:m}}tileIndex(u,m){for(let A=0;A<this.screenSlices.length;A++){const S=this.screenSlices[A].leftTopWidthHeight;if(u>S[0]&&m>S[1]&&u<S[0]+S[2]&&m<S[1]+S[3])return A}return-1}inRenderTile(u,m){const A=this.tileIndex(u,m);return A>=0&&this.screenSlices[A].axCorSag===4?A:-1}sliceScroll3D(u=0){if(u!==0){if(this.volumes.length>0&&this.scene.clipPlaneDepthAziElev[0]<1.8){const m=this.scene.clipPlaneDepthAziElev.slice();return u>0&&(m[0]=Math.min(1.5,m[0]+.025)),u<0&&(m[0]=Math.max(-1.5,m[0]-.025)),m[0]!==this.scene.clipPlaneDepthAziElev[0]?(this.scene.clipPlaneDepthAziElev=m,this.setClipPlane(this.scene.clipPlaneDepthAziElev)):void 0}u>0&&(this.scene.volScaleMultiplier=Math.min(2,this.scene.volScaleMultiplier*1.1)),u<0&&(this.scene.volScaleMultiplier=Math.max(.5,this.scene.volScaleMultiplier*.9)),this.drawScene()}}deleteThumbnail(){this.bmpTexture&&(this.gl.deleteTexture(this.bmpTexture),this.bmpTexture=null,this.thumbnailVisible=!1)}inGraphTile(u,m){if(this.graph.opacity<=0||this.volumes.length<1||this.volumes[0].nFrame4D<1||!this.graph.plotLTWH||this.graph.plotLTWH[2]<1||this.graph.plotLTWH[3]<1)return!1;const A=[(u-this.graph.LTWH[0])/this.graph.LTWH[2],(m-this.graph.LTWH[1])/this.graph.LTWH[3]];return A[0]>0&&A[1]>0&&A[0]<=1&&A[1]<=1}mouseClick(u,m,A=0,S=!0){if(u*=this.uiData.dpr,m*=this.uiData.dpr,this.canvas.focus(),this.thumbnailVisible){this.thumbnailVisible=!1,Promise.all([this.loadVolumes(this.deferredVolumes),this.loadMeshes(this.deferredMeshes)]).catch(v=>{throw v});return}if(this.inGraphTile(u,m)){if(!this.graph.plotLTWH)throw new Error("plotLTWH undefined");const v=[(u-this.graph.plotLTWH[0])/this.graph.plotLTWH[2],(m-this.graph.plotLTWH[1])/this.graph.plotLTWH[3]];if(v[0]>0&&v[1]>0&&v[0]<=1&&v[1]<=1){const w=Math.round(v[0]*(this.volumes[0].nFrame4D-1));this.setFrame4D(this.volumes[0].id,w);return}v[0]>.5&&v[1]>1&&this.loadDeferred4DVolumes(this.volumes[0].id).catch(w=>{throw w});return}if(this.inRenderTile(u,m)>=0){this.sliceScroll3D(A),this.drawScene();return}if(!(this.screenSlices.length<1||this.gl.canvas.height<1||this.gl.canvas.width<1))for(let v=0;v<this.screenSlices.length;v++){const w=this.screenSlices[v].axCorSag;if(this.drawPenAxCorSag>=0&&this.drawPenAxCorSag!==w||w>2)continue;const D=this.screenXY2TextureFrac(u,m,v,!1);if(!(D[0]<0)){if(!S){this.scene.crosshairPos[2-w]=A,this.drawScene();return}if(A!==0){let R=1;A<0&&(R=-1);const P=[0,0,0];P[2-w]=R,this.moveCrosshairInVox(P[0],P[1],P[2]),this.drawScene(),this.createOnLocationChange(w);return}if(this.opts.isForceMouseClickToVoxelCenters?this.scene.crosshairPos=clone$1(this.vox2frac(this.frac2vox(D))):this.scene.crosshairPos=clone$1(D),this.opts.drawingEnabled){const R=this.frac2vox(this.scene.crosshairPos);if(!isFinite(this.opts.penValue)||this.opts.penValue<0||Object.is(this.opts.penValue,-0)){isFinite(this.opts.penValue)?this.drawFloodFill(R,Math.abs(this.opts.penValue),this.opts.floodFillNeighbors):this.drawFloodFill(R,0,this.opts.penValue,this.opts.floodFillNeighbors);return}if(isNaN(this.drawPenLocation[0]))this.drawPenAxCorSag=w,this.drawPenFillPts=[],this.drawPt(...R,this.opts.penValue);else{if(R[0]===this.drawPenLocation[0]&&R[1]===this.drawPenLocation[1]&&R[2]===this.drawPenLocation[2])return;this.drawPenLine(R,this.drawPenLocation,this.opts.penValue)}this.drawPenLocation=R,this.opts.isFilledPen&&this.drawPenFillPts.push(R),this.refreshDrawing(!1)}this.drawScene(),this.createOnLocationChange(w);return}}}drawRuler(){let u=[],m=[];for(let R=0;R<this.screenSlices.length;R++)if(this.screenSlices[R].axCorSag!==4&&this.screenSlices[R].fovMM.length>1){m=this.screenSlices[R].leftTopWidthHeight,u=this.screenSlices[R].fovMM;break}if(m.length<4)return;const S=100/u[0]*m[2],v=m[0]+.5*m[2]-.5*S,w=m[1]+m[3]-2*this.opts.rulerWidth,D=[v,w,v+S,w];this.drawRuler10cm(D)}drawRuler10cm(u){if(!this.lineShader)throw new Error("lineShader undefined");this.gl.bindVertexArray(this.genericVAO),this.lineShader.use(this.gl),this.gl.uniform4fv(this.lineShader.uniforms.lineColor,this.opts.rulerColor),this.gl.uniform2fv(this.lineShader.uniforms.canvasWidthHeight,[this.gl.canvas.width,this.gl.canvas.height]),this.gl.uniform1f(this.lineShader.uniforms.thickness,this.opts.rulerWidth),this.gl.uniform4fv(this.lineShader.uniforms.startXYendXY,u),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4);const m=-.1*(u[0]-u[2]),A=u[1],S=A-2*this.opts.rulerWidth,v=A-4*this.opts.rulerWidth;for(let w=0;w<11;w++){const D=u[0]+w*m,R=[D,A,D,S];w%5===0&&(R[3]=v),this.gl.uniform4fv(this.lineShader.uniforms.startXYendXY,R),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4)}this.gl.bindVertexArray(this.unusedVAO)}screenXY2mm(u,m,A=-1){let S;for(let v=0;v<this.screenSlices.length;v++){let w=v;if(A>=0&&(w=A),this.screenSlices[w].axCorSag>2)continue;const R=this.screenSlices[w].leftTopWidthHeight;if(u<R[0]||m<R[1]||u>R[0]+R[2]||m>R[1]+R[3]||(S=this.screenXY2TextureFrac(u,m,w,!1),S[0]<0))continue;const P=this.frac2mm(S);return fromValues$1(P[0],P[1],P[2],w)}return fromValues$1(NaN,NaN,NaN,NaN)}dragForPanZoom(u){const m=this.screenXY2mm(u[2],u[3]);if(isNaN(m[0]))return;const A=this.screenXY2mm(u[0],u[1],m[3]);if(isNaN(A[0])||isNaN(m[0])||isNaN(m[3]))return;const S=create$1(),v=this.uiData.pan2DxyzmmAtMouseDown[3];sub(S,m,A),this.scene.pan2Dxyzmm[0]=this.uiData.pan2DxyzmmAtMouseDown[0]+v*S[0],this.scene.pan2Dxyzmm[1]=this.uiData.pan2DxyzmmAtMouseDown[1]+v*S[1],this.scene.pan2Dxyzmm[2]=this.uiData.pan2DxyzmmAtMouseDown[2]+v*S[2]}dragForCenterButton(u){this.dragForPanZoom(u)}dragForSlicer3D(u){let m=this.uiData.pan2DxyzmmAtMouseDown[3];const A=u[3]-u[1];m+=A*.01,m=Math.max(m,.1),m=Math.min(m,10);const v=this.scene.pan2Dxyzmm[3]-m;this.opts.yoke3Dto2DZoom&&(this.scene.volScaleMultiplier=m),this.scene.pan2Dxyzmm[3]=m;const w=this.frac2mm(this.scene.crosshairPos);this.scene.pan2Dxyzmm[0]+=v*w[0],this.scene.pan2Dxyzmm[1]+=v*w[1],this.scene.pan2Dxyzmm[2]+=v*w[2]}drawMeasurementTool(u){const m=this.gl;if(m.bindVertexArray(this.genericVAO),m.depthFunc(m.ALWAYS),m.enable(m.BLEND),m.blendFunc(m.SRC_ALPHA,m.ONE_MINUS_SRC_ALPHA),!this.lineShader)throw new Error("lineShader undefined");this.lineShader.use(this.gl),m.uniform4fv(this.lineShader.uniforms.lineColor,this.opts.rulerColor),m.uniform2fv(this.lineShader.uniforms.canvasWidthHeight,[m.canvas.width,m.canvas.height]),m.uniform1f(this.lineShader.uniforms.thickness,this.opts.rulerWidth),m.uniform4fv(this.lineShader.uniforms.startXYendXY,u),m.drawArrays(m.TRIANGLE_STRIP,0,4);const A=this.opts.rulerColor;A[3]=1,m.uniform4fv(this.lineShader.uniforms.lineColor,A);const S=this.opts.rulerWidth;m.uniform1f(this.lineShader.uniforms.thickness,S*2);let v=[u[0],u[1]-S,u[0],u[1]+S];m.uniform4fv(this.lineShader.uniforms.startXYendXY,v),m.drawArrays(m.TRIANGLE_STRIP,0,4),v=[u[2],u[3]-S,u[2],u[3]+S],m.uniform4fv(this.lineShader.uniforms.startXYendXY,v),m.drawArrays(m.TRIANGLE_STRIP,0,4);let w=this.canvasPos2frac([u[0],u[1]]),D=this.canvasPos2frac([u[2],u[3]]);if(w[0]>=0&&D[0]>=0){const R=this.frac2mm(w);w=fromValues$2(R[0],R[1],R[2]);const P=this.frac2mm(D);D=fromValues$2(P[0],P[1],P[2]);const L=create$2();sub$1(L,w,D);const d=len(L);let E=2;d>9&&(E=1),d>99&&(E=0);const e=d.toFixed(E);this.drawTextBetween(u,e,1,A)}m.bindVertexArray(this.unusedVAO)}drawRect(u,m=[1,0,0,-1]){if(m[3]<0&&(m=this.opts.crosshairColor),!this.rectShader)throw new Error("rectShader undefined");this.rectShader.use(this.gl),this.gl.enable(this.gl.BLEND),this.gl.uniform4fv(this.rectShader.uniforms.lineColor,m),this.gl.uniform2fv(this.rectShader.uniforms.canvasWidthHeight,[this.gl.canvas.width,this.gl.canvas.height]),this.gl.uniform4f(this.rectShader.uniforms.leftTopWidthHeight,u[0],u[1],u[2],u[3]),this.gl.bindVertexArray(this.genericVAO),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4),this.gl.bindVertexArray(this.unusedVAO)}drawCircle(u,m=this.opts.fontColor,A=1){if(!this.circleShader)throw new Error("circleShader undefined");this.circleShader.use(this.gl),this.gl.enable(this.gl.BLEND),this.gl.uniform4fv(this.circleShader.uniforms.circleColor,m),this.gl.uniform2fv(this.circleShader.uniforms.canvasWidthHeight,[this.gl.canvas.width,this.gl.canvas.height]),this.gl.uniform4f(this.circleShader.uniforms.leftTopWidthHeight,u[0],u[1],u[2],u[3]),this.gl.uniform1f(this.circleShader.uniforms.fillPercent,A),this.gl.uniform4fv(this.circleShader.uniforms.circleColor,m),this.gl.bindVertexArray(this.genericVAO),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4),this.gl.bindVertexArray(this.unusedVAO)}drawSelectionBox(u){this.drawRect(u,this.opts.selectionBoxColor)}effectiveCanvasHeight(){return this.gl.canvas.height-this.colorbarHeight}effectiveCanvasWidth(){return this.gl.canvas.width-this.getLegendPanelWidth()}getAllLabels(){const S=this.meshes.filter(w=>w.type==="connectome").flatMap(w=>w.nodes).map(w=>w.label).filter(w=>w!==void 0);return[...this.document.labels,...S]}getBulletMarginWidth(){let u=0;const m=this.getAllLabels();if(m.length===0)return 0;const A=m.length===1?m[0].style.bulletScale:m.reduce((w,D)=>w.style.bulletScale>D.style.bulletScale?w:D).style.bulletScale,S=m.length===1?m[0]:m.reduce((w,D)=>{const R=this.opts.textHeight*this.gl.canvas.height*w.style.textScale,P=this.opts.textHeight*this.gl.canvas.height*D.style.textScale;return this.textHeight(R,w.text)>this.textHeight(P,D.text)?w:D}),v=this.opts.textHeight*this.gl.canvas.height*S.style.textScale;return u=this.textHeight(v,S.text)*A,u+=v,u}getLegendPanelWidth(){const u=this.getAllLabels();if(!this.opts.showLegend||u.length===0)return 0;const A=this.opts.textHeight*this.gl.canvas.height*1;let S=0;const v=u.reduce((P,L)=>{const d=this.opts.textHeight*this.gl.canvas.height*P.style.textScale,E=this.opts.textHeight*this.gl.canvas.height*L.style.textScale;return this.textWidth(d,P.text)>this.textWidth(E,L.text)?P:L}),w=this.opts.textHeight*this.gl.canvas.height*v.style.textScale,D=this.textWidth(w,v.text),R=this.getBulletMarginWidth();return D&&(S=R+D,S+=A*2),S}getLegendPanelHeight(){const u=this.getAllLabels();let m=0;const S=this.opts.textHeight*this.gl.canvas.height*1;for(const v of u){const w=this.opts.textHeight*this.gl.canvas.height*v.style.textScale,D=this.textHeight(w,v.text);m+=D}return m&&(m+=S/2*(u.length+1)),m}reserveColorbarPanel(){let u=Math.max(this.opts.textHeight,.01);u=u*Math.min(this.gl.canvas.height,this.gl.canvas.width);const m=3*u,A=[0,this.gl.canvas.height-m,this.gl.canvas.width,m];return this.colorbarHeight=A[3]+1,A}drawColorbarCore(u=0,m=[0,0,0,0],A=!1,S=0,v=1,w){if(m[2]<=0||m[3]<=0)return;let D=Math.max(this.opts.textHeight,.01);D=D*Math.min(this.gl.canvas.height,this.gl.canvas.width);let R=D;const P=3*D;let L=D;if(m[3]<P){if(m[3]<3)return;R=1,L=m[3]-2}this.gl.disable(this.gl.DEPTH_TEST),this.colorbarHeight=m[3]+1;const d=[m[0]+R,m[1],m[2]-2*R,L],E=[d[0]-1,d[1]-1,d[2]+2,d[3]+2];if(this.drawRect(E,this.opts.crosshairColor),!this.colorbarShader)throw new Error("colorbarShader undefined");this.colorbarShader.use(this.gl),this.gl.activeTexture(TEXTURE1_COLORMAPS),this.gl.bindTexture(this.gl.TEXTURE_2D,this.colormapTexture),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_MIN_FILTER,this.gl.NEAREST),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_MAG_FILTER,this.gl.NEAREST);const e=u;if(this.gl.uniform1f(this.colorbarShader.uniforms.layer,e),this.gl.uniform2fv(this.colorbarShader.uniforms.canvasWidthHeight,[this.gl.canvas.width,this.gl.canvas.height]),this.gl.disable(this.gl.CULL_FACE),A){const J=[d[0]+d[2],d[1],-d[2],d[3]];this.gl.uniform4fv(this.colorbarShader.uniforms.leftTopWidthHeight,J)}else this.gl.uniform4fv(this.colorbarShader.uniforms.leftTopWidthHeight,d);this.gl.bindVertexArray(this.genericVAO),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4),this.gl.bindVertexArray(this.unusedVAO),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_MIN_FILTER,this.gl.LINEAR),this.gl.texParameteri(this.gl.TEXTURE_2D,this.gl.TEXTURE_MAG_FILTER,this.gl.LINEAR);let N=0;if(w&&v<0&&A?(N=v,v=0):w&&S>0&&(N=S,S=0),S===v||D<1)return;const q=Math.abs(v-S);let[z,e0]=tickSpacing(S,v);e0<S&&(e0+=z);function H(J){return J.toFixed(6).replace(/\.?0*$/,"")}let j=e0;const y=[0,d[1]+d[3]-D*.5,2,D*.75],Y=y[1]+y[3],G=1;for(;j<=v;){y[0]=d[0]+(j-S)/q*d[2],this.drawRect(y);const J=H(G*j);this.drawTextBelow([y[0],Y],J),j+=z}if(N!==0){const J=[d[0]+(N-S)/q*d[2],d[1]-d[3]*.25,2,d[3]*1.5];this.drawRect(J)}}drawColorbar(){const u=this.colormapLists,m=u.length;if(m<1)return;let A=0;for(let R=0;R<m;R++)u[R].visible&&A++;if(A<1)return;let S=this.reserveColorbarPanel(),v=Math.max(this.opts.textHeight,.01);v=v*Math.min(this.gl.canvas.height,this.gl.canvas.width);const w=3*v;let D=S[2]/A;(S[2]<=0||S[3]<=0)&&(D=this.gl.canvas.width/A,S=[0,this.gl.canvas.height-w,D,w]),S[2]=D;for(let R=0;R<m;R++)u[R].visible&&(this.drawColorbarCore(R,S,u[R].negative,u[R].min,u[R].max,u[R].alphaThreshold),S[0]+=D)}textWidth(u,m){if(!m)return 0;let A=0;const S=new TextEncoder().encode(m);for(let v=0;v<m.length;v++)A+=u*this.fontMets.mets[S[v]].xadv;return A}textHeight(u,m){if(!m)return 0;const A=new Set(Array.from(m)),S=new TextEncoder().encode(Array.from(A).join("")),w=Object.values(this.fontMets.mets).filter((D,R)=>S.includes(R)).reduce((D,R)=>D.lbwh[3]>R.lbwh[3]?D:R).lbwh[3];return u*w}drawChar(u,m,A){if(!this.fontShader)throw new Error("fontShader undefined");const S=this.fontMets.mets[A],v=u[0]+m*S.lbwh[0],w=-(m*S.lbwh[1]),D=m*S.lbwh[2],R=m*S.lbwh[3],P=u[1]+(w-R)+m;return this.gl.uniform4f(this.fontShader.uniforms.leftTopWidthHeight,v,P,D,R),this.gl.uniform4fv(this.fontShader.uniforms.uvLeftTopWidthHeight,S.uv_lbwh),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4),m*S.xadv}drawLoadingText(u){if(!this.canvas)throw new Error("canvas undefined");this.gl.viewport(0,0,this.gl.canvas.width,this.gl.canvas.height),this.gl.enable(this.gl.CULL_FACE),this.gl.enable(this.gl.BLEND),this.drawTextBelow([this.canvas.width/2,this.canvas.height/2],u,3)}drawText(u,m,A=1,S=null){if(this.opts.textHeight<=0)return;if(!this.fontShader)throw new Error("fontShader undefined");this.fontShader.use(this.gl);const v=this.opts.textHeight*Math.min(this.gl.canvas.height,this.gl.canvas.width)*A;this.gl.enable(this.gl.BLEND),this.gl.uniform2f(this.fontShader.uniforms.canvasWidthHeight,this.gl.canvas.width,this.gl.canvas.height),S===null&&(S=this.opts.fontColor),this.gl.uniform4fv(this.fontShader.uniforms.fontColor,S);let w=v/this.fontMets.size*this.fontMets.distanceRange;w=Math.max(w,1),this.gl.uniform1f(this.fontShader.uniforms.screenPxRange,w);const D=new TextEncoder().encode(m);this.gl.bindVertexArray(this.genericVAO);for(let R=0;R<m.length;R++)u[0]+=this.drawChar(u,v,D[R]);this.gl.bindVertexArray(this.unusedVAO)}drawTextRight(u,m,A=1,S=null){this.opts.textHeight<=0||(u[1]-=.5*this.opts.textHeight*this.gl.canvas.height,this.drawText(u,m,A,S))}drawTextLeft(u,m,A=1,S=null){if(this.opts.textHeight<=0)return;const v=this.opts.textHeight*this.gl.canvas.height*A;u[0]-=this.textWidth(v,m),u[1]-=.5*v,this.drawText(u,m,A,S)}drawTextRightBelow(u,m,A=1,S=null){this.opts.textHeight<=0||this.drawText(u,m,A,S)}drawTextBetween(u,m,A=1,S=null){if(this.opts.textHeight<=0)return;const v=[(u[0]+u[2])*.5,(u[1]+u[3])*.5],w=this.opts.textHeight*this.gl.canvas.height*A,D=this.textWidth(w,m);v[0]-=.5*D,v[1]-=.5*w;const R=[v[0]-1,v[1]-1,D+2,w+2];let P=S;P===null&&(P=this.opts.crosshairColor),P&&P[0]+P[1]+P[2]>.8?P=[0,0,0,.5]:P=[1,1,1,.5],this.drawRect(R,P),this.drawText(v,m,A,S)}drawTextBelow(u,m,A=1,S=null){if(this.opts.textHeight<=0)return;if(!this.canvas)throw new Error("canvas undefined");let v=this.opts.textHeight*this.gl.canvas.height*A,w=this.textWidth(v,m);w>this.canvas.width&&(A*=(this.canvas.width-2)/w,v=this.opts.textHeight*this.gl.canvas.height*A,w=this.textWidth(v,m)),u[0]-=.5*this.textWidth(v,m),u[0]=Math.max(u[0],1),u[0]=Math.min(u[0],this.canvas.width-w-1),this.drawText(u,m,A,S)}updateInterpolation(u,m=!1){let A=this.gl.LINEAR;!m&&this.opts.isNearestInterpolation&&(A=this.gl.NEAREST),u===0?this.gl.activeTexture(TEXTURE0_BACK_VOL):this.gl.activeTexture(TEXTURE2_OVERLAY_VOL),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_MIN_FILTER,A),this.gl.texParameteri(this.gl.TEXTURE_3D,this.gl.TEXTURE_MAG_FILTER,A)}setAtlasOutline(u){this.opts.atlasOutline=u,this.updateGLVolume(),this.drawScene()}setInterpolation(u){this.opts.isNearestInterpolation=u;const m=this.volumes.length;if(!(m<1)){for(let A=0;A<m;A++)this.updateInterpolation(A);this.drawScene()}}calculateMvpMatrix2D(u,m,A,S=1/0,v=0,w=0,D=0,R){this.gl.viewport(u[0],this.gl.canvas.height-(u[1]+u[3]),u[2],u[3]);let L=m[0],d=A[0],E=[L,m[1]],e=[d-L,A[1]-m[1]];R&&(E=[A[0],m[1]],e=[m[0]-A[0],A[1]-m[1]],L=-A[0],d=-m[0]);const N=2*Math.max(Math.abs(m[2]),Math.abs(A[2])),q=create$3();let z=.01,e0=N*8;if(S!==1/0){let J=R;D===0&&(w===0||w===180)&&(J=!J);let i0=N*1.8-v;J||(i0=N*1.8+v),z=i0-S,e0=i0+S}ortho(q,L,d,m[1],A[1],z,e0);const H=create$3();H[0]=-1;const j=fromValues$2(0,0,-N*1.8);translate(H,H,j),rotateX(H,H,deg2rad(270-D)),rotateZ(H,H,deg2rad(w-180));const y=create$3();invert(y,H);const Y=create$3();transpose(Y,y);const G=create$3();return multiply(G,q,H),{modelViewProjectionMatrix:G,modelMatrix:H,normalMatrix:Y,leftTopMM:E,fovMM:e}}swizzleVec3MM(u,m){return m===1?u=swizzleVec3(u,[0,2,1]):m===2&&(u=swizzleVec3(u,[1,2,0])),u}screenFieldOfViewVox(u=0){const m=clone$1(this.volumeObject3D.fieldOfViewDeObliqueMM);return this.swizzleVec3MM(m,u)}screenFieldOfViewMM(u=0,m=!1){if(!m&&!this.opts.isSliceMM)return this.screenFieldOfViewVox(u);const A=this.volumeObject3D.extentsMin,S=this.volumeObject3D.extentsMax;let v=fromValues$2(A[0],A[1],A[2]),w=fromValues$2(S[0],S[1],S[2]);v=this.swizzleVec3MM(v,u),w=this.swizzleVec3MM(w,u);const D=create$2();return subtract$1(D,w,v),D}screenFieldOfViewExtendedVox(u=0){const m=this.volumes[0].extentsMinOrtho,A=this.volumes[0].extentsMaxOrtho;let S=fromValues$2(m[0],m[1],m[2]),v=fromValues$2(A[0],A[1],A[2]);const w=create$3();S=this.swizzleVec3MM(S,u),v=this.swizzleVec3MM(v,u);const D=create$2();return subtract$1(D,v,S),{mnMM:S,mxMM:v,rotation:w,fovMM:D}}screenFieldOfViewExtendedMM(u=0){if(!this.volumeObject3D)throw new Error("volumeObject3D undefined");const m=this.volumeObject3D.extentsMin,A=this.volumeObject3D.extentsMax;let S=fromValues$2(m[0],m[1],m[2]),v=fromValues$2(A[0],A[1],A[2]);const w=create$3();S=this.swizzleVec3MM(S,u),v=this.swizzleVec3MM(v,u);const D=create$2();return subtract$1(D,v,S),{mnMM:S,mxMM:v,rotation:w,fovMM:D}}drawSliceOrientationText(u,m){this.gl.viewport(0,0,this.gl.canvas.width,this.gl.canvas.height);let A="S";m===0&&(A="A");let S=this.opts.isRadiologicalConvention?"R":"L";if(m===2&&(S=this.opts.sagittalNoseLeft?"A":"P"),this.opts.isCornerOrientationText){this.drawTextRightBelow([u[0],u[1]],S+A);return}this.drawTextBelow([u[0]+u[2]*.5,u[1]],A),this.drawTextRight([u[0],u[1]+u[3]*.5],S)}xyMM2xyzMM(u,m){let A=2;u===1&&(A=1),u===2&&(A=0);let S=[0,0,0],v=[1,1,0],w=[1,0,1];S[A]=m,v[A]=m,w[A]=m,S=this.frac2mm(S),v=this.frac2mm(v),w=this.frac2mm(w),S=this.swizzleVec3MM(fromValues$2(S[0],S[1],S[2]),u),v=this.swizzleVec3MM(fromValues$2(v[0],v[1],v[2]),u),w=this.swizzleVec3MM(fromValues$2(w[0],w[1],w[2]),u);const D=(v[0]-S[0])*(w[1]-S[1])-(w[0]-S[0])*(v[1]-S[1]);let R=(v[0]-S[0])*(w[2]-S[2])-(w[0]-S[0])*(v[2]-S[2]);R/=D;let P=(v[1]-S[1])*(w[2]-S[2])-(w[1]-S[1])*(v[2]-S[2]);P/=D;const L=[0,0,0,0,0];return L[0]=S[0],L[1]=S[1],L[2]=S[2],L[3]=P,L[4]=R,L}draw2D(u,m,A=NaN){let S=this.volumes[0].frac2mm.slice(),v=this.screenFieldOfViewExtendedMM(m),w=create$3();this.opts.isSliceMM||(S=this.volumes[0].frac2mmOrtho.slice(),w=clone$2(this.volumes[0].mm2ortho),v=this.screenFieldOfViewExtendedVox(m));let D=this.opts.isRadiologicalConvention&&m<2;A===1/0||A===-1/0?(D=A!==1/0,m===1&&(D=!D)):this.opts.sagittalNoseLeft&&m===2&&(D=!D);let R=0,P=0;m===2?P=D?90:-90:m===1?P=D?180:0:(P=D?180:0,R=D?-90:90);const L=this.gl;if(u[2]===0||u[3]===0){const e0=L.canvas.width/v.fovMM[0],H=L.canvas.height/v.fovMM[1],j=Math.min(e0,H),y=e0/j,Y=H/j;v.fovMM[0]*=y,v.fovMM[1]*=Y;let G=(v.mnMM[0]+v.mxMM[0])*.5;v.mnMM[0]=G-v.fovMM[0]*.5,v.mxMM[0]=G+v.fovMM[0]*.5,G=(v.mnMM[1]+v.mxMM[1])*.5,v.mnMM[1]=G-v.fovMM[1]*.5,v.mxMM[1]=G+v.fovMM[1]*.5,u=[0,0,L.canvas.width,L.canvas.height]}if(isNaN(A)){const e0=this.scene.pan2Dxyzmm,H=this.swizzleVec3MM(fromValues$2(e0[0],e0[1],e0[2]),m),j=this.scene.pan2Dxyzmm[3];v.mnMM[0]-=H[0],v.mxMM[0]-=H[0],v.mnMM[1]-=H[1],v.mxMM[1]-=H[1],v.mnMM[0]/=j,v.mxMM[0]/=j,v.mnMM[1]/=j,v.mxMM[1]/=j}let d=2;m===1&&(d=1),m===2&&(d=0);let E=this.scene.crosshairPos[d],e=this.frac2mm(this.scene.crosshairPos);!isNaN(A)&&A!==1/0&&A!==-1/0&&(e=this.frac2mm([.5,.5,.5]),e[d]=A,E=this.mm2frac(e)[d]);const N=e[d];L.clear(L.DEPTH_BUFFER_BIT);let q=this.calculateMvpMatrix2D(u,v.mnMM,v.mxMM,1/0,0,P,R,D);if(A===1/0||A===-1/0){const e0=u.slice();this.draw3D(u,q.modelViewProjectionMatrix,q.modelMatrix,q.normalMatrix,P,R);const H=this.screenSlices[this.screenSlices.length-1];H.leftTopWidthHeight=e0,H.axCorSag=m,H.sliceFrac=1/0,H.AxyzMxy=this.xyMM2xyzMM(m,E),H.leftTopMM=q.leftTopMM,H.fovMM=q.fovMM;return}L.enable(L.DEPTH_TEST),L.blendFunc(L.SRC_ALPHA,L.ONE_MINUS_SRC_ALPHA),L.disable(L.BLEND),L.depthFunc(L.GREATER),L.disable(L.CULL_FACE);let z=this.sliceMMShader;if(this.opts.isV1SliceShader&&(z=this.sliceV1Shader),!z)throw new Error("slice Shader undefined");if(z.use(this.gl),L.uniform1f(z.uniforms.overlayOutlineWidth,this.overlayOutlineWidth),L.uniform1f(z.uniforms.overlayAlphaShader,this.overlayAlphaShader),L.uniform1i(z.uniforms.isAlphaClipDark,this.isAlphaClipDark?1:0),L.uniform1i(z.uniforms.backgroundMasksOverlays,this.backgroundMasksOverlays),L.uniform1f(z.uniforms.drawOpacity,this.drawOpacity),L.enable(L.BLEND),L.blendFunc(L.SRC_ALPHA,L.ONE_MINUS_SRC_ALPHA),L.uniform1f(z.uniforms.opacity,this.volumes[0].opacity),L.uniform1i(z.uniforms.axCorSag,m),L.uniform1f(z.uniforms.slice,E),L.uniformMatrix4fv(z.uniforms.frac2mm,!1,S),L.uniformMatrix4fv(z.uniforms.mvpMtx,!1,q.modelViewProjectionMatrix.slice()),L.bindVertexArray(this.genericVAO),L.drawArrays(L.TRIANGLE_STRIP,0,4),L.bindVertexArray(this.unusedVAO),this.screenSlices.push({leftTopWidthHeight:u,axCorSag:m,sliceFrac:E,AxyzMxy:this.xyMM2xyzMM(m,E),leftTopMM:q.leftTopMM,screen2frac:[],fovMM:q.fovMM}),isNaN(A)&&this.drawCrosshairs3D(!0,1,q.modelViewProjectionMatrix,!0,this.opts.isSliceMM),this.opts.meshThicknessOn2D>0){this.opts.meshThicknessOn2D!==1/0&&(q=this.calculateMvpMatrix2D(u,v.mnMM,v.mxMM,this.opts.meshThicknessOn2D,N,P,R,D));const e0=clone$2(q.modelViewProjectionMatrix);multiply(e0,e0,w),this.drawMesh3D(!0,1,e0,q.modelMatrix,q.normalMatrix)}isNaN(A)&&this.drawCrosshairs3D(!1,.15,q.modelViewProjectionMatrix,!0,this.opts.isSliceMM),this.drawSliceOrientationText(u,m),this.readyForSync=!0}calculateMvpMatrix(u,m=[0,0,0,0],A,S){(m[2]===0||m[3]===0)&&(m=[0,0,this.gl.canvas.width,this.gl.canvas.height]);const v=m[2]/m[3];let w=this.furthestFromPivot;const D=this.pivot3D,R=create$3();w=.8*w/this.scene.volScaleMultiplier,v<1?ortho(R,-w,w,-w/v,w/v,w*.01,w*8):ortho(R,-w*v,w*v,-w,w,w*.01,w*8);const P=create$3();P[0]=-1;const L=fromValues$2(0,0,-w*1.8);translate(P,P,L),this.position&&translate(P,P,this.position),rotateX(P,P,deg2rad(270-S)),rotateZ(P,P,deg2rad(A-180)),translate(P,P,[-D[0],-D[1],-D[2]]);const d=create$3();invert(d,P);const E=create$3();transpose(E,d);const e=create$3();return multiply(e,R,P),[e,P,E]}calculateModelMatrix(u,m){if(!this.back)throw new Error("back undefined");const A=create$3();if(A[0]=-1,rotateX(A,A,deg2rad(270-m)),rotateZ(A,A,deg2rad(u-180)),this.back.obliqueRAS){const S=clone$2(this.back.obliqueRAS);multiply(A,A,S)}return A}calculateRayDirection(u,m){const A=this.calculateModelMatrix(u,m),S=fromValues$3(1,0,0,0,0,-1,0,0,0,0,-1,0,0,0,0,1),v=create$3();multiply(v,S,A);const w=create$3();invert(w,v);const D=fromValues$1(0,0,-1,1);transformMat4(D,D,w);const R=fromValues$2(D[0],D[1],D[2]);normalize$1(R,R);const P=5e-5;return Math.abs(R[0])<P&&(R[0]=P),Math.abs(R[1])<P&&(R[1]=P),Math.abs(R[2])<P&&(R[2]=P),R}sceneExtentsMinMax(u=!0){let m=fromValues$2(0,0,0),A=fromValues$2(0,0,0);if(this.volumes.length>0){if(!this.volumeObject3D)throw new Error("volumeObject3D undefined");m=fromValues$2(this.volumeObject3D.extentsMin[0],this.volumeObject3D.extentsMin[1],this.volumeObject3D.extentsMin[2]),A=fromValues$2(this.volumeObject3D.extentsMax[0],this.volumeObject3D.extentsMax[1],this.volumeObject3D.extentsMax[2]),u||(m=fromValues$2(this.volumes[0].extentsMinOrtho[0],this.volumes[0].extentsMinOrtho[1],this.volumes[0].extentsMinOrtho[2]),A=fromValues$2(this.volumes[0].extentsMaxOrtho[0],this.volumes[0].extentsMaxOrtho[1],this.volumes[0].extentsMaxOrtho[2]))}if(this.meshes.length>0){if(this.volumes.length<1){const v=this.meshes[0].extentsMin,w=this.meshes[0].extentsMax;m=fromValues$2(v[0],v[1],v[2]),A=fromValues$2(w[0],w[1],w[2])}for(let v=0;v<this.meshes.length;v++){const w=this.meshes[v].extentsMin,D=this.meshes[v].extentsMax,R=fromValues$2(w[0],w[1],w[2]);min(m,m,R);const P=fromValues$2(D[0],D[1],D[2]);max$2(A,A,P)}}const S=create$2();return subtract$1(S,A,m),[m,A,S]}setPivot3D(){const[u,m]=this.sceneExtentsMinMax(),A=create$2();add$1(A,u,m),scale$2(A,A,.5),this.pivot3D=[A[0],A[1],A[2]],subtract$1(A,m,u),this.extentsMin=u,this.extentsMax=m,this.furthestFromPivot=length$1(A)*.5}getMaxVols(){if(this.volumes.length<1)return 0;let u=0;for(let m=0;m<this.volumes.length;m++)u=Math.max(u,this.volumes[m].nFrame4D);return u}detectPartialllyLoaded4D(){if(this.volumes.length<1)return!1;for(let u=0;u<this.volumes.length;u++)if(this.volumes[u].nFrame4D<this.volumes[u].hdr.dims[4])return!0;return!1}drawGraph(){if(this.getMaxVols()<2)return;const u=this.graph;let m=0;if(this.graph.autoSizeMultiplanar&&this.opts.sliceType===3)for(let l0=0;l0<this.screenSlices.length;l0++){const c0=this.screenSlices[l0].axCorSag;if(c0===0&&(m=this.screenSlices[l0].leftTopWidthHeight[1]),c0!==2)continue;const d0=this.screenSlices[l0].leftTopWidthHeight.slice();d0[1]===m?(u.LTWH[0]=d0[0]+d0[2],u.LTWH[1]=d0[1]):(u.LTWH[0]=d0[0],u.LTWH[1]=d0[1]+d0[3]),u.LTWH[2]=d0[2],u.LTWH[3]=d0[2]}if(u.opacity<=0||u.LTWH[2]<=5||u.LTWH[3]<=5||u.LTWH[0]+u.LTWH[2]>this.gl.canvas.width||u.LTWH[1]+u.LTWH[3]>this.gl.canvas.height)return;u.backColor=[.15,.15,.15,u.opacity],u.lineColor=[1,1,1,1],this.opts.backColor[0]+this.opts.backColor[1]+this.opts.backColor[2]>1.5&&(u.backColor=[.95,.95,.95,u.opacity],u.lineColor=[0,0,0,1]),u.textColor=u.lineColor.slice(),u.lineThickness=4,u.lineAlpha=1,u.lines=[];const A=[];if(u.vols.length<1)this.volumes[0]!=null&&A.push(0);else for(let l0=0;l0<u.vols.length;l0++){const c0=u.vols[l0];this.volumes[c0]==null||this.volumes[c0].nFrame4D<2||A.push(c0)}if(A.length<1)return;const S=this.volumes[A[0]].nFrame4D;if(this.graph.selectedColumn=this.volumes[A[0]].frame4D,S<2){log.debug("Unable to generate a graph: Selected volume is 3D not 4D");return}for(let l0=0;l0<A.length;l0++){u.lines[l0]=[];const c0=this.frac2vox(this.scene.crosshairPos),d0=this.volumes[A[l0]];let S0=d0.nFrame4D;S0=Math.min(S0,S);for(let p0=0;p0<S0;p0++){const T0=d0.getValue(c0[0],c0[1],c0[2],p0);u.lines[l0].push(T0)}}u.lineRGB=[[1,0,0],[0,.7,0],[0,0,1],[1,1,0],[1,0,1],[0,1,1],[1,1,1],[0,0,0]];let v=u.lines[0][0],w=u.lines[0][0];for(let l0=0;l0<u.lines.length;l0++)for(let c0=0;c0<u.lines[l0].length;c0++){const d0=u.lines[l0][c0];v=Math.min(d0,v),w=Math.max(d0,w)}if(u.normalizeValues&&w>v){const l0=w-v;for(let c0=0;c0<u.lines.length;c0++)for(let d0=0;d0<u.lines[c0].length;d0++){const S0=u.lines[c0][d0];u.lines[c0][d0]=(S0-v)/l0}v=0,w=1}v>=w&&(w=v+1),this.drawRect(u.LTWH,u.backColor);const[D,R,P]=tickSpacing(v,w),L=Math.max(0,-1*Math.floor(Math.log(D)/Math.log(10)));v=Math.min(R,v),w=Math.max(P,w);function d(l0){return l0.toFixed(6).replace(/\.?0*$/,"")}const e=.07*(Math.min(u.LTWH[2],u.LTWH[3])/(this.fontMets.size*this.uiData.dpr));let N=this.opts.textHeight*this.gl.canvas.height*e;N<16&&(N=0);let q=0,z=R;if(N>0)for(;z<=w;){const l0=z.toFixed(L),c0=this.textWidth(N,l0);q=Math.max(c0,q),z+=D}const e0=.05,H=Math.abs(u.LTWH[2]),j=Math.abs(u.LTWH[3]),y=[u.LTWH[0]+e0*H+q,u.LTWH[1]+e0*j,u.LTWH[2]-q-2*e0*H,u.LTWH[3]-N-2*e0*j];this.graph.LTWH=u.LTWH,this.graph.plotLTWH=y,this.drawRect(y,this.opts.backColor);const Y=w-v,G=y[3]/Y,J=y[2]/(u.lines[0].length-1),i0=y[1]+y[3];z=R+.5*D;const Q=u.lineColor.slice();for(Q[3]=.25*u.lineColor[3];z<=w;){const l0=i0-(z-v)*G;this.drawLine([y[0],l0,y[0]+y[2],l0],.5*u.lineThickness,Q),z+=D}z=R;const f=.5*u.lineThickness;for(;z<=w;){const l0=i0-(z-v)*G;this.drawLine([y[0]-f,l0,y[0]+y[2]+u.lineThickness,l0],u.lineThickness,u.lineColor);const c0=z.toFixed(L);N>0&&this.drawTextLeft([y[0]-6,l0],c0,e,u.textColor),z+=D}let o0=1;for(;u.lines[0].length/o0>20;)o0*=5;for(let l0=0;l0<u.lines[0].length;l0+=o0){const c0=l0*J+y[0];let d0=u.lineThickness;if(l0%2===1)d0*=.5,this.drawLine([c0,y[1],c0,y[1]+y[3]],d0,Q);else{const S0=d(l0);N>0&&this.drawTextBelow([c0,2+y[1]+y[3]],S0,e,u.textColor),this.drawLine([c0,y[1],c0,y[1]+y[3]],d0,u.lineColor)}}for(let l0=0;l0<u.lines.length;l0++){let c0=[1,0,0,u.lineAlpha];l0<u.lineRGB.length&&(c0=[u.lineRGB[l0][0],u.lineRGB[l0][1],u.lineRGB[l0][2],u.lineAlpha]);for(let d0=1;d0<u.lines[l0].length;d0++){const S0=(d0-1)*J,p0=d0*J,T0=(u.lines[l0][d0-1]-v)*G,x0=(u.lines[l0][d0]-v)*G,M0=[y[0]+S0,y[1]+y[3]-T0,y[0]+p0,y[1]+y[3]-x0];this.drawLine(M0,u.lineThickness,c0)}}if(u.selectedColumn>=0&&u.selectedColumn<u.lines[0].length){const l0=u.selectedColumn*J+y[0];this.drawLine([l0,y[1],l0,y[1]+y[3]],u.lineThickness,[u.lineRGB[3][0],u.lineRGB[3][1],u.lineRGB[3][2],1])}this.detectPartialllyLoaded4D()&&this.drawTextBelow([y[0]+y[2],y[1]+y[3]+N*.5],"...",e,u.textColor)}depthPicker(u,m){if(!this.uiData.mouseDepthPicker)return;this.uiData.mouseDepthPicker=!1;const A=this.gl,S=this.mousePos[0]*u[2]/u[2],v=A.canvas.height-this.mousePos[1]*u[3]/u[3]-1,w=new Uint8Array(4);if(A.readPixels(S,v,1,1,A.RGBA,A.UNSIGNED_BYTE,w),this.selectedObjectId=w[3],this.selectedObjectId===this.VOLUME_ID){this.scene.crosshairPos=new Float32Array(w.slice(0,3)).map(E=>E/255);return}const D=unpackFloatFromVec4i(w);if(D>1)return;const R=(this.mousePos[0]-u[0])/u[2],P=(A.canvas.height-this.mousePos[1]-u[1])/u[3],L=unProject(R,P,D,m),d=this.mm2frac(L,0,!0);d[0]<0||d[0]>1||d[1]<0||d[1]>1||d[2]<0||d[2]>1||(this.scene.crosshairPos=this.mm2frac(L,0,!0))}drawImage3D(u,m,A){if(this.volumes.length===0)return;const S=this.gl,v=this.calculateRayDirection(m,A),w=this.volumeObject3D;if(w){S.enable(S.BLEND),S.blendFunc(S.SRC_ALPHA,S.ONE_MINUS_SRC_ALPHA),S.enable(S.CULL_FACE),S.cullFace(S.FRONT);let D=this.renderShader;if(this.uiData.mouseDepthPicker&&(D=this.pickingImageShader),D.use(this.gl),S.uniform1i(D.uniforms.backgroundMasksOverlays,this.backgroundMasksOverlays),this.gradientTextureAmount>0){S.activeTexture(TEXTURE6_GRADIENT),S.bindTexture(S.TEXTURE_3D,this.gradientTexture);const R=this.calculateModelMatrix(m,A),P=create$3();invert(P,R);const L=create$3();transpose(L,P),S.uniformMatrix4fv(D.uniforms.normMtx,!1,L)}this.drawBitmap&&this.drawBitmap.length>8?S.uniform2f(D.uniforms.renderDrawAmbientOcclusionXY,this.renderDrawAmbientOcclusion,this.drawOpacity):S.uniform2f(D.uniforms.renderDrawAmbientOcclusionXY,this.renderDrawAmbientOcclusion,0),S.uniformMatrix4fv(D.uniforms.mvpMtx,!1,u),S.uniformMatrix4fv(D.uniforms.matRAS,!1,this.back.matRAS),S.uniform3fv(D.uniforms.rayDir,v),this.gradientTextureAmount<0?S.uniform4fv(D.uniforms.clipPlane,[this.scene.crosshairPos[0],this.scene.crosshairPos[1],this.scene.crosshairPos[2],30]):S.uniform4fv(D.uniforms.clipPlane,this.scene.clipPlane),S.uniform1f(D.uniforms.drawOpacity,1),S.bindVertexArray(w.vao),S.drawElements(w.mode,w.indexCount,S.UNSIGNED_SHORT,0),S.bindVertexArray(this.unusedVAO)}}drawOrientationCube(u,m=0,A=0){if(!this.opts.isOrientCube)return;const S=.05*Math.min(u[2],u[3]);if(S<5)return;const v=this.gl;v.enable(v.CULL_FACE),v.cullFace(v.BACK),this.orientCubeShader.use(this.gl),v.bindVertexArray(this.orientCubeShaderVAO);const w=create$3(),D=create$3();ortho(D,0,v.canvas.width,0,v.canvas.height,-10*S,10*S);let R=0;u[1]===0&&(R=v.canvas.height-this.effectiveCanvasHeight()),translate(w,w,[1.8*S+u[0],R+1.8*S+u[1],0]),scale$3(w,w,[S,S,S]),rotateX(w,w,deg2rad(270-A)),rotateZ(w,w,deg2rad(-m));const P=create$3();multiply(P,D,w),v.uniformMatrix4fv(this.orientCubeShader.uniforms.u_matrix,!1,P),v.drawArrays(v.TRIANGLE_STRIP,0,168),v.bindVertexArray(this.unusedVAO),this.gl.disable(this.gl.CULL_FACE)}createOnLocationChange(u=NaN){const[m,A,S]=this.sceneExtentsMinMax(!0),v=Math.max(Math.max(S[0],S[1]),S[2]);function w(E){return Math.max(0,-Math.ceil(Math.log10(Math.abs(E))))}let D=w(v*.001);const R=this.frac2mm(this.scene.crosshairPos,0,!0);function P(E,e=0){return parseFloat(E.toFixed(e))}let L=P(R[0],D)+"×"+P(R[1],D)+"×"+P(R[2],D);if(this.volumes.length>0&&this.volumes[0].nFrame4D>0&&(L+="×"+P(this.volumes[0].frame4D)),this.volumes.length>0){let E=" = ";for(let q=0;q<this.volumes.length;q++){const z=this.volumes[q].mm2vox(R);let e0=this.volumes[q].getValue(z[0],z[1],z[2],this.volumes[q].frame4D);if(D=3,this.volumes[q].colormapLabel!==null){const H=Math.round(e0);H>=0&&H<this.volumes[q].colormapLabel.labels.length?E+=this.volumes[q].colormapLabel.labels[H]:E+="undefined("+P(e0,D)+")"}else E+=P(e0,D);this.volumes[q].imaginary&&(e0=this.volumes[q].getValue(z[0],z[1],z[2],this.volumes[q].frame4D,!0),e0>=0&&(E+="+"),E+=P(e0,D)),E+="   "}L+=E;const e=this.back.dimsRAS,N=e[1]*e[2]*e[3];if(this.drawBitmap&&this.drawBitmap.length===N){const q=this.frac2vox(this.scene.crosshairPos),z=q[0]+q[1]*e[1]+q[2]*e[1]*e[2];L+=" "+this.drawLut.labels[this.drawBitmap[z]]}}const d={mm:this.frac2mm(this.scene.crosshairPos,0,!0),axCorSag:u,vox:this.frac2vox(this.scene.crosshairPos),frac:this.scene.crosshairPos,xy:[this.mousePos[0],this.mousePos[1]],values:this.volumes.map(E=>{const e=this.frac2mm(this.scene.crosshairPos,0,!0),N=E.mm2vox(e),q=E.getValue(N[0],N[1],N[2],E.frame4D);return{name:E.name,value:q,id:E.id,mm:e,vox:N}}),string:L};this.onLocationChange(d)}addLabel(u,m,A){const S={textColor:this.opts.legendTextColor,textScale:1,textAlignment:"left",lineWidth:0,lineColor:this.opts.legendTextColor,lineTerminator:"none",bulletScale:0,bulletColor:this.opts.legendTextColor},v=m?{...S,...m}:{...S},w=new NVLabel3D(u,v,A);return this.document.labels.push(w),w}calculateScreenPoint(u,m,A){const S=create$1();return transformMat4(S,[...u,1],m),S[3]!==0&&(S[0]=(S[0]/S[3]+1)*.5*A[2],S[1]=(1-S[1]/S[3])*.5*A[3],S[2]/=S[3],S[0]+=A[0],S[1]+=A[1]),S}getLabelAtPoint(u){log.debug("screenPoint",u);const m=this.getLegendPanelHeight(),A=this.getLegendPanelWidth(),S=this.gl.canvas.width-A;let v=(this.canvas.height-m)/2;if(log.debug("panelrect",S,v,S+A,v+m),u[0]<S||u[1]<v||u[0]>S+A||u[1]>v+m)return null;const D=this.opts.textHeight*Math.min(this.gl.canvas.height,this.gl.canvas.width)*1,R=this.getAllLabels();for(const P of R){const L=this.opts.textHeight*this.gl.canvas.height*P.style.textScale,d=this.textHeight(L,P.text);if(u[1]>=v&&u[1]<=v+d+D/2)return P;v+=d,v+=D/2}return null}drawLabelLine(u,m,A,S,v=!1){const w=Array.isArray(u.points)&&Array.isArray(u.points[0])?u.points:[u.points];for(const D of w){const R=this.calculateScreenPoint(D,A,S);v?this.drawDottedLine([...m,R[0],R[1]],u.style.lineWidth,u.style.lineColor):this.draw3DLine(m,[R[0],R[1],R[2]],u.style.lineWidth,u.style.lineColor)}}draw3DLabel(u,m,A,S,v=0,w,D=!1){const R=u.text,P=m[0],L=m[1],d=this.opts.textHeight*Math.min(this.gl.canvas.height,this.gl.canvas.width)*1,E=this.textHeight(u.style.textScale,R)*d;if(u.style.lineWidth>0&&Array.isArray(u.points)&&this.drawLabelLine(u,[P,L+E],A,S,D),u.style.bulletScale){const N=u.style.bulletScale*E,q=E-N,z=L+q/2+N/2,e0=P+(v-N)/2;this.drawCircle([e0,z,N,N],u.style.bulletColor)}let e=P;if(u.style.textAlignment!=="left"){const N=this.textWidth(u.style.textScale,u.text)*d;if(u.style.textAlignment==="right")e=P+w-d*1.5-N;else{const q=w-(v||d);e+=(q-N)/2}}else e+=v;this.drawText([e,L],R,u.style.textScale,u.style.textColor)}draw3DLabels(u,m,A=!1){const S=this.getAllLabels();if(!this.opts.showLegend||S.length===0)return;if(!this.canvas)throw new Error("canvas undefined");const v=this.gl;v.disable(v.CULL_FACE),v.viewport(0,0,this.canvas.width,this.canvas.height);const D=this.opts.textHeight*Math.min(this.gl.canvas.height,this.gl.canvas.width)*1,R=this.getBulletMarginWidth(),P=this.getLegendPanelHeight(),L=this.getLegendPanelWidth(),d=v.canvas.width-L;let E=(this.canvas.height-P)/2;this.drawRect([v.canvas.width-L,E,L-D,P],this.opts.legendBackgroundColor);const e=v.getParameter(v.BLEND),N=v.getParameter(v.DEPTH_FUNC);A||(v.disable(v.BLEND),v.depthFunc(v.GREATER));for(const q of S){this.draw3DLabel(q,[d,E],u,m,R,L,A);const z=this.opts.textHeight*this.gl.canvas.height*q.style.textScale,e0=this.textHeight(z,q.text);E+=e0,E+=D/2}A||(v.depthFunc(N),e&&v.enable(v.BLEND))}draw3D(u=[0,0,0,0],m=null,A=null,S=null,v=null,w=0){const D=v!==null;this.setPivot3D(),D||(v=this.scene.renderAzimuth,w=this.scene.renderElevation);const R=this.gl;m===null&&([m,A,S]=this.calculateMvpMatrix(null,u,v,w));let P=[...u];if(u[2]===0||u[3]===0?(u=[0,0,R.canvas.width,R.canvas.height],P=[...u],this.screenSlices.push({leftTopWidthHeight:u,axCorSag:4,sliceFrac:0,AxyzMxy:[],leftTopMM:[],fovMM:[isRadiological(A),0]})):(this.screenSlices.push({leftTopWidthHeight:u.slice(),axCorSag:4,sliceFrac:0,AxyzMxy:[],leftTopMM:[],fovMM:[isRadiological(A),0]}),u[1]=R.canvas.height-u[3]-u[1]),R.enable(R.DEPTH_TEST),R.depthFunc(R.ALWAYS),R.depthMask(!0),R.clearDepth(0),this.draw3DLabels(m,P,!1),R.viewport(u[0],u[1],u[2],u[3]),this.volumes.length>0&&(this.updateInterpolation(0,!0),this.updateInterpolation(1,!0),this.drawImage3D(m,v,w)),this.updateInterpolation(0),this.updateInterpolation(1),D||this.drawCrosshairs3D(!0,1,m),this.drawMesh3D(!0,1,m,A,S),this.uiData.mouseDepthPicker){this.depthPicker(u,m),this.createOnLocationChange(),this.draw3D(u,m,A,S,v,w);return}this.opts.meshXRay>0&&this.drawMesh3D(!1,this.opts.meshXRay,m,A,S),this.draw3DLabels(m,P,!1),R.viewport(u[0],u[1],u[2],u[3]),D||this.drawCrosshairs3D(!1,.15,m),R.viewport(0,0,R.canvas.width,R.canvas.height),this.drawOrientationCube(u,v,w);const L="azimuth: "+this.scene.renderAzimuth.toFixed(0)+" elevation: "+this.scene.renderElevation.toFixed(0);return this.readyForSync=!0,this.sync(),this.draw3DLabels(m,P,!0),L}drawMesh3D(u=!0,m=1,A,S,v){if(this.meshes.length<1)return;const w=this.gl;A||([A,S,v]=this.calculateMvpMatrix(this.volumeObject3D,void 0,this.scene.renderAzimuth,this.scene.renderElevation)),w.enable(w.DEPTH_TEST),w.blendFunc(w.SRC_ALPHA,w.ONE_MINUS_SRC_ALPHA),w.disable(w.BLEND),w.depthFunc(w.GREATER),w.disable(w.CULL_FACE),u?(w.disable(w.BLEND),w.depthFunc(w.GREATER)):(w.enable(w.BLEND),w.depthFunc(w.ALWAYS),w.enable(w.CULL_FACE)),w.cullFace(w.BACK);let D=this.meshShaders[0].shader,R=!1;for(let P=0;P<this.meshes.length;P++)if(this.meshes[P].visible!==!1&&(D=this.meshShaders[this.meshes[P].meshShaderIndex].shader,this.uiData.mouseDepthPicker&&(D=this.pickingMeshShader),D.use(this.gl),w.uniformMatrix4fv(D.uniforms.mvpMtx,!1,A),w.uniformMatrix4fv(D.uniforms.normMtx,!1,v),w.uniform1f(D.uniforms.opacity,m),!(this.meshes[P].indexCount<3))){if(this.meshes[P].offsetPt0&&(this.meshes[P].fiberSides<3||this.meshes[P].fiberRadius<=0)){R=!0;continue}D.isMatcap&&(w.activeTexture(TEXTURE5_MATCAP),w.bindTexture(w.TEXTURE_2D,this.matCapTexture)),w.bindVertexArray(this.meshes[P].vao),w.drawElements(w.TRIANGLES,this.meshes[P].indexCount,w.UNSIGNED_INT,0),w.bindVertexArray(this.unusedVAO)}if(!R){w.enable(w.BLEND),w.depthFunc(w.ALWAYS);return}D=this.fiberShader,D.use(this.gl),w.uniformMatrix4fv(D.uniforms.mvpMtx,!1,A),w.uniform1f(D.uniforms.opacity,m);for(let P=0;P<this.meshes.length;P++)this.meshes[P].indexCount<3||this.meshes[P].offsetPt0&&(this.meshes[P].fiberSides>=3&&this.meshes[P].fiberRadius>0||(w.bindVertexArray(this.meshes[P].vaoFiber),w.drawElements(w.LINE_STRIP,this.meshes[P].indexCount,w.UNSIGNED_INT,0),w.bindVertexArray(this.unusedVAO)));w.enable(w.BLEND),w.depthFunc(w.ALWAYS),this.readyForSync=!0}drawCrosshairs3D(u=!0,m=1,A=null,S=!1,v=!0){if(!this.opts.show3Dcrosshair&&!S||this.opts.crosshairWidth<=0&&S)return;const w=this.gl,D=this.frac2mm(this.scene.crosshairPos,0,v);if(this.crosshairs3D===null||this.crosshairs3D.mm[0]!==D[0]||this.crosshairs3D.mm[1]!==D[1]||this.crosshairs3D.mm[2]!==D[2]){this.crosshairs3D!==null&&(w.deleteBuffer(this.crosshairs3D.indexBuffer),w.deleteBuffer(this.crosshairs3D.vertexBuffer));const[L,d,E]=this.sceneExtentsMinMax(v);let e=1;if(this.volumes.length>0){if(!this.back)throw new Error("back undefined");e=.5*Math.min(Math.min(this.back.pixDims[1],this.back.pixDims[2]),this.back.pixDims[3])}else(E[0]<50||E[0]>1e3)&&(e=E[0]*.02);e*=this.opts.crosshairWidth,this.crosshairs3D=NiivueObject3D.generateCrosshairs(this.gl,1,D,L,d,e,20,this.opts.crosshairGap),this.crosshairs3D.mm=D}if(!this.surfaceShader)throw new Error("surfaceShader undefined");const R=this.surfaceShader;R.use(this.gl),A==null&&([A]=this.calculateMvpMatrix(this.crosshairs3D,void 0,this.scene.renderAzimuth,this.scene.renderElevation)),w.uniformMatrix4fv(R.uniforms.mvpMtx,!1,A),w.bindBuffer(w.ELEMENT_ARRAY_BUFFER,this.crosshairs3D.indexBuffer),w.enable(w.DEPTH_TEST);const P=[...this.opts.crosshairColor];u?(w.disable(w.BLEND),w.depthFunc(w.GREATER)):(w.enable(w.BLEND),w.blendFunc(w.SRC_ALPHA,w.ONE_MINUS_SRC_ALPHA),w.depthFunc(w.ALWAYS)),P[3]=m,w.uniform4fv(R.uniforms.surfaceColor,P),w.bindVertexArray(this.crosshairs3D.vao),w.drawElements(w.TRIANGLES,this.crosshairs3D.indexCount,w.UNSIGNED_INT,0),w.bindVertexArray(this.unusedVAO)}mm2frac(u,m=0,A=!1){if(this.volumes.length<1){const S=fromValues$2(.1,.5,.5),[v,w,D]=this.sceneExtentsMinMax();return S[0]=(u[0]-v[0])/D[0],S[1]=(u[1]-v[1])/D[1],S[2]=(u[2]-v[2])/D[2],isFinite(S)||(isFinite(S[0])||(S[0]=.5),isFinite(S[1])||(S[1]=.5),isFinite(S[2])||(S[2]=.5),this.meshes.length<1&&log.error("mm2frac() not finite: objects not (yet) loaded.")),S}return this.volumes[m].convertMM2Frac(u,A||this.opts.isSliceMM)}vox2frac(u,m=0){return this.volumes[m].convertVox2Frac(u)}frac2vox(u,m=0){return this.volumes.length<=m?[0,0,0]:this.volumes[m].convertFrac2Vox(u)}moveCrosshairInVox(u,m,A){const S=this.frac2vox(this.scene.crosshairPos);S[0]+=u,S[1]+=m,S[2]+=A,S[0]=clamp(S[0],0,this.volumes[0].dimsRAS[1]-1),S[1]=clamp(S[1],0,this.volumes[0].dimsRAS[2]-1),S[2]=clamp(S[2],0,this.volumes[0].dimsRAS[3]-1),this.scene.crosshairPos=this.vox2frac(S),this.createOnLocationChange(),this.drawScene()}frac2mm(u,m=0,A=!1){const S=fromValues$1(u[0],u[1],u[2],1);if(this.volumes.length>0)return this.volumes[m].convertFrac2MM(u,A||this.opts.isSliceMM);{const[v,w]=this.sceneExtentsMinMax(),D=(R,P,L)=>R*(1-L)+P*L;S[0]=D(v[0],w[0],u[0]),S[1]=D(v[1],w[1],u[1]),S[2]=D(v[2],w[2],u[2])}return S}screenXY2TextureFrac(u,m,A,S=!0){const v=fromValues$2(-1,-1,-1),w=this.screenSlices[A].axCorSag;if(w>2)return v;const D=this.screenSlices[A].leftTopWidthHeight.slice();let R=!1;D[2]<0&&(R=!0,D[0]+=D[2],D[2]=-D[2]);let P=(u-D[0])/D[2];R&&(P=1-P);const L=1-(m-D[1])/D[3];if(P<0||P>1||L<0||L>1||this.screenSlices[A].AxyzMxy.length<4)return v;let d=fromValues$2(0,0,0);d[0]=this.screenSlices[A].leftTopMM[0]+P*this.screenSlices[A].fovMM[0],d[1]=this.screenSlices[A].leftTopMM[1]+L*this.screenSlices[A].fovMM[1];const E=this.screenSlices[A].AxyzMxy;d[2]=E[2]+E[4]*(d[1]-E[1])-E[3]*(d[0]-E[0]),w===1&&(d=swizzleVec3(d,[0,2,1])),w===2&&(d=swizzleVec3(d,[2,0,1]));const e=this.mm2frac(d);return S&&(e[0]<0||e[0]>1||e[1]<0||e[1]>1||e[2]<0||e[2]>1)?v:e}canvasPos2frac(u){for(let m=0;m<this.screenSlices.length;m++){const A=this.screenXY2TextureFrac(u[0],u[1],m);if(A[0]>=0)return A}return[-1,-1,-1]}scaleSlice(u,m,A=0,S=0){const v=this.effectiveCanvasWidth()-A,w=this.effectiveCanvasHeight()-S;let D=v/u;m*D>w&&(D=w/m);const R=u*D,P=m*D;return[(v-R)*.5,(w-P)*.5,R,P,D]}drawThumbnail(){if(!this.bmpShader)throw new Error("bmpShader undefined");this.bmpShader.use(this.gl),this.gl.uniform2f(this.bmpShader.uniforms.canvasWidthHeight,this.gl.canvas.width,this.gl.canvas.height);let u=this.gl.canvas.height,m=this.gl.canvas.height*this.bmpTextureWH;m>this.gl.canvas.width&&(u=this.gl.canvas.width/this.bmpTextureWH,m=this.gl.canvas.width),this.gl.uniform4f(this.bmpShader.uniforms.leftTopWidthHeight,0,0,m,u),this.gl.bindVertexArray(this.genericVAO),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4),this.gl.bindVertexArray(this.unusedVAO)}drawLine(u,m=1,A=[1,0,0,-1]){if(this.gl.bindVertexArray(this.genericVAO),!this.lineShader)throw new Error("lineShader undefined");this.lineShader.use(this.gl),A[3]<0&&(A=this.opts.crosshairColor),this.gl.uniform4fv(this.lineShader.uniforms.lineColor,A),this.gl.uniform2fv(this.lineShader.uniforms.canvasWidthHeight,[this.gl.canvas.width,this.gl.canvas.height]),this.gl.uniform1f(this.lineShader.uniforms.thickness,m),this.gl.uniform4fv(this.lineShader.uniforms.startXYendXY,u),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4),this.gl.bindVertexArray(this.unusedVAO)}draw3DLine(u,m,A=1,S=[1,0,0,-1]){if(this.gl.bindVertexArray(this.genericVAO),!this.line3DShader)throw new Error("line3DShader undefined");this.line3DShader.use(this.gl),S[3]<0&&(S=this.opts.crosshairColor),this.gl.uniform4fv(this.line3DShader.uniforms.lineColor,S),this.gl.uniform2fv(this.line3DShader.uniforms.canvasWidthHeight,[this.gl.canvas.width,this.gl.canvas.height]),this.gl.uniform1f(this.line3DShader.uniforms.thickness,A),this.gl.uniform2fv(this.line3DShader.uniforms.startXY,u),this.gl.uniform3fv(this.line3DShader.uniforms.endXYZ,m),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4),this.gl.bindVertexArray(this.unusedVAO)}drawDottedLine(u,m=1,A=[1,0,0,-1]){if(this.gl.bindVertexArray(this.genericVAO),!this.lineShader)throw new Error("lineShader undefined");this.lineShader.use(this.gl);const S=A[3]<0?[...this.opts.crosshairColor]:[...A];S[3]=.3;const v=fromValues(u[2]-u[0],u[3]-u[1]),w=length(v);normalize(v,v);const R=this.opts.textHeight*Math.min(this.gl.canvas.height,this.gl.canvas.width)*1;scale(v,v,R/2);const P=length(v);let L=Math.floor(w/P);w%P&&L++;const d=[u[0],u[1]];this.gl.uniform4fv(this.lineShader.uniforms.lineColor,S),this.gl.uniform2fv(this.lineShader.uniforms.canvasWidthHeight,[this.gl.canvas.width,this.gl.canvas.height]),this.gl.uniform1f(this.lineShader.uniforms.thickness,m);for(let E=0;E<L-1;E++){if(E%2){d[0]+=v[0],d[1]+=v[1];continue}const e=[d[0],d[1],d[0]+v[0],d[1]+v[1]];this.gl.uniform4fv(this.lineShader.uniforms.startXYendXY,e),this.gl.drawArrays(this.gl.TRIANGLE_STRIP,0,4),d[0]+=v[0],d[1]+=v[1]}this.gl.bindVertexArray(this.unusedVAO)}drawGraphLine(u,m=[1,0,0,.5],A=2){this.drawLine(u,A,m)}drawCrossLinesMM(u,m,A,S,v){if(u<0||this.screenSlices.length<=u)return;const w=this.screenSlices[u];let D=w.sliceFrac;D===1/0&&log.warn("Rendering approximate cross lines in world view mode"),D===1/0&&(D=.5);let P=S.slice(),L=v.slice();const d=Math.max(1,this.opts.crosshairWidth);m===1&&(P=A.slice()),m===2&&(P=A.slice(),L=S.slice());function E(e){const N=fromValues(0,0);return N[0]=w.leftTopWidthHeight[0]+(e[0]-w.leftTopMM[0])/w.fovMM[0]*w.leftTopWidthHeight[2],N[1]=w.leftTopWidthHeight[1]+w.leftTopWidthHeight[3]-(e[1]-w.leftTopMM[1])/w.fovMM[1]*w.leftTopWidthHeight[3],N}if(P.length>0&&m===0){const e=D,N=1;for(let q=0;q<P.length;q++){const z=this.frac2mm([.5,.5,.5]);z[N]=P[q];let e0=this.mm2frac(z);e0=e0[N];let H=this.frac2mm([0,e0,e]);H=swizzleVec3(H,[0,1,2]);let j=this.frac2mm([1,e0,e]);j=swizzleVec3(j,[0,1,2]),H=E(H),j=E(j),this.drawLine([H[0],H[1],j[0],j[1]],d)}}if(P.length>0&&m===1){const e=D,N=2;for(let q=0;q<P.length;q++){const z=this.frac2mm([.5,.5,.5]);z[N]=P[q];let e0=this.mm2frac(z);e0=e0[N];let H=this.frac2mm([0,e,e0]);H=swizzleVec3(H,[0,2,1]);let j=this.frac2mm([1,e,e0]);j=swizzleVec3(j,[0,2,1]),H=E(H),j=E(j),this.drawLine([H[0],H[1],j[0],j[1]],d)}}if(P.length>0&&m===2){const e=D,N=2;for(let q=0;q<P.length;q++){const z=this.frac2mm([.5,.5,.5]);z[N]=P[q];let e0=this.mm2frac(z);e0=e0[N];let H=this.frac2mm([e,0,e0]);H=swizzleVec3(H,[1,2,0]);let j=this.frac2mm([e,1,e0]);j=swizzleVec3(j,[1,2,0]),H=E(H),j=E(j),this.drawLine([H[0],H[1],j[0],j[1]],d)}}if(L.length>0&&m===0){const e=D,N=0;for(let q=0;q<L.length;q++){const z=this.frac2mm([.5,.5,.5]);z[N]=L[q];let e0=this.mm2frac(z);e0=e0[N];let H=this.frac2mm([e0,0,e]);H=swizzleVec3(H,[0,1,2]);let j=this.frac2mm([e0,1,e]);j=swizzleVec3(j,[0,1,2]),H=E(H),j=E(j),this.drawLine([H[0],H[1],j[0],j[1]],d)}}if(L.length>0&&m===1){const e=D,N=0;for(let q=0;q<L.length;q++){const z=this.frac2mm([.5,.5,.5]);z[N]=L[q];let e0=this.mm2frac(z);e0=e0[N];let H=this.frac2mm([e0,e,0]);H=swizzleVec3(H,[0,2,1]);let j=this.frac2mm([e0,e,1]);j=swizzleVec3(j,[0,2,1]),H=E(H),j=E(j),this.drawLine([H[0],H[1],j[0],j[1]],d)}}if(L.length>0&&m===2){const e=D,N=1;for(let q=0;q<L.length;q++){const z=this.frac2mm([.5,.5,.5]);z[N]=L[q];let e0=this.mm2frac(z);e0=e0[N];let H=this.frac2mm([e,e0,0]);H=swizzleVec3(H,[1,2,0]);let j=this.frac2mm([e,e0,1]);j=swizzleVec3(j,[1,2,0]),H=E(H),j=E(j),this.drawLine([H[0],H[1],j[0],j[1]],d)}}}drawCrossLines(u,m,A,S,v){if(u<0||this.screenSlices.length<=u)return;if(this.opts.isSliceMM)return this.drawCrossLinesMM(u,m,A,S,v);if(this.screenSlices[u].sliceFrac===1/0)return this.drawCrossLinesMM(u,m,A,S,v);const w=this.screenSlices[u];let D=S.slice(),R=v.slice();if(m===1&&(D=A.slice()),m===2&&(D=A.slice(),R=S.slice()),D.length>0){const P=w.leftTopWidthHeight.slice();let L=2;m===0&&(L=1);const d=this.frac2mm([.5,.5,.5]);for(let E=0;E<D.length;E++){d[L]=D[E];const e=this.mm2frac(d);this.drawRect([P[0],P[1]+P[3]-e[L]*P[3],P[2],1])}}if(R.length>0){const P=w.leftTopWidthHeight.slice(),L=w.fovMM[0]<0;let d=0;m===2&&(d=1);const E=this.frac2mm([.5,.5,.5]);for(let e=0;e<R.length;e++){E[d]=R[e];const N=this.mm2frac(E);L?this.drawRect([P[0]+(P[2]-N[d]*P[2]),P[1],1,P[3]]):this.drawRect([P[0]+N[d]*P[2],P[1],1,P[3]])}}}drawMosaic(u){if(this.volumes.length===0){log.debug("Unable to draw mosaic until voxel-based image is loaded");return}this.screenSlices=[];const m=this.screenFieldOfViewMM(0,!0),A=this.screenFieldOfViewMM(0);u=u.replaceAll(";"," ;").trim();const S=[],v=[],w=[],D=u.split(/\s+/);let R=1;const P=this.opts.textHeight;for(let L=0;L<2;L++){let d=!1,E=!1;d=!1;let e=0,N=0,q=0,z=0,e0=!1,H=0;for(let Y=0;Y<D.length;Y++){const G=D[Y];if(G.includes("X")){E=!0;continue}if(G.includes("L")){e0=!G.includes("-");continue}if(G.includes("V")||G.includes("H")){Y++;continue}G.includes("A")&&(H=0),G.includes("C")&&(H=1),G.includes("S")&&(H=2),G.includes("R")&&(d=!0),G.includes(";")&&(q+=e,z=Math.max(z,N),e=0,N=0);const J=parseFloat(G);if(isNaN(J))continue;let i0=0,Q=0,f=A;if(d&&(f=m),H===2?i0=f[1]:i0=f[0],H===0?Q=f[1]:Q=f[2],L===0)d||(H===0&&S.push(J),H===1&&v.push(J),H===2&&w.push(J));else{const o0=[R*N,R*q,R*i0,R*Q];if(this.opts.textHeight=e0?P:0,d){let l0=J<0?-1/0:1/0;Object.is(J,-0)&&(l0=-1/0),this.draw2D(o0,H,l0)}else this.draw2D(o0,H,J);E&&this.drawCrossLines(this.screenSlices.length-1,H,S,v,w),d=!1,E=!1}N+=i0,e=Math.max(e,Q)}if(q+=e,z=Math.max(z,N),z<=0||q<=0)break;const j=this.gl.canvas.width/z,y=this.effectiveCanvasHeight()/q;R=Math.min(j,y)}this.opts.textHeight=P}drawSceneCore(){if(!this.initialized)return;if(this.colorbarHeight=0,this.gl.clearColor(this.opts.backColor[0],this.opts.backColor[1],this.opts.backColor[2],this.opts.backColor[3]),this.gl.clear(this.gl.COLOR_BUFFER_BIT|this.gl.DEPTH_BUFFER_BIT),this.bmpTexture&&this.thumbnailVisible){this.drawThumbnail();return}let u="";if(this.volumes.length===0||typeof this.volumes[0].dims>"u"){if(this.meshes.length>0){this.screenSlices=[],this.opts.sliceType=4,this.draw3D(),this.opts.isColorbar&&this.drawColorbar();return}this.drawLoadingText(this.loadingText);return}if(this.back===null)return;if(this.uiData.isDragging&&this.scene.clipPlaneDepthAziElev[0]<1.8&&this.inRenderTile(this.uiData.dragStart[0],this.uiData.dragStart[1])>=0){const v=this.uiData.dragStart[0]-this.uiData.dragEnd[0],w=this.uiData.dragStart[1]-this.uiData.dragEnd[1],D=this.uiData.dragClipPlaneStartDepthAziElev.slice();if(D[1]-=v,D[1]=D[1]%360,D[2]+=w,D[1]!==this.scene.clipPlaneDepthAziElev[1]||D[2]!==this.scene.clipPlaneDepthAziElev[2])return this.scene.clipPlaneDepthAziElev=D,this.setClipPlane(this.scene.clipPlaneDepthAziElev)}if(this.sliceMosaicString.length<1&&this.opts.sliceType===4){this.opts.isColorbar&&this.reserveColorbarPanel(),this.screenSlices=[],this.draw3D(),this.opts.isColorbar&&this.drawColorbar();return}this.opts.isColorbar&&this.reserveColorbarPanel();const m=this.getMaxVols(),A=this.opts.sliceType===3&&m>1&&this.graph.autoSizeMultiplanar&&this.graph.opacity>0;if(this.sliceMosaicString.length>0)this.drawMosaic(this.sliceMosaicString);else if(this.gl.viewport(0,0,this.gl.canvas.width,this.gl.canvas.height),this.screenSlices=[],this.opts.sliceType===0)this.draw2D([0,0,0,0],0);else if(this.opts.sliceType===1)this.draw2D([0,0,0,0],1);else if(this.opts.sliceType===2)this.draw2D([0,0,0,0],2);else{const v=isFinite(this.drawPenLocation[0])&&this.opts.drawingEnabled,{volScale:w}=this.sliceScale();typeof this.opts.multiplanarPadPixels!="number"&&log.debug("multiplanarPadPixels must be numeric");const D=parseFloat(`${this.opts.multiplanarPadPixels}`),R=this.scaleSlice(w[0]+w[1],w[1]+w[2],D*1,D*1),P=Math.max(Math.max(w[1],w[2]),w[0]),L=this.scaleSlice(w[0]+w[0]+w[1],Math.max(w[1],w[2]),D*2),d=this.scaleSlice(w[0]+w[0]+w[1]+P,Math.max(w[1],w[2]),D*3),E=this.scaleSlice(P,w[1]+w[2]+w[2],0,D*2),e=this.scaleSlice(P,w[1]+w[2]+w[2]+P,0,D*3);let N=!v&&(m<2||!A),q=!1,z=!1,e0=!1;if(this.opts.multiplanarLayout===1?q=!0:this.opts.multiplanarLayout===2?z=!0:this.opts.multiplanarLayout===3?e0=!0:E[4]>L[4]&&E[4]>R[4]?q=!0:L[4]>R[4]?e0=!0:z=!0,q){let H=E;this.opts.multiplanarForceRender||e[4]>=E[4]?H=e:N=!1;const j=w[0]*H[4],y=w[1]*H[4],Y=w[2]*H[4],G=P*H[4];this.draw2D([H[0],H[1],j,y],0),this.draw2D([H[0],H[1]+y+D,j,Y],1),this.draw2D([H[0],H[1]+y+D+Y+D,y,Y],2),N&&this.draw3D([H[0],H[1]+y+Y+Y+D*3,G,G])}else if(e0){let H=L;this.opts.multiplanarForceRender||d[4]>=L[4]?H=d:N=!1;const j=w[0]*H[4],y=w[1]*H[4],Y=w[2]*H[4];this.draw2D([H[0],H[1],j,y],0),this.draw2D([H[0]+j+D,H[1],j,Y],1),this.draw2D([H[0]+j+j+D*2,H[1],y,Y],2),N&&this.draw3D([H[0]+j+j+y+D*3,H[1],H[3],H[3]])}else if(z){const H=R,j=w[0]*H[4],y=w[1]*H[4],Y=w[2]*H[4];this.draw2D([H[0],H[1]+Y+D,j,y],0),this.draw2D([H[0],H[1],j,Y],1),this.draw2D([H[0]+j+D,H[1],y,Y],2),N&&this.draw3D([H[0]+j+D,H[1]+Y+D,y,y])}}if(this.opts.isRuler&&this.drawRuler(),this.opts.isColorbar&&this.drawColorbar(),A&&this.drawGraph(),this.uiData.isDragging){if(this.uiData.mouseButtonCenterDown){this.dragForCenterButton([this.uiData.dragStart[0],this.uiData.dragStart[1],this.uiData.dragEnd[0],this.uiData.dragEnd[1]]);return}if(this.opts.dragMode===4){this.dragForSlicer3D([this.uiData.dragStart[0],this.uiData.dragStart[1],this.uiData.dragEnd[0],this.uiData.dragEnd[1]]);return}if(this.opts.dragMode===3){this.dragForPanZoom([this.uiData.dragStart[0],this.uiData.dragStart[1],this.uiData.dragEnd[0],this.uiData.dragEnd[1]]);return}if(this.inRenderTile(this.uiData.dragStart[0],this.uiData.dragStart[1])>=0)return;if(this.opts.dragMode===2){this.drawMeasurementTool([this.uiData.dragStart[0],this.uiData.dragStart[1],this.uiData.dragEnd[0],this.uiData.dragEnd[1]]);return}const v=Math.abs(this.uiData.dragStart[0]-this.uiData.dragEnd[0]),w=Math.abs(this.uiData.dragStart[1]-this.uiData.dragEnd[1]);this.drawSelectionBox([Math.min(this.uiData.dragStart[0],this.uiData.dragEnd[0]),Math.min(this.uiData.dragStart[1],this.uiData.dragEnd[1]),v,w])}const S=this.frac2mm([this.scene.crosshairPos[0],this.scene.crosshairPos[1],this.scene.crosshairPos[2]]);return u=S[0].toFixed(2)+"×"+S[1].toFixed(2)+"×"+S[2].toFixed(2),this.readyForSync=!0,this.sync(),u}drawScene(){if(this.isBusy){this.needsRefresh=!0;return}this.isBusy=!1,this.needsRefresh=!1;let u=this.drawSceneCore();return this._gl!==null&&this.gl.finish(),this.needsRefresh&&(u=this.drawScene()),u}get gl(){if(!this._gl)throw new Error("unable to get WebGL context. Maybe the browser doesn't support WebGL2.");return this._gl}set gl(u){this._gl=u}};async function main(){async function T(R){await v.loadFromFile(R[0]),v.setColormap(v.volumes[2].id,"red"),lesionSlider.oninput()}function u(R,P,L){let d=L-P,E=(R-P)/d;return Math.min(Math.max(E,0),1)}function m(){if(v.volumes.length!==3){window.alert('Please reload the page and open lesions with the "Open Lesion Map" button');return}let R=v.volumes[1].img,P=v.volumes[2].img,L=P.length;if(L!==R.length){window.alert("Lesion must precisely match mask");return}let d=0,E=0;const e=5,N=Math.floor(w.length/e);let q=new Float64Array(L);for(let y=0;y<L;y++)P[y]>0&&R[y]>0&&d++,R[y]>1&&(q[E++]=P[y]);if(N!==E){window.alert("PCA and mask have different number of elements");return}let z=new Float64Array(e);for(let y=0;y<E;y++)q[y]-=D[y];for(let y=0;y<e;y++){const Y=y*E;for(let G=0;G<E;G++)z[y]+=q[G]*w[G+Y]}for(let y=0;y<e;y++)z[y]=u(z[y],-51.9073,110.0535);let e0=parseFloat(cocNumber.value);e0=u(e0,-.024243014,.951938077);let H=d/1e3;H=u(H,0,21.625);const j=[z[0],z[1],z[2],z[3],z[4],e0,H];console.log(j)}openBtn.onclick=async function(){let R=document.createElement("input");R.style.display="none",R.type="file",document.body.appendChild(R),await v.removeVolume(v.volumes[2]),R.onchange=async function(P){await T(P.target.files)},R.click()},predictBtn.onclick=function(){window.alert("Outcome prediction: 0.70289")},aboutBtn.onclick=function(){window.alert("Drag and drop NIfTI images. Use pulldown menu to choose brainchop model")};function A(R){document.getElementById("intensity").innerHTML=R.string}const S={backColor:[.4,.4,.4,1],show3Dcrosshair:!0,onLocationChange:A};maskSlider.oninput=function(){v.setOpacity(1,this.value/255)},lesionSlider.oninput=function(){v.setOpacity(2,this.value/255)};const v=new Niivue(S);v.attachToCanvas(gl1);const w=(await v.loadFromUrl("./pca_values_coeff.nii.gz")).img,D=(await v.loadFromUrl("./pca_values_mu.nii.gz")).img;v.opts.dragMode=v.dragModes.pan,v.opts.multiplanarForceRender=!0,v.opts.yoke3Dto2DZoom=!0,v.opts.crosshairGap=11,v.setInterpolation(!0),await v.loadVolumes([{url:"./betsct1_unsmooth.nii.gz"},{url:"./mask.nii.gz",colormap:"blue"},{url:"./M2095_lesion.nii.gz",colormap:"red"}]),maskSlider.oninput(),lesionSlider.oninput(),m()}main();
diff --git a/index.html b/index.html
index 9b0b6dd..26a5b05 100644
--- a/index.html
+++ b/index.html
@@ -6,7 +6,7 @@
   <meta name="viewport" content="width=device-width, initial-scale=1.0" />
   <link rel="stylesheet" href="./niivue.css" />
   <title>Niivue Neglect Predictions</title>
-  <script type="module" crossorigin src="./assets/index-DgVXdzB4.js"></script>
+  <script type="module" crossorigin src="./assets/index-jh3uDHfm.js"></script>
 </head>
 
 <body>
diff --git a/mask.nii.gz b/mask.nii.gz
new file mode 100644
index 0000000..5700fda
Binary files /dev/null and b/mask.nii.gz differ
diff --git a/models_5x10_diff.mat b/models_5x10_diff.mat
new file mode 100644
index 0000000..6f28513
Binary files /dev/null and b/models_5x10_diff.mat differ
diff --git a/neglect_predict.m b/neglect_predict.m
new file mode 100644
index 0000000..829c65b
--- /dev/null
+++ b/neglect_predict.m
@@ -0,0 +1,100 @@
+function neglect_predict(fnm, acuteCoC)
+%Predict chronic recovery based on lesion map and acute center of
+%cancellation score.
+%Examples
+% neglect_predict; %use GUI;
+% neglect_predict('M2095_lesion.nii.gz', 0.65)
+    if nargin < 1
+        [p_file, p_path] = uigetfile('*.nii.gz;*.nii', 'Select lesion map');
+        if p_file==0
+            return
+        end
+        fnm = fullfile(p_path, p_file);
+    end
+    if nargin < 2
+        prompt = {'Enter acute CoC score (-1..1):'};
+        dlgtitle = 'Input';
+        fieldsize = [1];
+        definput = {'0.65'};
+        answer = inputdlg(prompt,dlgtitle,fieldsize,definput);
+        if length(answer) == 0
+            return
+        end
+        acuteCoC = str2double(answer{1});
+    end
+    lesion = niftiread(fnm);
+    mpath = fileparts(mfilename("fullpath"));
+    fnmMsk = fullfile(mpath, 'mask.nii.gz');
+    if ~exist(fnmMsk,'file')
+        error('Unable to find %s', fnmMsk)
+    end
+    maskVox012 = uint8(niftiread(fnmMsk));
+    maskROI = uint8(maskVox012 > 0);
+    maskVox = uint8(maskVox012 > 1);
+    ROI_volVox = nnz(maskROI & lesion);
+    ROI_volML = ROI_volVox / 1000; %convert voxels to ML
+    fprintf("%d lesioned voxels in ROI mask: %g ml\n", ROI_volVox, ROI_volML);
+    %PCA
+    fnmPCA = fullfile(mpath, 'pca_values_5x21220.mat');
+    if ~exist(fnmPCA,'file')
+        error('Unable to find %s', fnmPCA)
+    end
+    pca_val = load(fnmPCA).pca_val;
+    if size(pca_val.mu,2) ~= nnz(maskVox)
+        error('maskVox size does not match pca_val')
+    end
+    map_org = double(lesion(:));
+    mask_log = maskVox ~=0;
+    map = map_org(mask_log);
+    PC = (map'-pca_val.mu)*pca_val.coeff;
+    if numel(PC) ~= 5
+        error('scores_new should have 5 values')
+    end
+    %normalize values to range 0..1
+    %       -> PCs: min = -51.9073; max = 110.0535
+    %       -> CoC: min = -0.024243014; max = 0.951938077
+    %       -> ROI_vol: min = 0; max = 21.625
+    for i = 1:5
+        PC(i) = norm0to1(PC(i), -51.9073, 110.0535);
+    end
+   acuteCoC = norm0to1(acuteCoC, -0.024243014, 0.951938077);
+    ROI_volML = norm0to1(ROI_volML, 0, 21.625);
+    % input_vector = [PC1, PC2, PC3, PC4, PC5, CoC, ROI_vol];
+    % Exemplary patient data (with outcome = 0.8155):
+    input_vector = [PC(1), PC(2), PC(3), PC(4), PC(5), acuteCoC, ROI_volML];
+    fnmModel = fullfile(mpath, 'models_5x10_diff.mat');
+    if ~exist(fnmPCA,'file')
+        error('Unable to find %s', fnmPCA)
+    end
+    models = load(fnmModel).models;
+    % Blank prediction-array
+    predictions_mdls = zeros(numel(models),1);
+    
+    % For each single model (5-fold nested cross validation x 10 model repetitions)
+    for i = 1:numel(models)
+        % Extract information from the i-th model
+        support_vectors_i = models{i}.SVs;
+        coefficients_i = models{i}.sv_coef;
+        bias_i = -models{i}.rho;
+        % Define RBF kernel function
+        gamma_i = models{i}.Parameters(4);
+        rbf_kernel = @(x1, x2) exp(-gamma_i * sum((x1 - x2).^2));
+        % Calculate the kernel values
+        kernel_values_i = arrayfun(@(j) rbf_kernel(input_vector, support_vectors_i(j, :)), 1:size(support_vectors_i, 1));
+        % Calculate the prediction using the regression function
+        predictions_mdls(i) = sum(coefficients_i' .* kernel_values_i) + bias_i;
+        % Feature weights and bias term
+        % w = coefficients_i' * support_vectors_i; 
+        % b = bias_i;
+    end
+    % Calculate mean prediction
+    prediction_mean = mean(predictions_mdls);
+    disp(['Mean prediction: ' num2str(prediction_mean)]);
+end
+function ret = norm0to1(val, mn, mx)
+    %return normalized 0..1, linearly interpolated min..max
+    range = mx - mn;
+    ret = (val - mn)/range;
+    ret = min(max(ret,0),1);
+end
+
diff --git a/pca_values_5x21220.mat b/pca_values_5x21220.mat
new file mode 100644
index 0000000..8c838f7
Binary files /dev/null and b/pca_values_5x21220.mat differ
diff --git a/pca_values_coeff.nii.gz b/pca_values_coeff.nii.gz
new file mode 100644
index 0000000..a00910f
Binary files /dev/null and b/pca_values_coeff.nii.gz differ
diff --git a/pca_values_mu.nii.gz b/pca_values_mu.nii.gz
new file mode 100644
index 0000000..ef9ddb9
Binary files /dev/null and b/pca_values_mu.nii.gz differ