Ipyleaflet

Velocity map

import os

import xarray as xr

from ipyleaflet import Map, basemaps

from ipyleaflet.velocity import Velocity

if not os.path.exists("wind-global.nc"):

url = "https://github.com/benbovy/xvelmap/raw/master/notebooks/wind-global.nc"

import requests

r = requests.get(url)

wind_data = r.content

with open("wind-global.nc", "wb") as f:

f.write(wind_data)

center = [0, 0]

zoom = 2

m = Map(

center=center,

zoom=zoom,

interpolation="nearest",

basemap=basemaps.CartoDB.DarkMatter,

)

ds = xr.open_dataset("wind-global.nc")

display_options = {

"velocityType": "Global Wind",

"displayPosition": "bottomleft",

"displayEmptyString": "No wind data",

}

wind = Velocity(

data=ds,

zonal_speed="u_wind",

meridional_speed="v_wind",

latitude_dimension="lat",

longitude_dimension="lon",

velocity_scale=0.01,

max_velocity=20,

display_options=display_options,

)

m.add_layer(wind)

m

Split map

from ipyleaflet import Map, basemaps, basemap_to_tiles, SplitMapControl

m = Map(center=(42.6824, 365.581), zoom=5)

right_layer = basemap_to_tiles(basemaps.NASAGIBS.ModisTerraTrueColorCR, "2017-11-11")

left_layer = basemap_to_tiles(basemaps.NASAGIBS.ModisAquaBands721CR, "2017-11-11")

control = SplitMapControl(left_layer=left_layer, right_layer=right_layer)

m.add_control(control)

m