MainView.xaml.cs

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
using System.Windows;
using System.Windows.Controls;
using System.Windows.Data;
using System.Windows.Documents;
using System.Windows.Input;
using System.Windows.Media;
using System.Windows.Media.Imaging;
using System.Windows.Navigation;
using System.Windows.Shapes;
using VMS.TPS.Common.Model.API;
using VMS.TPS.Common.Model.Types;


namespace BED_Calc
{
    /// <summary>
    /// Interaction logic for MainView.xaml
    /// </summary>
    public partial class MainView : UserControl
    {
        // specify properties, including one for the open script context
        private int numberOfFractions;
        private string organName;
        private float alphaBetaRatio;
        public float bed;
        public float eqd2;
        private PlanSetup planSetup;

        private ScriptContext context;

        public MainView(ScriptContext context)
        {
            InitializeComponent();
            this.context = context; //Copy currently loaded patient data content
            planSetup = context.PlanSetup;
            numberOfFractions= (int)planSetup.NumberOfFractions; //extract number of fractions from open plancontext

        } 
            public MainView()
        {
            InitializeComponent();
        }

        private void Button_Click(object sender, RoutedEventArgs e)
        {
            //set organ name and alpha beta ratio on what was specified
            organName=OrganNameTB.Text;
            alphaBetaRatio = float.Parse(AlphaBetaRatioTB.Text);

            StructureSet structureSet = planSetup.StructureSet;
            if (structureSet == null)
                throw new ApplicationException("The selected plan does not reference a StructureSet.");

            //find structure if it exists, if it does, set target equal to that structure. throw exce
            Structure target = null;

            bool StructureEmpty = false;//boolean to check if empty structure
            foreach (var structure in structureSet.Structures)
                {   
                    if (structure.Id == organName) //if structure Id matches the f
                    {
                        if (structure.IsEmpty)
                        {
                            StructureEmpty = true;
                            MessageBox.Show(String.Format("Structure {0} is empty, try another structure", organName));
                            target = null;
                            break;
                        }
                        else
                        {
                            target = structure;
                            break;
                        }

                    }
                }
            //Message for debugging
            //MessageBox.Show(String.Format("StructureEmpty is {0}. targetnull is {1}", StructureEmpty, target == null));


            if (!StructureEmpty)
            {

                if (target == null)
                    MessageBox.Show(String.Format("I can't find a structure called {0}", organName));

                else if (DoseTypeTB.Text.Equals("Mean")) //if target is not null and structure is not empty then we can retriev our dose, we split into cases 
                {
                    DVHData dvhData = planSetup.GetDVHCumulativeData(target,
                                            DoseValuePresentation.Absolute,
                                            VolumePresentation.AbsoluteCm3, 0.1);
                    float dose = (float)dvhData.MeanDose.Dose;//compute mean dose            
                    float dosePerFraction = dose / numberOfFractions;
                    bed = dose * (1.0f + dosePerFraction / alphaBetaRatio);
                    eqd2 = dose * ((dosePerFraction + alphaBetaRatio) / (2.0f + alphaBetaRatio));


                    AbsoluteDoseTB.Text=dose.ToString("0.0");
                    BEDTB.Text = bed.ToString("0.0");
                    EQD2TB.Text = eqd2.ToString("0.0");

                    /*this is for debugging only
                    MessageBox.Show(String.Format("Organ={0}, MeanDose={1}," +
                        "Dose Per Fraction={2}, Number of Fractions= {3}", organName, dose,dosePerFraction, numberOfFractions));*/
                }

                else if (DoseTypeTB.Text.Equals("Max")) //if target is not null and structure is not empty then we can retriev our dose, we split into cases 
                {
                    DVHData dvhData = planSetup.GetDVHCumulativeData(target,
                                            DoseValuePresentation.Absolute,
                                            VolumePresentation.AbsoluteCm3, 0.1);
                    float dose = (float)dvhData.MaxDose.Dose;//compute mean dose            
                    float dosePerFraction = dose / numberOfFractions;
                    bed = dose * (1.0f + dosePerFraction / alphaBetaRatio);
                    eqd2 = dose * ((dosePerFraction + alphaBetaRatio) / (2.0f + alphaBetaRatio));

                    AbsoluteDoseTB.Text = dose.ToString("0.0");
                    BEDTB.Text = bed.ToString("0.0");
                    EQD2TB.Text = eqd2.ToString("0.0");

                    /*this is for debugging only
                    MessageBox.Show(String.Format("Organ={0}, MeanDose={1}," +
                        "Dose Per Fraction={2}, Number of Fractions= {3}", organName, dose,dosePerFraction, numberOfFractions));*/
                }
                else if (DoseTypeTB.Text.Equals("Min")) //if target is not null and structure is not empty then we can retriev our dose, we split into cases 
                {
                    DVHData dvhData = planSetup.GetDVHCumulativeData(target,
                                            DoseValuePresentation.Absolute,
                                            VolumePresentation.AbsoluteCm3, 0.1);
                    float dose = (float)dvhData.MinDose.Dose;//compute mean dose            
                    float dosePerFraction = dose / numberOfFractions;
                    bed = dose * (1.0f + dosePerFraction / alphaBetaRatio);
                    eqd2 = dose * ((dosePerFraction + alphaBetaRatio) / (2.0f + alphaBetaRatio));

                    AbsoluteDoseTB.Text = dose.ToString("0.0");
                    BEDTB.Text = bed.ToString("0.0");
                    EQD2TB.Text = eqd2.ToString("0.0");

                    /*this is for debugging only
                    MessageBox.Show(String.Format("Organ={0}, MeanDose={1}," +
                        "Dose Per Fraction={2}, Number of Fractions= {3}", organName, dose,dosePerFraction, numberOfFractions));*/
                }
                else if (DoseTypeTB.Text.EndsWith("cc") && DoseTypeTB.Text.StartsWith("D"))
                {

                    string volString = FindStringBetween("D", "cc", DoseTypeTB.Text);


                    bool success = float.TryParse(volString, out float volFloat);

                    if (success)
                    {
                        //MessageBox.Show(String.Format("Debug: Sucessfully converted absolute volume to float, value is {0}cc", volFloat));
                    }
                    else
                    {
                        MessageBox.Show("Incorrect Dose Format, your dose must have format \"DNcc\" where N is the volume in cc. N must be a number");
                    }
                    //at this point, we have sucessfully converted the volume in cc to a float from the entered string. Let's find the dose that corresponds to this
                    DVHData dvhData = planSetup.GetDVHCumulativeData(target, //extract absolute dvh
                        DoseValuePresentation.Absolute,
                        VolumePresentation.AbsoluteCm3, 0.1);
                    float dose = (float)DoseAtVolume(dvhData, volFloat).Dose;
                    float dosePerFraction = dose / numberOfFractions;
                    bed = dose * (1.0f + dosePerFraction / alphaBetaRatio);
                    eqd2 = dose * ((dosePerFraction + alphaBetaRatio) / (2.0f + alphaBetaRatio));

                    AbsoluteDoseTB.Text = dose.ToString("0.0");
                    BEDTB.Text = bed.ToString("0.0");
                    EQD2TB.Text = eqd2.ToString("0.0");


                    //MessageBox.Show("Debug: volume is " + volString+"cc");
                }
                else if (DoseTypeTB.Text.EndsWith("%") && DoseTypeTB.Text.StartsWith("D"))
                {
                    //since this is the % case, we convert from % volume to cc by multiplying through. This is built into the Dose at volume function if you extract your dvh properly




                    string volString = FindStringBetween("D", "%", DoseTypeTB.Text);
                    
                    bool success = float.TryParse(volString, out float volFloat);

                    if (success)
                    {
                        //MessageBox.Show(String.Format("Debug: Sucessfully converted volume % to float, value is {0}%", volFloat));
                    }
                    else
                    {
                        MessageBox.Show("Incorrect Dose Format, your dose must have format \"DN%\" where N is the volume in %. N must be a number");
                    }

                    DVHData dvhData = planSetup.GetDVHCumulativeData(target,
                        DoseValuePresentation.Absolute,
                        VolumePresentation.Relative,0.1);

                    float dose = (float)DoseAtVolume(dvhData, volFloat).Dose;
                    float dosePerFraction = dose / numberOfFractions;
                    bed = dose * (1.0f + dosePerFraction / alphaBetaRatio);
                    eqd2 = dose * ((dosePerFraction + alphaBetaRatio) / (2.0f + alphaBetaRatio));

                    AbsoluteDoseTB.Text = dose.ToString("0.0");
                    BEDTB.Text = bed.ToString("0.0");
                    EQD2TB.Text = eqd2.ToString("0.0");


                    //MessageBox.Show("Debug: volume is "+volString+"%");
                    //MessageBox.Show("Still working on retriving % doses");
                }
                else
                {
                    MessageBox.Show("Invalid Dose type, please enter either enter either \"Mean\", \"Max\", \"Min\", \"DNcc\" or \"DN%\" into the Input Dose box.");
                }
            }


        }
        private static string FindStringBetween(string firstString, string secondString, string wholeString)
        {
            int pFrom = wholeString.IndexOf(firstString) + firstString.Length;
            int pTo = wholeString.LastIndexOf(secondString);

            return wholeString.Substring(pFrom, pTo - pFrom);
        }

        private static DoseValue DoseAtVolume(DVHData dvhData, double volume) //volume lookup
        {
            if (dvhData == null || dvhData.CurveData.Count() == 0) //if the DVH is empty then return undefined dose
                return DoseValue.UndefinedDose();
            double absVolume = dvhData.CurveData[0].VolumeUnit == "%" ? volume * dvhData.Volume * 0.01 : volume;
            if (volume < 0.0 || absVolume > dvhData.Volume)
                return DoseValue.UndefinedDose();

            DVHPoint[] hist = dvhData.CurveData;
            for (int i = 0; i < hist.Length; i++)
            {
                if (hist[i].Volume < volume)
                    return hist[i].DoseValue;
            }
            return DoseValue.UndefinedDose();
        }

    }
}