@fideus-labs/ngff-zarr 0.2.0 → 0.2.1

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@fideus-labs/ngff-zarr",
-  "version": "0.2.0",
+  "version": "0.2.1",
   "description": "TypeScript implementation of ngff-zarr for reading and writing OME-Zarr files",
   "keywords": [
     "ome-zarr",
@@ -26,11 +26,35 @@
   "types": "./esm/mod.d.ts",
   "exports": {
     ".": {
+      "types": "./esm/mod.d.ts",
+      "browser": "./esm/browser-mod.js",
       "import": "./esm/mod.js",
       "require": "./script/mod.js",
-      "types": "./esm/mod.d.ts"
+      "default": "./esm/mod.js"
+    },
+    "./browser": {
+      "types": "./types/browser-mod.d.ts",
+      "import": "./esm/browser-mod.js",
+      "require": "./script/browser-mod.js",
+      "default": "./esm/browser-mod.js"
+    },
+    "./methods/itkwasm-browser": {
+      "types": "./types/methods/itkwasm-browser.d.ts",
+      "import": "./esm/methods/itkwasm-browser.js",
+      "require": "./script/methods/itkwasm-browser.js",
+      "default": "./esm/methods/itkwasm-browser.js"
+    },
+    "./methods/itkwasm-node": {
+      "types": "./types/methods/itkwasm-node.d.ts",
+      "import": "./esm/methods/itkwasm-node.js",
+      "require": "./script/methods/itkwasm-node.js",
+      "default": "./esm/methods/itkwasm-node.js"
     }
   },
+  "browser": {
+    "./esm/methods/itkwasm-node.js": "./esm/methods/itkwasm-browser.js",
+    "./script/methods/itkwasm-node.js": "./script/methods/itkwasm-browser.js"
+  },
   "files": [
     "esm/",
     "script/",
@@ -40,7 +64,7 @@
   ],
   "dependencies": {
     "@itk-wasm/downsample": "^1.8.1",
-    "itk-wasm": "^1.0.0-b.195",
+    "itk-wasm": "^1.0.0-b.196",
     "p-queue": "^8.1.0",
     "@zarrita/storage": "^0.1.1",
     "zod": "^4.0.2",

package/script/browser-mod.d.ts

@@ -0,0 +1,14 @@
+export * from "./types/units.js";
+export * from "./types/methods.js";
+export * from "./types/array_interface.js";
+export * from "./types/zarr_metadata.js";
+export * from "./types/ngff_image.js";
+export * from "./types/multiscales.js";
+export * from "./schemas/units.js";
+export * from "./schemas/methods.js";
+export * from "./schemas/zarr_metadata.js";
+export * from "./schemas/ngff_image.js";
+export * from "./schemas/multiscales.js";
+export { isValidDimension, isValidUnit, validateMetadata, } from "./utils/validation.js";
+export { createAxis, createDataset, createMetadata, createMultiscales, createNgffImage, } from "./utils/factory.js";
+//# sourceMappingURL=browser-mod.d.ts.map

package/script/browser-mod.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"browser-mod.d.ts","sourceRoot":"","sources":["../src/browser-mod.ts"],"names":[],"mappings":"AASA,cAAc,kBAAkB,CAAC;AACjC,cAAc,oBAAoB,CAAC;AACnC,cAAc,4BAA4B,CAAC;AAC3C,cAAc,0BAA0B,CAAC;AACzC,cAAc,uBAAuB,CAAC;AACtC,cAAc,wBAAwB,CAAC;AAEvC,cAAc,oBAAoB,CAAC;AACnC,cAAc,sBAAsB,CAAC;AACrC,cAAc,4BAA4B,CAAC;AAC3C,cAAc,yBAAyB,CAAC;AACxC,cAAc,0BAA0B,CAAC;AAEzC,OAAO,EACL,gBAAgB,EAChB,WAAW,EACX,gBAAgB,GACjB,MAAM,uBAAuB,CAAC;AAC/B,OAAO,EACL,UAAU,EACV,aAAa,EACb,cAAc,EACd,iBAAiB,EACjB,eAAe,GAChB,MAAM,oBAAoB,CAAC"}

package/script/browser-mod.js

@@ -0,0 +1,48 @@
+"use strict";
+var __createBinding = (this && this.__createBinding) || (Object.create ? (function(o, m, k, k2) {
+    if (k2 === undefined) k2 = k;
+    var desc = Object.getOwnPropertyDescriptor(m, k);
+    if (!desc || ("get" in desc ? !m.__esModule : desc.writable || desc.configurable)) {
+      desc = { enumerable: true, get: function() { return m[k]; } };
+    }
+    Object.defineProperty(o, k2, desc);
+}) : (function(o, m, k, k2) {
+    if (k2 === undefined) k2 = k;
+    o[k2] = m[k];
+}));
+var __exportStar = (this && this.__exportStar) || function(m, exports) {
+    for (var p in m) if (p !== "default" && !Object.prototype.hasOwnProperty.call(exports, p)) __createBinding(exports, m, p);
+};
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.createNgffImage = exports.createMultiscales = exports.createMetadata = exports.createDataset = exports.createAxis = exports.validateMetadata = exports.isValidUnit = exports.isValidDimension = void 0;
+// SPDX-FileCopyrightText: Copyright (c) Fideus Labs LLC
+// SPDX-License-Identifier: MIT
+// Browser-compatible module exports
+// This module excludes I/O functionality (from_ngff_zarr, to_ngff_zarr)
+// because those modules depend on Node.js/Deno-specific filesystem APIs
+// that are not available in browser environments.
+//
+// For browser use cases, the schemas, types, and validation utilities
+// are the most commonly needed functionality.
+__exportStar(require("./types/units.js"), exports);
+__exportStar(require("./types/methods.js"), exports);
+__exportStar(require("./types/array_interface.js"), exports);
+__exportStar(require("./types/zarr_metadata.js"), exports);
+__exportStar(require("./types/ngff_image.js"), exports);
+__exportStar(require("./types/multiscales.js"), exports);
+__exportStar(require("./schemas/units.js"), exports);
+__exportStar(require("./schemas/methods.js"), exports);
+__exportStar(require("./schemas/zarr_metadata.js"), exports);
+__exportStar(require("./schemas/ngff_image.js"), exports);
+__exportStar(require("./schemas/multiscales.js"), exports);
+var validation_js_1 = require("./utils/validation.js");
+Object.defineProperty(exports, "isValidDimension", { enumerable: true, get: function () { return validation_js_1.isValidDimension; } });
+Object.defineProperty(exports, "isValidUnit", { enumerable: true, get: function () { return validation_js_1.isValidUnit; } });
+Object.defineProperty(exports, "validateMetadata", { enumerable: true, get: function () { return validation_js_1.validateMetadata; } });
+var factory_js_1 = require("./utils/factory.js");
+Object.defineProperty(exports, "createAxis", { enumerable: true, get: function () { return factory_js_1.createAxis; } });
+Object.defineProperty(exports, "createDataset", { enumerable: true, get: function () { return factory_js_1.createDataset; } });
+Object.defineProperty(exports, "createMetadata", { enumerable: true, get: function () { return factory_js_1.createMetadata; } });
+Object.defineProperty(exports, "createMultiscales", { enumerable: true, get: function () { return factory_js_1.createMultiscales; } });
+Object.defineProperty(exports, "createNgffImage", { enumerable: true, get: function () { return factory_js_1.createNgffImage; } });
+// Note: Excluding I/O modules for browser compatibility

package/script/methods/itkwasm-browser.d.ts

@@ -0,0 +1,6 @@
+import { NgffImage } from "../types/ngff_image.js";
+/**
+ * Main downsampling function for ITK-Wasm (browser version)
+ */
+export declare function downsampleItkWasm(ngffImage: NgffImage, scaleFactors: (Record<string, number> | number)[], smoothing: "gaussian" | "bin_shrink" | "label_image"): Promise<NgffImage[]>;
+//# sourceMappingURL=itkwasm-browser.d.ts.map

package/script/methods/itkwasm-browser.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"itkwasm-browser.d.ts","sourceRoot":"","sources":["../../src/methods/itkwasm-browser.ts"],"names":[],"mappings":"AAcA,OAAO,EAAE,SAAS,EAAE,MAAM,wBAAwB,CAAC;AA2fnD;;GAEG;AACH,wBAAsB,iBAAiB,CACrC,SAAS,EAAE,SAAS,EACpB,YAAY,EAAE,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,GAAG,MAAM,CAAC,EAAE,EACjD,SAAS,EAAE,UAAU,GAAG,YAAY,GAAG,aAAa,GACnD,OAAO,CAAC,SAAS,EAAE,CAAC,CAsHtB"}

