Index of /scripts/goesr

Name	Last modified	Size

Parent Directory		-
strip_header	2020-05-04 00:09	43
composite_goesr	2026-10-16 12:32	10K

Compositing GOES-16 Tiles

It starts with compositing the tiles. I have written a simple Python script using the netCDF module that creates a single netCDF image file from the tiles. Here are the LDM entries:

NOTHER  TI([RSTU]...) KNES ([0-3][0-9])([0-2][0-9])([0-5][0-9]) (...)
        PIPE    -close /home/ldm/scripts/strip_header /data/noaaport/satraw/(\2:yyyy)(\2:mm)\2_\3\4_g16_\1_\5.nc
# Composite vis-blue
NOTHER  TI(R[ES]01) KNES ([0-3][0-9])([0-2][0-9])([0-5][0-9]) (PAA)
        EXEC    /home/ldm/scripts/composite_goes16 /data/noaaport/satraw/(\2:yyyy)(\2:mm)\2_\3\4_g16_\1_\5.nc
# Composite vis-red
NOTHER  TI(R[ES]02) KNES ([0-3][0-9])([0-2][0-9])([0-5][0-9]) (PAA)
        EXEC    /home/ldm/scripts/composite_goes16 /data/noaaport/satraw/(\2:yyyy)(\2:mm)\2_\3\4_g16_\1_\5.nc
# Composite wv
NOTHER  TI(R[ES]09) KNES ([0-3][0-9])([0-2][0-9])([0-5][0-9]) (PAA)
        EXEC    /home/ldm/scripts/composite_goes16 /data/noaaport/satraw/(\2:yyyy)(\2:mm)\2_\3\4_g16_\1_\5.nc
# Composite ir
NOTHER  TI(R[ES]14) KNES ([0-3][0-9])([0-2][0-9])([0-5][0-9]) (PAA)
        EXEC    /home/ldm/scripts/composite_goes16 /data/noaaport/satraw/(\2:yyyy)(\2:mm)\2_\3\4_g16_\1_\5.nc

The strip header script removes the WMO header from the files so they are true netCDF files. The netCDF library can't handle a netCDF that's embedded inside other data.

Once it's written to disk, you can start the composite process. This Python script determines the number of tiles from the PAA (first tile) and then slowly adds tiles into the composite file as they come in. When all the tiles are received, it writes the full file. It will also time out after about 5 minutes and write what it has seen so far.

You can see from this that I'm writing out the full disk and CONUS sectors and compositing the two vis, IR and WV images.

The satplot program reads the netCDF file and converts it into a 8 bit WXP image for processing. WXP has the concept of a valued image which means there is, for IR images, a temperature value attached to each pixel value 0-255. You can then pass a color_fill file to satplot to enhance the image.

satplot -if=sat_ir_g16 -cof=sat_mb.cfl

The file type can be specified in the name convention file:

sat_vib_g16   %S/g16/%Y%m%d/%Y%m%d_%h%n_g16_us_vib.nc -    850000,cdf min=45
sat_vis_g16   %S/g16/%Y%m%d/%Y%m%d_%h%n_g16_us_vir.nc -    850000,cdf min=45
sat_ir_g16    %S/g16/%Y%m%d/%Y%m%d_%h%n_g16_us_ir.nc -     850000,cdf min=45
sat_wv_g16    %S/g16/%Y%m%d/%Y%m%d_%h%n_g16_us_wv.nc -     850000,cdf min=45

Otherwise, you'll have to add a -in=cdf to the command line. Here are a few samples: