@fideus-labs/fidnii 0.1.0

package/dist/ResolutionSelector.js

@@ -0,0 +1,224 @@
+// SPDX-FileCopyrightText: Copyright (c) Fideus Labs LLC
+// SPDX-License-Identifier: MIT
+import { clipPlanesToPixelRegion } from "./ClipPlanes.js";
+/**
+ * Select the appropriate resolution level based on pixel budget and clip planes.
+ *
+ * The selection process:
+ * 1. Starts from the highest resolution (level 0)
+ * 2. Finds the highest resolution that fits within maxPixels
+ * 3. Considers the clipped region size, not full volume
+ *
+ * @param multiscales - The OME-Zarr multiscales data
+ * @param maxPixels - Maximum number of pixels to use
+ * @param clipPlanes - Current clip planes in world space
+ * @param volumeBounds - Full volume bounds in world space
+ * @param viewportBounds - Optional viewport bounds (for viewport-aware mode)
+ * @returns The selected resolution level and buffer dimensions
+ */
+export function selectResolution(multiscales, maxPixels, clipPlanes, volumeBounds, viewportBounds) {
+    const images = multiscales.images;
+    // Try each resolution from highest to lowest
+    for (let i = 0; i < images.length; i++) {
+        const image = images[i];
+        const region = clipPlanesToPixelRegion(clipPlanes, volumeBounds, image, viewportBounds);
+        const alignedRegion = alignRegionToChunks(region, image);
+        const dimensions = [
+            alignedRegion.end[0] - alignedRegion.start[0],
+            alignedRegion.end[1] - alignedRegion.start[1],
+            alignedRegion.end[2] - alignedRegion.start[2],
+        ];
+        const pixelCount = dimensions[0] * dimensions[1] * dimensions[2];
+        if (pixelCount <= maxPixels) {
+            return {
+                levelIndex: i,
+                dimensions,
+                pixelCount,
+            };
+        }
+    }
+    // Fall back to lowest resolution
+    const lowestImage = images[images.length - 1];
+    const region = clipPlanesToPixelRegion(clipPlanes, volumeBounds, lowestImage, viewportBounds);
+    const alignedRegion = alignRegionToChunks(region, lowestImage);
+    const dimensions = [
+        alignedRegion.end[0] - alignedRegion.start[0],
+        alignedRegion.end[1] - alignedRegion.start[1],
+        alignedRegion.end[2] - alignedRegion.start[2],
+    ];
+    return {
+        levelIndex: images.length - 1,
+        dimensions,
+        pixelCount: dimensions[0] * dimensions[1] * dimensions[2],
+    };
+}
+/**
+ * Get the chunk shape for a volume.
+ *
+ * @param ngffImage - The NgffImage to get chunk shape from
+ * @returns Chunk shape as [z, y, x]
+ */
+export function getChunkShape(ngffImage) {
+    const chunks = ngffImage.data.chunks;
+    const dims = ngffImage.dims;
+    // Find z, y, x indices in dims
+    const zIdx = dims.indexOf("z");
+    const yIdx = dims.indexOf("y");
+    const xIdx = dims.indexOf("x");
+    if (zIdx === -1 || yIdx === -1 || xIdx === -1) {
+        // Fallback: assume last 3 dimensions are z, y, x
+        const n = chunks.length;
+        return [chunks[n - 3] || 1, chunks[n - 2] || 1, chunks[n - 1] || 1];
+    }
+    return [chunks[zIdx], chunks[yIdx], chunks[xIdx]];
+}
+/**
+ * Get the shape of a volume as [z, y, x].
+ *
+ * @param ngffImage - The NgffImage
+ * @returns Shape as [z, y, x]
+ */
+export function getVolumeShape(ngffImage) {
+    const shape = ngffImage.data.shape;
+    const dims = ngffImage.dims;
+    // Find z, y, x indices in dims
+    const zIdx = dims.indexOf("z");
+    const yIdx = dims.indexOf("y");
+    const xIdx = dims.indexOf("x");
+    if (zIdx === -1 || yIdx === -1 || xIdx === -1) {
+        // Fallback: assume last 3 dimensions are z, y, x
+        const n = shape.length;
+        return [shape[n - 3] || 1, shape[n - 2] || 1, shape[n - 1] || 1];
+    }
+    return [shape[zIdx], shape[yIdx], shape[xIdx]];
+}
+/**
+ * Align a pixel region to chunk boundaries.
+ * Expands the region to include complete chunks.
+ *
+ * @param region - The pixel region to align
+ * @param ngffImage - The NgffImage (for chunk shape)
+ * @returns Aligned region
+ */
+export function alignRegionToChunks(region, ngffImage) {
+    const chunkShape = getChunkShape(ngffImage);
+    const volumeShape = getVolumeShape(ngffImage);
+    // Align start down to chunk boundary
+    const alignedStart = [
+        Math.floor(region.start[0] / chunkShape[0]) * chunkShape[0],
+        Math.floor(region.start[1] / chunkShape[1]) * chunkShape[1],
+        Math.floor(region.start[2] / chunkShape[2]) * chunkShape[2],
+    ];
+    // Align end up to chunk boundary (but don't exceed volume size)
+    const alignedEnd = [
+        Math.min(Math.ceil(region.end[0] / chunkShape[0]) * chunkShape[0], volumeShape[0]),
+        Math.min(Math.ceil(region.end[1] / chunkShape[1]) * chunkShape[1], volumeShape[1]),
+        Math.min(Math.ceil(region.end[2] / chunkShape[2]) * chunkShape[2], volumeShape[2]),
+    ];
+    return {
+        start: alignedStart,
+        end: alignedEnd,
+    };
+}
+/**
+ * Calculate the middle resolution level index.
+ *
+ * @param multiscales - The OME-Zarr multiscales data
+ * @returns Middle resolution level index
+ */
+export function getMiddleResolutionIndex(multiscales) {
+    return Math.floor(multiscales.images.length / 2);
+}
+/**
+ * Calculate upsample factor from one resolution level to another.
+ *
+ * @param multiscales - The OME-Zarr multiscales data
+ * @param fromLevel - Source resolution level
+ * @param toLevel - Target resolution level (should be higher resolution, i.e., lower index)
+ * @returns Upsample factor for each dimension [z, y, x]
+ */
+export function calculateUpsampleFactor(multiscales, fromLevel, toLevel) {
+    const fromImage = multiscales.images[fromLevel];
+    const toImage = multiscales.images[toLevel];
+    const fromShape = getVolumeShape(fromImage);
+    const toShape = getVolumeShape(toImage);
+    return [
+        toShape[0] / fromShape[0],
+        toShape[1] / fromShape[1],
+        toShape[2] / fromShape[2],
+    ];
+}
+/**
+ * Get dimensions for the full volume at a given resolution level.
+ */
+export function getFullVolumeDimensions(multiscales, levelIndex) {
+    return getVolumeShape(multiscales.images[levelIndex]);
+}
+/**
+ * Select the appropriate resolution level for a 2D slice view.
+ *
+ * Unlike `selectResolution` which counts all 3D pixels (z*y*x), this function
+ * counts only the 2D in-plane pixels (e.g., x*y for axial), ignoring the
+ * orthogonal axis. This allows much higher resolution for 2D slice views.
+ *
+ * The slab dimensions returned include one chunk of thickness in the
+ * orthogonal direction (needed for zarr fetching efficiency). The pixel
+ * budget accounts for the full 3D slab volume (in-plane area multiplied
+ * by the chunk thickness along the orthogonal axis).
+ *
+ * @param multiscales - The OME-Zarr multiscales data
+ * @param maxPixels - Maximum number of voxels for the slab (in-plane area × chunk depth)
+ * @param clipPlanes - Current clip planes in world space
+ * @param volumeBounds - Full volume bounds in world space
+ * @param orthogonalAxis - The axis perpendicular to the slice plane (0=Z, 1=Y, 2=X)
+ * @param viewportBounds - Optional viewport bounds (for viewport-aware mode)
+ * @returns The selected resolution level and slab dimensions
+ */
+export function select2DResolution(multiscales, maxPixels, clipPlanes, volumeBounds, orthogonalAxis, viewportBounds) {
+    const images = multiscales.images;
+    // Try each resolution from highest to lowest
+    for (let i = 0; i < images.length; i++) {
+        const image = images[i];
+        const region = clipPlanesToPixelRegion(clipPlanes, volumeBounds, image, viewportBounds);
+        const alignedRegion = alignRegionToChunks(region, image);
+        const dimensions = [
+            alignedRegion.end[0] - alignedRegion.start[0],
+            alignedRegion.end[1] - alignedRegion.start[1],
+            alignedRegion.end[2] - alignedRegion.start[2],
+        ];
+        // Count the full slab volume: in-plane area × chunk depth along the
+        // orthogonal axis. The slab is always one chunk thick, so we use the
+        // chunk shape rather than the full volume extent in that axis.
+        const chunkShape = getChunkShape(image);
+        const slabDepth = chunkShape[orthogonalAxis];
+        const inPlaneAxes = [0, 1, 2].filter((a) => a !== orthogonalAxis);
+        const slabVoxelCount = dimensions[inPlaneAxes[0]] *
+            dimensions[inPlaneAxes[1]] * slabDepth;
+        if (slabVoxelCount <= maxPixels) {
+            return {
+                levelIndex: i,
+                dimensions,
+                pixelCount: slabVoxelCount,
+            };
+        }
+    }
+    // Fall back to lowest resolution
+    const lowestImage = images[images.length - 1];
+    const region = clipPlanesToPixelRegion(clipPlanes, volumeBounds, lowestImage, viewportBounds);
+    const alignedRegion = alignRegionToChunks(region, lowestImage);
+    const dimensions = [
+        alignedRegion.end[0] - alignedRegion.start[0],
+        alignedRegion.end[1] - alignedRegion.start[1],
+        alignedRegion.end[2] - alignedRegion.start[2],
+    ];
+    const chunkShape = getChunkShape(lowestImage);
+    const slabDepth = chunkShape[orthogonalAxis];
+    const inPlaneAxes = [0, 1, 2].filter((a) => a !== orthogonalAxis);
+    return {
+        levelIndex: images.length - 1,
+        dimensions,
+        pixelCount: dimensions[inPlaneAxes[0]] * dimensions[inPlaneAxes[1]] *
+            slabDepth,
+    };
+}
+//# sourceMappingURL=ResolutionSelector.js.map

