diff --git a/scikits/image/morphology/tests/test_watershed.py b/scikits/image/morphology/tests/test_watershed.py
index ca7c43bf..265a25e2 100644
--- a/scikits/image/morphology/tests/test_watershed.py
+++ b/scikits/image/morphology/tests/test_watershed.py
@@ -388,11 +388,9 @@ class TestWatershed(unittest.TestCase):
         before = time.clock() 
         out = watershed(image, markers, self.eight)
         elapsed = time.clock() - before
-        print "Fast watershed ran a megapixel image in %f seconds"%(elapsed)
         before = time.clock()
         out = scipy.ndimage.watershed_ift(image.astype(np.uint16), markers, self.eight)
         elapsed = time.clock() - before
-        print "Scipy watershed ran a megapixel image in %f seconds"%(elapsed)
 
 
 class TestIsLocalMaximum(unittest.TestCase):
diff --git a/scikits/image/transform/tests/test_project.py b/scikits/image/transform/tests/test_project.py
index 363ff2b7..f1756f99 100644
--- a/scikits/image/transform/tests/test_project.py
+++ b/scikits/image/transform/tests/test_project.py
@@ -39,8 +39,6 @@ def test_fast_homography():
     H[:2, 2] = [tx, ty]
 
     for mode in ('constant', 'mirror', 'wrap'):
-        print 'Transform mode:', mode
-
         p0 = homography(img, H, mode=mode, order=1)
         p1 = fast_homography(img, H, mode=mode)
         p1 = np.round(p1)
@@ -54,7 +52,6 @@ def test_fast_homography():
         ## plt.show()
 
         d = np.mean(np.abs(p0 - p1))
-        print "delta=", d
         assert d < 0.2