How to compare contours?

Kshitij Singh answered . 2024-12-21 16:44:49

1. Create a Grid: Use meshgrid to create a grid of x and y coordinates corresponding to your dose_data matrix.
2. Generate Binary Masks: Create binary masks where the pixels within the 11 Gy isodose line and within the tumour circle are 1, and all other pixels are 0.
3. Calculate the Intersection: Find the intersection of the binary masks to see which pixels are within both the isodose line and the circle.
4. Calculate the DSC: The DSC is calculated as 2 * (|A ∩ B|) / (|A| + |B|), where A and B are the binary masks of the two regions, and |A ∩ B| is the count of pixels in the intersection.
5. Calculate the Percentage Area Covered: The percentage of the circle covered by the 11 Gy dose contour is the number of pixels in the intersection divided by the number of pixels in the circle mask.

% Assume col_scale and row_scale are vectors that scale the matrix indices to your physical coordinates.

[x, y] = meshgrid(col_scale, row_scale);

% Binary mask for the 11 Gy isodose line
isodose_11Gy_mask = dose_data >= 11;

% Binary mask for the circle (tumour)
% This is a logical operation to create a binary mask of the circle.
circle_mask = ((x.^2 + y.^2) <= 5^2);

% Find the intersection of the isodose line and the circle
intersection_mask = isodose_11Gy_mask & circle_mask;

% Calculate Dice Similarity Coefficient (DSC)
DSC = 2 * nnz(intersection_mask) / (nnz(isodose_11Gy_mask) + nnz(circle_mask));

% Calculate the percentage area of circle covered by the 11 Gy dose contour
percent_coverage = nnz(intersection_mask) / nnz(circle_mask) * 100;

% Display the results
disp(['DSC: ', num2str(DSC)]);
disp(['Percentage of Circle Covered by 11 Gy Contour: ', num2str(percent_coverage), '%']);

Not satisfied with the answer ?? ASK NOW