viewinline 0.2.3__tar.gz → 0.3.1__tar.gz

{viewinline-0.2.3 → viewinline-0.3.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: viewinline
-Version: 0.2.3
+Version: 0.3.1
 Summary: Quick look geospatial viewer for the terminal, with inline image previews
 Project-URL: Homepage, https://github.com/nkeikon/viewinline
 Project-URL: Repository, https://github.com/nkeikon/viewinline
@@ -37,7 +37,7 @@ Description-Content-Type: text/markdown
 [![Python version](https://img.shields.io/badge/python-%3E%3D3.9-blue.svg)](https://pypi.org/project/viewinline/)
 
 **Quick-look geospatial viewer for compatible terminals.**  
-Displays rasters, vectors, and tabular data directly in the terminal with no GUI and no temporary files.
+Displays rasters, vectors, and tabular data directly in the terminal with no GUI.
 
 <p align="center">
   <a href="viewinline_gif1.gif"><img src="viewinline_gif1.gif" width="49%"></a>
@@ -65,13 +65,17 @@ pip install viewinline
 viewinline path/to/file.tif
 viewinline R.tif G.tif B.tif                 # RGB composite (also works with --rgbfiles)
 viewinline path/to/multiband.tif --rgb 3 2 1
-viewinline path/to/folder --gallery 4x3      # show image gallery (e.g. 4x3 grid)
+viewinline path/to/folder --gallery 4x3      # show image gallery of all files in the folder (e.g. 4x3 grid)
+viewinline path/to/hyperspectral.tif --bands 10-50 # show image gallery of selected bands (also works with --bands 11,15,30,45; --gallery 5x5)
 
 # NetCDF and HDF
 viewinline file.nc                           # list variables
 viewinline file.nc --subset 2                # display variable 2
 viewinline file.nc --subset 1 --band 10      # variable 1, timestep 10 --band or --timestep
 viewinline temp.nc --subset 1 --colormap plasma --vmin 273 --vmax 310
+viewinline hyperspectral.nc --subset 1 --reduce NumberOfScanlines  # override auto-detected axis
+viewinline hyperspectral.nc --subset 22 --band 50 # show image gallery of selected bands
+viewinline hyperspectral.nc --subset 22 --bands 10-54 --gallery 5x11
 
 # Vectors
 viewinline path/to/vector.geojson
@@ -113,6 +117,7 @@ Install chafa once (it's a system binary, available across all conda/virtualenv
 brew install chafa          # macOS
 sudo apt install chafa      # Debian/Ubuntu
 sudo dnf install chafa      # Fedora
+scoop install chafa         # Windows
 ```
 Without chafa, terminals outside the native list above show an info message instead of an image.
 You can also force the chafa path on any terminal by setting `INLINE_VIEWER_ENGINE=chafa`.
@@ -166,10 +171,13 @@ You can also force the chafa path on any terminal by setting `INLINE_VIEWER_ENGI
 
 **Gallery view**
 - Display all images in a folder with `--gallery 4x4`
+- Display selected bands of a single raster as a grid with `--bands 101-120` or `--bands 11,12,45,55`. Works with GeoTIFF and NetCDF files.
 
 **NetCDF/HDF notes:**
 - viewinline lists only variables that can be displayed as 2D or 3D arrays
-- Variables with additional dimensions (e.g., vertical levels) may be listed but will fail to display with a clear error message
+- 3D variables with time or known spatial dimensions are auto-handled (slices along the non-spatial axis)
+- For 3D variables with non-standard dimensions (e.g., hyperspectral cubes like PICARD), viewinline auto-detects the band axis by smallest dimension. Use `--reduce DIM_NAME` to override.
+- Variables with 4+ dimensions are not supported
 - For a complete variable list, use `ncdump -h file.nc` or `viewtif`
 
 ## Dependencies
@@ -198,7 +206,7 @@ You can also force the chafa path on any terminal by setting `INLINE_VIEWER_ENGI
 
 **Note on HDF support:**
 - **HDF5** (.h5, .hdf5): Supported via rasterio if GDAL has HDF5 support (most installations)
-- **HDF4** (.hdf): Requires GDAL compiled with HDF4 support (common in MODIS data)
+- **HDF4** (.hdf): Requires GDAL compiled with HDF4 support (the legacy format used by MODIS and older NASA products)
 - **NetCDF** (.nc): Supported via rasterio (uses GDAL's NetCDF driver)
 
 ## Available options
@@ -208,8 +216,10 @@ General:
 
 Raster:
   --band BAND           Band number to display (single raster), or slice number for NetCDF. (default: 1)
+  --bands BANDS         Display multiple bands as a grid. Accepts ranges (30-40), lists (3,4,5), or mixed (1,5,10-15).
   --timestep INTEGER    Alias for --band when working with NetCDF files.
   --subset INTEGER      Variable index for NetCDF/HDF files (e.g., --subset 1).
+  --reduce DIM_NAME     For 3D NetCDF variables, specify which dimension to use as the band/slider axis.  Auto-detected if omitted.
   --colormap            Apply colormap to single-band rasters. Flag without the color scheme → 'terrain'.
   --rgb R G B           Three band numbers for RGB display (e.g., --rgb 4 3 2). Overrides default 1 2 3.
   --rgbfiles R G B      Three single-band rasters for RGB composite. Can also provide as positional arguments.

{viewinline-0.2.3 → viewinline-0.3.1}/README.md

@@ -4,7 +4,7 @@
 [![Python version](https://img.shields.io/badge/python-%3E%3D3.9-blue.svg)](https://pypi.org/project/viewinline/)
 
 **Quick-look geospatial viewer for compatible terminals.**  
-Displays rasters, vectors, and tabular data directly in the terminal with no GUI and no temporary files.
+Displays rasters, vectors, and tabular data directly in the terminal with no GUI.
 
 <p align="center">
   <a href="viewinline_gif1.gif"><img src="viewinline_gif1.gif" width="49%"></a>
@@ -32,13 +32,17 @@ pip install viewinline
 viewinline path/to/file.tif
 viewinline R.tif G.tif B.tif                 # RGB composite (also works with --rgbfiles)
 viewinline path/to/multiband.tif --rgb 3 2 1
-viewinline path/to/folder --gallery 4x3      # show image gallery (e.g. 4x3 grid)
+viewinline path/to/folder --gallery 4x3      # show image gallery of all files in the folder (e.g. 4x3 grid)
+viewinline path/to/hyperspectral.tif --bands 10-50 # show image gallery of selected bands (also works with --bands 11,15,30,45; --gallery 5x5)
 
 # NetCDF and HDF
 viewinline file.nc                           # list variables
 viewinline file.nc --subset 2                # display variable 2
 viewinline file.nc --subset 1 --band 10      # variable 1, timestep 10 --band or --timestep
 viewinline temp.nc --subset 1 --colormap plasma --vmin 273 --vmax 310
+viewinline hyperspectral.nc --subset 1 --reduce NumberOfScanlines  # override auto-detected axis
+viewinline hyperspectral.nc --subset 22 --band 50 # show image gallery of selected bands
+viewinline hyperspectral.nc --subset 22 --bands 10-54 --gallery 5x11
 
 # Vectors
 viewinline path/to/vector.geojson
@@ -80,6 +84,7 @@ Install chafa once (it's a system binary, available across all conda/virtualenv
 brew install chafa          # macOS
 sudo apt install chafa      # Debian/Ubuntu
 sudo dnf install chafa      # Fedora
+scoop install chafa         # Windows
 ```
 Without chafa, terminals outside the native list above show an info message instead of an image.
 You can also force the chafa path on any terminal by setting `INLINE_VIEWER_ENGINE=chafa`.
@@ -133,10 +138,13 @@ You can also force the chafa path on any terminal by setting `INLINE_VIEWER_ENGI
 
 **Gallery view**
 - Display all images in a folder with `--gallery 4x4`
+- Display selected bands of a single raster as a grid with `--bands 101-120` or `--bands 11,12,45,55`. Works with GeoTIFF and NetCDF files.
 
 **NetCDF/HDF notes:**
 - viewinline lists only variables that can be displayed as 2D or 3D arrays
-- Variables with additional dimensions (e.g., vertical levels) may be listed but will fail to display with a clear error message
+- 3D variables with time or known spatial dimensions are auto-handled (slices along the non-spatial axis)
+- For 3D variables with non-standard dimensions (e.g., hyperspectral cubes like PICARD), viewinline auto-detects the band axis by smallest dimension. Use `--reduce DIM_NAME` to override.
+- Variables with 4+ dimensions are not supported
 - For a complete variable list, use `ncdump -h file.nc` or `viewtif`
 
 ## Dependencies
@@ -165,7 +173,7 @@ You can also force the chafa path on any terminal by setting `INLINE_VIEWER_ENGI
 
 **Note on HDF support:**
 - **HDF5** (.h5, .hdf5): Supported via rasterio if GDAL has HDF5 support (most installations)
-- **HDF4** (.hdf): Requires GDAL compiled with HDF4 support (common in MODIS data)
+- **HDF4** (.hdf): Requires GDAL compiled with HDF4 support (the legacy format used by MODIS and older NASA products)
 - **NetCDF** (.nc): Supported via rasterio (uses GDAL's NetCDF driver)
 
 ## Available options
@@ -175,8 +183,10 @@ General:
 
 Raster:
   --band BAND           Band number to display (single raster), or slice number for NetCDF. (default: 1)
+  --bands BANDS         Display multiple bands as a grid. Accepts ranges (30-40), lists (3,4,5), or mixed (1,5,10-15).
   --timestep INTEGER    Alias for --band when working with NetCDF files.
   --subset INTEGER      Variable index for NetCDF/HDF files (e.g., --subset 1).
+  --reduce DIM_NAME     For 3D NetCDF variables, specify which dimension to use as the band/slider axis.  Auto-detected if omitted.
   --colormap            Apply colormap to single-band rasters. Flag without the color scheme → 'terrain'.
   --rgb R G B           Three band numbers for RGB display (e.g., --rgb 4 3 2). Overrides default 1 2 3.
   --rgbfiles R G B      Three single-band rasters for RGB composite. Can also provide as positional arguments.

{viewinline-0.2.3 → viewinline-0.3.1}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "hatchling.build"
 
 [project]
 name = "viewinline"
-version = "0.2.3"
+version = "0.3.1"
 description = "Quick look geospatial viewer for the terminal, with inline image previews"
 readme = "README.md"
 license = { text = "Apache-2.0" }

{viewinline-0.2.3 → viewinline-0.3.1}/src/viewinline/viewinline.py

@@ -28,12 +28,19 @@ from matplotlib import colormaps
 import matplotlib as mpl
 import subprocess
 
+try:
+    import netCDF4
+    HAS_NETCDF4 = True
+except ImportError:
+    HAS_NETCDF4 = False
+
 import warnings
 
 warnings.filterwarnings("ignore", message="More than one layer found", category=UserWarning)
 warnings.filterwarnings("ignore", message="Dataset has no geotransform", category=UserWarning)
+warnings.filterwarnings("ignore", message="invalid scale_factor or add_offset attribute", category=UserWarning)
 
-__version__ = "0.2.3"
+__version__ = "0.3.1"
 
 AVAILABLE_COLORMAPS = [
     "viridis", "inferno", "magma", "plasma",
@@ -184,14 +191,7 @@ def show_inline_image(image_array: np.ndarray, display_scale = None, is_vector:
 
     if _TERMINAL_SUPPORTS_IMAGES:
         sys.stdout.write(f"\033]1337;File=inline=1;width={width_pct}%:{encoded}\a\n")
-    else:
-        # if is_chafa_available():
-        #     chafa_output = subprocess.check_output(
-        #         ["chafa", "-"],  
-        #         input=image_bytes
-        #     ).decode()
-            
-        #     sys.stdout.write(f"\n{chafa_output}\a\n")
+    else: 
         if is_chafa_available():
             # Inside tmux, force chafa to use block-art symbols instead of
             # graphics protocols. Tmux mangles kitty graphics and sixel
@@ -227,9 +227,9 @@ def show_image_auto(img: np.ndarray, display_scale=None, is_vector: bool = False
     try:
         show_inline_image(img, display_scale, is_vector)
         if _TERMINAL_SUPPORTS_IMAGES:
-            print("[VIEW] Image rendered")
+            print("[VIEW] Inline render complete")
         elif is_chafa_available():
-            print("[VIEW] Rendered via chafa")
+            print("[VIEW] Inline render complete via chafa")
         # If neither path applies, show_inline_image already printed the info message
     except Exception as e:
         print(f"[ERROR] Failed to render image: {e}")
@@ -248,6 +248,24 @@ def resize_to_terminal(img: np.ndarray) -> tuple[np.ndarray, float]:
     pil_img = ImageOps.contain(pil_img, (new_w, new_h))
     return np.array(pil_img), scale
 
+def parse_bands(s: str) -> list[int]:
+    """Parse --bands argument: accepts ranges (30-40), lists (3,4,5), or mixed (1,3,10-15)."""
+    bands = []
+    for part in s.split(","):
+        part = part.strip()
+        if "-" in part:
+            try:
+                start, end = map(int, part.split("-", 1))
+                bands.extend(range(start, end + 1))
+            except ValueError:
+                print(f"[WARN] Could not parse band range: {part}")
+        else:
+            try:
+                bands.append(int(part))
+            except ValueError:
+                print(f"[WARN] Could not parse band: {part}")
+    return sorted(set(bands))
+
 # ---------------------------------------------------------------------
 # CSV handling
 # ---------------------------------------------------------------------
@@ -262,7 +280,6 @@ def load_csv_to_df(path: str) -> pd.DataFrame:
         print(f"[ERROR] Failed to read CSV: {e}")
         return pd.DataFrame()
 
-
 # =============================================================
 # Preview
 # =============================================================
@@ -643,6 +660,196 @@ def render_simple_image(filepath: str, args) -> None:
     except Exception as e:
         print(f"[ERROR] Failed to load image: {e}")
 
+def render_netcdf_via_netcdf4(path, args):
+    """Read a NetCDF file via netCDF4 (bypassing GDAL). Handles hierarchical
+    groups and hyperspectral cubes where GDAL aborts or interprets axes wrong.
+    """
+    if not HAS_NETCDF4:
+        print("[ERROR] netCDF4 not installed. Install with:")
+        print("        pip install netCDF4")
+        print("        or:  pip install viewinline[netcdf]")
+        return
+
+    try:
+        nc = netCDF4.Dataset(path)
+    except Exception as e:
+        print(f"[ERROR] Could not open NetCDF file: {e}")
+        return
+
+    # Recursively collect (path, variable) pairs across all groups
+    def collect_vars(group, prefix=""):
+        out = []
+        for name, var in group.variables.items():
+            full_name = f"{prefix}{name}"
+            out.append((full_name, var))
+        for sub_name, sub in group.groups.items():
+            out.extend(collect_vars(sub, f"{prefix}{sub_name}/"))
+        return out
+
+    all_vars = collect_vars(nc)
+
+    if not all_vars:
+        print("[ERROR] No variables found in file.")
+        nc.close()
+        return
+
+    # If no --subset, list all variables and exit
+    if not args.subset:
+        print(f"Found {len(all_vars)} variables in {os.path.basename(path)}:")
+        for i, (name, var) in enumerate(all_vars, 1):
+            shape_str = "x".join(str(s) for s in var.shape)
+            print(f"  [{i}] {name}  ({shape_str}, {var.dtype})")
+        print(f"\nUse --subset <N> to display a specific variable.")
+        nc.close()
+        return
+
+    # Validate --subset
+    if args.subset < 1 or args.subset > len(all_vars):
+        print(f"[ERROR] --subset must be between 1 and {len(all_vars)}")
+        nc.close()
+        return
+
+    var_name, var = all_vars[args.subset - 1]
+    print(f"[INFO] Displaying variable {args.subset}: {var_name}")
+    print(f"[DATA] Shape: {var.shape}  dtype: {var.dtype}  dims: {var.dimensions}")
+
+    # Detect dimensionality and read the right slice
+    if var.ndim == 2:
+        data = np.asarray(var[:, :], dtype=np.float64)
+        slice_info = "2D variable"
+
+    elif var.ndim == 3:
+        spatial_dims = {'lat', 'lon', 'latitude', 'longitude', 'y', 'x'}
+        
+        spectral_axis = None
+        
+        # 1. User override via --reduce
+        if args.reduce_dim is not None:
+            if args.reduce_dim in var.dimensions:
+                spectral_axis = list(var.dimensions).index(args.reduce_dim)
+                print(f"[INFO] Using user-specified --reduce '{args.reduce_dim}'")
+            else:
+                print(f"[ERROR] --reduce '{args.reduce_dim}' is not a dimension of this variable.")
+                print(f"[INFO] Available dimensions: {list(var.dimensions)}")
+                nc.close()
+                return
+        
+        # 2. Standard convention: reduce along the non-spatial dim
+        if spectral_axis is None:
+            has_standard_spatial = any(d in spatial_dims for d in var.dimensions)
+            if has_standard_spatial:
+                for i, dim_name in enumerate(var.dimensions):
+                    if dim_name not in spatial_dims:
+                        spectral_axis = i
+                        break
+        
+        # 3. Fallback heuristic: smallest dim is typically the band axis
+        if spectral_axis is None:
+            sizes = [(i, var.shape[i]) for i in range(3)]
+            spectral_axis = min(sizes, key=lambda x: x[1])[0]
+            print(f"[INFO] Non-standard dimensions detected: {list(var.dimensions)}")
+            print(f"[INFO] Reducing along '{var.dimensions[spectral_axis]}' (size {var.shape[spectral_axis]}, assumed band/spectral axis)")
+            print(f"[INFO] If this is not correct, use --reduce DIM_NAME to override.")
+        
+        # Slice along chosen axis
+        band_count = var.shape[spectral_axis]
+        band_num = args.band if args.band is not None else 1
+        band_idx = max(0, min(band_num - 1, band_count - 1))
+
+        # BANDS GALLERY for NetCDF
+        if getattr(args, "bands", None):
+            band_list = parse_bands(args.bands)
+            print(f"[DEBUG] band_count={band_count}, band_list={band_list}")
+            valid_bands = [b for b in band_list if 1 <= b <= band_count]
+            if not valid_bands:
+                print(f"[ERROR] No valid bands. Variable has {band_count} bands along '{var.dimensions[spectral_axis]}'.")
+                nc.close()
+                return
+            slices = []
+            for b in valid_bands:
+                slicer = [slice(None)] * 3
+                slicer[spectral_axis] = b - 1
+                slices.append(np.asarray(var[tuple(slicer)], dtype=np.float64))
+            nc.close()
+            colormap = args.colormap if args.colormap else "viridis"
+            render_bands_gallery(np.stack(slices, axis=0), valid_bands, band_count,
+                                display_scale=getattr(args, "display", None),
+                                colormap=colormap)
+            return
+
+        slicer = [slice(None)] * 3
+        slicer[spectral_axis] = band_idx
+        data = np.asarray(var[tuple(slicer)], dtype=np.float64)
+        slice_info = f"slice along axis {spectral_axis} ({var.dimensions[spectral_axis]}), band {band_idx + 1} of {band_count}"
+
+    else:
+        print(f"[ERROR] viewinline only supports 2D or 3D variables. This one is {var.ndim}D.")
+        nc.close()
+        return
+
+    print(f"[DATA] {slice_info}")
+    # Apply fill value
+    fill = getattr(var, '_FillValue', None)
+    if fill is not None:
+        data = np.where(data == fill, np.nan, data)
+    # Flip vertically if data is stored south-to-north so north appears at top.
+    # Determine which dims remain after slicing — for the 2D result, figure out
+    # which axis (0 or 1) corresponds to latitude, and check that dim's coord values.
+    if var.ndim == 2:
+        remaining_dims = list(var.dimensions)
+    elif var.ndim == 3:
+        if spectral_axis is not None:
+            remaining_dims = [d for i, d in enumerate(var.dimensions) if i != spectral_axis]
+        else:
+            remaining_dims = list(var.dimensions[1:])  # axis 0 was reduced
+    else:
+        remaining_dims = []
+    lat_names = {'lat', 'latitude', 'y'}
+    for axis_in_2d, dim_name in enumerate(remaining_dims):
+        if dim_name in lat_names and dim_name in nc.variables:
+            lat_vals = nc[dim_name][:]
+            if len(lat_vals) > 1 and lat_vals[0] < lat_vals[-1]:
+                data = np.flip(data, axis=axis_in_2d)
+                print(f"[INFO] Flipped along '{dim_name}' for display (data stored south-to-north).")
+            break
+    nc.close()
+
+    # Normalize and display
+    band_u8 = normalize_to_uint8(data, vmin=args.vmin, vmax=args.vmax,
+                                 nodata=args.nodata)
+
+    if args.colormap:
+        cmap = colormaps[args.colormap]
+        colored = cmap(band_u8 / 255.0)
+        img = (colored[:, :, :3] * 255).astype(np.uint8)
+        print(f"[INFO] Applying colormap: {args.colormap}")
+    else:
+        img = np.stack([band_u8] * 3, axis=-1)
+        print("[INFO] Displaying grayscale")
+
+    # Resize to terminal
+    H, W = img.shape[:2]
+    if args.display:
+        new_w, new_h = max(1, int(W * args.display)), max(1, int(H * args.display))
+        img = np.array(Image.fromarray(img).resize((new_w, new_h), Image.BILINEAR))
+        print(f"[VIEW] Manual resize ×{args.display:.2f} → {new_w}×{new_h}px")
+    else:
+            max_dim = 2000
+            if max(img.shape[:2]) > max_dim:
+                scale = max_dim / max(img.shape[:2])
+                new_w = int(img.shape[1] * scale)
+                new_h = int(img.shape[0] * scale)
+                img = np.array(Image.fromarray(img).resize((new_w, new_h), Image.BILINEAR))
+                print(f"[VIEW] Downsampled from {W}×{H}px to {new_w}×{new_h}px (display scale={scale:.2f})")
+                print(f"[INFO] Use --display 1 for full resolution.")
+            else:
+            # (matches the width_pct logic in show_inline_image)
+                display_pct = args.display if args.display is not None else 0.33
+            
+                print(f"[VIEW] Rendered image size → {img.shape[1]}×{img.shape[0]}px (size={display_pct:.2f})")  
+       
+    show_image_auto(img, getattr(args, "display", None), is_vector=False)
+
 def render_raster(paths: list[str], args) -> None:
     try:
         import rasterio
@@ -655,8 +862,13 @@ def render_raster(paths: list[str], args) -> None:
 
         if len(paths) == 1:
             path = paths[0]
+
+            if path.lower().endswith('.nc'):
+                render_netcdf_via_netcdf4(path, args)
+                return
             
             # Handle NetCDF/HDF with subdatasets
+        
             if path.lower().endswith(('.nc', '.hdf', '.hdf5', '.h5')):
                 try:
                 
@@ -732,7 +944,8 @@ def render_raster(paths: list[str], args) -> None:
                         resampling=rasterio.enums.Resampling.bilinear
                     )
 
-                    print(f"[PROC] Downsampled for preview → {out_w}×{out_h}px (scale={scale:.3f})")
+                    print(f"[VIEW] Downsampled for preview → {out_w}×{out_h}px (display scale={scale:.3f})")
+                    print(f"[INFO] Use --display 1 for full resolution.")
                 else:
                     data = ds.read()
 
@@ -743,9 +956,18 @@ def render_raster(paths: list[str], args) -> None:
                 else:
                     print(f"[INFO] Multi-band raster detected ({band_count} bands)")
 
-            # MULTI BAND RGB (skip for NetCDF - treat as slices/timesteps, not RGB)
-            if band_count >= 3 and not paths[0].lower().endswith('.nc'):
+            # BANDS GALLERY
+            if getattr(args, "bands", None):
+                band_list = parse_bands(args.bands)
+                colormap = args.colormap if args.colormap else "viridis"
+                render_bands_gallery(data, band_list, band_count,
+                                    grid=getattr(args, "gallery", None),
+                                    display_scale=getattr(args, "display", None),
+                                    colormap=colormap)
+                return
 
+            user_specified_band = args.band is not None
+            if band_count >= 3 and not paths[0].lower().endswith('.nc') and not user_specified_band and getattr(args, 'rgb', None):
 
                 if getattr(args, "rgb", None):
                     try:
@@ -774,7 +996,8 @@ def render_raster(paths: list[str], args) -> None:
             # SINGLE BAND
             else:
 
-                band_idx = max(0, min(args.band - 1, band_count - 1))
+                band_num = args.band if args.band is not None else 1
+                band_idx = max(0, min(band_num - 1, band_count - 1))
                 # print(f"[INFO] Displaying band {band_idx + 1} of {band_count}")
                 raw_band = data[band_idx].astype(float)
 
@@ -851,7 +1074,7 @@ def render_raster(paths: list[str], args) -> None:
             print(f"[ERROR] Cannot display this variable.")
             print("[INFO] viewinline only supports 2D or 3D NetCDF variables")
         else:
-            print(f"[ERROR] Raster rendering failed: {e}")
+            print(f"[ERROR] Inline render failed: {e}")
 
 
 def render_gallery(folder: str, grid: str = "4x4", display_scale=None, is_vector=False) -> None:
@@ -866,12 +1089,11 @@ def render_gallery(folder: str, grid: str = "4x4", display_scale=None, is_vector
             cols, rows = 4, 4
         nmax = cols * rows
 
-        # Collect files
-        exts = (".png", ".jpg", ".jpeg", ".tif", ".tiff")
+        # Collect all files, let the loader decide what's valid
         files = [os.path.join(folder, f) for f in sorted(os.listdir(folder))
-                 if f.lower().endswith(exts)]
+                 if os.path.isfile(os.path.join(folder, f))]
         if not files:
-            print(f"[WARN] No image/raster files found in {folder}")
+            print(f"[WARN] No files found in {folder}")
             return
 
         files = files[:nmax]
@@ -884,8 +1106,12 @@ def render_gallery(folder: str, grid: str = "4x4", display_scale=None, is_vector
                 ext = os.path.splitext(f)[1].lower()
                 if ext in [".tif", ".tiff"]:
                     import rasterio
+                    from rasterio.enums import Resampling
                     with rasterio.open(f) as ds:
-                        arr = ds.read()
+                        arr = ds.read(
+                            out_shape=(ds.count, thumb_size[1], thumb_size[0]),
+                            resampling=Resampling.nearest
+                        )
                         if arr.shape[0] >= 3:
                             rgb = np.stack([normalize_to_uint8(arr[i]) for i in range(3)], axis=-1)
                         else:
@@ -896,8 +1122,9 @@ def render_gallery(folder: str, grid: str = "4x4", display_scale=None, is_vector
                     img = Image.open(f).convert("RGB")
                 img.thumbnail(thumb_size)
                 thumbs.append(img)
+
             except Exception as e:
-                print(f"[WARN] Skipped {os.path.basename(f)} ({e})")
+                print(f"[SKIP] {os.path.basename(f)} — {e}")
 
         if not thumbs:
             print("[WARN] No valid images loaded.")
@@ -925,6 +1152,87 @@ def render_gallery(folder: str, grid: str = "4x4", display_scale=None, is_vector
     except Exception as e:
         print(f"[ERROR] Failed to render gallery: {e}")
 
+def render_bands_gallery(data: np.ndarray, band_list: list[int], band_count: int,
+                          grid: str = None, display_scale=None, colormap: str = "viridis") -> None:
+    """Render multiple bands from a single raster as a grid of thumbnails."""
+    import math
+    from PIL import ImageDraw
+
+    # Validate bands
+    valid_bands = [b for b in band_list if 1 <= b <= band_count]
+    if not valid_bands:
+        print(f"[ERROR] No valid bands. File has {band_count} bands.")
+        return
+    skipped = set(band_list) - set(valid_bands)
+    if skipped:
+        print(f"[WARN] Skipped out-of-range bands: {sorted(skipped)}")
+
+    n = len(valid_bands)
+
+    # Auto grid or user-specified
+    if grid:
+        try:
+            cols, rows = map(int, grid.lower().split("x"))
+        except Exception:
+            cols = math.ceil(math.sqrt(n))
+            rows = math.ceil(n / cols)
+    else:
+        cols = math.ceil(math.sqrt(n))
+        rows = math.ceil(n / cols)
+
+    # Use fixed longest dimension to preserve aspect ratio
+    max_thumb = 256
+    label_height = 14  # space below each thumbnail for band label
+    margin = 8
+    cmap = colormaps[colormap]
+
+    thumbs = []
+    thumb_w = thumb_h = max_thumb  # will be updated from first image
+
+    for i, b in enumerate(valid_bands):
+        arr = data[i].astype(float)
+        normalized = normalize_to_uint8(arr)
+        colored = cmap(normalized / 255.0)
+        rgb = (colored[:, :, :3] * 255).astype(np.uint8)
+        img = Image.fromarray(rgb)
+        img.thumbnail((max_thumb, max_thumb), Image.LANCZOS)
+        if i == 0:
+            thumb_w, thumb_h = img.size
+        thumbs.append(img)
+
+    # Build canvas with extra height per row for labels
+    cols = min(cols, n)
+    rows = math.ceil(n / cols)
+    cell_w = thumb_w + margin
+    cell_h = thumb_h + label_height + margin
+    canvas_w = cols * cell_w + margin
+    canvas_h = rows * cell_h + margin
+    canvas = Image.new("RGB", (canvas_w, canvas_h), (220, 220, 220))
+    draw = ImageDraw.Draw(canvas)
+    from PIL import ImageFont
+    try:
+        font = ImageFont.truetype("/System/Library/Fonts/Helvetica.ttc", 10)
+    except Exception:
+        try:
+            font = ImageFont.truetype("/usr/share/fonts/truetype/dejavu/DejaVuSans.ttf", 10)
+        except Exception:
+            font = ImageFont.load_default()
+
+    for i, img in enumerate(thumbs):
+        r, c = divmod(i, cols)
+        x = margin + c * cell_w
+        y = margin + r * cell_h
+        canvas.paste(img, (x, y))
+        # Draw label on canvas background below the thumbnail
+        label = f"B{valid_bands[i]}"
+        lx = x + (thumb_w - len(label) * 6) // 2
+        ly = y + thumb_h + 2
+        draw.text((lx, ly), label, fill=(0, 0, 0), font=font)
+
+    print(f"[INFO] Displaying {n} bands ({cols}×{rows} grid, colormap: {colormap})")
+    # print(f"[DEBUG] canvas size: {canvas_w}×{canvas_h}px")
+    show_image_auto(np.array(canvas), display_scale)
+
 # ---------------------------------------------------------------------
 # Vector handling
 # ---------------------------------------------------------------------
@@ -956,13 +1264,6 @@ def render_vector(path, args):
     except Exception as e:
         print(f"[WARN] Could not list layers: {e}")
 
-    # try:
-    #     gdf = gpd.read_file(path, layer=getattr(args, "layer", None))
-    #     print(f"[DATA] Vector loaded: {os.path.basename(path)} ({len(gdf)} features)")
-    # except Exception as e:
-    #     print(f"[ERROR] Failed to read vector: {e}")
-    #     return
-
     try:
         # Use read_parquet for parquet/geoparquet files
         if path.lower().endswith(('.parquet', '.geoparquet')):
@@ -1304,7 +1605,7 @@ def main() -> None:
 
     # Raster options
     parser.add_argument(
-        "--band", type=int, default=1,
+        "--band", type=int, default=None,
         help="Band number to display (single raster case), or slice number for NetCDF."
     )
     parser.add_argument(
@@ -1320,6 +1621,10 @@ def main() -> None:
         "--rgb", nargs=3, type=int, metavar=('R', 'G', 'B'), default=None,
         help="Three band numbers for RGB display (e.g., --rgb 4 3 2). Overrides default 1 2 3."
     )
+    parser.add_argument(
+        "--rgbfiles", nargs=3, type=str, metavar=('R', 'G', 'B'),
+        help="Three single-band rasters for RGB composite (e.g., --rgbfiles R.tif G.tif B.tif). Can also provide as positional arguments without the flag."
+    )
     parser.add_argument(
         "--vmin", type=float, default=None,
         help="Minimum pixel value for raster display scaling."
@@ -1341,8 +1646,14 @@ def main() -> None:
         help="Variable index for NetCDF files (e.g. --subset 1)."
     )
     parser.add_argument(
-        "--rgbfiles", nargs=3, type=str, metavar=('R', 'G', 'B'),
-        help="Three single-band rasters for RGB composite (e.g., --rgbfiles R.tif G.tif B.tif). Can also provide as positional arguments without the flag."
+        "--reduce", dest="reduce_dim", type=str, default=None,
+        metavar="DIM_NAME",
+        help="For 3D NetCDF variables, specify which dimension to use as the band axis (auto-detected if omitted)."
+    )
+    parser.add_argument(
+        "--bands",
+        type=str,
+        help="Display multiple bands as a grid. Accepts ranges (30-40), lists (3,4,5), or mixed (1,3,10-15)."
     )
 
     # CSV options
@@ -1426,7 +1737,7 @@ def main() -> None:
     parser.add_argument(
     "--table", action="store_true",
     help="Display vector/parquet file as tabular data instead of rendering geometry."
-)
+    ) 
 
     parser.add_argument("--version", action="version", version=f"%(prog)s {__version__}")