ashaw · June 26, 2014 20:48 · Jun 26, 2014
diff --git a/charlie-loyd.rake b/charlie-loyd.rake
@@ -0,0 +1,26 @@
+# https://www.mapbox.com/blog/processing-landsat-8/
+task :landsat_2014 => :environment do
+  dirs = Dir["#{Rails.root.to_s}/db/initial/raster/*"].reject {|q| q =~ /tar\.gz/ }
+  dirs.each do |dir|
+    scene_id = dir.split("/")[-1]
+    [4,3,2].each do |band|
+      `gdalwarp -t_srs EPSG:3857 #{dir}/#{scene_id}_B#{band}.TIF #{dir}/#{scene_id}_B#{band}-projected.tif`
+    end
+    `convert -combine #{dir}/#{scene_id}_B{4,3,2}-projected.tif #{dir}/#{scene_id}_RGB-projected.tif && \
+     convert -channel B -gamma 0.925 -channel R -gamma 1.03 -channel RGB -sigmoidal-contrast 50x16% #{dir}/#{scene_id}_RGB-projected.tif #{dir}/#{scene_id}_RGB-projected-corrected.tif && \
+     convert -depth 8 #{dir}/#{scene_id}_RGB-projected-corrected.tif  #{dir}/#{scene_id}_RGB-projected-corrected-8bit.tif && \
+     listgeo -tfw #{dir}/#{scene_id}_B4-projected.tif && \
+     mv #{dir}/#{scene_id}_B4-projected.tfw #{dir}/#{scene_id}_RGB-projected-corrected-8bit.tfw && \
+     gdal_edit.py -a_srs EPSG:3857 #{dir}/#{scene_id}_RGB-projected-corrected-8bit.tif && \
+     gdal_translate -a_nodata 0 #{dir}/#{scene_id}_RGB-projected-corrected-8bit.tif #{dir}/#{scene_id}_RGB-projected-corrected-8bit-nodata.tif`
+  end
+end
+
+task :landsat_2014_merge => :environment do
+  dirs = Dir["#{Rails.root.to_s}/db/initial/raster/*"].reject {|q| q =~ /\..*$/ }
+  scenes = dirs.map do |q|
+    scene_id = q.split("/")[-1]
+    corrected = "#{q}/#{scene_id}_RGB-projected-corrected-8bit-nodata.tif"
+  end
+  `cd #{Rails.root.to_s}/db/initial/raster/ && gdalwarp --config GDAL_CACHEMAX 3000 -wm 3000 #{scenes.join(" ")} #{Rails.root.to_s}/db/initial/raster/merged.tif`
+end