diff --git a/docs/source/changelog.rst b/docs/source/changelog.rst
index 281d24ead..9878b3bed 100644
--- a/docs/source/changelog.rst
+++ b/docs/source/changelog.rst
@@ -3,6 +3,28 @@
 Changelog
 =========
 
+v 2.2.7
+-------
+
+Winston-Lutz
+^^^^^^^^^^^^
+
+* A small change was made to the Winston-Lutz BB finding algorithm to be more robust and use less custom code. The output from WL analyses should be within 0.1mm of previous values.
+* A section was added to the documentation to describe how images are classified and the analysis of output from the .results() method.
+
+Bug Fixes
+^^^^^^^^^
+
+* `#187 <https://github.com/jrkerns/pylinac/issues/187>`_ Scipy's imresize function has been deprecated. Functionality was converted to use skimage.transform.resize().
+* `#185 <https://github.com/jrkerns/pylinac/issues/185>`_ Winston-Lutz PDF generation had an artifact causing catastrophic failure.
+* `#183 <https://github.com/jrkerns/pylinac/issues/183>`_ The Bakai fomula of the gamma calculation had an operational inconsistency such that dose-to-agreement other than 1% would give incorrect values of the gamma value.
+* `#190 <https://github.com/jrkerns/pylinac/issues/190>`_ The Catphan module had an inconsistency in the rMTF/spatial resolution determination. Some line pair regions would be detected for some phantoms and not for others. This was caused by the different CatPhan models having slighly different rotations of the CTP528 module. Pylinac now has model-specific boundaries.
+* `#192 <https://github.com/jrkerns/pylinac/issues/192>`_ The FlatSym plot would conflate the vertical and horizontal lines shown on the analyzed image. Analysis is unaffected, only the depiction of position.
+* `#194 <https://github.com/jrkerns/pylinac/issues/194>`_ The Leeds low contrast ROI color on the analyzed image was not consistent with the contrast plots. ROI color is now based on the pass/fail of the contrast constant, not the contrast.
+* `#196 <https://github.com/jrkerns/pylinac/issues/196>`_ Winston-Lutz images with a dense BB and low photon energy could cause BB detection to fail. A better BB-finding algorithm has been implemented.
+* `#197 <https://github.com/jrkerns/pylinac/issues/197>`_ EPID RMS deviation would return 0 for the .results() method always. This now calculates correctly.
+
+
 V 2.2.6
 -------
 
diff --git a/pylinac/core/profile.py b/pylinac/core/profile.py
index cd75bb383..9a315cfb9 100644
--- a/pylinac/core/profile.py
+++ b/pylinac/core/profile.py
@@ -531,7 +531,7 @@ def __init__(self, values: Union[np.ndarray, Sequence]):
         self.peaks = []
         self.valleys = []
 
-    def plot(self, show_peaks: bool=True):
+    def plot(self, show_peaks: bool=True, ax=None):
         """Plot the profile.
 
         Parameters
@@ -540,7 +540,8 @@ def plot(self, show_peaks: bool=True):
             Whether to plot the peak locations as well. Will not show if a peak search has
             not yet been done.
         """
-        fig, ax = plt.subplots()
+        if ax is None:
+            fig, ax = plt.subplots()
         ax.plot(self.values)
         if show_peaks:
             peaks_x = [peak.idx for peak in self.peaks]
diff --git a/setup.py b/setup.py
index c1400bb4a..cfe2854ba 100644
--- a/setup.py
+++ b/setup.py
@@ -1,6 +1,6 @@
 from setuptools import setup, find_packages
 
-__version__ = '2.2.6'
+__version__ = '2.2.7'
 
 
 setup(
@@ -14,11 +14,11 @@
     author='James Kerns',
     author_email='jkerns100@gmail.com',
     description='A toolkit for performing TG-142 QA-related tasks on a linear accelerator',
-    install_requires=["numpy >= 1.11",
-                      "scipy >= 0.17",
+    install_requires=["numpy >= 1.14, <= 1.17",
+                      "scipy >= 1.0",
                       "pydicom >= 1.0",
                       "matplotlib >= 1.4",
-                      "scikit-image >= 0.12",
+                      "scikit-image >= 0.12, <=0.14",
                       "scikit-learn >= 0.18",
                       "Pillow >= 4.0",
                       "tqdm >= 3.8",