i$ddlmZddlmZmZmZmZddlmZddl m Z ddl m Z m Z ddlmZddlmZddlmZdd lmZmZdd lmZmZesd gZesd d gZGd d eZy))Sequence)AnyListOptionalUnion)Tensor)Literal)!_spatial_distortion_index_compute _spatial_distortion_index_update)Metric)rank_zero_warn) dim_zero_cat)_MATPLOTLIB_AVAILABLE_TORCHVISION_AVAILABLE)_AX_TYPE_PLOT_OUT_TYPEzSpatialDistortionIndex.plotSpatialDistortionIndexc ,eZdZUdZdZeed<dZeed<dZeed<dZ e ed<d Z e ed <e e ed <e e ed <e e ed <e e ed< ddedededdeddf fd Zd e deee fddfdZde fdZ ddeee ee fdeedefdZxZS)ra Compute Spatial Distortion Index (SpatialDistortionIndex_) also now as D_s. The metric is used to compare the spatial distortion between two images. A value of 0 indicates no distortion (optimal value) and corresponds to the case where the high resolution panchromatic image is equal to the low resolution panchromatic image. The metric is defined as: .. math:: D_s = \\sqrt[q]{\frac{1}{L}\\sum_{l=1}^L|Q(\\hat{G_l}, P) - Q(\tilde{G}, \tilde{P})|^q} where :math:`Q` is the universal image quality index (see this :class:`~torchmetrics.image.UniversalImageQualityIndex` for more info), :math:`\\hat{G_l}` is the l-th band of the high resolution multispectral image, :math:`\tilde{G}` is the high resolution panchromatic image, :math:`P` is the high resolution panchromatic image, :math:`\tilde{P}` is the low resolution panchromatic image, :math:`L` is the number of bands and :math:`q` is the order of the norm applied on the difference. As input to ``forward`` and ``update`` the metric accepts the following input - ``preds`` (:class:`~torch.Tensor`): High resolution multispectral image of shape ``(N,C,H,W)``. - ``target`` (:class:`~Dict`): A dictionary containing the following keys: - ``ms`` (:class:`~torch.Tensor`): Low resolution multispectral image of shape ``(N,C,H',W')``. - ``pan`` (:class:`~torch.Tensor`): High resolution panchromatic image of shape ``(N,C,H,W)``. - ``pan_lr`` (:class:`~torch.Tensor`): Low resolution panchromatic image of shape ``(N,C,H',W')``. where H and W must be multiple of H' and W'. As output of `forward` and `compute` the metric returns the following output - ``sdi`` (:class:`~torch.Tensor`): if ``reduction!='none'`` returns float scalar tensor with average SDI value over sample else returns tensor of shape ``(N,)`` with SDI values per sample Args: norm_order: Order of the norm applied on the difference. window_size: Window size of the filter applied to degrade the high resolution panchromatic image. reduction: a method to reduce metric score over labels. - ``'elementwise_mean'``: takes the mean (default) - ``'sum'``: takes the sum - ``'none'``: no reduction will be applied kwargs: Additional keyword arguments, see :ref:`Metric kwargs` for more info. Example: >>> from torch import rand >>> from torchmetrics.image import SpatialDistortionIndex >>> preds = rand([16, 3, 32, 32]) >>> target = { ... 'ms': rand([16, 3, 16, 16]), ... 'pan': rand([16, 3, 32, 32]), ... } >>> sdi = SpatialDistortionIndex() >>> sdi(preds, target) tensor(0.0090) Fhigher_is_betterTis_differentiablefull_state_updategplot_lower_boundg?plot_upper_boundpredsmspanpan_lr norm_order window_size reductionelementwise_meansumnonekwargsreturnNc t|di|tdt|tr|dkrt d|d||_t|tr|dkrt d|d||_d}||vrt d|d|||_|jd gd |jd gd |jd gd |jdgd y)NzMetric `SpatialDistortionIndex` will save all targets and predictions in buffer. For large datasets this may lead to large memory footprint.rz@Expected `norm_order` to be a positive integer. Got norm_order: .zBExpected `window_size` to be a positive integer. Got window_size: r!z(Expected argument `reduction` be one of z but got rcat)defaultdist_reduce_fxrrr) super__init__r isinstanceint ValueErrorrrr add_state)selfrrr r%allowed_reductions __class__s k/Volumes/fast/ai/experiments/voice-extract-mac/.venv/lib/python3.12/site-packages/torchmetrics/image/d_s.pyr.zSpatialDistortionIndex.__init__fs "6" * *c*jAo_`j_kklmn n$+s+{a/?abmannopq q&@ . .GHZG[[dendopq q" w5A tR> ub? xEBtargetcd|vrtd|jdd|vrtd|jd|d}|d}|jd}t||||\}}}}|jj ||j j |d|jj |dd|vr|jj |dyy)aUpdate state with preds and target. Args: preds: High resolution multispectral image. target: A dictionary containing the following keys: - ``'ms'``: low resolution multispectral image. - ``'pan'``: high resolution panchromatic image. - ``'pan_lr'``: (optional) low resolution panchromatic image. Raises: ValueError: If ``target`` doesn't have ``ms`` and ``pan``. rz0Expected `target` to have key `ms`. Got target: r(rz1Expected `target` to have key `pan`. Got target: rN) r1keysgetr rappendrrr)r3rr8rrrs r6updatezSpatialDistortionIndex.updates v OPVP[P[P]^_`a a  PQWQ\Q\Q^P__`ab b D\UmH%!A%SRX!Yr3 %  vd|$ u & v  KK  vh/ 0 r7c tt|j}t|j}t|j}t |j dkDrt|j nd}||d}|j |d|init|||||j|j|jS)z-Compute and returns spatial distortion index.rN)rrr) rrrrlenrr=r rrr )r3rrrrr8s r6computezSpatialDistortionIndex.computesTZZ( $'' "488$.1$++.>.Bdkk*3' F,>x(BG0 2sFDOOT5E5Et~~  r7valaxc&|j||S)aPlot a single or multiple values from the metric. Args: val: Either a single result from calling `metric.forward` or `metric.compute` or a list of these results. If no value is provided, will automatically call `metric.compute` and plot that result. ax: An matplotlib axis object. If provided will add plot to that axis Returns: Figure and Axes object Raises: ModuleNotFoundError: If `matplotlib` is not installed .. plot:: :scale: 75 >>> # Example plotting a single value >>> from torch import rand >>> from torchmetrics.image import SpatialDistortionIndex >>> preds = rand([16, 3, 32, 32]) >>> target = { ... 'ms': rand([16, 3, 16, 16]), ... 'pan': rand([16, 3, 32, 32]), ... } >>> metric = SpatialDistortionIndex() >>> metric.update(preds, target) >>> fig_, ax_ = metric.plot() .. plot:: :scale: 75 >>> # Example plotting multiple values >>> from torch import rand >>> from torchmetrics.image import SpatialDistortionIndex >>> preds = rand([16, 3, 32, 32]) >>> target = { ... 'ms': rand([16, 3, 16, 16]), ... 'pan': rand([16, 3, 32, 32]), ... } >>> metric = SpatialDistortionIndex() >>> values = [ ] >>> for _ in range(10): ... values.append(metric(preds, target)) >>> fig_, ax_ = metric.plot(values) )_plot)r3rArBs r6plotzSpatialDistortionIndex.plotsdzz#r""r7)r")NN)__name__ __module__ __qualname____doc__rbool__annotations__rrrfloatrrrr0r rr.dictstrr=r@rrrrrrE __classcell__)r5s@r6rr#s#5n#d""t"#t#!e!!e! < V  f L@R CCC<= C  C  C:1F1Df,=1$1<    _c2#E&(6*:":;<2#IQRZI[2# 2#r7N)collections.abcrtypingrrrrtorchrtyping_extensionsr !torchmetrics.functional.image.d_sr r torchmetrics.metricr torchmetrics.utilitiesr torchmetrics.utilities.datartorchmetrics.utilities.importsrrtorchmetrics.utilities.plotrr__doctest_skip__rr,r7r6r]sQ%--%q&14X@5602OP|#V|#r7