Automated computation of nerve fibre inclinations from 3D polarised light imaging measurements of brain tissue

The method 3D polarised light imaging (3D-PLI) measures the birefringence of histological brain sections to determine the spatial course of nerve fibres (myelinated axons). While the in-plane fibre directions can be determined with high accuracy, the computation of the out-of-plane fibre inclinations is more challenging because they are derived from the amplitude of the birefringence signals, which depends e.g. on the amount of nerve fibres. One possibility to improve the accuracy is to consider the average transmitted light intensity (transmittance weighting). The current procedure requires effortful manual adjustment of parameters and anatomical knowledge. Here, we introduce an automated, optimised computation of the fibre inclinations, allowing for a much faster, reproducible determination of fibre orientations in 3D-PLI. Depending on the degree of myelination, the algorithm uses different models (transmittance-weighted, unweighted, or a linear combination), allowing to account for regionally specific behaviour. As the algorithm is parallelised and GPU optimised, it can be applied to large data sets. Moreover, it only uses images from standard 3D-PLI measurements without tilting, and can therefore be applied to existing data sets from previous measurements. The functionality is demonstrated on unstained coronal and sagittal histological sections of vervet monkey and rat brains.