package/dist/ResolutionSelector.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ResolutionSelector.js","sourceRoot":"","sources":["../src/ResolutionSelector.ts"],"names":[],"mappings":"AAAA,wDAAwD;AACxD,+BAA+B;AAS/B,OAAO,EAAE,uBAAuB,EAAE,MAAM,iBAAiB,CAAC;AAQ1D;;;;;;;;;;;;;;GAcG;AACH,MAAM,UAAU,gBAAgB,CAC9B,WAAwB,EACxB,SAAiB,EACjB,UAAsB,EACtB,YAA0B,EAC1B,cAA6B;IAE7B,MAAM,MAAM,GAAG,WAAW,CAAC,MAAM,CAAC;IAElC,6CAA6C;IAC7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACvC,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QACxB,MAAM,MAAM,GAAG,uBAAuB,CACpC,UAAU,EACV,YAAY,EACZ,KAAK,EACL,cAAc,CACf,CAAC;QACF,MAAM,aAAa,GAAG,mBAAmB,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;QAEzD,MAAM,UAAU,GAA6B;YAC3C,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC;YAC7C,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC;YAC7C,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC;SAC9C,CAAC;QAEF,MAAM,UAAU,GAAG,UAAU,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;QAEjE,IAAI,UAAU,IAAI,SAAS,EAAE,CAAC;YAC5B,OAAO;gBACL,UAAU,EAAE,CAAC;gBACb,UAAU;gBACV,UAAU;aACX,CAAC;QACJ,CAAC;IACH,CAAC;IAED,iCAAiC;IACjC,MAAM,WAAW,GAAG,MAAM,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAC9C,MAAM,MAAM,GAAG,uBAAuB,CACpC,UAAU,EACV,YAAY,EACZ,WAAW,EACX,cAAc,CACf,CAAC;IACF,MAAM,aAAa,GAAG,mBAAmB,CAAC,MAAM,EAAE,WAAW,CAAC,CAAC;IAE/D,MAAM,UAAU,GAA6B;QAC3C,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC;QAC7C,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC;QAC7C,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC;KAC9C,CAAC;IAEF,OAAO;QACL,UAAU,EAAE,MAAM,CAAC,MAAM,GAAG,CAAC;QAC7B,UAAU;QACV,UAAU,EAAE,UAAU,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC;KAC1D,CAAC;AACJ,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,aAAa,CAAC,SAAoB;IAChD,MAAM,MAAM,GAAG,SAAS,CAAC,IAAI,CAAC,MAAM,CAAC;IACrC,MAAM,IAAI,GAAG,SAAS,CAAC,IAAI,CAAC;IAE5B,+BAA+B;IAC/B,MAAM,IAAI,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;IAC/B,MAAM,IAAI,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;IAC/B,MAAM,IAAI,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;IAE/B,IAAI,IAAI,KAAK,CAAC,CAAC,IAAI,IAAI,KAAK,CAAC,CAAC,IAAI,IAAI,KAAK,CAAC,CAAC,EAAE,CAAC;QAC9C,iDAAiD;QACjD,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;QACxB,OAAO,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;IACtE,CAAC;IAED,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,EAAE,MAAM,CAAC,IAAI,CAAC,EAAE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC;AACpD,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,cAAc,CAAC,SAAoB;IACjD,MAAM,KAAK,GAAG,SAAS,CAAC,IAAI,CAAC,KAAK,CAAC;IACnC,MAAM,IAAI,GAAG,SAAS,CAAC,IAAI,CAAC;IAE5B,+BAA+B;IAC/B,MAAM,IAAI,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;IAC/B,MAAM,IAAI,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;IAC/B,MAAM,IAAI,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;IAE/B,IAAI,IAAI,KAAK,CAAC,CAAC,IAAI,IAAI,KAAK,CAAC,CAAC,IAAI,IAAI,KAAK,CAAC,CAAC,EAAE,CAAC;QAC9C,iDAAiD;QACjD,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC;QACvB,OAAO,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;IACnE,CAAC;IAED,OAAO,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC;AACjD,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,mBAAmB,CACjC,MAAmB,EACnB,SAAoB;IAEpB,MAAM,UAAU,GAAG,aAAa,CAAC,SAAS,CAAC,CAAC;IAC5C,MAAM,WAAW,GAAG,cAAc,CAAC,SAAS,CAAC,CAAC;IAE9C,qCAAqC;IACrC,MAAM,YAAY,GAA6B;QAC7C,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC;QAC3D,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC;QAC3D,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC;KAC5D,CAAC;IAEF,gEAAgE;IAChE,MAAM,UAAU,GAA6B;QAC3C,IAAI,CAAC,GAAG,CACN,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,EACxD,WAAW,CAAC,CAAC,CAAC,CACf;QACD,IAAI,CAAC,GAAG,CACN,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,EACxD,WAAW,CAAC,CAAC,CAAC,CACf;QACD,IAAI,CAAC,GAAG,CACN,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,EACxD,WAAW,CAAC,CAAC,CAAC,CACf;KACF,CAAC;IAEF,OAAO;QACL,KAAK,EAAE,YAAY;QACnB,GAAG,EAAE,UAAU;KAChB,CAAC;AACJ,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,wBAAwB,CAAC,WAAwB;IAC/D,OAAO,IAAI,CAAC,KAAK,CAAC,WAAW,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;AACnD,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,uBAAuB,CACrC,WAAwB,EACxB,SAAiB,EACjB,OAAe;IAEf,MAAM,SAAS,GAAG,WAAW,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC;IAChD,MAAM,OAAO,GAAG,WAAW,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC;IAE5C,MAAM,SAAS,GAAG,cAAc,CAAC,SAAS,CAAC,CAAC;IAC5C,MAAM,OAAO,GAAG,cAAc,CAAC,OAAO,CAAC,CAAC;IAExC,OAAO;QACL,OAAO,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC;QACzB,OAAO,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC;QACzB,OAAO,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC;KAC1B,CAAC;AACJ,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,uBAAuB,CACrC,WAAwB,EACxB,UAAkB;IAElB,OAAO,cAAc,CAAC,WAAW,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC;AACxD,CAAC;AAED;;;;;;;;;;;;;;;;;;;GAmBG;AACH,MAAM,UAAU,kBAAkB,CAChC,WAAwB,EACxB,SAAiB,EACjB,UAAsB,EACtB,YAA0B,EAC1B,cAA8B,EAC9B,cAA6B;IAE7B,MAAM,MAAM,GAAG,WAAW,CAAC,MAAM,CAAC;IAElC,6CAA6C;IAC7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACvC,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QACxB,MAAM,MAAM,GAAG,uBAAuB,CACpC,UAAU,EACV,YAAY,EACZ,KAAK,EACL,cAAc,CACf,CAAC;QACF,MAAM,aAAa,GAAG,mBAAmB,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;QAEzD,MAAM,UAAU,GAA6B;YAC3C,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC;YAC7C,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC;YAC7C,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC;SAC9C,CAAC;QAEF,oEAAoE;QACpE,qEAAqE;QACrE,+DAA+D;QAC/D,MAAM,UAAU,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC;QACxC,MAAM,SAAS,GAAG,UAAU,CAAC,cAAc,CAAC,CAAC;QAC7C,MAAM,WAAW,GAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAW,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CACpD,CAAC,KAAK,cAAc,CACrB,CAAC;QACF,MAAM,cAAc,GAAG,UAAU,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC;YAC/C,UAAU,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC;QAEzC,IAAI,cAAc,IAAI,SAAS,EAAE,CAAC;YAChC,OAAO;gBACL,UAAU,EAAE,CAAC;gBACb,UAAU;gBACV,UAAU,EAAE,cAAc;aAC3B,CAAC;QACJ,CAAC;IACH,CAAC;IAED,iCAAiC;IACjC,MAAM,WAAW,GAAG,MAAM,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAC9C,MAAM,MAAM,GAAG,uBAAuB,CACpC,UAAU,EACV,YAAY,EACZ,WAAW,EACX,cAAc,CACf,CAAC;IACF,MAAM,aAAa,GAAG,mBAAmB,CAAC,MAAM,EAAE,WAAW,CAAC,CAAC;IAE/D,MAAM,UAAU,GAA6B;QAC3C,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC;QAC7C,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC;QAC7C,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC;KAC9C,CAAC;IAEF,MAAM,UAAU,GAAG,aAAa,CAAC,WAAW,CAAC,CAAC;IAC9C,MAAM,SAAS,GAAG,UAAU,CAAC,cAAc,CAAC,CAAC;IAC7C,MAAM,WAAW,GAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAW,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,KAAK,cAAc,CAAC,CAAC;IAE7E,OAAO;QACL,UAAU,EAAE,MAAM,CAAC,MAAM,GAAG,CAAC;QAC7B,UAAU;QACV,UAAU,EAAE,UAAU,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC;YACjE,SAAS;KACZ,CAAC;AACJ,CAAC"}

