diff --git a/console/nii_dicom.cpp b/console/nii_dicom.cpp
index 0be768d7..13605a85 100644
--- a/console/nii_dicom.cpp
+++ b/console/nii_dicom.cpp
@@ -5982,6 +5982,7 @@ const uint32_t kEffectiveTE = 0x0018 + uint32_t(0x9082 << 16); //FD
 			if (d.manufacturer != kMANUFACTURER_SIEMENS)
 				break;
 			dcmMultiFloat(lLength, (char *)&buffer[lPos], 3, &d.CSA.tablePos[0]); //slice position
+			d.CSA.tablePos[3] = -d.CSA.tablePos[3]; //reverse Z polarity, issue 726
 			d.CSA.tablePos[0] = 1.0; //set
 			break;
 		}
@@ -6794,6 +6795,7 @@ const uint32_t kEffectiveTE = 0x0018 + uint32_t(0x9082 << 16); //FD
 			if (d.manufacturer != kMANUFACTURER_SIEMENS)
 				break;
 			dcmMultiFloat(lLength, (char *)&buffer[lPos], 3, &d.CSA.tablePos[0]); //slice position
+			d.CSA.tablePos[3] = -d.CSA.tablePos[3]; //reverse Z polarity, issue 726
 			d.CSA.tablePos[0] = 1.0; //set
 			break;
 		}
@@ -7138,7 +7140,7 @@ const uint32_t kEffectiveTE = 0x0018 + uint32_t(0x9082 << 16); //FD
 		case kMRStackTablePosLong: //FL
 			if (d.manufacturer != kMANUFACTURER_PHILIPS)
 				break;
-			d.CSA.tablePos[3] = dcmFloat(lLength, &buffer[lPos], d.isLittleEndian);
+			d.CSA.tablePos[3] = - dcmFloat(lLength, &buffer[lPos], d.isLittleEndian);
 			d.CSA.tablePos[0] = 1.0;
 			break;
 		case kMRImageDiffBValueNumber: