diff --git a/metpy/calc/tests/test_calc_tools.py b/metpy/calc/tests/test_calc_tools.py
index ac222ac782d..b9f46f41fbf 100644
--- a/metpy/calc/tests/test_calc_tools.py
+++ b/metpy/calc/tests/test_calc_tools.py
@@ -158,6 +158,21 @@ def test_reduce_point_density(thin_point_data, radius, truth):
     assert_array_equal(reduce_point_density(thin_point_data, radius=radius), truth)
 
 
+@pytest.mark.parametrize('radius, truth',
+                         [(2.0, np.array([1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+                                          0, 0, 0, 0, 0, 0, 0, 0, 0, 0], dtype=np.bool)),
+                          (1.0, np.array([1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+                                          0, 0, 0, 0, 0, 0, 0, 0, 0, 0], dtype=np.bool)),
+                          (0.3, np.array([1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+                                          0, 0, 0, 0, 0, 0, 0, 0, 0, 0], dtype=np.bool)),
+                          (0.1, np.array([1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+                                          0, 0, 0, 0, 0, 0, 0, 0, 1, 0], dtype=np.bool))
+                          ])
+def test_reduce_point_density_units(thin_point_data, radius, truth):
+    r"""Test that reduce_point_density works with units."""
+    assert_array_equal(reduce_point_density(thin_point_data, radius=radius * units.dam), truth)
+
+
 @pytest.mark.parametrize('radius, truth',
                          [(2.0, np.array([0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
                                           0, 0, 0, 0, 0, 0, 0, 0, 0, 1], dtype=np.bool)),
diff --git a/metpy/calc/tools.py b/metpy/calc/tools.py
index f3997465f39..1981ad51dbb 100644
--- a/metpy/calc/tools.py
+++ b/metpy/calc/tools.py
@@ -259,14 +259,15 @@ def reduce_point_density(points, radius, priority=None):
     data), returning a mask that can be used to select the points from one or more arrays
     (e.g. arrays of temperature and dew point). The points selected can be controlled by
     providing an array of ``priority`` values (e.g. rainfall totals to ensure that
-    stations with higher precipitation remain in the mask).
+    stations with higher precipitation remain in the mask). The radius can be specified
+    as a `pint.Quantity` with units. If none are provided, meters are assumed.
 
     Parameters
     ----------
     points : (N, K) array-like
         N locations of the points in K dimensional space
-    radius : float
-        minimum radius allowed between points
+    radius : `pint.Quantity` or float
+        Minimum radius allowed between points. If units are not provided, meters is assumed.
     priority : (N, K) array-like, optional
         If given, this should have the same shape as ``points``; these values will
         be used to control selection priority for points.
@@ -285,6 +286,10 @@ def reduce_point_density(points, radius, priority=None):
     array([False,  True, False])
 
     """
+    # Handle a radius with units. Assume meters if units are not specified
+    if hasattr(radius, 'units'):
+        radius = radius.to('m').m
+
     # Handle 1D input
     if points.ndim < 2:
         points = points.reshape(-1, 1)