Hi GEONETCasters!

The following GOES-16 RGB Composites (GeoTIFF format) were added to the GEONETCast-Americas broadcast today, June 4th, 2019, 11:20 UTC:

• Day Land Cloud
• Natural True Color
• Airmass
• Day Microphysics
• Night Microphysics
• Day Cloud Phase Distinction
• Cloud Phase
• Day Convection
• Dust

Please find below a sample of each RGB, and their Quick Guides:

— DAY LAND CLOUD RGB —

Naming convention: G16_DLCREG_YYYYMMDDHHMN.tif

Quick Guide (click on the image to access):

— NATURAL TRUE COLOR RGB —

Naming convention: G16_NTCREG_YYYYMMDDHHMN.tif

Quick Guide (click on the image to access):

— AIRMASS RGB —

Naming convention: G16_ARMREG_YYYYMMDDHHMN.tif

Quick Guide (click on the image to access):

— DAY MICROPHYSICS RGB —

Naming convention: G16_DMPREG_YYYYMMDDHHMN.tif

Quick Guide (click on the image to access):

— NIGHT MICROPHYSICS RGB —

Naming convention: G16_NMPREG_YYYYMMDDHHMN.tif

Quick Guide (click on the image to access):

— DAY CLOUD PHASE DISTINCTION RGB —

Naming convention: G16_DCPREG_YYYYMMDDHHMN.tif

Quick Guide (click on the image to access):

— CLOUD PHASE RGB —

Naming convention: G16_CLPREG_YYYYMMDDHHMN.tif

Quick Guide (click on the image to access):

— DAY CONVECTION RGB —

Naming convention: G16_CONREG_YYYYMMDDHHMN.tif

Quick Guide (click on the image to access):

— DUST RGB —

Naming convention: G16_DSTREG_YYYYMMDDHHMN.tif

Quick Guide (click on the image to access):

— GNC-A RECEIVE FOLDER —

The RGB’s may be found on the new GOES-R-RGB-Composites GNC-A ingestion folder:

— NAMING CONVENTION AND SECTORS —

The RGB’s naming convention is as follows:

G16_RGBREG_YYYYMMDDHHMN.tif

Where:

YYYYMMDDHHMN: Year, Month, Day, Hour and Minutes (UTC)

RGB: 

• DLC: Day Land Cloud
• NTC: Natural True Color
• ARM: Airmass
• DMP: Day Microphysics
• NMP: Night Microphysics
• DCP: Day Cloud Phase Distinction
• CLP: Cloud Phase
• CON: Day Convection
• DST: Dust

In order to fit the RGB’s in the GNC-A available space, for each time frame, the RGB’s are divided in 8 parts: A Low Resolution Full Disk and 7 higher resolution regions (3 km), as below:

REG: 

• FDK: Full Disk (Low Resolution)
• S01: Sector 01
• S02: Sector 02
• S03: Sector 03
• S04: Sector 04
• S05: Sector 05
• S06: Sector 06
• S07: Sector 07
• S08: Sector 08

The table below shows the extent of each sector:

The image below shows the 8 parts of each RGB:

The RGB’s are broadcasted for the following minutes: 00, 20, 30 and 50, each hour.

— JOINING THE SECTORS AND SUBSECT A GIVEN REGION —

Users may join the 8 regions into a single GeoTIFF and subsect it for a given region using a single GDAL instruction, called “gdalwarp”:

gdalwarp FDK.tif* S01.tif* …… S07.tif output.tif -te min_lon min_lat max_lon max_lat -overwrite  

Below, an example instruction to join all the Natural True Color RGB sectors, and cut it to the following region:

• min lon: -120
• min lat: 0
• max lon: -75
• max lat: 35

gdalwarp G16_NTCFDK_201905311400.tif* G16_NTCS01_201905311400.tif* G16_NTCS02_201905311400.tif* G16_NTCS03_201905311400.tif* G16_NTCS04_201905311400.tif* G16_NTCS05_201905311400.tif* G16_NTCS06_201905311400.tif* G16_NTCS07_201905311400.tif* G16_NTC_201905311400.tif -te -120 0 -75 35 -overwrite

The result GeoTIFF, called G16_NTC_201905311400.tif, is seen below:

Alternatively, you may use the following Python code:

# Required libraries
from osgeo import gdal, osr, ogr  # Import GDAL
import os 			  # Miscellaneous operating system interfaces

# Desired Extent
extent = [-120, 0, -75, 35]

# Define KM_PER_DEGREE
KM_PER_DEGREE = 111.32

# Calculate the total number of degrees in lat and lon
deg_lon = extent[2] - extent[0]
deg_lat = extent[3] - extent[1]

# Calculate the number of pixels
resolution = 2
width = (KM_PER_DEGREE * deg_lon) /  resolution
height = (KM_PER_DEGREE * deg_lat) /  resolution

# Create the mosaic
grid = gdal.BuildVRT('tmp.vrt', ['s8.tif', 's1.tif', 's2.tif', 's3.tif', 's4.tif', 's5.tif', 's6.tif', 's7.tif'])
ds = gdal.Translate('output.tif', 'tmp.vrt', projWin = [extent[0], extent[3], extent[2], extent[1]], width = width, height = height)
grid = None
ds = None