diff --git a/src/toolbox_scs/routines/boz.py b/src/toolbox_scs/routines/boz.py
index f118fb80e1b5cd49545d3f4b8429becc6aceea4f..5e47f03de15d01794df3d5a80aa6cc25443e49ab 100644
--- a/src/toolbox_scs/routines/boz.py
+++ b/src/toolbox_scs/routines/boz.py
@@ -451,6 +451,23 @@ def inspect_rois(data_mean, rois, threshold=None, allrois=False):
                        vmax=np.percentile(data_mean[:, :256], 99))
     main_ax.set_aspect('equal')
 
+    from matplotlib.patches import Rectangle
+    roi = rois['n']
+    main_ax.add_patch(Rectangle((roi['xl'], 128-roi['yh']),
+        roi['xh'] - roi['xl'],
+        roi['yh'] - roi['yl'],
+        alpha=0.3, color='b'))
+    roi = rois['0']
+    main_ax.add_patch(Rectangle((roi['xl'], 128-roi['yh']),
+        roi['xh'] - roi['xl'],
+        roi['yh'] - roi['yl'],
+        alpha=0.3, color='g'))
+    roi = rois['p']
+    main_ax.add_patch(Rectangle((roi['xl'], 128-roi['yh']),
+        roi['xh'] - roi['xl'],
+        roi['yh'] - roi['yl'],
+        alpha=0.3, color='r'))
+
     x.plot(Xs, pX)
     x.invert_yaxis()
     if threshold is not None: