monai-weekly 1.5.dev2506__py3-none-any.whl → 1.5.dev2507__py3-none-any.whl

tests/transforms/test_sobel_gradientd.py

@@ -0,0 +1,212 @@
+# Copyright (c) MONAI Consortium
+# 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 __future__ import annotations
+
+import unittest
+
+import torch
+from parameterized import parameterized
+
+from monai.transforms import SobelGradientsd
+from tests.test_utils import assert_allclose
+
+IMAGE = torch.zeros(1, 16, 16, dtype=torch.float32)
+IMAGE[0, 8, :] = 1
+
+# Output with reflect padding
+OUTPUT_3x3 = torch.zeros(2, 16, 16, dtype=torch.float32)
+OUTPUT_3x3[0, 7, :] = 0.5
+OUTPUT_3x3[0, 9, :] = -0.5
+
+# Output with zero padding
+OUTPUT_3x3_ZERO_PAD = OUTPUT_3x3.clone()
+OUTPUT_3x3_ZERO_PAD[1, 7, 0] = OUTPUT_3x3_ZERO_PAD[1, 9, 0] = 0.125
+OUTPUT_3x3_ZERO_PAD[1, 8, 0] = 0.25
+OUTPUT_3x3_ZERO_PAD[1, 7, -1] = OUTPUT_3x3_ZERO_PAD[1, 9, -1] = -0.125
+OUTPUT_3x3_ZERO_PAD[1, 8, -1] = -0.25
+OUTPUT_3x3_ZERO_PAD[0, 7, 0] = OUTPUT_3x3_ZERO_PAD[0, 7, -1] = 3.0 / 8.0
+OUTPUT_3x3_ZERO_PAD[0, 9, 0] = OUTPUT_3x3_ZERO_PAD[0, 9, -1] = -3.0 / 8.0
+
+TEST_CASE_0 = [{"image": IMAGE}, {"keys": "image", "kernel_size": 3, "dtype": torch.float32}, {"image": OUTPUT_3x3}]
+TEST_CASE_1 = [{"image": IMAGE}, {"keys": "image", "kernel_size": 3, "dtype": torch.float64}, {"image": OUTPUT_3x3}]
+TEST_CASE_2 = [
+    {"image": IMAGE},
+    {"keys": "image", "kernel_size": 3, "dtype": torch.float32, "new_key_prefix": "sobel_"},
+    {"sobel_image": OUTPUT_3x3},
+]
+TEST_CASE_3 = [
+    {"image": IMAGE},
+    {"keys": "image", "kernel_size": 3, "spatial_axes": 0, "dtype": torch.float32},
+    {"image": OUTPUT_3x3[0][None, ...]},
+]
+TEST_CASE_4 = [
+    {"image": IMAGE},
+    {"keys": "image", "kernel_size": 3, "spatial_axes": 1, "dtype": torch.float32},
+    {"image": OUTPUT_3x3[1][None, ...]},
+]
+TEST_CASE_5 = [
+    {"image": IMAGE},
+    {"keys": "image", "kernel_size": 3, "spatial_axes": [1], "dtype": torch.float32},
+    {"image": OUTPUT_3x3[1][None, ...]},
+]
+TEST_CASE_6 = [
+    {"image": IMAGE},
+    {"keys": "image", "kernel_size": 3, "spatial_axes": [0, 1], "normalize_kernels": True, "dtype": torch.float32},
+    {"image": OUTPUT_3x3},
+]
+TEST_CASE_7 = [
+    {"image": IMAGE},
+    {"keys": "image", "kernel_size": 3, "spatial_axes": (0, 1), "padding_mode": "reflect", "dtype": torch.float32},
+    {"image": OUTPUT_3x3},
+]
+TEST_CASE_8 = [
+    {"image": IMAGE},
+    {"keys": "image", "kernel_size": 3, "spatial_axes": (0, 1), "padding_mode": "zeros", "dtype": torch.float32},
+    {"image": OUTPUT_3x3_ZERO_PAD},
+]
+TEST_CASE_9 = [  # Non-normalized kernels
+    {"image": IMAGE},
+    {"keys": "image", "kernel_size": 3, "spatial_axes": (0, 1), "normalize_kernels": False, "dtype": torch.float32},
+    {"image": OUTPUT_3x3 * 8.0},
+]
+TEST_CASE_10 = [  # Normalized gradients and normalized kernels
+    {"image": IMAGE},
+    {
+        "keys": "image",
+        "kernel_size": 3,
+        "spatial_axes": (0, 1),
+        "normalize_kernels": True,
+        "normalize_gradients": True,
+        "dtype": torch.float32,
+    },
+    {"image": torch.cat([OUTPUT_3x3[0:1] + 0.5, OUTPUT_3x3[1:2]])},
+]
+TEST_CASE_11 = [  # Normalized gradients but non-normalized kernels
+    {"image": IMAGE},
+    {
+        "keys": "image",
+        "kernel_size": 3,
+        "spatial_axes": (0, 1),
+        "normalize_kernels": False,
+        "normalize_gradients": True,
+        "dtype": torch.float32,
+    },
+    {"image": torch.cat([OUTPUT_3x3[0:1] + 0.5, OUTPUT_3x3[1:2]])},
+]
+
+TEST_CASE_KERNEL_0 = [
+    {"keys": "image", "kernel_size": 3, "dtype": torch.float64},
+    (torch.tensor([-0.5, 0.0, 0.5], dtype=torch.float64), torch.tensor([0.25, 0.5, 0.25], dtype=torch.float64)),
+]
+TEST_CASE_KERNEL_1 = [
+    {"keys": "image", "kernel_size": 5, "dtype": torch.float64},
+    (
+        torch.tensor([-0.1250, -0.2500, 0.0000, 0.2500, 0.1250], dtype=torch.float64),
+        torch.tensor([0.0625, 0.2500, 0.3750, 0.2500, 0.0625], dtype=torch.float64),
+    ),
+]
+TEST_CASE_KERNEL_2 = [
+    {"keys": "image", "kernel_size": 7, "dtype": torch.float64},
+    (
+        torch.tensor([-0.03125, -0.125, -0.15625, 0.0, 0.15625, 0.125, 0.03125], dtype=torch.float64),
+        torch.tensor([0.015625, 0.09375, 0.234375, 0.3125, 0.234375, 0.09375, 0.015625], dtype=torch.float64),
+    ),
+]
+TEST_CASE_KERNEL_NON_NORMALIZED_0 = [
+    {"keys": "image", "kernel_size": 3, "normalize_kernels": False, "dtype": torch.float64},
+    (torch.tensor([-1.0, 0.0, 1.0], dtype=torch.float64), torch.tensor([1.0, 2.0, 1.0], dtype=torch.float64)),
+]
+TEST_CASE_KERNEL_NON_NORMALIZED_1 = [
+    {"keys": "image", "kernel_size": 5, "normalize_kernels": False, "dtype": torch.float64},
+    (
+        torch.tensor([-1.0, -2.0, 0.0, 2.0, 1.0], dtype=torch.float64),
+        torch.tensor([1.0, 4.0, 6.0, 4.0, 1.0], dtype=torch.float64),
+    ),
+]
+TEST_CASE_KERNEL_NON_NORMALIZED_2 = [
+    {"keys": "image", "kernel_size": 7, "normalize_kernels": False, "dtype": torch.float64},
+    (
+        torch.tensor([-1.0, -4.0, -5.0, 0.0, 5.0, 4.0, 1.0], dtype=torch.float64),
+        torch.tensor([1.0, 6.0, 15.0, 20.0, 15.0, 6.0, 1.0], dtype=torch.float64),
+    ),
+]
+TEST_CASE_ERROR_0 = [{"image": IMAGE}, {"keys": "image", "kernel_size": 1}]  # kernel size less than 3
+TEST_CASE_ERROR_1 = [{"image": IMAGE}, {"keys": "image", "kernel_size": 4}]  # even kernel size
+TEST_CASE_ERROR_2 = [{"image": IMAGE}, {"keys": "image", "spatial_axes": "horizontal"}]  # wrong type direction
+TEST_CASE_ERROR_3 = [{"image": IMAGE}, {"keys": "image", "spatial_axes": 3}]  # wrong direction
+TEST_CASE_ERROR_4 = [{"image": IMAGE}, {"keys": "image", "spatial_axes": [3]}]  # wrong direction in a list
+TEST_CASE_ERROR_5 = [
+    {"image": IMAGE},
+    {"keys": "image", "spatial_axes": [0, 4]},
+]  # correct and wrong direction in a list
+
+
+class SobelGradientTests(unittest.TestCase):
+    backend = None
+
+    @parameterized.expand(
+        [
+            TEST_CASE_0,
+            TEST_CASE_1,
+            TEST_CASE_2,
+            TEST_CASE_3,
+            TEST_CASE_4,
+            TEST_CASE_5,
+            TEST_CASE_6,
+            TEST_CASE_7,
+            TEST_CASE_8,
+            TEST_CASE_9,
+            TEST_CASE_10,
+            TEST_CASE_11,
+        ]
+    )
+    def test_sobel_gradients(self, image_dict, arguments, expected_grad):
+        sobel = SobelGradientsd(**arguments)
+        grad = sobel(image_dict)
+        key = "image" if "new_key_prefix" not in arguments else arguments["new_key_prefix"] + arguments["keys"]
+        assert_allclose(grad[key], expected_grad[key])
+
+    @parameterized.expand(
+        [
+            TEST_CASE_KERNEL_0,
+            TEST_CASE_KERNEL_1,
+            TEST_CASE_KERNEL_2,
+            TEST_CASE_KERNEL_NON_NORMALIZED_0,
+            TEST_CASE_KERNEL_NON_NORMALIZED_1,
+            TEST_CASE_KERNEL_NON_NORMALIZED_2,
+        ]
+    )
+    def test_sobel_kernels(self, arguments, expected_kernels):
+        sobel = SobelGradientsd(**arguments)
+        self.assertEqual(sobel.kernel_diff.dtype, expected_kernels[0].dtype)
+        self.assertEqual(sobel.kernel_smooth.dtype, expected_kernels[0].dtype)
+        assert_allclose(sobel.kernel_diff, expected_kernels[0])
+        assert_allclose(sobel.kernel_smooth, expected_kernels[1])
+
+    @parameterized.expand(
+        [
+            TEST_CASE_ERROR_0,
+            TEST_CASE_ERROR_1,
+            TEST_CASE_ERROR_2,
+            TEST_CASE_ERROR_3,
+            TEST_CASE_ERROR_4,
+            TEST_CASE_ERROR_5,
+        ]
+    )
+    def test_sobel_gradients_error(self, image_dict, arguments):
+        with self.assertRaises(ValueError):
+            sobel = SobelGradientsd(**arguments)
+            sobel(image_dict)
+
+
+if __name__ == "__main__":
+    unittest.main()

