diff --git a/docs/tutorial/figures/kymographs/background.png b/docs/tutorial/figures/kymographs/background.png
new file mode 100644
index 000000000..e789086ef
--- /dev/null
+++ b/docs/tutorial/figures/kymographs/background.png
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:e24221d1969438ac81f50ae67363399811f05b71da4eacd4ad308f1186f856f1
+size 11316
diff --git a/docs/tutorial/figures/kymographs/kymograph_selection.png b/docs/tutorial/figures/kymographs/kymograph_selection.png
new file mode 100644
index 000000000..32b09199d
--- /dev/null
+++ b/docs/tutorial/figures/kymographs/kymograph_selection.png
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:9f56509037b53fa3aa5b5ba390cc0eda3921209013f57854ce6fcc6c4937b660
+size 13137
diff --git a/docs/tutorial/figures/kymographs/photon_counts.png b/docs/tutorial/figures/kymographs/photon_counts.png
new file mode 100644
index 000000000..41eefe1f0
--- /dev/null
+++ b/docs/tutorial/figures/kymographs/photon_counts.png
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:979e81baee4f80b54e50e9993a38ed19ab16009c3a2c6bf16f63fd0f43346abb
+size 45412
diff --git a/docs/tutorial/kymographs.rst b/docs/tutorial/kymographs.rst
index c8ecb6e28..46a85a8af 100644
--- a/docs/tutorial/kymographs.rst
+++ b/docs/tutorial/kymographs.rst
@@ -140,6 +140,50 @@ This returns a new (but flipped) :class:`~lumicks.pylake.kymo.Kymo`::
 
 .. image:: figures/kymographs/kymo_flipped.png
 
+Extracting photon counts
+------------------------
+
+Photon counts can be extracted from the raw image, for example for stoichiometry analysis.
+First, select a region for which the photon counts have to be determined::
+
+    kymo_selection = kymo["88.5s":"93s"].crop_by_distance(21.4, 22.3)
+
+    plt.figure()
+    kymo_selection.plot("rgb", adjustment=lk.ColorAdjustment(0, 98, mode="percentile"))
+
+Get the raw data (photon counts for each pixel) for the selected region::
+
+    selection_green = kymo_selection.get_image(channel="green")
+
+.. image:: figures/kymographs/kymograph_selection.png
+
+If the background has to be subtracted, choose a region without binding events as background::
+
+    background =  kymo["92.1s":"100s"].crop_by_distance(21.4, 22.3)
+    
+    plt.figure()
+    background.plot("rgb", adjustment=lk.ColorAdjustment(0, 98, mode="percentile"))
+
+.. image:: figures/kymographs/background.png
+
+Sum over the vertical axis (along the width of the binding event) to obtain the total photon counts over the selected regions.
+Since every pixel has an average background signal of `background_mean`, we need to sum the background over the same region to obtain the total background::
+
+    summed_photon_counts = np.sum(selection_green, axis=0)
+    background_mean = np.mean(background.get_image(channel="green"))
+    summed_background = background_mean * selection_green.shape[0]
+    time = np.arange(len(summed_photon_counts)) * kymo.line_time_seconds
+
+    plt.figure()
+    plt.plot(time, summed_photon_counts, label = "Binding event")
+    plt.axhline(summed_background, color='k', linestyle='--', label = "background")
+    plt.xlabel("time (s)")
+    plt.ylabel("Photon counts")
+    plt.title("Photon counts along binding event")
+    plt.legend()
+
+.. image:: figures/kymographs/photon_counts.png
+
 Calibrating to base pairs
 -------------------------