diff --git a/lib/matplotlib/tests/baseline_images/test_axes/aitoff_proj.png b/lib/matplotlib/tests/baseline_images/test_axes/aitoff_proj.png
new file mode 100644
index 000000000000..989e2bd5961a
Binary files /dev/null and b/lib/matplotlib/tests/baseline_images/test_axes/aitoff_proj.png differ
diff --git a/lib/matplotlib/tests/test_axes.py b/lib/matplotlib/tests/test_axes.py
index c36070442a4c..66f33ca413f0 100644
--- a/lib/matplotlib/tests/test_axes.py
+++ b/lib/matplotlib/tests/test_axes.py
@@ -797,6 +797,23 @@ def test_polar_rlim_zero():
     assert ax.get_ylim()[0] == 0
 
 
+@image_comparison(['aitoff_proj'], extensions=["png"],
+                  remove_text=True, style='mpl20')
+def test_aitoff_proj():
+    """
+    Test aitoff projection ref.:
+    https://github.com/matplotlib/matplotlib/pull/14451
+    """
+    x = np.linspace(-np.pi, np.pi, 20)
+    y = np.linspace(-np.pi / 2, np.pi / 2, 20)
+    X, Y = np.meshgrid(x, y)
+
+    fig, ax = plt.subplots(figsize=(8, 4.2),
+                           subplot_kw=dict(projection="aitoff"))
+    ax.grid()
+    ax.plot(X.flat, Y.flat, 'o', markersize=4)
+
+
 @image_comparison(['axvspan_epoch'])
 def test_axvspan_epoch():
     from datetime import datetime