diff --git a/CHANGELOG.rst b/CHANGELOG.rst
index 1421d22..e0578b5 100644
--- a/CHANGELOG.rst
+++ b/CHANGELOG.rst
@@ -41,6 +41,8 @@ Fixed
   out.
 - nss_simulator.py, when asked to produce an output set of maps, uses a zero nodata value
   rather than whatever nodata value was present in the input burial depth map.
+- anaglyph.py needed some minor changes to align with upcoming changes in the
+  scikit-image architecture.
 
 
 0.6.0 (2023-09-25)
diff --git a/src/vipersci/vis/anaglyph.py b/src/vipersci/vis/anaglyph.py
index 57f8893..5d4730e 100755
--- a/src/vipersci/vis/anaglyph.py
+++ b/src/vipersci/vis/anaglyph.py
@@ -165,8 +165,8 @@ def correlate_and_shift(left: NDArray, right: NDArray) -> NDArray:
     """
     shift = tuple(
         phase_cross_correlation(
-            left, right, return_error=False, normalization=None
-        ).astype(int)
+            left, right, return_error="always", normalization=None
+        )[0].astype(int)
     )
     logger.info(f"Right image shift: {shift}")
     return np.roll(right, shift, (0, 1))