Week 2 - Hannah Sheffield
This week, I shadowed Dr. Allison Levy and Dr. Carolyn Eisen for breast imaging. When shadowing Dr. Levy, I saw her diagnose patients based on their mammograms, ultrasound, and MRI images. She taught me some new terminology such as birads, tomosynthesis, and hypoechoic images. With Dr. Eisen, I shadowed a mammogram biopsy where I learned about the necessary procedures used for obtaining a tissue sample. I was intrigued with how I was able to see how tissue samples were obtained, and were given to a pathologist for further analysis.
As for research, I delve into researching our problem on shifting the delayed arterial input function (AIF - concentration of the contrast agent in the artery throughout time) to align with the reference AIF. To specify, we have a map with a delayed response (called delay map) in capturing the contrast agent concentration, where the TK parameters undergo a loss throughout each epoch loop. To understand the problem on a deeper level, I needed to understand what was happening in the code. I did this by mathematically deriving the gradients of the TK parameters throughout the epoch loops. This gave me a deeper understanding in why a loss occurs between each epoch iteration, and we found that vp (blood plasma volume) is the AIF curve, ktrans (volume transfer constant - rate of contrast agent leaking from the blood vessel to the tissue) is the integral of the AIF, and the delay map is the AIF's derivative (or change in slope) throughout time. Through this, we were able to pin-point the two issues. The first is that the delay map has a very negligible temporal change compared to the TK parameters, which leads to an imbalance between the two metrics. Second is that all variables have very negligible values, which prohibits the epoch loop from capturing the values and iteratively shifting the delayed AIF with the reference AIF.
To counteract these issues, Ari and I decided to explore options to resolve the issue. First, we artificially changed the learning rate of the delay map to check if there is a response. In changing the learning rate from 1e-2 to 1000, there was an appropriate response where delay map shifted. We then explored normalization techniques to see if the map changes there rather than arbitrarily changing the learning rate. We looked into a gradient normalization technique called "GradNorm," which takes the magnitude of the gradients of the delay learning rate, ktrans, and vp so that the magnitude imbalance is resolved. We found a working code that fixes this issue, and attempted to integrate the GradNorms function into our working code. However, this attempt was unsuccessful (when the learning rate was at its original value), and did not shift the delayed AIF. To this end, we have decided to move forward to the next stage of the project, which is running the working code through all cases, and developing a two compartment model pipeline. We hope to revisit this issue at a later time.
Comments
Post a Comment