diff --git a/plantcv/geospatial/analyze/coverage.py b/plantcv/geospatial/analyze/coverage.py
index 7b6f5ae..45b0fbe 100644
--- a/plantcv/geospatial/analyze/coverage.py
+++ b/plantcv/geospatial/analyze/coverage.py
@@ -85,8 +85,9 @@ def coverage(img, bin_mask, geojson):
     plt.title("Shapefile on GeoTIFF")
     plt.xlabel("Longitude")
     plt.ylabel("Latitude")
-    # Show the plot
-    plt.show()
-    _debug(visual=bounds.boundary.plot(),
+    # Store the plot
+    plotting_img = plt.gcf()
+
+    _debug(visual=plotting_img,
            filename=os.path.join(params.debug_outdir, f"{params.device}_analyze_coverage.png"))
-    return plt.
\ No newline at end of file
+    return plotting_img