Пример #1
0
        private static double CalculateRatio(List <ScatterData> datas, List <int> waveLengths)
        {
            var R11 = datas[0].DetectedPhotons1 / datas[0].DetectedPhotons2;
            var R22 = datas[1].DetectedPhotons1 / datas[1].DetectedPhotons2;
            var R33 = datas[2].DetectedPhotons1 / datas[2].DetectedPhotons2;

            var R12 = R11 / R22;
            var R13 = R11 / R33;

            var alpha12 = Math.Log(R12);
            var alpha13 = Math.Log(R13);

            var alpha1 = Coefficients.ObtainAbsorptionCoefficients((uint)waveLengths[0]);
            var alpha2 = Coefficients.ObtainAbsorptionCoefficients((uint)waveLengths[1]);
            var alpha3 = Coefficients.ObtainAbsorptionCoefficients((uint)waveLengths[2]);

            var alphaD1 = alpha1.AbsorptionBlood;
            var alphaD2 = alpha2.AbsorptionBlood;
            var alphaD3 = alpha3.AbsorptionBlood;

            var alphaO1 = alpha1.AbsorptionOxygenatedBlood;
            var alphaO2 = alpha2.AbsorptionOxygenatedBlood;
            var alphaO3 = alpha3.AbsorptionOxygenatedBlood;

            var topLine    = alpha12 * (alphaD3 + alphaD1) - alpha13 * (alphaD1 + alphaD2);
            var bottomLine = alpha13 * (alpha1.AbsorptionOxygenatedBlood - alpha2.AbsorptionOxygenatedBlood +
                                        alpha1.AbsorptionBlood - alpha2.AbsorptionBlood) -
                             alpha12 * (alpha1.AbsorptionOxygenatedBlood + alpha3.AbsorptionOxygenatedBlood -
                                        alpha1.AbsorptionBlood - alpha3.AbsorptionBlood);

            var ratio = Math.Sqrt(Math.Pow(topLine / bottomLine, 2));

            return(ratio);
        }
Пример #2
0
        public void AbsorptionCoefficientsTest()
        {
            const uint testWavelength1 = 660;
            const uint testWavelength2 = 940;

            var validReturnData1 = new AbsorptionCoefficients
            {
                AbsorptionBlood           = 3200,
                AbsorptionOxygenatedBlood = 316
            };

            var validReturnData2 = new AbsorptionCoefficients
            {
                AbsorptionBlood           = 802,
                AbsorptionOxygenatedBlood = 1214
            };

            var coefficients = coefficientsClass.ObtainAbsorptionCoefficients(testWavelength1);

            Assert.Equal(validReturnData1.AbsorptionBlood, coefficients.AbsorptionBlood);
            Assert.Equal(validReturnData1.AbsorptionOxygenatedBlood, coefficients.AbsorptionOxygenatedBlood);

            coefficients = coefficientsClass.ObtainAbsorptionCoefficients(testWavelength2);

            Assert.Equal(validReturnData2.AbsorptionBlood, coefficients.AbsorptionBlood);
            Assert.Equal(validReturnData2.AbsorptionOxygenatedBlood, coefficients.AbsorptionOxygenatedBlood);
        }
Пример #3
0
        static void Main(string[] args)
        {
            var absorption             = Coefficients.ObtainAbsorptionCoefficients(660);
            var scatteringCoefficients = Coefficients.ObtainScatteringCoefficients(660);

            var data400 = Scatter.Scatterlight(450, 0.1, 0.3, 0.05, 0.05, 1, 3, 0.150, 0.9);
            var data660 = Scatter.Scatterlight(660, 0.1, 0.3, 0.05, 0.05, 1, 3, 0.150, 0.9);
            var data900 = Scatter.Scatterlight(900, 0.1, 0.3, 0.05, 0.05, 1, 3, 0.150, 0.9);

            var datas = new List <ScatterData>
            {
                data400,
                data660,
                data900
            };
            var wavelengths = new List <int> {
                450, 660, 900
            };

            CalculateRatio(datas, wavelengths);
        }