Bulk download Ocean Color images from ESA and NASA platforms.

getOC is a python utility simplifying bulk download of Ocean Color images from NASA Ocean Data and ESA COPERNICUS Dataspace APIs. getOC takes as input a csv file with the list of positions and dates of interest (ex: test.csv) and it downloads images for the selected instrument, processing level, and product type. Depending on the options selected an EARTHDATA or a [COPERNICUS dataspace] (https://identity.dataspace.copernicus.eu/auth/realms/CDSE/login-actions/registration?client_id=cdse-public&tab_id=G8kjxyKCxI8) account is required. A strong internet connection is recommended.

Synthax for ESA satellite collections can be found here (example for Sentinel-3 collection)

DEPRECATED: Previous getOC versions used the Creodias API but the Creodias v2 API required 2 step authentications preventing automatic downloads. ESA satellites (+Landsat over europe) are now downloaded using the new Copernicus API

If you only need to download images in a given time frame without specifying a position, the File Search utility from NASA might be a better tool. The NASA utility also provides a few wget examples on this Ocean Color Forum Post.

Resolutions are indicative for Level 1 and 2 data, Level 3 is only available at 4 or 9 km for NASA satellites.

Typical call from bash:

python getOC.py -i VIIRSJ1 -l L2 <csv_filename> -u <earthdata-username> -w --box 60

• -i INSTRUMENT, --instrument=INSTRUMENT: specify instrument, available options are:

  • SeaWiFS
  • MODIS-Aqua
  • MODIS-Terra
  • OCTS
  • CZCS
  • MERIS
  • VIIRSN
  • VIIRSJ1
  • HICO
  • OLCI
  • SLSTR
  • SYN
  • SRAL
  • MSI
  • GOCI
  • L5-TM
  • L7-ETM
  • L8-OLI-TIRS

• -l LEVEL, --level=LEVEL: specify processing level, available options are listed below (c.f. table above for level compatible with the instrument selected).

  • MODIS-Aqua MODIS-Terra VIIRSN VIIRSJ1 SeaWiFS OCTS CZCS:
    • L1A
    • L2
    • L3b
    • L3m
  • OLCI:
    • L1-EFR full resolution or L1_ERR reduced resolution
    • L2-WFR full resolution or L2_WRR reduced resolution [Water]
    • L2-LFR full resolution or L2_LRR reduced resolution [Land]
  • SLSTR:
    • L1-RBT
    • L2_WST [Water temperature]
    • L2_LST [Land temprature]
    • L2_AOD [Aerosol optical depth]
    • L2_FRP [Fire radiative power]
  • SRAL:
    • L1-SRA
    • L1-SRA-A
    • L1-SRA-BS
    • L2-WAT [Water altimetry]
    • L2-LAN [Land altimetry]
    • L2-LAN-HY
    • L2-LAN-LY
    • L2-LAN-SI
  • SYN:
    • L2-SYN [Reflectance and aerosol parameter over land]
    • L2-V10 [Vegetation-Like product (~VGT-S10), 10 day synthesis surface reflectance and NDVI]
    • L2-VG1 [Vegetation-Like product (~VGT-S1), 1 day synthesis surface reflectance and NDVI]
    • L2-VGP [Vegetation-Like product (~VGT-P), TOA reflectance]
    • L2-AOD [Aerosol optical depth over land and water 4.5x4.5 km resolution]
  • MERIS:
    • L1
    • L2
    • L3b
    • L3m
  • SAR:
    • L0-RAW
    • L1-GRD
    • L1-GRD-COG
    • L1-SLC
    • L2-OCN
    • L2-CARD-BS
    • L2-CARD-COH12
  • GOCI:
    • L1B
    • L2
  • MSI:
    • L1C
    • L2A

• -u USERNAME, --username=USERNAME: specify username to login to Copernicus (for OLCI, SLSTR, SLAR, SAR, MSI, L5-TM, L7-ETM and L8-OLI-TIRS instruments) or EarthData (for any other sensor). Password will be prompted when the script is executed. (TO DO: add option to use a credential file credentials.ini)

  • Earthdata login for NASA satellites
  • Copernicus dataspace login for ESA satellites

• -w, --write-image-name: getOC first query an api to retrieve the list of images to download. The output of that query can be written to a csv file. getOC can then be restarted directly from that file saving that query time.

• -r, --read-image-list: getOC loads the list previously queried and printed, to avoid querying twice the same list

• -q, --quiet: Quiet please ! getOC does not output any information relative to the download and querying of the points of interest.

• -p product
  Specify the product type to download:

  • OC [default]
  • IOP
  • SST

Options specific to a level:

• Level 1 & 2:
  • --box=BOUNDING_BOX_SIZE, --bounding-box-size=BOUNDING_BOX_SIZE: Define the size of the bounding box around the point of interest in nautical miles. getOC downloads all images that intersect with this box (must be > 0).
• Level 2 & 3:
  • -p PRODUCT, --product=PRODUCT: Specify the product type to download (only for NASA satellites) (c.f. table above for level compatible with the instrument selected):
    • Level 2 (NASA):
      • OC (default)
      • IOP (deprecated)
      • SST
    • Level 3 (NASA):
      • CHL
      • POC
      • not tested: GSM, IOP, KD, LAND, PAR, PIC, QAA, RRS, and ZLEE
• Level 3: at this level the world's ocean is downloaded, getOC ignores the latitude and longitudes in the input csv file.
  • -b BINNING_PERIOD, --binning-period=BINNING_PERIOD: specify binning period (only for L3), available options are:
    • DAY (default)
    • 8D
    • MO
    • YR
  • --res=SRESOL, --spatial-resolution=SRESOL: specify spatial resolution (only for L3), available options are:
    • 4km (default)
    • 9km

Options specific to a downloading platform:

• ESA Copernicus downloads only:

  • -c CLOUD_COVER_INTERVAL, --cloud-cover=CLOUD_COVER_INTERVAL: Cloud cover interval to download (in %):
    • -c [0 100]: download all images [default]
    • -c [0 80]: download all images with less than 80% cloud cover
    • -c [0 20]: download all images with less than 20% cloud cover
    • ...

• NASA Ocean Color Level 1&2 Browser:

  • -d QUERY_DELAY, --delay=QUERY_DELAY: Delay between queries only needed to query L1L2_browser. (default=1 second)

Instruments specificity:

• VIIRS: GEO files required to process L1A files from that sensor are downloaded automatically when the level L1A is selected for this instrument.
• Level synthax for ESA satellites are specific to sensors providing different options: Use L1-EFR for OLCI full resolution data and L1-ERR for low resolution (cf table and level paragraph)

Level 1:

python getOC.py -i MODIS-Aqua -l L1A test.csv -u <earthdata-username> -w --box 60
python getOC.py -i GOCI -l L1B test.csv -u <copernicus-username> -w --box 60 
python getOC.py -i VIIRSN -l L1A test.csv -u <earthdata-username> -w --box 60
python getOC.py -i VIIRSJ1 -l L1A test.csv -u <earthdata-username> -w --box 60
python getOC.py -i MSI -l L1C test.csv -u <copernicus-username> -w --box 60
python getOC.py -i OLCI -l L1-EFR test.csv -u <copernicus-username> -w --box 60
python getOC.py -i SLSTR -l L1-WST test.csv -u <copernicus-username> -w --box 60
python getOC.py -i SRAL -l L1-RBT test.csv -u <copernicus-username> -w --box 60
python getOC.py -i L5-TM -l L1T test.csv -u <copernicus-username> -w --box 60

Level 2:

python getOC.py -i MODIS-Aqua -l L2 test.csv -u <earthdata-username> -p OC -w --box 60
python getOC.py -i GOCI -l L2 test.csv -u <copernicus-username> -w --box 60
python getOC.py -i VIIRSN -l L2 test.csv -u <earthdata-username> -p OC -w --box 60
python getOC.py -i VIIRSJ1 -l L2 test.csv -u <earthdata-username> -p OC -w --box 60
python getOC.py -i MSI -l L2A test.csv -u <copernicus-username> -w --box 60 
python getOC.py -i OLCI -l L2-WRR test.csv -u <copernicus-username> -w --box 60
python getOC.py -i SLSTR -l L2-WST test.csv -u <copernicus-username> -w --box 60
python getOC.py -i SYN -l L2-AOD test.csv -u <copernicus-username> -w --box 60
python getOC.py -i L8-OLI-TIRS -l L2 test.csv -u <copernicus-username> -w --box 60

Level 3:

python getOC.py -i MODIS-Aqua -l L3b test.csv <earthdata-username> -p CHL -b 8D --res 4km -w
python getOC.py -i MODIS-Aqua -l L3m test.csv <earthdata-username> -p CHL -b YR --res 9km -w