diff --git a/ehtim/image.py b/ehtim/image.py
index 90e98664..3dbf1933 100644
--- a/ehtim/image.py
+++ b/ehtim/image.py
@@ -3485,7 +3485,7 @@ def display(self, pol=None, cfun=False, interp='gaussian',
             else:
                 areaunit = 'beam'
                 beamarea = (2.0 * np.pi * beamparams[0] * beamparams[1] / (8.0 * np.log(2)))
-                factor = beamarea / (self.psize**2)
+                factor *= beamarea / (self.psize**2)
                 if power != 0:
                     areaunit = areaunit + (r' ($10^{{' + str(power) + '}}$ K)')