diff --git a/lightly/transforms/__init__.py b/lightly/transforms/__init__.py
index 949fbe905..2cdb45963 100644
--- a/lightly/transforms/__init__.py
+++ b/lightly/transforms/__init__.py
@@ -18,6 +18,7 @@
 from lightly.transforms.dino_transform import DINOTransform, DINOViewTransform
 from lightly.transforms.fast_siam_transform import FastSiamTransform
 from lightly.transforms.gaussian_blur import GaussianBlur
+from lightly.transforms.gaussian_mixture_masks_transform import GaussianMixtureMasks
 from lightly.transforms.irfft2d_transform import IRFFT2DTransform
 from lightly.transforms.jigsaw import Jigsaw
 from lightly.transforms.mae_transform import MAETransform
diff --git a/lightly/transforms/gaussian_mixture_masks_transform.py b/lightly/transforms/gaussian_mixture_masks_transform.py
new file mode 100644
index 000000000..4ab98aeaf
--- /dev/null
+++ b/lightly/transforms/gaussian_mixture_masks_transform.py
@@ -0,0 +1,103 @@
+from typing import Tuple
+
+import torch
+import torch.fft
+from torch import Tensor
+
+from lightly.transforms.irfft2d_transform import IRFFT2DTransform
+from lightly.transforms.rfft2d_transform import RFFT2DTransform
+
+
+class GaussianMixtureMasks:
+    """Applies a Gaussian Mixture Mask in the Fourier domain to RGB images.
+
+    The mask is created using random Gaussian kernels, which are applied in
+    the frequency domain via RFFT2D, and then the IRFFT2D is used to return
+    to the spatial domain. The transformation is applied to each RGB channel separately.
+
+    Attributes:
+        num_gaussians: Number of Gaussian kernels to generate in the mixture mask.
+        std_range: Tuple containing the minimum and maximum standard deviation for the Gaussians.
+    """
+
+    def __init__(self, num_gaussians: int = 20, std_range: Tuple[int, int] = (10, 15)):
+        """Initializes GaussianMixtureMasks with the given parameters.
+
+        Args:
+            num_gaussians: Number of Gaussian kernels to generate in the mixture mask.
+            std_range: Tuple containing the minimum and maximum standard deviation for the Gaussians.
+        """
+
+        self.rfft2d_transform = RFFT2DTransform()
+        self.num_gaussians = num_gaussians
+        self.std_range = std_range
+
+    def gaussian_kernel(
+        self, size: Tuple[int, int], sigma: Tensor, center: Tensor
+    ) -> Tensor:
+        """Generates a 2D Gaussian kernel.
+
+        Args:
+            size: Tuple specifying the dimensions of the Gaussian kernel (C, H, W).
+            sigma: Tensor specifying the standard deviation of the Gaussian.
+            center: Tensor specifying the center of the Gaussian kernel.
+
+        Returns:
+            Tensor: A 2D Gaussian kernel.
+        """
+        u, v = torch.meshgrid(torch.arange(0, size[0]), torch.arange(0, size[1]))
+        u0, v0 = center
+        gaussian = torch.exp(
+            -((u - u0) ** 2 / (2 * sigma[0] ** 2) + (v - v0) ** 2 / (2 * sigma[1] ** 2))
+        )
+
+        return gaussian
+
+    def apply_gaussian_mixture_mask(
+        self, image_channel: Tensor, num_gaussians: int, std: Tuple[int, int]
+    ) -> Tensor:
+        """Applies the Gaussian mixture mask to a single channel in the frequency domain.
+
+        Args:
+            image_channel: Tensor representing a single channel of the image.
+            num_gaussians: Number of Gaussian kernels to generate in the mask.
+            std: Tuple specifying the standard deviation range for the Gaussians.
+
+        Returns:
+            Tensor: Image after applying the Gaussian mixture mask.
+        """
+        image_size = image_channel[0].shape
+
+        self.irfft2d_transform = IRFFT2DTransform((image_size[0], image_size[1]))
+        f_transform = self.rfft2d_transform(image_channel)
+
+        size = f_transform[0].shape
+
+        mask = torch.ones(size)
+
+        for _ in range(num_gaussians):
+            u0 = torch.randint(0, size[0], (1,))
+            v0 = torch.randint(0, size[1], (1,))
+            center = torch.tensor((u0, v0))
+            sigma = torch.rand(2) * 5 + 10
+
+            g_kernel = self.gaussian_kernel((size[0], size[1]), sigma, center)
+            mask -= g_kernel
+
+        filtered_f_transform = f_transform * mask
+        filtered_image = self.irfft2d_transform(filtered_f_transform).abs()
+        return filtered_image
+
+    def __call__(self, image_tensor: Tensor) -> Tensor:
+        """Applies the Gaussian mixture mask transformation to the input image.
+
+        Args:
+            image_tensor: Tensor representing an RGB image of shape (C, H, W).
+
+        Returns:
+            Tensor: The transformed image after applying the Gaussian mixture mask.
+        """
+        transformed_channel: Tensor = self.apply_gaussian_mixture_mask(
+            image_tensor, self.num_gaussians, self.std_range
+        )
+        return transformed_channel
diff --git a/tests/transforms/test_gaussian_mixture_masks.py b/tests/transforms/test_gaussian_mixture_masks.py
new file mode 100644
index 000000000..ae2fd9fe8
--- /dev/null
+++ b/tests/transforms/test_gaussian_mixture_masks.py
@@ -0,0 +1,10 @@
+import torch
+
+from lightly.transforms import GaussianMixtureMasks
+
+
+def test() -> None:
+    transform = GaussianMixtureMasks(20, (10, 15))
+    image = torch.rand(3, 32, 32)
+    output = transform(image)
+    assert output.shape == (3, 32, 32)