Triangle Society for Neuroscience 2025 Conference
Matthew McCormick, PhD1, Taylor Hanayik, PhD2, and Chris Rorden, PhD3
1Fideus Labs, Research Triangle Park, NC 27709, matt@fideus.io, 2QuantCo virtual diagnostics team, 3University of South Carolina
Matthew McCormick, PhD1, Taylor Hanayik, PhD2, and Chris Rorden, PhD3
1Fideus Labs, Research Triangle Park, NC 27709, matt@fideus.io, 2QuantCo virtual diagnostics team, 3University of South Carolina
A novel normalized logarithmic gradient magnitude approach, implemented in NiiVue, enhances visualization of internal structures while maintaining real-time performance across devices.
Introduction
-
Neuroimaging Visualization Challenge:
- Visualizing internal brain structures is critical for understanding complex neuroanatomy, functional relationships, and pathology detection.
- Our goal is to provide advanced visualization techniques through an intuitive web interface accessible to researchers without specialized hardware or software installations.
-
NiiVue Capabilities:
- NiiVue [1] is a web-based visualization tool for neuroimaging that runs on any operating system and web device (phone, tablet, computer).
- Uniquely displays multiple neuroimaging datatypes simultaneously: voxels, meshes, tractography streamlines, statistical maps, and connectomes.
- Gradient Opacity Role: Traditional volume rendering often obscures internal structures behind high-intensity regions; gradient opacity allows visualization through tissue based on structural boundaries. [2]
Methods
-
Normalized Logarithmic Gradient Approach:
- Applies logarithmic transformation to compress wide dynamic range of gradient magnitudes
- Normalizes values to range for intuitive user control
$$o = \left( \log_2\left(G_x^2 + G_y^2 + G_z^2 + \frac{1}{255^2 \cdot 8}\right) + \log_2\left(\frac{1}{255^2 \cdot 8}\right) \right)^{8.0}$$Where final exponentiation to 8.0 adds translucency to rendering.
-
Gradient Magnitude Calculation:
- First-order gradients computed using 3D Sobel-approximation operators
- Optionally uses second-order gradients based on Taylor Series expansion for enhanced detail [3]: \(f^1_0 = \frac{f_{-2} - 8 f_{-1} + 8 f_1 - f_2}{ 12 h }\)
-
Real-Time, GPGPU Performance [4]:
- Pre-computes gradients using WebGL2 fragment shaders for parallel processing
- Leverages GPU texture interpolation hardware for optimized sampling
- Encodes 3D gradient orientation and magnitude in only 4 bytes total per voxel
- One-time preprocessing cost of ~100ms for typical brain volumes
- Minimal impact on frame rate during interactive viewing
Results
Gradient Opacity 0.5: Ventricles and major sulci become visible.
Gradient Opacity 1.0: Complete visualization of internal structures.
Discussion
- Power-Law Distribution: Natural structures exhibit gradient magnitudes that follow power-law distributions following Zipf's law [6].
-
A logarithmic transformation is particularly
well-suited for neuroimaging data as it:
- Compresses the wide dynamic range of gradient magnitudes
- Emphasizes subtle structural boundaries within similar intensity regions
- Produces more perceptually linear visual results
- Ultrasound Comparison: Similar logarithmic transformations are standard in B-mode ultrasound imaging, where power-law distributions are well-documented in tissue interactions [5].
-
User Interface Considerations:
- Single intuitive parameter (0-1 scale) controls the effect across all modalities
- Consistent behavior regardless of image resolution or intensity range
- Preserves original color mapping and intensity relationships
-
Future Directions: Implementation designed to
handle datasets from microscopic to macroscopic scales:
- Will support terabyte-scale OME-Zarr [7] imaging datasets in future work
- Maintain consistent visualization from subcellular to whole-brain volumes [6]
Figure 8: OME-Zarr lightsheet human cortex volume [8]
rendered at multiple scales.
References
- NiiVue/niivue: a WebGL2 based medical image viewer. GitHub. (2025). https://github.com/niivue/niivue
- Engel, K., Hadwiger, M., Kniss, J. M., & Rezk-Salama, C. (2006). Local volume illumination. In Real-time volume graphics (pp. 47-52). Eurographics Association. Retrieved April 1, 2025, from https://webdocs.cs.ualberta.ca/~pierreb/Visualization2006/Real-Time-Volume-Rendering.pdf
- Khan, IR and Ohba, Ryoji. Closed-form expressions for the finite difference approximations of first and higher derivatives based on Taylor series. Journal of Computational and Applied Mathematics, 107:179-193, 1999. https://doi.org/10.1016/S0377-0427(99)00088-6.
- Ikits, M., Kniss, J., Lefohn, A., & Hansen, C. (2004). Volume rendering techniques. In R. Fernando (Ed.), GPU Gems: Programming Techniques, Tips, and Tricks for Real-Time Graphics (pp. 667-690). Addison-Wesley. Retrieved April 1, 2025, from https://developer.nvidia.com/gpugems/gpugems/part-vi-beyond-triangles/chapter-39-volume-rendering-techniques
- Parker, K.J. (2022). Power laws prevail in ultrasound-tissue interactions. PMC, 9118335.
- Roman Sabin and Bertolotti Francesco (2022). A master equation for power laws. R. Soc. Open Sci.9220531. https://royalsocietypublishing.org/doi/10.1098/rsos.220531
- Moore, J., Basurto-Lozada, D., Besson, S. et al. OME-Zarr: a cloud-optimized bioimaging file format with international community support. Histochem Cell Biol 160, 223-251 (2023). https://doi.org/10.1007/s00418-023-02209-1
- Kamentsky, Lee; Marx, Slayton; Park, Juhyuk; Su-Arcaro, Clover; Moukheiber, Mira; Zhao, Victor (2023) Light sheet imaging of the human brain (Version draft) [Data set]. DANDI archive. https://doi.org/10.80507/dandi.123456/0.123456.1234