It doesn't mention the constraint that the coordinate systems must be equal, otherwise it raises a ValueError.

And the matching-coord-sys constraint also keeps making me wonder if that should be reflected in the function name. Perhaps by using "rectilinear" in the name?