tests/transforms/test_spacing.py

@@ -0,0 +1,381 @@
+# Copyright (c) MONAI Consortium
+# 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 __future__ import annotations
+
+import unittest
+
+import numpy as np
+import torch
+from parameterized import parameterized
+
+from monai.config import USE_COMPILED
+from monai.data.meta_obj import set_track_meta
+from monai.data.meta_tensor import MetaTensor
+from monai.data.utils import affine_to_spacing
+from monai.transforms import Spacing
+from monai.utils import fall_back_tuple
+from tests.lazy_transforms_utils import test_resampler_lazy
+from tests.test_utils import TEST_DEVICES, TEST_NDARRAYS_ALL, assert_allclose, skip_if_quick
+
+TESTS: list[list] = []
+for device in TEST_DEVICES:
+    TESTS.append(
+        [
+            {"pixdim": (1.0, 1.5), "padding_mode": "zeros", "dtype": float},
+            torch.arange(4).reshape((1, 2, 2)) + 1.0,  # data
+            torch.eye(4),
+            {},
+            torch.tensor([[[1.0, 1.0], [3.0, 2.0]]]),
+            *device,
+        ]
+    )
+    TESTS.append(
+        [
+            {"pixdim": 1.0, "padding_mode": "zeros", "dtype": float},
+            torch.ones((1, 2, 1, 2)),  # data
+            torch.eye(4),
+            {},
+            torch.tensor([[[[1.0, 1.0]], [[1.0, 1.0]]]]),
+            *device,
+        ]
+    )
+    TESTS.append(
+        [
+            {"pixdim": 2.0, "padding_mode": "zeros", "dtype": float},
+            torch.arange(4).reshape((1, 2, 2)) + 1.0,  # data
+            torch.eye(4),
+            {},
+            torch.tensor([[[1.0, 0.0], [0.0, 0.0]]]),
+            *device,
+        ]
+    )
+    TESTS.append(
+        [
+            {"pixdim": (1.0, 1.0, 1.0), "padding_mode": "zeros", "dtype": float},
+            torch.ones((1, 2, 1, 2)),  # data
+            torch.eye(4),
+            {},
+            torch.tensor([[[[1.0, 1.0]], [[1.0, 1.0]]]]),
+            *device,
+        ]
+    )
+    TESTS.append(
+        [
+            {"pixdim": (1.0, 0.2, 1.5), "diagonal": False, "padding_mode": "zeros", "align_corners": True},
+            torch.ones((1, 2, 1, 2)),  # data
+            torch.tensor([[2, 1, 0, 4], [-1, -3, 0, 5], [0, 0, 2.0, 5], [0, 0, 0, 1]]),
+            {},
+            (
+                torch.tensor([[[[0.75, 0.75]], [[0.75, 0.75]], [[0.75, 0.75]]]])
+                if USE_COMPILED
+                else torch.tensor([[[[0.95527864, 0.95527864]], [[1.0, 1.0]], [[1.0, 1.0]]]])
+            ),
+            *device,
+        ]
+    )
+    TESTS.append(
+        [
+            {"pixdim": (3.0, 1.0), "padding_mode": "zeros"},
+            torch.arange(24).reshape((2, 3, 4)),  # data
+            torch.as_tensor(np.diag([-3.0, 0.2, 1.5, 1])),
+            {},
+            torch.tensor([[[0, 0], [4, 0], [8, 0]], [[12, 0], [16, 0], [20, 0]]]),
+            *device,
+        ]
+    )
+    TESTS.append(
+        [
+            {"pixdim": (3.0, 1.0), "padding_mode": "zeros"},
+            torch.arange(24).reshape((2, 3, 4)),  # data
+            torch.eye(4),
+            {},
+            torch.tensor([[[0, 1, 2, 3], [0, 0, 0, 0]], [[12, 13, 14, 15], [0, 0, 0, 0]]]),
+            *device,
+        ]
+    )
+    TESTS.append(
+        [
+            {"pixdim": (1.0, 1.0), "align_corners": True},
+            torch.arange(24).reshape((2, 3, 4)),  # data
+            torch.eye(4),
+            {},
+            torch.tensor(
+                [[[0, 1, 2, 3], [4, 5, 6, 7], [8, 9, 10, 11]], [[12, 13, 14, 15], [16, 17, 18, 19], [20, 21, 22, 23]]]
+            ),
+            *device,
+        ]
+    )
+    TESTS.append(
+        [
+            {"pixdim": (4.0, 5.0, 6.0)},
+            torch.arange(24).reshape((1, 2, 3, 4)),  # data
+            torch.tensor([[-4, 0, 0, 4], [0, 5, 0, -5], [0, 0, 6, -6], [0, 0, 0, 1]]),
+            {},
+            torch.arange(24).reshape((1, 2, 3, 4)),  # data
+            *device,
+        ]
+    )
+    TESTS.append(
+        [
+            {"pixdim": (4.0, 5.0, 6.0), "diagonal": True},
+            torch.arange(24).reshape((1, 2, 3, 4)),  # data
+            torch.tensor([[-4, 0, 0, 4], [0, 5, 0, 0], [0, 0, 6, 0], [0, 0, 0, 1]]),
+            {},
+            torch.tensor(
+                [[[[12, 13, 14, 15], [16, 17, 18, 19], [20, 21, 22, 23]], [[0, 1, 2, 3], [4, 5, 6, 7], [8, 9, 10, 11]]]]
+            ),
+            *device,
+        ]
+    )
+    TESTS.append(
+        [
+            {"pixdim": (4.0, 5.0, 6.0), "padding_mode": "border", "diagonal": True},
+            torch.arange(24).reshape((1, 2, 3, 4)),  # data
+            torch.tensor([[-4, 0, 0, -4], [0, 5, 0, 0], [0, 0, 6, 0], [0, 0, 0, 1]]),
+            {},
+            torch.tensor(
+                [[[[12, 13, 14, 15], [16, 17, 18, 19], [20, 21, 22, 23]], [[0, 1, 2, 3], [4, 5, 6, 7], [8, 9, 10, 11]]]]
+            ),
+            *device,
+        ]
+    )
+    TESTS.append(
+        [
+            {"pixdim": (4.0, 5.0, 6.0), "padding_mode": "border", "diagonal": True},
+            torch.arange(24).reshape((1, 2, 3, 4)),  # data
+            torch.tensor([[-4, 0, 0, -4], [0, 5, 0, 0], [0, 0, 6, 0], [0, 0, 0, 1]]),
+            {"mode": "nearest"},
+            torch.tensor(
+                [[[[12, 13, 14, 15], [16, 17, 18, 19], [20, 21, 22, 23]], [[0, 1, 2, 3], [4, 5, 6, 7], [8, 9, 10, 11]]]]
+            ),
+            *device,
+        ]
+    )
+    TESTS.append(
+        [
+            {"pixdim": (1.9, 4.0), "padding_mode": "zeros", "diagonal": True},
+            torch.arange(24).reshape((1, 4, 6)),  # data
+            torch.tensor([[-4, 0, 0, -4], [0, 5, 0, 0], [0, 0, 6, 0], [0, 0, 0, 1]]),
+            {"mode": "nearest"},
+            torch.tensor(
+                [
+                    [
+                        [18.0, 19.0, 20.0, 20.0, 21.0, 22.0, 23.0],
+                        [18.0, 19.0, 20.0, 20.0, 21.0, 22.0, 23.0],
+                        [12.0, 13.0, 14.0, 14.0, 15.0, 16.0, 17.0],
+                        [12.0, 13.0, 14.0, 14.0, 15.0, 16.0, 17.0],
+                        [6.0, 7.0, 8.0, 8.0, 9.0, 10.0, 11.0],
+                        [6.0, 7.0, 8.0, 8.0, 9.0, 10.0, 11.0],
+                        [0.0, 1.0, 2.0, 2.0, 3.0, 4.0, 5.0],
+                    ]
+                ]
+            ),
+            *device,
+        ]
+    )
+    TESTS.append(
+        [
+            {"pixdim": (5.0, 3.0), "padding_mode": "border", "diagonal": True, "dtype": torch.float32},
+            torch.arange(24).reshape((1, 4, 6)),  # data
+            torch.tensor([[-4, 0, 0, 0], [0, 5, 0, 0], [0, 0, 6, 0], [0, 0, 0, 1]]),
+            {"mode": "bilinear"},
+            torch.tensor(
+                [
+                    [
+                        [18.0, 18.6, 19.2, 19.8, 20.400002, 21.0, 21.6, 22.2, 22.8],
+                        [10.5, 11.1, 11.700001, 12.299999, 12.900001, 13.5, 14.1, 14.700001, 15.3],
+                        [3.0, 3.6000001, 4.2000003, 4.8, 5.4000006, 6.0, 6.6000004, 7.200001, 7.8],
+                    ]
+                ]
+            ),
+            *device,
+        ]
+    )
+    TESTS.append(
+        [
+            {"pixdim": (5.0, 3.0), "padding_mode": "zeros", "diagonal": True, "dtype": torch.float32},
+            torch.arange(24).reshape((1, 4, 6)),  # data
+            torch.tensor([[-4, 0, 0, 0], [0, 5, 0, 0], [0, 0, 6, 0], [0, 0, 0, 1]]),
+            {"mode": "bilinear"},
+            torch.tensor(
+                [
+                    [
+                        [18.0000, 18.6000, 19.2000, 19.8000, 20.4000, 21.0000, 21.6000, 22.2000, 22.8000],
+                        [10.5000, 11.1000, 11.7000, 12.3000, 12.9000, 13.5000, 14.1000, 14.7000, 15.3000],
+                        [3.0000, 3.6000, 4.2000, 4.8000, 5.4000, 6.0000, 6.6000, 7.2000, 7.8000],
+                    ]
+                ]
+            ),
+            *device,
+        ]
+    )
+    TESTS.append(
+        [
+            {"pixdim": [-1, -1, 0.5], "padding_mode": "zeros", "dtype": float},
+            torch.ones((1, 2, 1, 2)),  # data
+            torch.eye(4),
+            {},
+            torch.tensor([[[[1.0, 1.0, 1.0]], [[1.0, 1.0, 1.0]]]]),
+            *device,
+        ]
+    )
+    TESTS.append(  # 5D input
+        [
+            {"pixdim": [-1, -1, 0.5], "padding_mode": "zeros", "dtype": float, "align_corners": True},
+            torch.ones((1, 2, 2, 2, 1)),  # data
+            torch.eye(4),
+            {},
+            torch.ones((1, 2, 2, 3, 1)),
+            *device,
+        ]
+    )
+    TESTS.append(  # 5D input
+        [
+            {"pixdim": 0.5, "padding_mode": "constant", "mode": "nearest", "scale_extent": True},
+            torch.ones((1, 368, 336, 368)),  # data
+            torch.tensor(
+                [
+                    [0.41, 0.005, 0.008, -79.7],
+                    [-0.0049, 0.592, 0.0664, -57.4],
+                    [-0.0073, -0.0972, 0.404, -32.1],
+                    [0.0, 0.0, 0.0, 1.0],
+                ]
+            ),
+            {},
+            torch.ones((1, 302, 403, 301)),
+            *device,
+        ]
+    )
+
+TESTS_TORCH = []
+for track_meta in (False, True):
+    for p in TEST_NDARRAYS_ALL:
+        TESTS_TORCH.append([[1.2, 1.3, 0.9], p(torch.zeros((1, 3, 4, 5))), track_meta])
+
+TEST_INVERSE = []
+for d in TEST_DEVICES:
+    for recompute in (False, True):
+        for align in (False, True):
+            for scale_extent in (False, True):
+                TEST_INVERSE.append([*d, recompute, align, scale_extent])
+
+
+@skip_if_quick
+class TestSpacingCase(unittest.TestCase):
+    @parameterized.expand(TESTS)
+    def test_spacing(
+        self,
+        init_param: dict,
+        img: torch.Tensor,
+        affine: torch.Tensor,
+        data_param: dict,
+        expected_output: torch.Tensor,
+        device: torch.device,
+    ):
+        img = MetaTensor(img, affine=affine).to(device)
+        tr = Spacing(**init_param)
+        call_param = data_param.copy()
+        call_param["data_array"] = img
+        res: MetaTensor = tr(**call_param)  # type: ignore
+        self.assertEqual(img.device, res.device)
+
+        test_resampler_lazy(tr, res, init_param=init_param, call_param=call_param)
+
+        assert_allclose(res, expected_output, atol=1e-1, rtol=1e-1)
+        sr = min(len(res.shape) - 1, 3)
+        if isinstance(init_param["pixdim"], float):
+            init_param["pixdim"] = [init_param["pixdim"]] * sr
+        init_pixdim = init_param["pixdim"][:sr]
+        norm = affine_to_spacing(res.affine, sr).cpu().numpy()
+        assert_allclose(fall_back_tuple(init_pixdim, norm), norm, type_test=False)  # type: ignore
+
+    @parameterized.expand(TESTS_TORCH)
+    def test_spacing_torch(self, pixdim, img, track_meta: bool):
+        set_track_meta(track_meta)
+        init_param = {"pixdim": pixdim}
+        tr = Spacing(**init_param)
+        call_param = {"data_array": img}
+        res = tr(**call_param)
+
+        if track_meta:
+            self.assertIsInstance(res, MetaTensor)
+            new_spacing = affine_to_spacing(res.affine, 3)  # type: ignore
+            assert_allclose(new_spacing, pixdim, type_test=False)
+            self.assertNotEqual(img.shape, res.shape)
+            test_resampler_lazy(tr, res, init_param=init_param, call_param=call_param)
+        else:
+            self.assertIsInstance(res, torch.Tensor)
+            self.assertNotIsInstance(res, MetaTensor)
+            self.assertNotEqual(img.shape, res.shape)
+        set_track_meta(True)
+
+    @parameterized.expand(TEST_INVERSE)
+    def test_inverse(self, device, recompute, align, scale_extent):
+        img_t = torch.rand((1, 10, 9, 8), dtype=torch.float32, device=device)
+        affine = torch.tensor(
+            [[0, 0, -1, 0], [1, 0, 0, 0], [0, 1, 0, 0], [0, 0, 0, 1]], dtype=torch.float32, device="cpu"
+        )
+        meta = {"fname": "somewhere"}
+        img = MetaTensor(img_t, affine=affine, meta=meta)
+        tr = Spacing(pixdim=[1.1, 1.2, 0.9], recompute_affine=recompute, align_corners=align, scale_extent=scale_extent)
+        # check that image and affine have changed
+        img = tr(img)
+        self.assertNotEqual(img.shape, img_t.shape)
+        l2_norm_affine = ((affine - img.affine) ** 2).sum() ** 0.5
+        self.assertGreater(l2_norm_affine, 5e-2)
+        # check that with inverse, image affine are back to how they were
+        img = tr.inverse(img)
+        self.assertEqual(img.applied_operations, [])
+        self.assertEqual(img.shape, img_t.shape)
+        l2_norm_affine = ((affine - img.affine) ** 2).sum() ** 0.5
+        self.assertLess(l2_norm_affine, 5e-2)
+
+    @parameterized.expand(TEST_INVERSE)
+    def test_inverse_mn_mx(self, device, recompute, align, scale_extent):
+        img_t = torch.rand((1, 10, 9, 8), dtype=torch.float32, device=device)
+        affine = torch.tensor(
+            [[0, 0, -1, 0], [1, 0, 0, 0], [0, 1, 0, 0], [0, 0, 0, 1]], dtype=torch.float32, device="cpu"
+        )
+        img = MetaTensor(img_t, affine=affine, meta={"fname": "somewhere"})
+        choices = [(None, None), [1.2, None], [None, 0.7], [0.7, 0.9]]
+        idx = np.random.choice(range(len(choices)), size=1)[0]
+        tr = Spacing(
+            pixdim=[1.1, 1.2, 0.9],
+            recompute_affine=recompute,
+            align_corners=align,
+            scale_extent=scale_extent,
+            min_pixdim=[0.9, None, choices[idx][0]],
+            max_pixdim=[1.1, 1.1, choices[idx][1]],
+        )
+        img_out = tr(img)
+        if isinstance(img_out, MetaTensor):
+            assert_allclose(
+                img_out.pixdim, [1.0, 1.125, 0.888889] if recompute else [1.0, 1.2, 0.9], type_test=False, rtol=1e-4
+            )
+        img_out = tr.inverse(img_out)
+        self.assertEqual(img_out.applied_operations, [])
+        self.assertEqual(img_out.shape, img_t.shape)
+        self.assertLess(((affine - img_out.affine) ** 2).sum() ** 0.5, 5e-2)
+
+    def test_property_no_change(self):
+        affine = torch.tensor(
+            [[1, 0, 0, 0], [0, 1, 0, 0], [0, 0, 1, 0], [0, 0, 0, 1]], dtype=torch.float32, device="cpu"
+        )
+        affine[:3] *= 1 - 1e-4  # make sure it's not exactly target but close to the target
+        img = MetaTensor(torch.rand((1, 10, 9, 8), dtype=torch.float32), affine=affine, meta={"fname": "somewhere"})
+        tr = Spacing(pixdim=[1.0, 1.0, 1.0])
+        tr(img)
+        assert_allclose(tr.pixdim, [1.0, 1.0, 1.0], type_test=False)
+
+
+if __name__ == "__main__":
+    unittest.main()