package/dist/ViewportBounds.d.ts

@@ -0,0 +1,50 @@
+import type { Niivue } from "@niivue/niivue";
+import { SLICE_TYPE } from "@niivue/niivue";
+import type { VolumeBounds } from "./types.js";
+/**
+ * Intersect two axis-aligned bounding boxes.
+ * Returns the overlapping region, clamped so max >= min.
+ */
+export declare function intersectBounds(a: VolumeBounds, b: VolumeBounds): VolumeBounds;
+/**
+ * Compute the world-space axis-aligned bounding box of the visible region
+ * in a 3D render view.
+ *
+ * NiiVue uses an orthographic projection whose extents depend on
+ * `volScaleMultiplier`, `renderAzimuth`, and `renderElevation`.
+ * We build the 8 corners of the ortho frustum, rotate them by the inverse
+ * of the view rotation, and take the AABB.
+ *
+ * @param nv - The Niivue instance
+ * @param volumeBounds - Full volume bounds to intersect with
+ * @returns Viewport bounds in world space, intersected with volumeBounds
+ */
+export declare function computeViewportBounds3D(nv: Niivue, volumeBounds: VolumeBounds): VolumeBounds;
+/**
+ * Compute the world-space bounding box of the visible region in a 2D slice
+ * view, accounting for pan and zoom.
+ *
+ * NiiVue's 2D slice renderer applies pan (`pan2Dxyzmm[0..2]`) and zoom
+ * (`pan2Dxyzmm[3]`) to the base field of view. We replicate this math to
+ * determine what mm range is visible on screen.
+ *
+ * @param nv - The Niivue instance
+ * @param sliceType - Current 2D slice type (AXIAL, CORONAL, SAGITTAL)
+ * @param volumeBounds - Full volume bounds to intersect with
+ * @param normalizationScale - If the NVImage affine was normalized (multiplied
+ *   by 1/maxVoxelSize to avoid NiiVue precision issues), pass that scale here
+ *   so we can convert the NiiVue mm-space FOV back to physical world
+ *   coordinates. Pass 1.0 if no normalization was applied.
+ * @returns Viewport bounds in world space, intersected with volumeBounds
+ */
+export declare function computeViewportBounds2D(nv: Niivue, sliceType: SLICE_TYPE, volumeBounds: VolumeBounds, normalizationScale?: number): VolumeBounds;
+/**
+ * Check if two VolumeBounds are approximately equal.
+ *
+ * @param a - First bounds
+ * @param b - Second bounds
+ * @param tolerance - Relative tolerance (default: 0.01 = 1%)
+ * @returns True if bounds are within tolerance
+ */
+export declare function boundsApproxEqual(a: VolumeBounds, b: VolumeBounds, tolerance?: number): boolean;
+//# sourceMappingURL=ViewportBounds.d.ts.map

