Simple combine script with rasterio

.gitignore

@@ -6,3 +6,4 @@ config*
 token*
 local*
 data/
+output/

landsat/ReadMe.md

@@ -36,3 +36,10 @@ cd CLASS-Examples/landsat/
 bash get-index.sh
 cat data/index.csv | python3 search.py | bash download.sh
 ```
+
+## Analysis
+```
+sudo apt-get install python3-rasterio
+install -dv output/
+python3 combine.py
+```

landsat/combine.py

@@ -0,0 +1,17 @@
+#!/usr/bin/python3
+import os
+import rasterio
+
+# Open the first directory, could walk entire tree
+for dirname, dirs, files in os.walk('data'):
+    source = dirs[0]
+    break
+
+with rasterio.open("data/%s/%s_B2.TIF" % (source, source)) as band2:
+    meta = band2.meta
+
+meta.update(count = 3, driver='PNG')
+with rasterio.open('output/result.png', 'w+', **meta) as output:
+    for i in range(1, 4):
+        print(i)
+        output.write_band(i,rasterio.open("data/%s/%s_B%d.TIF" % (source, source, i+3)).read(1)) # Red