Methods, systems and computer readable storage media storing instructions for image-guided treatment planning and assessment

1. A method for determining a treatment response of a treatment for a region of interest, comprising:

registering at least a first image data and a second image data of a region of interest to generate a registered image, the first image data corresponding to a first time point, the second image data corresponding to a second time point, the first time point and the second time point being different with respect to a session and/or treatment stage;

determining a treatment response, by a computer having a processor and a memory, based on at least the registered image, treatment plan information associated with the treatment, or a combination thereof, the determining including:

processing the registered image to determine a difference image;

processing at least the difference image to determine one or more signal regions, each signal region corresponding to one or more voxels that are grouped together based on location and signal value, the one or more signal regions including a positive signal region and/or a negative signal region; and

determining a treatment response for each signal region based on the signal value associated with the each region, the positive signal region indicating a treatment-resistant region and the negative signal region indicating a respondent region; and

outputting treatment response results, the treatment response results including the treatment response.

2. The method according to claim 1 , wherein the processing at least the difference image to determine one or more signal regions include:

determining an estimate signal region of one or more voxels based on the signal value and location with respect to a treatment site; and

processing each voxel of the estimate signal region to determine the signal region, the processing including comparing each voxel to a threshold range associated with the signal value and estimated signal region.

3. The method according to claim 2 , wherein the processing each voxel of the estimate signal region to determine the signal region includes:

determining to include a voxel in the signal region when a signal value of the voxel is within the threshold range and is consistent with a signal value of one or more voxels neighboring the each voxel, and

determining to exclude a voxel from the signal region when a value of the voxel is outside the threshold, the value of the voxel is inconsistent with a signal value of one or more voxels neighboring the each voxel, or a combination thereof.

4. The method according to claim 1 , wherein the treatment stage includes pre-treatment, pre-treatment planning, post-treatment planning, post-treatment, or a combination thereof.

5. The method according to claim 1 , wherein the image data for each time point includes PET, MR, CT, SPECT-CT, or a combination thereof.

6. The method according to claim 1 , further comprising:

correlating the treatment response for each region with the treatment plan information;

wherein the treatment response results includes the registered image with the one or more signal regions with the associated treatment response and the treatment plan information,

wherein the treatment plan information includes planning dose, planning contours, delivered dose, or a combination thereof.

7. The method according to claim 1 , further comprising:

determining a dose for each image; and

correlating the dose with the one or more signal regions.

8. The method according to claim 1 , further comprising:

determining a dose for each image; and

determining a sum distribution file based on the registered image.

9. The method according to claim 1 , wherein the registration is deformable registration.

10. A method for determining a treatment response of a treatment for a region of interest, comprising:

registering at least a first image data and a second image data of a region of interest to generate a registered image, the first image data corresponding to a first time point, the second image data corresponding to a second time point, the first time point and the second time point being different with respect to a session and/or treatment stage;

determining a dose for each of voxel of the first image data and the second image data;

determining a treatment response, by a computer having a processor and a memory, based on at least the registered image, the determining including;

determining a sum dose distribution using the registered image by adding the dose associated with each voxel of the first image data to each corresponding voxel of the second image data;

classifying each voxel of the registered image into one or more signal regions based on signal value and/or location;

correlating the sum dose distribution to the one or more signal regions; and

determining the treatment response based on the sum dose distribution; and

outputting treatment response results that include the treatment response.

11. The method according to claim 10 , further comprising:

generating a dose volume histogram from the dose, the sum dose distribution, or a combination thereof; and

wherein the treatment response results includes the dose volume histogram.

12. The method according to claim 10 ,

wherein the classifying each voxel of the registered image includes:

processing the registered image to determine a difference image;

processing at least the difference image to determine the one or more signal regions, each signal region corresponding to one or more voxels that are grouped together based on the location and the signal value, the one or more signal regions including a positive signal region and/or a negative signal region; and

determining the treatment response associated with the signal value of the one or more signal regions, the positive signal region indicating a treatment-resistant region and the negative signal region indicating a respondent region.

13. The method according to claim 12 , wherein the processing at least the difference image to determine one or more signal regions include:

determining an estimate signal region of one or more voxels based on the signal value and location with respect to a treatment site; and

processing each voxel of the estimate signal region to determine the signal region, the processing including comparing each voxel to a threshold range associated with the signal value and estimated signal region.

14. The method according to claim 13 , wherein the processing each voxel of the estimate signal region to determine the signal region includes:

determining to include a voxel in the signal region when a signal value of the voxel is within the threshold range and is consistent with a signal value of one or more voxels neighboring the each voxel, and

determining to exclude a voxel from the signal region when a value of the voxel is outside the threshold, the value of the voxel is inconsistent with a signal value of one or more voxels neighboring the each voxel, or a combination thereof.

15. The method according to claim 10 , wherein the treatment stage includes pre-treatment, pre-treatment planning, post-treatment planning, post-treatment, or a combination thereof.

16. The method according to claim 10 , wherein the image data for each time point includes PET, MR, CT, SPECT-CT, or a combination thereof.

17. The method according to claim 10 , wherein the registration is deformable registration.

18. The method according to claim 10 , wherein the sum dose distribution is determined from PET data.

19. A non-transitory computer-readable storage medium storing instructions for determining a treatment response of a treatment for a region of interest, the instructions comprising:

registering at least a first image data and a second image data of a region of interest to generate a registered image, the first image data corresponding to a first time point, the second image data corresponding to a second time point, the first time point and the second time point being different with respect to a session and/or treatment stage;

determining a dose for each of voxel of the first image data and the second image data;

determining a treatment response, by a computer having a processor and a memory, based on at least the registered image, the determining including;

determining a sum dose distribution using the registered image by adding the dose associated with each voxel of the first image data to each corresponding voxel of the second image data;

classifying each voxel of the registered image into one or more signal regions based on signal value and/or location;

correlating the sum dose distribution to the one or more signal regions; and

determining the treatment response based on the sum dose distribution; and

outputting treatment response results that include the treatment response.

20. The computer-readable storage medium according to claim 19 , the instructions further comprising:

generating a dose volume histogram from the dose, the sum dose distribution, or a combination thereof; and

wherein the treatment response results includes the dose volume histogram.