package/script/methods/itkwasm-browser.js

@@ -0,0 +1,488 @@
+"use strict";
+// SPDX-FileCopyrightText: Copyright (c) Fideus Labs LLC
+// SPDX-License-Identifier: MIT
+var __createBinding = (this && this.__createBinding) || (Object.create ? (function(o, m, k, k2) {
+    if (k2 === undefined) k2 = k;
+    var desc = Object.getOwnPropertyDescriptor(m, k);
+    if (!desc || ("get" in desc ? !m.__esModule : desc.writable || desc.configurable)) {
+      desc = { enumerable: true, get: function() { return m[k]; } };
+    }
+    Object.defineProperty(o, k2, desc);
+}) : (function(o, m, k, k2) {
+    if (k2 === undefined) k2 = k;
+    o[k2] = m[k];
+}));
+var __setModuleDefault = (this && this.__setModuleDefault) || (Object.create ? (function(o, v) {
+    Object.defineProperty(o, "default", { enumerable: true, value: v });
+}) : function(o, v) {
+    o["default"] = v;
+});
+var __importStar = (this && this.__importStar) || function (mod) {
+    if (mod && mod.__esModule) return mod;
+    var result = {};
+    if (mod != null) for (var k in mod) if (k !== "default" && Object.prototype.hasOwnProperty.call(mod, k)) __createBinding(result, mod, k);
+    __setModuleDefault(result, mod);
+    return result;
+};
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.downsampleItkWasm = downsampleItkWasm;
+/**
+ * Browser-compatible ITK-Wasm downsampling support
+ * Uses WebWorker-based implementations from @itk-wasm/downsample
+ */
+const downsample_1 = require("@itk-wasm/downsample");
+const zarr = __importStar(require("zarrita"));
+const ngff_image_js_1 = require("../types/ngff_image.js");
+const itkwasm_shared_js_1 = require("./itkwasm-shared.js");
+/**
+ * Perform Gaussian downsampling using ITK-Wasm (browser version)
+ */
+async function downsampleGaussian(image, dimFactors, spatialDims) {
+    // Handle time dimension by processing each time slice independently
+    if (image.dims.includes("t")) {
+        const tDimIndex = image.dims.indexOf("t");
+        const tSize = image.data.shape[tDimIndex];
+        const newDims = image.dims.filter((dim) => dim !== "t");
+        // Downsample each time slice
+        const downsampledSlices = [];
+        for (let t = 0; t < tSize; t++) {
+            // Extract time slice
+            const selection = new Array(image.data.shape.length).fill(null);
+            selection[tDimIndex] = t;
+            const sliceData = await zarr.get(image.data, selection);
+            // Create temporary zarr array for this slice
+            const sliceStore = new Map();
+            const sliceRoot = zarr.root(sliceStore);
+            const sliceShape = image.data.shape.filter((_, i) => i !== tDimIndex);
+            const sliceChunkShape = sliceShape.map((s) => Math.min(s, 256));
+            const sliceArray = await zarr.create(sliceRoot.resolve("slice"), {
+                shape: sliceShape,
+                chunk_shape: sliceChunkShape,
+                data_type: image.data.dtype,
+                fill_value: 0,
+            });
+            const fullSelection = new Array(sliceShape.length).fill(null);
+            await zarr.set(sliceArray, fullSelection, sliceData);
+            // Create NgffImage for this slice (without 't' dimension)
+            const sliceImage = new ngff_image_js_1.NgffImage({
+                data: sliceArray,
+                dims: newDims,
+                scale: Object.fromEntries(Object.entries(image.scale).filter(([dim]) => dim !== "t")),
+                translation: Object.fromEntries(Object.entries(image.translation).filter(([dim]) => dim !== "t")),
+                name: image.name,
+                axesUnits: image.axesUnits
+                    ? Object.fromEntries(Object.entries(image.axesUnits).filter(([dim]) => dim !== "t"))
+                    : undefined,
+                computedCallbacks: image.computedCallbacks,
+            });
+            // Recursively downsample this slice (without 't', so no infinite loop)
+            const downsampledSlice = await downsampleGaussian(sliceImage, dimFactors, spatialDims);
+            downsampledSlices.push(downsampledSlice.data);
+        }
+        // Combine downsampled slices back into a single array with 't' dimension
+        const firstSlice = downsampledSlices[0];
+        const combinedShape = [...image.data.shape];
+        combinedShape[tDimIndex] = tSize;
+        // Update spatial dimensions based on downsampled size
+        for (let i = 0; i < image.dims.length; i++) {
+            if (i !== tDimIndex) {
+                const sliceIndex = i < tDimIndex ? i : i - 1;
+                combinedShape[i] = firstSlice.shape[sliceIndex];
+            }
+        }
+        // Create combined array
+        const combinedStore = new Map();
+        const combinedRoot = zarr.root(combinedStore);
+        const combinedArray = await zarr.create(combinedRoot.resolve("combined"), {
+            shape: combinedShape,
+            chunk_shape: combinedShape.map((s) => Math.min(s, 256)),
+            data_type: image.data.dtype,
+            fill_value: 0,
+        });
+        // Copy each downsampled slice into the combined array
+        for (let t = 0; t < tSize; t++) {
+            const sliceData = await zarr.get(downsampledSlices[t]);
+            const targetSelection = new Array(combinedShape.length).fill(null);
+            targetSelection[tDimIndex] = t;
+            await zarr.set(combinedArray, targetSelection, sliceData);
+        }
+        // Compute new metadata (time dimension unchanged, spatial dimensions downsampled)
+        const [translation, scale] = (0, itkwasm_shared_js_1.nextScaleMetadata)(image, dimFactors, spatialDims);
+        return new ngff_image_js_1.NgffImage({
+            data: combinedArray,
+            dims: image.dims,
+            scale: { ...image.scale, ...scale },
+            translation: { ...image.translation, ...translation },
+            name: image.name,
+            axesUnits: image.axesUnits,
+            computedCallbacks: image.computedCallbacks,
+        });
+    }
+    const isVector = image.dims.includes("c");
+    // Convert to ITK-Wasm format
+    const itkImage = await (0, itkwasm_shared_js_1.zarrToItkImage)(image.data, image.dims, isVector);
+    // Prepare shrink factors - need to be for ALL dimensions in ITK order (reversed)
+    const shrinkFactors = [];
+    for (let i = image.dims.length - 1; i >= 0; i--) {
+        const dim = image.dims[i];
+        if (itkwasm_shared_js_1.SPATIAL_DIMS.includes(dim)) {
+            shrinkFactors.push(dimFactors[dim] || 1);
+        }
+    }
+    // Use all zeros for cropRadius
+    const cropRadius = new Array(shrinkFactors.length).fill(0);
+    // Perform downsampling using browser-compatible function
+    const { downsampled } = await (0, downsample_1.downsample)(itkImage, {
+        shrinkFactors,
+        cropRadius: cropRadius,
+    });
+    // Compute new metadata
+    const [translation, scale] = (0, itkwasm_shared_js_1.nextScaleMetadata)(image, dimFactors, spatialDims);
+    // Convert back to zarr array in a new in-memory store
+    const store = new Map();
+    const chunkShape = downsampled.size.map((s) => Math.min(s, 256)).reverse();
+    const array = await (0, itkwasm_shared_js_1.itkImageToZarr)(downsampled, store, "image", chunkShape, image.dims);
+    return new ngff_image_js_1.NgffImage({
+        data: array,
+        dims: image.dims,
+        scale,
+        translation,
+        name: image.name,
+        axesUnits: image.axesUnits,
+        computedCallbacks: image.computedCallbacks,
+    });
+}
+/**
+ * Perform bin shrink downsampling using ITK-Wasm (browser version)
+ */
+async function downsampleBinShrinkImpl(image, dimFactors, spatialDims) {
+    // Handle time dimension by processing each time slice independently
+    if (image.dims.includes("t")) {
+        const tDimIndex = image.dims.indexOf("t");
+        const tSize = image.data.shape[tDimIndex];
+        const newDims = image.dims.filter((dim) => dim !== "t");
+        // Downsample each time slice
+        const downsampledSlices = [];
+        for (let t = 0; t < tSize; t++) {
+            // Extract time slice
+            const selection = new Array(image.data.shape.length).fill(null);
+            selection[tDimIndex] = t;
+            const sliceData = await zarr.get(image.data, selection);
+            // Create temporary zarr array for this slice
+            const sliceStore = new Map();
+            const sliceRoot = zarr.root(sliceStore);
+            const sliceShape = image.data.shape.filter((_, i) => i !== tDimIndex);
+            const sliceChunkShape = sliceShape.map((s) => Math.min(s, 256));
+            const sliceArray = await zarr.create(sliceRoot.resolve("slice"), {
+                shape: sliceShape,
+                chunk_shape: sliceChunkShape,
+                data_type: image.data.dtype,
+                fill_value: 0,
+            });
+            const fullSelection = new Array(sliceShape.length).fill(null);
+            await zarr.set(sliceArray, fullSelection, sliceData);
+            // Create NgffImage for this slice (without 't' dimension)
+            const sliceImage = new ngff_image_js_1.NgffImage({
+                data: sliceArray,
+                dims: newDims,
+                scale: Object.fromEntries(Object.entries(image.scale).filter(([dim]) => dim !== "t")),
+                translation: Object.fromEntries(Object.entries(image.translation).filter(([dim]) => dim !== "t")),
+                name: image.name,
+                axesUnits: image.axesUnits
+                    ? Object.fromEntries(Object.entries(image.axesUnits).filter(([dim]) => dim !== "t"))
+                    : undefined,
+                computedCallbacks: image.computedCallbacks,
+            });
+            // Recursively downsample this slice
+            const downsampledSlice = await downsampleBinShrinkImpl(sliceImage, dimFactors, spatialDims);
+            downsampledSlices.push(downsampledSlice.data);
+        }
+        // Combine downsampled slices back into a single array with 't' dimension
+        const firstSlice = downsampledSlices[0];
+        const combinedShape = [...image.data.shape];
+        combinedShape[tDimIndex] = tSize;
+        // Update spatial dimensions based on downsampled size
+        for (let i = 0; i < image.dims.length; i++) {
+            if (i !== tDimIndex) {
+                const sliceIndex = i < tDimIndex ? i : i - 1;
+                combinedShape[i] = firstSlice.shape[sliceIndex];
+            }
+        }
+        // Create combined array
+        const combinedStore = new Map();
+        const combinedRoot = zarr.root(combinedStore);
+        const combinedArray = await zarr.create(combinedRoot.resolve("combined"), {
+            shape: combinedShape,
+            chunk_shape: combinedShape.map((s) => Math.min(s, 256)),
+            data_type: image.data.dtype,
+            fill_value: 0,
+        });
+        // Copy each downsampled slice into the combined array
+        for (let t = 0; t < tSize; t++) {
+            const sliceData = await zarr.get(downsampledSlices[t]);
+            const targetSelection = new Array(combinedShape.length).fill(null);
+            targetSelection[tDimIndex] = t;
+            await zarr.set(combinedArray, targetSelection, sliceData);
+        }
+        // Compute new metadata
+        const [translation, scale] = (0, itkwasm_shared_js_1.nextScaleMetadata)(image, dimFactors, spatialDims);
+        return new ngff_image_js_1.NgffImage({
+            data: combinedArray,
+            dims: image.dims,
+            scale: { ...image.scale, ...scale },
+            translation: { ...image.translation, ...translation },
+            name: image.name,
+            axesUnits: image.axesUnits,
+            computedCallbacks: image.computedCallbacks,
+        });
+    }
+    const isVector = image.dims.includes("c");
+    // Convert to ITK-Wasm format
+    const itkImage = await (0, itkwasm_shared_js_1.zarrToItkImage)(image.data, image.dims, isVector);
+    // Prepare shrink factors - only for spatial dimensions in ITK order (reversed)
+    const shrinkFactors = [];
+    for (let i = image.dims.length - 1; i >= 0; i--) {
+        const dim = image.dims[i];
+        if (itkwasm_shared_js_1.SPATIAL_DIMS.includes(dim)) {
+            shrinkFactors.push(dimFactors[dim] || 1);
+        }
+    }
+    // Perform downsampling using browser-compatible function
+    const { downsampled } = await (0, downsample_1.downsampleBinShrink)(itkImage, {
+        shrinkFactors,
+    });
+    // Compute new metadata
+    const [translation, scale] = (0, itkwasm_shared_js_1.nextScaleMetadata)(image, dimFactors, spatialDims);
+    // Convert back to zarr array in a new in-memory store
+    const store = new Map();
+    const chunkShape = downsampled.size.map((s) => Math.min(s, 256)).reverse();
+    const array = await (0, itkwasm_shared_js_1.itkImageToZarr)(downsampled, store, "image", chunkShape, image.dims);
+    return new ngff_image_js_1.NgffImage({
+        data: array,
+        dims: image.dims,
+        scale,
+        translation,
+        name: image.name,
+        axesUnits: image.axesUnits,
+        computedCallbacks: image.computedCallbacks,
+    });
+}
+/**
+ * Perform label image downsampling using ITK-Wasm (browser version)
+ */
+async function downsampleLabelImageImpl(image, dimFactors, spatialDims) {
+    // Handle time dimension by processing each time slice independently
+    if (image.dims.includes("t")) {
+        const tDimIndex = image.dims.indexOf("t");
+        const tSize = image.data.shape[tDimIndex];
+        const newDims = image.dims.filter((dim) => dim !== "t");
+        // Downsample each time slice
+        const downsampledSlices = [];
+        for (let t = 0; t < tSize; t++) {
+            // Extract time slice
+            const selection = new Array(image.data.shape.length).fill(null);
+            selection[tDimIndex] = t;
+            const sliceData = await zarr.get(image.data, selection);
+            // Create temporary zarr array for this slice
+            const sliceStore = new Map();
+            const sliceRoot = zarr.root(sliceStore);
+            const sliceShape = image.data.shape.filter((_, i) => i !== tDimIndex);
+            const sliceChunkShape = sliceShape.map((s) => Math.min(s, 256));
+            const sliceArray = await zarr.create(sliceRoot.resolve("slice"), {
+                shape: sliceShape,
+                chunk_shape: sliceChunkShape,
+                data_type: image.data.dtype,
+                fill_value: 0,
+            });
+            const fullSelection = new Array(sliceShape.length).fill(null);
+            await zarr.set(sliceArray, fullSelection, sliceData);
+            // Create NgffImage for this slice (without 't' dimension)
+            const sliceImage = new ngff_image_js_1.NgffImage({
+                data: sliceArray,
+                dims: newDims,
+                scale: Object.fromEntries(Object.entries(image.scale).filter(([dim]) => dim !== "t")),
+                translation: Object.fromEntries(Object.entries(image.translation).filter(([dim]) => dim !== "t")),
+                name: image.name,
+                axesUnits: image.axesUnits
+                    ? Object.fromEntries(Object.entries(image.axesUnits).filter(([dim]) => dim !== "t"))
+                    : undefined,
+                computedCallbacks: image.computedCallbacks,
+            });
+            // Recursively downsample this slice
+            const downsampledSlice = await downsampleLabelImageImpl(sliceImage, dimFactors, spatialDims);
+            downsampledSlices.push(downsampledSlice.data);
+        }
+        // Combine downsampled slices back into a single array with 't' dimension
+        const firstSlice = downsampledSlices[0];
+        const combinedShape = [...image.data.shape];
+        combinedShape[tDimIndex] = tSize;
+        // Update spatial dimensions based on downsampled size
+        for (let i = 0; i < image.dims.length; i++) {
+            if (i !== tDimIndex) {
+                const sliceIndex = i < tDimIndex ? i : i - 1;
+                combinedShape[i] = firstSlice.shape[sliceIndex];
+            }
+        }
+        // Create combined array
+        const combinedStore = new Map();
+        const combinedRoot = zarr.root(combinedStore);
+        const combinedArray = await zarr.create(combinedRoot.resolve("combined"), {
+            shape: combinedShape,
+            chunk_shape: combinedShape.map((s) => Math.min(s, 256)),
+            data_type: image.data.dtype,
+            fill_value: 0,
+        });
+        // Copy each downsampled slice into the combined array
+        for (let t = 0; t < tSize; t++) {
+            const sliceData = await zarr.get(downsampledSlices[t]);
+            const targetSelection = new Array(combinedShape.length).fill(null);
+            targetSelection[tDimIndex] = t;
+            await zarr.set(combinedArray, targetSelection, sliceData);
+        }
+        // Compute new metadata
+        const [translation, scale] = (0, itkwasm_shared_js_1.nextScaleMetadata)(image, dimFactors, spatialDims);
+        return new ngff_image_js_1.NgffImage({
+            data: combinedArray,
+            dims: image.dims,
+            scale: { ...image.scale, ...scale },
+            translation: { ...image.translation, ...translation },
+            name: image.name,
+            axesUnits: image.axesUnits,
+            computedCallbacks: image.computedCallbacks,
+        });
+    }
+    const isVector = image.dims.includes("c");
+    // Convert to ITK-Wasm format
+    const itkImage = await (0, itkwasm_shared_js_1.zarrToItkImage)(image.data, image.dims, isVector);
+    // Prepare shrink factors - need to be for ALL dimensions in ITK order (reversed)
+    const shrinkFactors = [];
+    for (let i = image.dims.length - 1; i >= 0; i--) {
+        const dim = image.dims[i];
+        if (itkwasm_shared_js_1.SPATIAL_DIMS.includes(dim)) {
+            shrinkFactors.push(dimFactors[dim] || 1);
+        }
+        else {
+            shrinkFactors.push(1); // Non-spatial dimensions don't shrink
+        }
+    }
+    // Use all zeros for cropRadius
+    const cropRadius = new Array(shrinkFactors.length).fill(0);
+    // Perform downsampling using browser-compatible function
+    const { downsampled } = await (0, downsample_1.downsampleLabelImage)(itkImage, {
+        shrinkFactors,
+        cropRadius: cropRadius,
+    });
+    // Compute new metadata
+    const [translation, scale] = (0, itkwasm_shared_js_1.nextScaleMetadata)(image, dimFactors, spatialDims);
+    // Convert back to zarr array in a new in-memory store
+    const store = new Map();
+    const chunkShape = downsampled.size.map((s) => Math.min(s, 256)).reverse();
+    const array = await (0, itkwasm_shared_js_1.itkImageToZarr)(downsampled, store, "image", chunkShape, image.dims);
+    return new ngff_image_js_1.NgffImage({
+        data: array,
+        dims: image.dims,
+        scale,
+        translation,
+        name: image.name,
+        axesUnits: image.axesUnits,
+        computedCallbacks: image.computedCallbacks,
+    });
+}
+/**
+ * Main downsampling function for ITK-Wasm (browser version)
+ */
+async function downsampleItkWasm(ngffImage, scaleFactors, smoothing) {
+    const multiscales = [ngffImage];
+    const dims = ngffImage.dims;
+    const spatialDims = dims.filter((dim) => itkwasm_shared_js_1.SPATIAL_DIMS.includes(dim));
+    let previousImage = ngffImage;
+    let previousDimFactors = {};
+    for (const dim of dims)
+        previousDimFactors[dim] = 1;
+    for (let i = 0; i < scaleFactors.length; i++) {
+        const scaleFactor = scaleFactors[i];
+        let sourceImage;
+        let sourceDimFactors;
+        if (smoothing === "bin_shrink") {
+            // Purely incremental: scaleFactor is the shrink for this step
+            sourceImage = previousImage;
+            sourceDimFactors = {};
+            if (typeof scaleFactor === "number") {
+                for (const dim of spatialDims)
+                    sourceDimFactors[dim] = scaleFactor;
+            }
+            else {
+                for (const dim of spatialDims) {
+                    sourceDimFactors[dim] = scaleFactor[dim] || 1;
+                }
+            }
+            for (const dim of dims) {
+                if (!(dim in sourceDimFactors))
+                    sourceDimFactors[dim] = 1;
+            }
+        }
+        else {
+            // Hybrid absolute strategy
+            const dimFactors = (0, itkwasm_shared_js_1.dimScaleFactors)(dims, scaleFactor, previousDimFactors, ngffImage, previousImage);
+            let canDownsampleIncrementally = true;
+            for (const dim of Object.keys(dimFactors)) {
+                const dimIndex = ngffImage.dims.indexOf(dim);
+                if (dimIndex >= 0) {
+                    const originalSize = ngffImage.data.shape[dimIndex];
+                    const targetSize = Math.floor(originalSize /
+                        (typeof scaleFactor === "number"
+                            ? scaleFactor
+                            : scaleFactor[dim]));
+                    const prevDimIndex = previousImage.dims.indexOf(dim);
+                    const previousSize = previousImage.data.shape[prevDimIndex];
+                    if (Math.floor(previousSize / dimFactors[dim]) !== targetSize) {
+                        canDownsampleIncrementally = false;
+                        break;
+                    }
+                }
+            }
+            if (canDownsampleIncrementally) {
+                sourceImage = previousImage;
+                sourceDimFactors = dimFactors;
+            }
+            else {
+                sourceImage = ngffImage;
+                const originalDimFactors = {};
+                for (const dim of dims)
+                    originalDimFactors[dim] = 1;
+                sourceDimFactors = (0, itkwasm_shared_js_1.dimScaleFactors)(dims, scaleFactor, originalDimFactors);
+            }
+        }
+        let downsampled;
+        if (smoothing === "gaussian") {
+            downsampled = await downsampleGaussian(sourceImage, sourceDimFactors, spatialDims);
+        }
+        else if (smoothing === "bin_shrink") {
+            downsampled = await downsampleBinShrinkImpl(sourceImage, sourceDimFactors, spatialDims);
+        }
+        else if (smoothing === "label_image") {
+            downsampled = await downsampleLabelImageImpl(sourceImage, sourceDimFactors, spatialDims);
+        }
+        else {
+            throw new Error(`Unknown smoothing method: ${smoothing}`);
+        }
+        multiscales.push(downsampled);
+        previousImage = downsampled;
+        if (smoothing === "bin_shrink") {
+            if (typeof scaleFactor === "number") {
+                for (const dim of spatialDims) {
+                    previousDimFactors[dim] *= scaleFactor;
+                }
+            }
+            else {
+                for (const dim of spatialDims) {
+                    previousDimFactors[dim] *= scaleFactor[dim] || 1;
+                }
+            }
+        }
+        else {
+            previousDimFactors = (0, itkwasm_shared_js_1.updatePreviousDimFactors)(scaleFactor, spatialDims, previousDimFactors);
+        }
+    }
+    return multiscales;
+}

package/script/methods/itkwasm-node.d.ts

@@ -0,0 +1,6 @@
+import { NgffImage } from "../types/ngff_image.js";
+/**
+ * Main downsampling function for ITK-Wasm (browser version)
+ */
+export declare function downsampleItkWasm(ngffImage: NgffImage, scaleFactors: (Record<string, number> | number)[], smoothing: "gaussian" | "bin_shrink" | "label_image"): Promise<NgffImage[]>;
+//# sourceMappingURL=itkwasm-node.d.ts.map

package/script/methods/itkwasm-node.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"itkwasm-node.d.ts","sourceRoot":"","sources":["../../src/methods/itkwasm-node.ts"],"names":[],"mappings":"AAcA,OAAO,EAAE,SAAS,EAAE,MAAM,wBAAwB,CAAC;AA2fnD;;GAEG;AACH,wBAAsB,iBAAiB,CACrC,SAAS,EAAE,SAAS,EACpB,YAAY,EAAE,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,GAAG,MAAM,CAAC,EAAE,EACjD,SAAS,EAAE,UAAU,GAAG,YAAY,GAAG,aAAa,GACnD,OAAO,CAAC,SAAS,EAAE,CAAC,CAsHtB"}