code: def getVoxelsScalar(self, vmin=None, vmax=None):
                  
          nx, ny, nz = self.image.GetDimensions()
          
          voxdata = np.zeros([nx, ny, nz])
          lsx = range(nx)
          lsy = range(ny)
          lsz = range(nz)
          renorm = vmin is not None and vmax is not None
          for i in lsx:
              for j in lsy:
                  for k in lsz:
                      s = self.image.GetScalarComponentAsFloat(i, j, k, 0)
                      if renorm:
                          s = (s - vmin) / (vmax - vmin)
                      voxdata[i, j, k] = s
          return voxdata
docs: Return voxel content as a ``numpy.array``.
          
          :param float vmin: rescale scalar content to match `vmin` and `vmax` range. 
          :param float vmax: rescale scalar content to match `vmin` and `vmax` range.