diff --git a/coastsat/SDS_shoreline.py b/coastsat/SDS_shoreline.py
index de893fc0..0491438c 100644
--- a/coastsat/SDS_shoreline.py
+++ b/coastsat/SDS_shoreline.py
@@ -118,7 +118,7 @@ def extract_shorelines(metadata, settings):
         str_new = ''
         if not sklearn.__version__[:4] == '0.20':
             str_new = '_new'
-        if satname in ['L5','L7','L8']:
+        if satname in ['L5','L7','L8','L9']:
             pixel_size = 15
             if settings['sand_color'] == 'dark':
                 clf = joblib.load(os.path.join(filepath_models, 'NN_4classes_Landsat_dark%s.pkl'%str_new))