Enhance bursts transparency calculation in plot_scenes()

pygmtsar/pygmtsar/Stack.py

@@ -133,10 +133,13 @@ def plot_scenes(self, dem='auto', image=None, alpha=None, caption='Estimated Sce
         gdf = self.to_dataframe()
         cmap = matplotlib.colormaps[cmap]
         colors = dict([(v, cmap(k)) for k, v in enumerate(gdf.index.unique())])
-        # too small an alpha becomes invisible
-        gdf_alpha = 0.5/len(gdf)
-        gdf_alpha = gdf_alpha if gdf_alpha>=0.002 else 0.002
-        gdf.reset_index().plot(color=[colors[k] for k in gdf.index], alpha=gdf_alpha, edgecolor='black', ax=plt.gca())
+
+        # Calculate overlaps including self-overlap
+        overlap_count = [sum(1 for geom2 in gdf.geometry if geom1.intersects(geom2)) for geom1 in gdf.geometry]
+        # define transparency for the calculated overlaps
+        gdf_alpha = 1/max(overlap_count)
+        # apply minimum transparency threshold
+        gdf.reset_index().plot(color=[colors[k] for k in gdf.index], alpha=max(gdf_alpha, 0.002), edgecolor='black', ax=plt.gca())
         self.plot_AOI(**kwargs)
         self.plot_POI(**kwargs)
         if aspect is not None: