diff --git a/examples/cartesian_reconstruction.py b/examples/cartesian_reconstruction.py
index 90e4c6f8..c02b0270 100644
--- a/examples/cartesian_reconstruction.py
+++ b/examples/cartesian_reconstruction.py
@@ -11,8 +11,6 @@
 # acquired on a Cartesian grid but with partial echo and partial Fourier acceleration. Finally, we will reconstruct a
 # Cartesian scan with regular undersampling using iterative SENSE.
 
-# %% [markdown]
-# ## Import MRpro and download data
 
 # %% tags=["hide-cell"]
 # Get the raw data from zenodo
@@ -30,7 +28,8 @@
     print(f.name)
 
 # %% [markdown]
-# We have three different scans obtained from the same object with the same FOV and resolution:
+# We have three different scans obtained from the same object with the same FOV and resolution, saved as ISMRMRD
+# raw data files (*.mrd):
 #
 # - cart_t1.mrd is a fully sampled Cartesian acquisition
 #
@@ -42,7 +41,7 @@
 # %% [markdown]
 # ## Read in raw data and explore header
 #
-# To read in an ISMRMRD raw data file (*.mrd), we can simply pass on the file name to a `KData` object.
+# To read in an ISMRMRD raw data file (*.mrd), we can simply pass on the file name to a {py:class}`KData` object.
 # Additionally, we need to provide information about the trajectory. In MRpro, this is done using trajectory
 # calculators. These are functions that calculate the trajectory based on the acquisition information and additional
 # parameters provided to the calculators (e.g. the angular step for a radial acquisition).
@@ -58,7 +57,7 @@
 
 # %% [markdown]
 # Now we can explore this data object.
-# Simply calling ``print(kdata)``` gives us a basic overview of the KData object.
+# Simply calling ``print(kdata)`` gives us a basic overview of the `KData` object.
 
 # %%
 print(kdata)
@@ -168,7 +167,6 @@
 # This is because MRpro saves the trajectory in the most efficient way.
 # To get the full trajectory as a tensor, we can just call as_tensor().
 
-
 # %%
 # Plot the fully sampled trajectory (in blue)
 plt.plot(kdata.traj.as_tensor()[2, 0, 0, :, :].flatten(), kdata.traj.as_tensor()[1, 0, 0, :, :].flatten(), 'ob')