o F,h&@slddlZddlZddlmZddlmZddlmZddl mm Z GdddeZ Gdddej ZdS)N)Metric)cKDTreecsFeZdZdZd deffdd Zdejdejfd d Zd d Z Z S) Sensitivityz Computes Sensitivity (Recall) = TP / (TP + FN) Args: threshold (float, optional): Threshold for binarizing predictions. Default is 0.5. ?F thresholdcsPtj|d||_|jdtjdtjddd|jdtjdtjddddS)N)dist_sync_on_steptrue_positivesr)dtypesum)defaultdist_reduce_fxfalse_negatives)super__init__r add_statetorchtensorfloat32)selfrr __class__,/home/xmoravc/AI_project/utils/my_metrics.pyrs zSensitivity.__init__predstargetcCsb||jk}|}t|dk|dk@}t|dk|dk@}|j|7_|j|7_dS)z Update metric states with new batch of predictions and targets. Args: preds (torch.Tensor): Model predictions (logits or probabilities). target (torch.Tensor): Ground truth labels (binary). rN)rfloatrr rr )rrrtpfnrrrupdates  zSensitivity.updatecCs|j|j|jdS)zN Computes final sensitivity score across accumulated batches. g|=)rr )rrrrcompute.szSensitivity.compute)rF) __name__ __module__ __qualname____doc__rrrTensorrr __classcell__rrrrr s  rcs:eZdZdZd fdd ZddZdd Zd d ZZS) HausdorffDistancez Computes the Hausdorff Distance (HD) and the 95th percentile Hausdorff Distance (HD95) for 3D binary segmentation masks in PyTorch. _rcs&t||_||_td|_dS)z Initializes the HausdorffDistance class. :param percentile: Percentile for the HD computation (default is 95). :param threshold: Binarization threshold for segmentation masks. cpuN)rr percentilerrdevice)rr*rrrrr;s zHausdorffDistance.__init__cCs ||_|S)zD Ensures the metric is moved to the correct device. r+)rr+rrrtoGszHausdorffDistance.toc Cs||j}||j}||jk}||jk}tj|dd}tj|dd}|dks4|dkrHtjtd|jdtjtd|jdfS|||}|||}t ||g}t |}t ||j d} || fS)aE Computes the Hausdorff Distance (HD) and the 95th percentile Hausdorff Distance (HD95). :param seg1: Binary 3D PyTorch tensor (segmentation mask 1) :param seg2: Binary 3D PyTorch tensor (segmentation mask 2) :return: (HD, HD95) as a tuple of torch.Tensors (floats or NaN if invalid) F)as_tuplernanr,gY@) r-r+rrrnonzeronumelr _min_distancecatmaxquantiler*) rseg1seg2Z seg1_pointsZ seg2_pointsZ dists_1_to_2Z dists_2_to_1Z all_distshdhd95rrrforwardps  (   zHausdorffDistance.forwardcCs`|jddks|jddkrtjtdg|jdStj||dd}tj|dd\}}|S)a! Computes the minimum distance from each point in points1 to the closest point in points2. :param points1: Tensor of shape (N, 3) - voxel positions :param points2: Tensor of shape (M, 3) - voxel positions :return: Tensor of distances of shape (N,) rinfr,)pr)dim)shaperrrr+cdistmin)rZpoints1Zpoints2distsZ min_dists_rrrr2s zHausdorffDistance._min_distance)r(r) r!r"r#r$rr-r:r2r&rrrrr'5s  )&r')numpynpr torchmetricsr scipy.spatialrtorch.nnnnZtorch.nn.functional functionalFrModuler'rrrrs   ,