Stack: vertical image intensity correction

stack
#1

hello,

Problem: in a 3D stack, the mean intensity varies along the z - axis according to a certain function, e.g. like a Gaussian. How can I compensate this trend for the entire stack?

#2

Good day Andreas,

what about looping all slices and multiplying every slice by the inverse of the corresponding attenuation?

Regards

Herbie

PS:
Here is an example Macro that does this for the “MRI Stack”:

run("MRI Stack (528K)");
run("32-bit");
setSlice(1);
getRawStatistics(N, refMn);
for ( i=2; i<=nSlices; i++ ) { 
   setSlice(i);
   getRawStatistics(N, mn);
   run("Multiply...", "value="+d2s(refMn/mn,9)+" slice");
   setResult("Original", nResults, d2s(mn,3));
   getRawStatistics(N, mn);
   setResult("Compensated", nResults-1, d2s(mn,3));
}
setSlice(1);
exit();

Paste the above macro code to an empty macro window (Plugins >> New >> Macro) and run it.
(You need an open internet connection to load the demo image.)