This project has moved. For the latest updates, please go here.

Two layer model for brain

Sep 24, 2012 at 10:48 PM


I have created a two layer model using the Monte Carlo solver panel with the following thicknesses, mua, and musp:

I run the simulation and get reflectance as a function of rho.  I import this into Matlab and do a Hankel Transform to get reflectance as a function of frequency.  Then I use spatial frequency domain imaging code, assuming the reflectance is coming from a homogeneous medium, and fit for the perceived mua and musp. 

I'm am hypothesizing this two layer model to explain real data that I have acquired where the scattering does not follow the power law and total hemoglobin is lower and oxygen saturation is higher than expected.  

Am I doing everything correctly?  Do you have any other suggestions to help me explain my data?  -zander

Sep 26, 2012 at 7:26 PM

Any thoughts?

Sep 26, 2012 at 7:43 PM

Hi Xander,  I didn't reply earlier because I don't have much personal experience with your question.  We *do* have an ROfFx detector defined in the Monte Carlo code, but I didn't code it so I'm not sure the state of its completion.  Independent of using that detector, what you outline sounds fine to me.  I'd be interested in finding out what recovered optical properties you determine.

Sep 26, 2012 at 8:24 PM

Hey Carole- If you look at the image link (, the recovered optical properties are in the table on the right as the "fitted" optical properties and hemoglobin.

Sep 27, 2012 at 3:07 PM

Those results look promising.  Would you like to describe your process and results to our Modeling and Computation group that meets on Friday mornings sometime?

Sep 27, 2012 at 4:31 PM

sure, I can come by next Friday (10/5).  What time is it?