diff --git a/docs/reference/utils/index.rst b/docs/reference/utils/index.rst
index 7f872059b..af803d7be 100644
--- a/docs/reference/utils/index.rst
+++ b/docs/reference/utils/index.rst
@@ -8,3 +8,4 @@ Utils
    colors
    math
    visualization
+   tensor_dataclass
\ No newline at end of file
diff --git a/docs/reference/utils/tensor_dataclass.rst b/docs/reference/utils/tensor_dataclass.rst
new file mode 100644
index 000000000..fb6b28f4f
--- /dev/null
+++ b/docs/reference/utils/tensor_dataclass.rst
@@ -0,0 +1,8 @@
+.. _tensor_dataclass:
+
+TensorDataclass
+=================
+
+.. automodule:: pyrad.utils.tensor_dataclass
+   :members:
+   :show-inheritance:
diff --git a/pyrad/utils/tensor_dataclass.py b/pyrad/utils/tensor_dataclass.py
new file mode 100644
index 000000000..3afbab2ba
--- /dev/null
+++ b/pyrad/utils/tensor_dataclass.py
@@ -0,0 +1,148 @@
+# Copyright 2022 The Plenoptix Team. All rights reserved.
+#
+# 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.
+
+"""Tensor dataclass"""
+
+import dataclasses
+from typing import Tuple
+
+import numpy as np
+import torch
+
+
+class TensorDataclass:
+    """Data class of tensor with the same size batch. Allows indexing and standard tensor ops.
+
+    Example:
+
+    .. code-block:: python
+
+        @dataclass
+        class TestTensorDataclass(TensorDataclass):
+            a: torch.Tensor
+            b: torch.Tensor
+            c: torch.Tensor = None
+
+        # Create a new tensor dataclass with batch size of [2,3,4]
+        test = TestTensorDataclass(a=torch.ones((2, 3, 4, 2)), b=torch.ones((4, 3)))
+
+        test.shape  # [2, 3, 4]
+        test.a.shape  # [2, 3, 4, 2]
+        test.b.shape  # [2, 3, 4, 3]
+
+        test.reshape((6,4)).shape  # [6, 4]
+        test.flatten().shape  # [24,]
+
+        test[..., 0].shape  # [2, 3]
+        test[:, 0, :].shape  # [2, 4]
+    """
+
+    _shape: tuple
+
+    def __post_init__(self) -> None:
+        if not dataclasses.is_dataclass(self):
+            raise TypeError("TensorDataclass must be a dataclass")
+
+        fields = dataclasses.asdict(self)
+        batch_shapes = [v.shape[:-1] for v in fields.values() if v is not None]
+        if len(batch_shapes) == 0:
+            raise ValueError("TensorDataclass must have at least one tensor")
+        batch_shape = torch.broadcast_shapes(*batch_shapes)
+
+        for k, v in fields.items():
+            if v is not None:
+                self.__setattr__(k, torch.broadcast_to(v, (*batch_shape, v.shape[-1])))
+
+        self.__setattr__("_shape", batch_shape)
+
+    def __getitem__(self, indices) -> "TensorDataclass":
+        if isinstance(indices, int):
+            indices = (indices,)
+        return self._apply_fn_to_fields(lambda x: x[indices + (slice(None),)])
+
+    def __len__(self) -> int:
+        return self.shape[0]
+
+    def __bool__(self) -> bool:
+        if len(self) == 0:
+            raise ValueError(
+                f"The truth value of {self.__class__.__name__} when `len(x) == 0` "
+                "is ambiguous. Use `len(x)` or `x is not None`."
+            )
+        return True
+
+    @property
+    def shape(self) -> tuple:
+        """Returns the batch shape of the tensor dataclass."""
+        return self._shape
+
+    @property
+    def size(self) -> int:
+        """Returns the number of elements in the tensor dataclass batch dimension."""
+        if len(self._shape) == 0:
+            return 1
+        return int(np.prod(self._shape))
+
+    @property
+    def ndim(self) -> int:
+        """Returns the number of dimensions of the tensor dataclass."""
+        return len(self._shape)
+
+    def reshape(self, shape: Tuple[int, ...]) -> "TensorDataclass":
+        """Returns a new TensorDataclass with the same data but with a new shape.
+
+        Args:
+            shape (Tuple[int]): The new shape of the tensor dataclass.
+
+        Returns:
+            TensorDataclass: A new TensorDataclass with the same data but with a new shape.
+        """
+        if isinstance(shape, int):
+            shape = (shape,)
+        return self._apply_fn_to_fields(lambda x: x.reshape((*shape, x.shape[-1])))
+
+    def flatten(self) -> "TensorDataclass":
+        """Returns a new TensorDataclass with flattened batch dimensions
+
+        Returns:
+            TensorDataclass: A new TensorDataclass with the same data but with a new shape.
+        """
+        return self.reshape((-1,))
+
+    def broadcast_to(self, shape: Tuple[int]) -> "TensorDataclass":
+        """Returns a new TensorDataclass broadcast to new shape.
+
+        Args:
+            shape (Tuple[int]): The new shape of the tensor dataclass.
+
+        Returns:
+            TensorDataclass: A new TensorDataclass with the same data but with a new shape.
+        """
+
+        return self._apply_fn_to_fields(lambda x: x.broadcast_to((*shape, x.shape[-1])))
+
+    def _apply_fn_to_fields(self, fn: callable) -> "TensorDataclass":
+        """Applies a function to all fields of the tensor dataclass.
+
+        Args:
+            fn (callable): The function to apply to all fields.
+
+        Returns:
+            TensorDataclass: A new TensorDataclass with the same data but with a new shape.
+        """
+
+        fields = dataclasses.asdict(self)
+        fields = {k: fn(v) for k, v in fields.items() if v is not None}
+
+        return dataclasses.replace(self, **fields)
diff --git a/tests/utils/test_tensor_dataclass.py b/tests/utils/test_tensor_dataclass.py
new file mode 100644
index 000000000..e2dcff02d
--- /dev/null
+++ b/tests/utils/test_tensor_dataclass.py
@@ -0,0 +1,111 @@
+"""
+Test tensor dataclass
+"""
+from dataclasses import dataclass
+from typing import Any
+import pytest
+import torch
+
+from pyrad.utils.tensor_dataclass import TensorDataclass
+
+
+@dataclass
+class TestTensorDataclass(TensorDataclass):
+    """Dummy dataclass"""
+
+    a: torch.Tensor
+    b: torch.Tensor
+    c: torch.Tensor = None
+
+
+def test_init():
+    """Test that dataclass is properly initialized"""
+
+    @dataclass
+    class Dummy(TensorDataclass):
+        """Dummy dataclass"""
+
+        dummy_vals: Any = None
+
+    dummy = Dummy(dummy_vals=torch.ones(1))
+    with pytest.raises(ValueError):
+        dummy = Dummy()
+
+
+def test_broadcasting():
+    """Test broadcasting during init"""
+
+    a = torch.ones((4, 6, 3))
+    b = torch.ones((6, 2))
+    tensor_dataclass = TestTensorDataclass(a=a, b=b)
+    assert tensor_dataclass.b.shape == (4, 6, 2)
+
+    a = torch.ones((4, 6, 3))
+    b = torch.ones((2))
+    tensor_dataclass = TestTensorDataclass(a=a, b=b)
+    assert tensor_dataclass.b.shape == (4, 6, 2)
+
+    # Invalid broadcasting
+    a = torch.ones((4, 6, 3))
+    b = torch.ones((3, 2))
+    with pytest.raises(RuntimeError):
+        tensor_dataclass = TestTensorDataclass(a=a, b=b)
+
+
+def test_tensor_ops():
+    """Test get shape"""
+
+    a = torch.ones((4, 6, 3))
+    b = torch.ones((6, 2))
+    tensor_dataclass = TestTensorDataclass(a=a, b=b)
+
+    assert tensor_dataclass.shape == (4, 6)
+    assert tensor_dataclass.a.shape == (4, 6, 3)
+    assert tensor_dataclass.b.shape == (4, 6, 2)
+    assert tensor_dataclass.size == 24
+    assert tensor_dataclass.ndim == 2
+    assert len(tensor_dataclass) == 4
+
+    reshaped = tensor_dataclass.reshape((2, 12))
+    assert reshaped.shape == (2, 12)
+    assert reshaped.a.shape == (2, 12, 3)
+    assert reshaped.b.shape == (2, 12, 2)
+
+    flattened = tensor_dataclass.flatten()
+    assert flattened.shape == (24,)
+    assert flattened.a.shape == (24, 3)
+    assert flattened.b.shape == (24, 2)
+
+    # Test indexing operations
+    assert tensor_dataclass[:, 1].shape == (4,)
+    assert tensor_dataclass[:, 1].a.shape == (4, 3)
+    assert tensor_dataclass[..., 1].shape == (4,)
+    assert tensor_dataclass[..., 1].a.shape == (4, 3)
+    assert tensor_dataclass[0].shape == (6,)
+    assert tensor_dataclass[0].a.shape == (6, 3)
+    assert tensor_dataclass[0, ...].shape == (6,)
+    assert tensor_dataclass[0, ...].a.shape == (6, 3)
+
+    tensor_dataclass = TestTensorDataclass(a=torch.ones((2, 3, 4, 5)), b=torch.ones((4, 5)))
+    assert tensor_dataclass[0, ...].shape == (3, 4)
+    assert tensor_dataclass[0, ...].a.shape == (3, 4, 5)
+    assert tensor_dataclass[0, ..., 0].shape == (3,)
+    assert tensor_dataclass[0, ..., 0].a.shape == (3, 5)
+    assert tensor_dataclass[..., 0].shape == (2, 3)
+    assert tensor_dataclass[..., 0].a.shape == (2, 3, 5)
+
+
+def test_iter():
+    """Test iterating over tensor dataclass"""
+    tensor_dataclass = TestTensorDataclass(a=torch.ones((3, 4, 5)), b=torch.ones((3, 4, 5)))
+    for batch in tensor_dataclass:
+        assert batch.shape == (4,)
+        assert batch.a.shape == (4, 5)
+        assert batch.b.shape == (4, 5)
+
+
+if __name__ == "__main__":
+    test_init()
+    test_broadcasting()
+    test_tensor_ops()
+    test_iter()