diff --git a/vviewer/Image.py b/vviewer/Image.py index 016d7c44a4949178ec316466b6b701c1b5f801c6..f1254dc7482bf56a7a12ada8289bb5f8323b860d 100644 --- a/vviewer/Image.py +++ b/vviewer/Image.py @@ -340,7 +340,7 @@ class Image(object): return self.image.affine # probably not used? - def getUsedAffine(self): + def getAffineUsed(self): return self.affine_res_inv def getVoxelCoords(self, m): @@ -356,6 +356,7 @@ class Image(object): yet!") def setUnresampled(self): + # TODO: setze self.affine_res_inv = np.eye(4)? self.image_res = self.image.get_data() def setHistogram(self): diff --git a/vviewer/vviewer.py b/vviewer/vviewer.py index 4d38b17a51905cea91fe4e57a133cc1b5b083f6a..f72b4c50a7ef2d857da7d4dc2732024000312b92 100755 --- a/vviewer/vviewer.py +++ b/vviewer/vviewer.py @@ -1480,10 +1480,7 @@ class vviewer(QtGui.QMainWindow): index = self.imagelist.currentRow() if self.voxel_coord: if index >= 0: - img_affine = self.images[index].getAffine() - back_affine_inv = np.dot( - np.linalg.inv(self.affine), img_affine) - m = self.applyTransform(x, y, z, back_affine_inv) + m = self.applyTransform(x, y, z, self.images[index].getAffineUsed()) else: QtGui.QMessageBox.warning(self, "Warning", "Error: No image selected to display voxel \