new script

Author: nlsschim

scripts/split_nd2_videos.py

@@ -3,15 +3,16 @@
 import skimage.io as sio
 from pathlib import Path
 
-def process_nd2_to_quadrants(nd2_file):
-    """Convert ND2 to TIF and split into 4 quadrants (512x512 each)"""
+
+def process_nd2_to_sixteenths(nd2_file):
+    """Convert ND2 to TIF and split into 16 tiles (256x256 each)"""
     # Read ND2 file
     with ND2Reader(nd2_file) as images:
         img = np.array(images)
 
     # Pad to 1024x1024 if needed
     oshape = img.shape
-    padded = np.zeros((oshape[0], 1024, 1024), dtype=img.dtype)
+    padded = np.zeros((oshape[0], 2044, 2048), dtype=img.dtype)
     padded[:, :oshape[1], :oshape[2]] = img
 
     # Get original file's directory and name
@@ -23,22 +24,22 @@ def process_nd2_to_quadrants(nd2_file):
     full_tif = output_dir / f"{output_name}.tif"
     sio.imsave(str(full_tif), padded)
 
-    # Split into quadrants and save in original directory
-    quadrant_files = []
-    for row in range(2):
-        for col in range(2):
-            quadrant = padded[:, row*512:(row+1)*512, col*512:(col+1)*512]
+    # Split into 16ths (4x4 grid of 256x256 tiles) and save in original directory
+    tile_files = []
+    for row in range(4):
+        for col in range(4):
+            tile = padded[:, row*256:(row+1)*256, col*256:(col+1)*256]
             filename = output_dir / f"{output_name}_{row}_{col}.tif"
-            sio.imsave(str(filename), quadrant)
-            quadrant_files.append(str(filename))
+            sio.imsave(str(filename), tile)
+            tile_files.append(str(filename))
 
-    return quadrant_files
+    return tile_files
 
 # Usage
 if __name__ == "__main__":
     import sys
     if len(sys.argv) > 1:
-        files = process_nd2_to_quadrants(sys.argv[1])
-        print(f"Created: {', '.join(files)}")
+        files = process_nd2_to_sixteenths(sys.argv[1])
+        print(f"Created {len(files)} tiles: {', '.join(files)}")
     else:
         print("Usage: python script.py <file.nd2>")