VIIRS Nighttime Lights 2012 processing

README.md

This script processes VIIRS Nighttime Lights 2012 GeoTIFFs from the Earth Observation Group at NOAA National Geophysical Data Center to prepare them for rendering in TileMill and uploading to MapBox Hosting.

Read Lights of the Night on MapBox to learn more about NPP the functions of this script.

process.sh

#!/bin/bash
FILES=('FILE2' 'FILE1');
ADDO="2 4 8 16 32 64 128 256 512 1024 2048"
LUTS=('15:0,50:255' '20:0,75:255' '5:0,200:255');
RELPATH="PATH/TO/GEOTIFFS"
for FILE in ${FILES[@]}; do
  gdal_translate \
     -scale 0 50 \
     -ot Byte \
     $FILE \
     ${FILE%.tif}_scale.tif
  gdalwarp \
     -wm 2000 \
     -multi \
     -srcnodata 0 \
     -dstnodata 0 \
     -r lanczos \
     -co COMPRESS=LZW \
     -wo OPTIMIZE_SIZE=TRUE\
     -t_srs EPSG:3785 \
     ${FILE%.tif}_scale.tif \
     ${FILE%.tif}_3785.tif
  gdaladdo \
     -r cubic \
     --config COMPRESS_OVERVIEW LZW \
     ${FILE%.tif}_3785.tif \
     $ADDO 
  rm ${FILE%.tif}_scale.tif
done
gdalbuildvrt \
   -srcnodata 0 \
   -vrtnodata 0 \
   red.vrt \
   $RELPATH/nasa/*_3785.tif
gdalbuildvrt \
   -srcnodata 0 \
   -vrtnodata 0 \
   blue.vrt \
   $RELPATH/nasa/*_3785.tif
gdalbuildvrt \
   -srcnodata 0 \
   -vrtnodata 0 \
   green.vrt \
   $RELPATH/nasa/*_3785.tif
gdalbuildvrt \
   -separate \
   -srcnodata "0 0 0" \
   -vrtnodata "0 0 0" \
   -addalpha rgb.vrt \
   $RELPATH/nasa/red.vrt \
   $RELPATH/nasa/green.vrt \
   $RELPATH/nasa/blue.vrt 
cat rgb.vrt | sed -E "8s_(<SourceBand>1</SourceBand>)_\1<LUT>${LUTS[0]}</LUT>_1;19s_(<SourceBand>1</SourceBand>)_\1<LUT>${LUTS[1]}</LUT>_1;30s_(<SourceBand>1</SourceBand>)_\1<LUT>${LUTS[2]}</LUT>_1" |
   sed -E "s/SimpleSource/ComplexSource/g;s_(<VRTRasterBand dataType=\"Byte\" band=\"1\">)_\1\
    <ColorInterp>Red</ColorInterp>_g;s_(<VRTRasterBand dataType=\"Byte\" band=\"2\">)_\1\
    <ColorInterp>Green</ColorInterp>_g;s_(<VRTRasterBand dataType=\"Byte\" band=\"3\">)_\1\
    <ColorInterp>Blue</ColorInterp>_g;" > tmp.vrt
mv -f tmp.vrt rgb.vrt
gdalwarp \
   -dstalpha \
   -srcnodata "0 0 0" \
   -dstnodata "0 0 0" \
   -co COMPRESS=LZW \
   -wm 2000 \
   -multi \
   rgb.vrt \
   rgb.tif
 gdaladdo \
   -ro \
   -r gauss \
   rgb.tif \
   $ADDO

New global mosaics are available here: http://ngdc.noaa.gov/eog/viirs/download_monthly.html

@fractalien the link is also broken :/

New link: https://ngdc.noaa.gov/eog/viirs/download_dnb_composites.html

What a waste of my time to write the damn code