diff --git a/README.rst b/README.rst
index 77c31313..88fa5873 100644
--- a/README.rst
+++ b/README.rst
@@ -218,13 +218,15 @@ land cover types using ``cover_statistics``, and actual ET with ``ssebopeta_byge
stats = gh.cover_statistics(lulc.cover)
eta = gh.ssebopeta_bygeom(geometry, dates=("2005-10-01", "2005-10-05"))
-
+.. image:: https://raw.githubusercontent.com/cheginit/pygeohydro/master/docs/_static/lulc.png
+ :target: https://raw.githubusercontent.com/cheginit/pygeohydro/master/docs/_static/lulc.png
+ :width: 200
+ :height: 200
+
+.. image:: https://raw.githubusercontent.com/cheginit/pygeohydro/master/docs/_static/eta.png
+ :target: https://raw.githubusercontent.com/cheginit/pygeohydro/master/docs/_static/eta.png
+ :width: 200
+ :height: 200
Additionally, we can all the US dams data using ``get_nid`` and ``get_nid_codes``: