Changeset 75313af in sasview for src/sas

Timestamp:

Dec 12, 2018 1:50:57 PM (6 years ago)

Author:

Paul Kienzle <pkienzle@…>

Branches:

master, magnetic_scatt, release-4.2.2, ticket-1009, ticket-1094-headless, ticket-1242-2d-resolution, ticket-1249

Children:

1342f6a

Parents:

175b83f

Message:

fix 2D log-scale plotting; maybe broke vmin/vmax limits

File:

• src/sas/sasgui/plottools/PlotPanel.py (modified) (2 diffs)

src/sas/sasgui/plottools/PlotPanel.py

@@ -1434 +1434 @@
 
         """
+        # TODO: include mask info in plotter
         self.data = data
         self.qx_data = qx_data
@@ -1451 +1452 @@
         else:
             output = copy.deepcopy(self.data)
-        # check scale
+        # rescale data if necessary
         if self.scale == 'log_{10}':
-            try:
-                if  self.zmin_2D <= 0  and len(output[output > 0]) > 0:
-                    zmin_temp = self.zmin_2D
-                    output[output > 0] = np.log10(output[output > 0])
-                    #In log scale Negative values are not correct in general
-                    #output[output<=0] = math.log(np.min(output[output>0]))
-                elif self.zmin_2D <= 0:
-                    zmin_temp = self.zmin_2D
-                    output[output > 0] = np.zeros(len(output))
-                    output[output <= 0] = -32
-                else:
-                    zmin_temp = self.zmin_2D
-                    output[output > 0] = np.log10(output[output > 0])
-                    #In log scale Negative values are not correct in general
-                    #output[output<=0] = math.log(np.min(output[output>0]))
-            except:
-                #Too many problems in 2D plot with scale
-                pass
-
-        else:
-            zmin_temp = self.zmin_2D
+            with np.errstate(all='ignore'):
+                output = np.log10(output)
+            index = np.isfinite(output)
+            if not index.all():
+                cutoff = (np.min(output[index]) - np.log10(2)
+                          if index.any() else 0.)
+                output[~index] = cutoff
+        # TODO: fix handling of zmin_2D/zmax_2D in _onToggleScale
+        # For now, use default vmin/vmax from data
+        #vmin, vmax = self.zmin_2D, self.zmax_2D
+        vmin, vmax = None, None
         self.cmap = cmap
         if self.dimension != 3:
             #Re-adjust colorbar
             self.subplot.figure.subplots_adjust(left=0.2, right=.8, bottom=.2)
-
             im = self.subplot.imshow(output, interpolation='nearest',
                                      origin='lower',
-                                     vmin=zmin_temp, vmax=self.zmax_2D,
+                                     vmin=vmin, vmax=vmax,
                                      cmap=self.cmap,
                                      extent=(self.xmin_2D, self.xmax_2D,