package/dist/ViewportBounds.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ViewportBounds.d.ts","sourceRoot":"","sources":["../src/ViewportBounds.ts"],"names":[],"mappings":"AAGA,OAAO,KAAK,EAAE,MAAM,EAAE,MAAM,gBAAgB,CAAC;AAC7C,OAAO,EAAE,UAAU,EAAE,MAAM,gBAAgB,CAAC;AAC5C,OAAO,KAAK,EAAE,YAAY,EAAE,MAAM,YAAY,CAAC;AAE/C;;;GAGG;AACH,wBAAgB,eAAe,CAC7B,CAAC,EAAE,YAAY,EACf,CAAC,EAAE,YAAY,GACd,YAAY,CAwBd;AAED;;;;;;;;;;;;GAYG;AACH,wBAAgB,uBAAuB,CACrC,EAAE,EAAE,MAAM,EACV,YAAY,EAAE,YAAY,GACzB,YAAY,CA2Gd;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,uBAAuB,CACrC,EAAE,EAAE,MAAM,EACV,SAAS,EAAE,UAAU,EACrB,YAAY,EAAE,YAAY,EAC1B,kBAAkB,GAAE,MAAY,GAC/B,YAAY,CA6Jd;AAED;;;;;;;GAOG;AACH,wBAAgB,iBAAiB,CAC/B,CAAC,EAAE,YAAY,EACf,CAAC,EAAE,YAAY,EACf,SAAS,GAAE,MAAa,GACvB,OAAO,CAUT"}

package/dist/ViewportBounds.js

