Exemplo n.º 1
0
 def get_volume(self, offset, offset_inner=None):
     if offset_inner == None:
         offset_inner = -pow(10., 38)
     cnt = 0.
     scl = Scaler(self._voxels)
     vo = scl.scale(self._voxels.shape[0]*2, self._voxels.shape[1]*2, self._voxels.shape[2]*2).value
     print vo.shape
     for val in vo.flatten():
         if (val < offset) and (val > offset_inner):
             cnt += 1.
     return cnt*(self.spc/2)**3
Exemplo n.º 2
0
 def get_voxels(self, refine=1):
     scaler = Scaler(self._voxels)
     dims = (self._voxels.shape[0]*refine,self._voxels.shape[1]*refine,self._voxels.shape[2]*refine)
     return scaler.scale(*dims).value