We consider a class of maps from the unit square to itself preserving a contracting foliation and inducing a one-dimensional map having an absolutely continuous invariant measure. We show how the physical measure of those systems can be rigorously approximated with an explicitly given bound on the error with respect to the Wasserstein distance. We present a rigorous implementation of our algorithm using interval arithmetics, and the result of the computation on a non-trivial example of a Lorenz-like two-dimensional map and its attractor, obtaining a statement on its local dimension.