diff --git a/src/Math/testing/VoxImageTest.cpp b/src/Math/testing/VoxImageTest.cpp
index ad0c156..01ca09c 100644
--- a/src/Math/testing/VoxImageTest.cpp
+++ b/src/Math/testing/VoxImageTest.cpp
@@ -89,17 +89,31 @@ int main() {
   }
 
   {
-    VoxImage<TestVoxel> img(Vector3i(4, 4, 4));
+    const int size = 100;
+    VoxImage<TestVoxel> img(Vector3i(size, size, size));
     img.InitVoxels({0.f, 0});
-    for (int i = 0; i < 4; i++) {
-      for (int j = 0; j < 4; j++) {
-        for (int k = 0; k < 4; k++) {
+    for (int i = 0; i < size; i++) {
+      for (int j = 0; j < size; j++) {
+        for (int k = 0; k < size; k++) {
           img[Vector3i(i, j, k)] = {static_cast<float>(i + j + k), 0};
         }
       }
     }
-    img.ExportToVti("./vti_saved.vti", 0, 1);
-    // img.ImportFromVtkXml("./test_vox_image.vti");
+    img.ExportToVti("./vti_saved.vti", 0, true); // compressed
+
+    VoxImage<TestVoxel> imgR(Vector3i(0, 0, 0));
+    imgR.ImportFromVti("./vti_saved.vti", 0);
+    TEST1(img.GetDims() == imgR.GetDims());
+    TEST1(img.GetSpacing() == imgR.GetSpacing());
+    TEST1(img.GetOrigin() == imgR.GetOrigin());
+    TEST1(img.GetPosition() == imgR.GetPosition());
+    for (int i = 0; i < size; i++) {
+      for (int j = 0; j < size; j++) {
+        for (int k = 0; k < size; k++) {
+          TEST1(img[Vector3i(i, j, k)].Value == imgR[Vector3i(i, j, k)].Value);
+        }
+      }
+    }
   }
 
   {