@@ -0,0 +1,325 @@
+// SPDX-FileCopyrightText: Copyright (c) Fideus Labs LLC
+// SPDX-License-Identifier: MIT
+import { SLICE_TYPE } from "@niivue/niivue";
+/**
+ * Intersect two axis-aligned bounding boxes.
+ * Returns the overlapping region, clamped so max >= min.
+ */
+export function intersectBounds(a, b) {
+    const min = [
+        Math.max(a.min[0], b.min[0]),
+        Math.max(a.min[1], b.min[1]),
+        Math.max(a.min[2], b.min[2]),
+    ];
+    const max = [
+        Math.min(a.max[0], b.max[0]),
+        Math.min(a.max[1], b.max[1]),
+        Math.min(a.max[2], b.max[2]),
+    ];
+    // Ensure valid bounds (max >= min)
+    return {
+        min: [
+            Math.min(min[0], max[0]),
+            Math.min(min[1], max[1]),
+            Math.min(min[2], max[2]),
+        ],
+        max: [
+            Math.max(min[0], max[0]),
+            Math.max(min[1], max[1]),
+            Math.max(min[2], max[2]),
+        ],
+    };
+}
+/**
+ * Compute the world-space axis-aligned bounding box of the visible region
+ * in a 3D render view.
+ *
+ * NiiVue uses an orthographic projection whose extents depend on
+ * `volScaleMultiplier`, `renderAzimuth`, and `renderElevation`.
+ * We build the 8 corners of the ortho frustum, rotate them by the inverse
+ * of the view rotation, and take the AABB.
+ *
+ * @param nv - The Niivue instance
+ * @param volumeBounds - Full volume bounds to intersect with
+ * @returns Viewport bounds in world space, intersected with volumeBounds
+ */
+export function computeViewportBounds3D(nv, volumeBounds) {
+    // Get scene extents: [min, max, range]
+    const extents = nv.sceneExtentsMinMax(true);
+    const mn = extents[0]; // vec3
+    const mx = extents[1]; // vec3
+    const range = extents[2]; // vec3
+    // Pivot = center of scene
+    const pivotX = (mn[0] + mx[0]) * 0.5;
+    const pivotY = (mn[1] + mx[1]) * 0.5;
+    const pivotZ = (mn[2] + mx[2]) * 0.5;
+    // furthestFromPivot = half-diagonal of bounding box
+    const furthest = Math.sqrt(range[0] * range[0] + range[1] * range[1] + range[2] * range[2]) * 0.5;
+    // NiiVue's orthographic scale (matches calculateMvpMatrix)
+    const scale = (0.8 * furthest) / (nv.scene.volScaleMultiplier || 1);
+    // Canvas aspect ratio
+    const canvas = nv.canvas;
+    const canvasW = canvas?.width ?? 1;
+    const canvasH = canvas?.height ?? 1;
+    const whratio = canvasW / canvasH;
+    // Ortho extents in view space (before rotation)
+    let halfW, halfH;
+    if (whratio < 1) {
+        // Portrait
+        halfW = scale;
+        halfH = scale / whratio;
+    }
+    else {
+        // Landscape
+        halfW = scale * whratio;
+        halfH = scale;
+    }
+    // For viewport-aware resolution, we need the world-space extent that is
+    // visible on screen. Rather than rotating a full 3D frustum box (whose depth
+    // dominates the AABB after rotation), we compute the extent differently:
+    //
+    // Project the view-space axes onto world space and accumulate the half-extents
+    // contributed by each view axis. The view X axis (with extent halfW) and view
+    // Y axis (with extent halfH) determine what is visible on screen. The depth
+    // axis does NOT constrain visibility — we see through the full depth of the
+    // volume — so we include the full volume extent along the view Z axis.
+    // Build the inverse of NiiVue's view rotation.
+    // NiiVue applies: rotateX(270 - elevation) then rotateZ(azimuth - 180)
+    // Also mirrors X (modelMatrix[0] = -1).
+    // We need the inverse rotation to go from view space back to world space.
+    const azimuth = nv.scene.renderAzimuth ?? 0;
+    const elevation = nv.scene.renderElevation ?? 0;
+    const azRad = ((azimuth - 180) * Math.PI) / 180;
+    const elRad = ((270 - elevation) * Math.PI) / 180;
+    const cosAz = Math.cos(azRad);
+    const sinAz = Math.sin(azRad);
+    const cosEl = Math.cos(elRad);
+    const sinEl = Math.sin(elRad);
+    // Compute inverse rotation of unit view-space axes to world space.
+    // Inverse: un-mirror X, then rotateZ(-azRad), then rotateX(-elRad).
+    //
+    // View X axis (1, 0, 0) after un-mirror: (-1, 0, 0)
+    const viewXinWorld = [
+        -cosAz, // wx
+        sinAz * cosEl, // wy
+        -sinAz * sinEl, // wz
+    ];
+    // View Y axis (0, 1, 0) — no mirror effect
+    const viewYinWorld = [
+        sinAz, // wx
+        cosAz * cosEl, // wy
+        -cosAz * sinEl, // wz
+    ];
+    // View Z axis (0, 0, 1) — no mirror effect
+    const viewZinWorld = [
+        0, // wx
+        sinEl, // wy
+        cosEl, // wz
+    ];
+    // For each world axis, the visible half-extent is:
+    //   halfW * |viewXinWorld[axis]| + halfH * |viewYinWorld[axis]|
+    //   + furthest * |viewZinWorld[axis]|  (full volume depth along view Z)
+    const worldHalfExtent = [0, 0, 0];
+    for (let axis = 0; axis < 3; axis++) {
+        worldHalfExtent[axis] = halfW * Math.abs(viewXinWorld[axis]) +
+            halfH * Math.abs(viewYinWorld[axis]) +
+            furthest * Math.abs(viewZinWorld[axis]);
+    }
+    const frustumBounds = {
+        min: [
+            pivotX - worldHalfExtent[0],
+            pivotY - worldHalfExtent[1],
+            pivotZ - worldHalfExtent[2],
+        ],
+        max: [
+            pivotX + worldHalfExtent[0],
+            pivotY + worldHalfExtent[1],
+            pivotZ + worldHalfExtent[2],
+        ],
+    };
+    return intersectBounds(frustumBounds, volumeBounds);
+}
+/**
+ * Compute the world-space bounding box of the visible region in a 2D slice
+ * view, accounting for pan and zoom.
+ *
+ * NiiVue's 2D slice renderer applies pan (`pan2Dxyzmm[0..2]`) and zoom
+ * (`pan2Dxyzmm[3]`) to the base field of view. We replicate this math to
+ * determine what mm range is visible on screen.
+ *
+ * @param nv - The Niivue instance
+ * @param sliceType - Current 2D slice type (AXIAL, CORONAL, SAGITTAL)
+ * @param volumeBounds - Full volume bounds to intersect with
+ * @param normalizationScale - If the NVImage affine was normalized (multiplied
+ *   by 1/maxVoxelSize to avoid NiiVue precision issues), pass that scale here
+ *   so we can convert the NiiVue mm-space FOV back to physical world
+ *   coordinates. Pass 1.0 if no normalization was applied.
+ * @returns Viewport bounds in world space, intersected with volumeBounds
+ */
+export function computeViewportBounds2D(nv, sliceType, volumeBounds, normalizationScale = 1.0) {
+    // Compute the base field of view from the FULL volume bounds (in normalized
+    // mm space), then swizzle to screen axes for this slice orientation.
+    //
+    // IMPORTANT: We intentionally do NOT use nv.screenFieldOfViewExtendedMM()
+    // because that returns the extents of the *currently loaded* NVImage (the
+    // slab). After a viewport-aware reload shrinks the slab, the next call to
+    // screenFieldOfViewExtendedMM() would return a smaller FOV, creating a
+    // feedback loop that progressively shrinks the slab to nothing.
+    //
+    // Instead, we derive the base FOV from the constant full-volume bounds,
+    // scaled to normalized mm space (matching the slab NVImage's affine).
+    // We then apply the same swizzle that NiiVue uses, giving us a stable
+    // base FOV that doesn't depend on the current slab geometry.
+    const normMin = [
+        volumeBounds.min[0] * normalizationScale,
+        volumeBounds.min[1] * normalizationScale,
+        volumeBounds.min[2] * normalizationScale,
+    ];
+    const normMax = [
+        volumeBounds.max[0] * normalizationScale,
+        volumeBounds.max[1] * normalizationScale,
+        volumeBounds.max[2] * normalizationScale,
+    ];
+    // Swizzle to screen axes (same mapping as NiiVue's swizzleVec3MM):
+    //   AXIAL:    screen X = mm X, screen Y = mm Y
+    //   CORONAL:  screen X = mm X, screen Y = mm Z
+    //   SAGITTAL: screen X = mm Y, screen Y = mm Z
+    let mnMM0, mxMM0, mnMM1, mxMM1;
+    switch (sliceType) {
+        case SLICE_TYPE.CORONAL:
+            mnMM0 = normMin[0];
+            mxMM0 = normMax[0]; // screen X = mm X
+            mnMM1 = normMin[2];
+            mxMM1 = normMax[2]; // screen Y = mm Z
+            break;
+        case SLICE_TYPE.SAGITTAL:
+            mnMM0 = normMin[1];
+            mxMM0 = normMax[1]; // screen X = mm Y
+            mnMM1 = normMin[2];
+            mxMM1 = normMax[2]; // screen Y = mm Z
+            break;
+        default: // AXIAL
+            mnMM0 = normMin[0];
+            mxMM0 = normMax[0]; // screen X = mm X
+            mnMM1 = normMin[1];
+            mxMM1 = normMax[1]; // screen Y = mm Y
+            break;
+    }
+    // Account for canvas aspect ratio stretching (matches draw2DMain logic)
+    // NiiVue stretches the FOV to fill the canvas while preserving aspect ratio
+    const canvas = nv.canvas;
+    if (canvas) {
+        const canvasW = canvas.width || 1;
+        const canvasH = canvas.height || 1;
+        const fovW = Math.abs(mxMM0 - mnMM0);
+        const fovH = Math.abs(mxMM1 - mnMM1);
+        if (fovW > 0 && fovH > 0) {
+            const canvasAspect = canvasW / canvasH;
+            const fovAspect = fovW / fovH;
+            if (canvasAspect > fovAspect) {
+                // Canvas is wider than FOV: expand X
+                const midX = (mnMM0 + mxMM0) * 0.5;
+                const newHalfW = (fovH * canvasAspect) * 0.5;
+                mnMM0 = midX - newHalfW;
+                mxMM0 = midX + newHalfW;
+            }
+            else {
+                // Canvas is taller than FOV: expand Y
+                const midY = (mnMM1 + mxMM1) * 0.5;
+                const newHalfH = (fovW / canvasAspect) * 0.5;
+                mnMM1 = midY - newHalfH;
+                mxMM1 = midY + newHalfH;
+            }
+        }
+    }
+    // Apply pan and zoom (matching NiiVue's draw2DMain logic)
+    const pan = nv.scene.pan2Dxyzmm; // vec4: [panX, panY, panZ, zoom]
+    // Swizzle the pan vector to match the current orientation
+    const panSwizzled = nv.swizzleVec3MM([pan[0], pan[1], pan[2]], sliceType);
+    const zoom = pan[3] || 1;
+    // Apply pan: shift visible window
+    mnMM0 -= panSwizzled[0];
+    mxMM0 -= panSwizzled[0];
+    mnMM1 -= panSwizzled[1];
+    mxMM1 -= panSwizzled[1];
+    // Apply zoom: divide by zoom factor (zoom > 1 = zoomed in = smaller FOV)
+    mnMM0 /= zoom;
+    mxMM0 /= zoom;
+    mnMM1 /= zoom;
+    mxMM1 /= zoom;
+    // Convert from NiiVue's mm space back to physical world coordinates.
+    // If the slab affine was normalized (multiplied by normalizationScale),
+    // NiiVue's mm values are world * normalizationScale. Dividing by the
+    // normalization scale recovers physical coordinates.
+    if (normalizationScale !== 1.0 && normalizationScale > 0) {
+        const invNorm = 1.0 / normalizationScale;
+        mnMM0 *= invNorm;
+        mxMM0 *= invNorm;
+        mnMM1 *= invNorm;
+        mxMM1 *= invNorm;
+    }
+    // Now un-swizzle back to RAS world coordinates.
+    // The swizzle maps depend on the slice orientation:
+    // AXIAL:    screen X = R/L (world X), screen Y = A/P (world Y)
+    // CORONAL:  screen X = R/L (world X), screen Y = S/I (world Z)
+    // SAGITTAL: screen X = A/P (world Y), screen Y = S/I (world Z)
+    //
+    // The orthogonal axis (depth) is left as full volume extent.
+    const result = {
+        min: [...volumeBounds.min],
+        max: [...volumeBounds.max],
+    };
+    // Ensure min < max for each swizzled axis
+    const visMin0 = Math.min(mnMM0, mxMM0);
+    const visMax0 = Math.max(mnMM0, mxMM0);
+    const visMin1 = Math.min(mnMM1, mxMM1);
+    const visMax1 = Math.max(mnMM1, mxMM1);
+    switch (sliceType) {
+        case SLICE_TYPE.AXIAL:
+            // screen X = world X (R/L), screen Y = world Y (A/P)
+            result.min[0] = visMin0;
+            result.max[0] = visMax0;
+            result.min[1] = visMin1;
+            result.max[1] = visMax1;
+            // Z (S/I) = full extent (orthogonal axis)
+            break;
+        case SLICE_TYPE.CORONAL:
+            // screen X = world X (R/L), screen Y = world Z (S/I)
+            result.min[0] = visMin0;
+            result.max[0] = visMax0;
+            result.min[2] = visMin1;
+            result.max[2] = visMax1;
+            // Y (A/P) = full extent (orthogonal axis)
+            break;
+        case SLICE_TYPE.SAGITTAL:
+            // screen X = world Y (A/P), screen Y = world Z (S/I)
+            result.min[1] = visMin0;
+            result.max[1] = visMax0;
+            result.min[2] = visMin1;
+            result.max[2] = visMax1;
+            // X (R/L) = full extent (orthogonal axis)
+            break;
+    }
+    return intersectBounds(result, volumeBounds);
+}
+/**
+ * Check if two VolumeBounds are approximately equal.
+ *
+ * @param a - First bounds
+ * @param b - Second bounds
+ * @param tolerance - Relative tolerance (default: 0.01 = 1%)
+ * @returns True if bounds are within tolerance
+ */
+export function boundsApproxEqual(a, b, tolerance = 0.01) {
+    for (let i = 0; i < 3; i++) {
+        const rangeA = a.max[i] - a.min[i];
+        const rangeB = b.max[i] - b.min[i];
+        const maxRange = Math.max(Math.abs(rangeA), Math.abs(rangeB), 1e-10);
+        if (Math.abs(a.min[i] - b.min[i]) / maxRange > tolerance)
+            return false;
+        if (Math.abs(a.max[i] - b.max[i]) / maxRange > tolerance)
+            return false;
+    }
+    return true;
+}
+//# sourceMappingURL=ViewportBounds.js.map

package/dist/ViewportBounds.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ViewportBounds.js","sourceRoot":"","sources":["../src/ViewportBounds.ts"],"names":[],"mappings":"AAAA,wDAAwD;AACxD,+BAA+B;AAG/B,OAAO,EAAE,UAAU,EAAE,MAAM,gBAAgB,CAAC;AAG5C;;;GAGG;AACH,MAAM,UAAU,eAAe,CAC7B,CAAe,EACf,CAAe;IAEf,MAAM,GAAG,GAA6B;QACpC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC5B,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC5B,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;KAC7B,CAAC;IACF,MAAM,GAAG,GAA6B;QACpC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC5B,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC5B,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;KAC7B,CAAC;IACF,mCAAmC;IACnC,OAAO;QACL,GAAG,EAAE;YACH,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;YACxB,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;YACxB,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;SACzB;QACD,GAAG,EAAE;YACH,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;YACxB,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;YACxB,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;SACzB;KACF,CAAC;AACJ,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,uBAAuB,CACrC,EAAU,EACV,YAA0B;IAE1B,uCAAuC;IACvC,MAAM,OAAO,GAAG,EAAE,CAAC,kBAAkB,CAAC,IAAI,CAAC,CAAC;IAC5C,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO;IAC9B,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO;IAC9B,MAAM,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO;IAEjC,0BAA0B;IAC1B,MAAM,MAAM,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IACrC,MAAM,MAAM,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IACrC,MAAM,MAAM,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAErC,oDAAoD;IACpD,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CACxB,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAChE,GAAG,GAAG,CAAC;IAER,2DAA2D;IAC3D,MAAM,KAAK,GAAG,CAAC,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,kBAAkB,IAAI,CAAC,CAAC,CAAC;IAEpE,sBAAsB;IACtB,MAAM,MAAM,GAAG,EAAE,CAAC,MAAM,CAAC;IACzB,MAAM,OAAO,GAAG,MAAM,EAAE,KAAK,IAAI,CAAC,CAAC;IACnC,MAAM,OAAO,GAAG,MAAM,EAAE,MAAM,IAAI,CAAC,CAAC;IACpC,MAAM,OAAO,GAAG,OAAO,GAAG,OAAO,CAAC;IAElC,gDAAgD;IAChD,IAAI,KAAa,EAAE,KAAa,CAAC;IACjC,IAAI,OAAO,GAAG,CAAC,EAAE,CAAC;QAChB,WAAW;QACX,KAAK,GAAG,KAAK,CAAC;QACd,KAAK,GAAG,KAAK,GAAG,OAAO,CAAC;IAC1B,CAAC;SAAM,CAAC;QACN,YAAY;QACZ,KAAK,GAAG,KAAK,GAAG,OAAO,CAAC;QACxB,KAAK,GAAG,KAAK,CAAC;IAChB,CAAC;IACD,wEAAwE;IACxE,6EAA6E;IAC7E,yEAAyE;IACzE,EAAE;IACF,+EAA+E;IAC/E,8EAA8E;IAC9E,4EAA4E;IAC5E,4EAA4E;IAC5E,uEAAuE;IAEvE,+CAA+C;IAC/C,uEAAuE;IACvE,wCAAwC;IACxC,0EAA0E;IAC1E,MAAM,OAAO,GAAG,EAAE,CAAC,KAAK,CAAC,aAAa,IAAI,CAAC,CAAC;IAC5C,MAAM,SAAS,GAAG,EAAE,CAAC,KAAK,CAAC,eAAe,IAAI,CAAC,CAAC;IAChD,MAAM,KAAK,GAAG,CAAC,CAAC,OAAO,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;IAChD,MAAM,KAAK,GAAG,CAAC,CAAC,GAAG,GAAG,SAAS,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;IAElD,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;IAC9B,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;IAC9B,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;IAC9B,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;IAE9B,mEAAmE;IACnE,oEAAoE;IACpE,EAAE;IACF,oDAAoD;IACpD,MAAM,YAAY,GAA6B;QAC7C,CAAC,KAAK,EAAE,KAAK;QACb,KAAK,GAAG,KAAK,EAAE,KAAK;QACpB,CAAC,KAAK,GAAG,KAAK,EAAE,KAAK;KACtB,CAAC;IACF,2CAA2C;IAC3C,MAAM,YAAY,GAA6B;QAC7C,KAAK,EAAE,KAAK;QACZ,KAAK,GAAG,KAAK,EAAE,KAAK;QACpB,CAAC,KAAK,GAAG,KAAK,EAAE,KAAK;KACtB,CAAC;IACF,2CAA2C;IAC3C,MAAM,YAAY,GAA6B;QAC7C,CAAC,EAAE,KAAK;QACR,KAAK,EAAE,KAAK;QACZ,KAAK,EAAE,KAAK;KACb,CAAC;IAEF,mDAAmD;IACnD,gEAAgE;IAChE,wEAAwE;IACxE,MAAM,eAAe,GAA6B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5D,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,IAAI,EAAE,EAAE,CAAC;QACpC,eAAe,CAAC,IAAI,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,IAAI,CAAC,CAAC;YAC1D,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,IAAI,CAAC,CAAC;YACpC,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC;IAC5C,CAAC;IAED,MAAM,aAAa,GAAiB;QAClC,GAAG,EAAE;YACH,MAAM,GAAG,eAAe,CAAC,CAAC,CAAC;YAC3B,MAAM,GAAG,eAAe,CAAC,CAAC,CAAC;YAC3B,MAAM,GAAG,eAAe,CAAC,CAAC,CAAC;SAC5B;QACD,GAAG,EAAE;YACH,MAAM,GAAG,eAAe,CAAC,CAAC,CAAC;YAC3B,MAAM,GAAG,eAAe,CAAC,CAAC,CAAC;YAC3B,MAAM,GAAG,eAAe,CAAC,CAAC,CAAC;SAC5B;KACF,CAAC;IAEF,OAAO,eAAe,CAAC,aAAa,EAAE,YAAY,CAAC,CAAC;AACtD,CAAC;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,MAAM,UAAU,uBAAuB,CACrC,EAAU,EACV,SAAqB,EACrB,YAA0B,EAC1B,qBAA6B,GAAG;IAEhC,4EAA4E;IAC5E,qEAAqE;IACrE,EAAE;IACF,0EAA0E;IAC1E,0EAA0E;IAC1E,0EAA0E;IAC1E,uEAAuE;IACvE,gEAAgE;IAChE,EAAE;IACF,wEAAwE;IACxE,sEAAsE;IACtE,sEAAsE;IACtE,6DAA6D;IAC7D,MAAM,OAAO,GAA6B;QACxC,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,kBAAkB;QACxC,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,kBAAkB;QACxC,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,kBAAkB;KACzC,CAAC;IACF,MAAM,OAAO,GAA6B;QACxC,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,kBAAkB;QACxC,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,kBAAkB;QACxC,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,kBAAkB;KACzC,CAAC;IAEF,mEAAmE;IACnE,+CAA+C;IAC/C,+CAA+C;IAC/C,+CAA+C;IAC/C,IAAI,KAAa,EAAE,KAAa,EAAE,KAAa,EAAE,KAAa,CAAC;IAC/D,QAAQ,SAAS,EAAE,CAAC;QAClB,KAAK,UAAU,CAAC,OAAO;YACrB,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;YACnB,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,kBAAkB;YACtC,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;YACnB,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,kBAAkB;YACtC,MAAM;QACR,KAAK,UAAU,CAAC,QAAQ;YACtB,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;YACnB,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,kBAAkB;YACtC,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;YACnB,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,kBAAkB;YACtC,MAAM;QACR,SAAS,QAAQ;YACf,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;YACnB,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,kBAAkB;YACtC,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;YACnB,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,kBAAkB;YACtC,MAAM;IACV,CAAC;IAED,wEAAwE;IACxE,4EAA4E;IAC5E,MAAM,MAAM,GAAG,EAAE,CAAC,MAAM,CAAC;IACzB,IAAI,MAAM,EAAE,CAAC;QACX,MAAM,OAAO,GAAG,MAAM,CAAC,KAAK,IAAI,CAAC,CAAC;QAClC,MAAM,OAAO,GAAG,MAAM,CAAC,MAAM,IAAI,CAAC,CAAC;QACnC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC;QACrC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC;QACrC,IAAI,IAAI,GAAG,CAAC,IAAI,IAAI,GAAG,CAAC,EAAE,CAAC;YACzB,MAAM,YAAY,GAAG,OAAO,GAAG,OAAO,CAAC;YACvC,MAAM,SAAS,GAAG,IAAI,GAAG,IAAI,CAAC;YAC9B,IAAI,YAAY,GAAG,SAAS,EAAE,CAAC;gBAC7B,qCAAqC;gBACrC,MAAM,IAAI,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,GAAG,CAAC;gBACnC,MAAM,QAAQ,GAAG,CAAC,IAAI,GAAG,YAAY,CAAC,GAAG,GAAG,CAAC;gBAC7C,KAAK,GAAG,IAAI,GAAG,QAAQ,CAAC;gBACxB,KAAK,GAAG,IAAI,GAAG,QAAQ,CAAC;YAC1B,CAAC;iBAAM,CAAC;gBACN,sCAAsC;gBACtC,MAAM,IAAI,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,GAAG,CAAC;gBACnC,MAAM,QAAQ,GAAG,CAAC,IAAI,GAAG,YAAY,CAAC,GAAG,GAAG,CAAC;gBAC7C,KAAK,GAAG,IAAI,GAAG,QAAQ,CAAC;gBACxB,KAAK,GAAG,IAAI,GAAG,QAAQ,CAAC;YAC1B,CAAC;QACH,CAAC;IACH,CAAC;IAED,0DAA0D;IAC1D,MAAM,GAAG,GAAG,EAAE,CAAC,KAAK,CAAC,UAAU,CAAC,CAAC,iCAAiC;IAClE,0DAA0D;IAC1D,MAAM,WAAW,GAAG,EAAE,CAAC,aAAa,CAClC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAwC,EAC/D,SAAS,CACV,CAAC;IACF,MAAM,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAEzB,kCAAkC;IAClC,KAAK,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC;IACxB,KAAK,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC;IACxB,KAAK,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC;IACxB,KAAK,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC;IAExB,yEAAyE;IACzE,KAAK,IAAI,IAAI,CAAC;IACd,KAAK,IAAI,IAAI,CAAC;IACd,KAAK,IAAI,IAAI,CAAC;IACd,KAAK,IAAI,IAAI,CAAC;IAEd,qEAAqE;IACrE,wEAAwE;IACxE,qEAAqE;IACrE,qDAAqD;IACrD,IAAI,kBAAkB,KAAK,GAAG,IAAI,kBAAkB,GAAG,CAAC,EAAE,CAAC;QACzD,MAAM,OAAO,GAAG,GAAG,GAAG,kBAAkB,CAAC;QACzC,KAAK,IAAI,OAAO,CAAC;QACjB,KAAK,IAAI,OAAO,CAAC;QACjB,KAAK,IAAI,OAAO,CAAC;QACjB,KAAK,IAAI,OAAO,CAAC;IACnB,CAAC;IAED,gDAAgD;IAChD,oDAAoD;IACpD,+DAA+D;IAC/D,+DAA+D;IAC/D,+DAA+D;IAC/D,EAAE;IACF,6DAA6D;IAC7D,MAAM,MAAM,GAAiB;QAC3B,GAAG,EAAE,CAAC,GAAG,YAAY,CAAC,GAAG,CAAC;QAC1B,GAAG,EAAE,CAAC,GAAG,YAAY,CAAC,GAAG,CAAC;KAC3B,CAAC;IAEF,0CAA0C;IAC1C,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;IACvC,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;IACvC,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;IACvC,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;IAEvC,QAAQ,SAAS,EAAE,CAAC;QAClB,KAAK,UAAU,CAAC,KAAK;YACnB,qDAAqD;YACrD,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YACxB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YACxB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YACxB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YACxB,0CAA0C;YAC1C,MAAM;QACR,KAAK,UAAU,CAAC,OAAO;YACrB,qDAAqD;YACrD,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YACxB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YACxB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YACxB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YACxB,0CAA0C;YAC1C,MAAM;QACR,KAAK,UAAU,CAAC,QAAQ;YACtB,qDAAqD;YACrD,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YACxB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YACxB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YACxB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YACxB,0CAA0C;YAC1C,MAAM;IACV,CAAC;IAED,OAAO,eAAe,CAAC,MAAM,EAAE,YAAY,CAAC,CAAC;AAC/C,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,iBAAiB,CAC/B,CAAe,EACf,CAAe,EACf,YAAoB,IAAI;IAExB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACnC,MAAM,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACnC,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,KAAK,CAAC,CAAC;QAErE,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,SAAS;YAAE,OAAO,KAAK,CAAC;QACvE,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,SAAS;YAAE,OAAO,KAAK,CAAC;IACzE,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC"}