-
Notifications
You must be signed in to change notification settings - Fork 387
/
rase.py
100 lines (80 loc) · 3.92 KB
/
rase.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
# Copyright The PyTorch Lightning team.
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
# http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.
from typing import Tuple
import torch
from torch import Tensor
from torchmetrics.functional.image.helper import _uniform_filter
from torchmetrics.functional.image.rmse_sw import _rmse_sw_compute, _rmse_sw_update
def _rase_update(
preds: Tensor, target: Tensor, window_size: int, rmse_map: Tensor, target_sum: Tensor, total_images: Tensor
) -> Tuple[Tensor, Tensor, Tensor]:
"""Calculate the sum of RMSE map values for the batch of examples and update intermediate states.
Args:
preds: Deformed image
target: Ground truth image
window_size: Sliding window used for RMSE calculation
rmse_map: Sum of RMSE map values over all examples
target_sum: target...
total_images: Total number of images
Return:
Intermediate state of RMSE map
Updated total number of already processed images
"""
_, rmse_map, total_images = _rmse_sw_update(
preds, target, window_size, rmse_val_sum=None, rmse_map=rmse_map, total_images=total_images
)
target_sum += torch.sum(_uniform_filter(target, window_size) / (window_size**2), dim=0)
return rmse_map, target_sum, total_images
def _rase_compute(rmse_map: Tensor, target_sum: Tensor, total_images: Tensor, window_size: int) -> Tensor:
"""Compute RASE.
Args:
rmse_map: Sum of RMSE map values over all examples
target_sum: target...
total_images: Total number of images.
window_size: Sliding window used for rmse calculation
Return:
Relative Average Spectral Error (RASE)
"""
_, rmse_map = _rmse_sw_compute(rmse_val_sum=None, rmse_map=rmse_map, total_images=total_images)
target_mean = target_sum / total_images
target_mean = target_mean.mean(0) # mean over image channels
rase_map = 100 / target_mean * torch.sqrt(torch.mean(rmse_map**2, 0))
crop_slide = round(window_size / 2)
return torch.mean(rase_map[crop_slide:-crop_slide, crop_slide:-crop_slide])
def relative_average_spectral_error(preds: Tensor, target: Tensor, window_size: int = 8) -> Tensor:
"""Compute Relative Average Spectral Error (RASE) (RelativeAverageSpectralError_).
Args:
preds: Deformed image
target: Ground truth image
window_size: Sliding window used for rmse calculation
Return:
Relative Average Spectral Error (RASE)
Example:
>>> from torchmetrics.functional.image import relative_average_spectral_error
>>> g = torch.manual_seed(22)
>>> preds = torch.rand(4, 3, 16, 16)
>>> target = torch.rand(4, 3, 16, 16)
>>> relative_average_spectral_error(preds, target)
tensor(5114.6641)
Raises:
ValueError: If ``window_size`` is not a positive integer.
"""
if not isinstance(window_size, int) or isinstance(window_size, int) and window_size < 1:
raise ValueError("Argument `window_size` is expected to be a positive integer.")
img_shape = target.shape[1:] # [num_channels, width, height]
rmse_map = torch.zeros(img_shape, dtype=target.dtype, device=target.device)
target_sum = torch.zeros(img_shape, dtype=target.dtype, device=target.device)
total_images = torch.tensor(0.0, device=target.device)
rmse_map, target_sum, total_images = _rase_update(preds, target, window_size, rmse_map, target_sum, total_images)
return _rase_compute(rmse_map, target_sum, total_images, window_size)