diff --git a/skimage/transform/__init__.py b/skimage/transform/__init__.py
index 4f00a076..c6408309 100644
--- a/skimage/transform/__init__.py
+++ b/skimage/transform/__init__.py
@@ -5,6 +5,7 @@ from .radon_transform import radon, iradon, iradon_sart
 from .finite_radon_transform import frt2, ifrt2
 from .integral import integral_image, integrate
 from ._geometric import (warp, warp_coords, estimate_transform,
+                         matrix_transform,
                          SimilarityTransform, AffineTransform,
                          ProjectiveTransform, PolynomialTransform,
                          PiecewiseAffineTransform)
@@ -30,6 +31,7 @@ __all__ = ['hough_circle',
            'warp',
            'warp_coords',
            'estimate_transform',
+           'matrix_transform',
            'SimilarityTransform',
            'AffineTransform',
            'ProjectiveTransform',
diff --git a/skimage/transform/tests/test_geometric.py b/skimage/transform/tests/test_geometric.py
index 774d0888..06c4f652 100644
--- a/skimage/transform/tests/test_geometric.py
+++ b/skimage/transform/tests/test_geometric.py
@@ -3,7 +3,7 @@ from numpy.testing import (assert_equal, assert_array_almost_equal,
                            assert_raises)
 from skimage.transform._geometric import _stackcopy
 from skimage.transform._geometric import GeometricTransform
-from skimage.transform import (estimate_transform,
+from skimage.transform import (estimate_transform, matrix_transform,
                                SimilarityTransform, AffineTransform,
                                ProjectiveTransform, PolynomialTransform,
                                PiecewiseAffineTransform)
@@ -47,6 +47,11 @@ def test_estimate_transform():
                   SRC[:2, :], DST[:2, :])
 
 
+def test_matrix_transform():
+    tform = AffineTransform(scale=(0.1, 0.5), rotation=2)
+    assert_equal(tform(SRC), matrix_transform(SRC, tform._matrix))
+
+
 def test_similarity_estimation():
     # exact solution
     tform = estimate_transform('similarity', SRC[:2, :], DST[:2, :])