81 goes timelapse

Uncomment the following line to install geemap if needed.

This notebook is adapted from Justin Braaten's JavaScript example. Credits to Justin Braaten. See also his Medium post - GOES in Earth Engine.

In [ ]:

import os
import ee
import geemap

import os import ee import geemap

In [ ]:

# geemap.update_package()

# geemap.update_package()

In [ ]:

m = geemap.ee_initialize()

m = geemap.ee_initialize()

In [ ]:

region = ee.Geometry.Polygon(
    [
        [
            [-159.5954379282731, 60.40883060191719],
            [-159.5954379282731, 24.517881970830725],
            [-114.2438754282731, 24.517881970830725],
            [-114.2438754282731, 60.40883060191719],
        ]
    ],
    None,
    False,
)

region = ee.Geometry.Polygon( [ [ [-159.5954379282731, 60.40883060191719], [-159.5954379282731, 24.517881970830725], [-114.2438754282731, 24.517881970830725], [-114.2438754282731, 60.40883060191719], ] ], None, False, )

In [ ]:

start_date = "2021-10-24T14:00:00"
end_date = "2021-10-25T01:00:00"
data = "GOES-17"
scan = "full_disk"

start_date = "2021-10-24T14:00:00" end_date = "2021-10-25T01:00:00" data = "GOES-17" scan = "full_disk"

In [ ]:

col = geemap.goes_timeseries(start_date, end_date, data, scan, region)

col = geemap.goes_timeseries(start_date, end_date, data, scan, region)

In [ ]:

visParams = {
    'bands': ['CMI_C02', 'CMI_GREEN', 'CMI_C01'],
    'min': 0,
    'max': 0.8,
    'dimensions': 700,
    'framesPerSecond': 9,
    'region': region,
    'crs': col.first().projection(),
}

visParams = { 'bands': ['CMI_C02', 'CMI_GREEN', 'CMI_C01'], 'min': 0, 'max': 0.8, 'dimensions': 700, 'framesPerSecond': 9, 'region': region, 'crs': col.first().projection(), }

In [ ]:

out_dir = os.path.expanduser("~/Downloads")
out_gif = os.path.join(out_dir, "goes_timelapse.gif")
if not os.path.exists(out_dir):
    os.makedirs(out_dir)

out_dir = os.path.expanduser("~/Downloads") out_gif = os.path.join(out_dir, "goes_timelapse.gif") if not os.path.exists(out_dir): os.makedirs(out_dir)

In [ ]:

geemap.download_ee_video(col, visParams, out_gif)

geemap.download_ee_video(col, visParams, out_gif)

In [ ]:

timestamps = geemap.image_dates(col, date_format='YYYY-MM-dd HH:mm').getInfo()

timestamps = geemap.image_dates(col, date_format='YYYY-MM-dd HH:mm').getInfo()

In [ ]:

geemap.add_text_to_gif(
    out_gif,
    out_gif,
    xy=('3%', '3%'),
    text_sequence=timestamps,
    font_size=20,
    font_color='#ffffff',
)

geemap.add_text_to_gif( out_gif, out_gif, xy=('3%', '3%'), text_sequence=timestamps, font_size=20, font_color='#ffffff', )

All the above steps can be reduced to just one line of code using the goes_timelapse() function.

In [ ]:

timelapse = geemap.goes_timelapse(
    region, out_gif, start_date, end_date, data, scan, framesPerSecond=10
)
geemap.show_image(timelapse)

timelapse = geemap.goes_timelapse( region, out_gif, start_date, end_date, data, scan, framesPerSecond=10 ) geemap.show_image(timelapse)