def create_wildfire_propagation(area,
ignition: ty.Sequence[TimedPoint],
wind_speed=1.,
wind_direction=3 * np.pi / 4,
world=THE_WORLD) -> FirePropagation:
"""Reference wildfire propagation"""
fire_env = Environment(area, wind_speed, wind_direction, world)
fire_prop = FirePropagation(fire_env)
for ig in ignition:
fire_prop.set_ignition_point(ig)
propagation_end_time = np.inf
fire_prop.propagate(propagation_end_time)
return fire_prop
fire_prop = FirePropagation(fire_env)
# now = 0
now = datetime.datetime(2020, 1, 1, 0, 0).timestamp()
ignitionpoint1 = TimedPoint(482000, 6212000.0, now)
ignitionpoint2 = TimedPoint(482000 - 250, 6212000.0 + 1800, now)
ignitionpoint3 = TimedPoint(482000 + 2000, 6212000.0 + 2000, now)
fire_prop.set_ignition_point(ignitionpoint1)
fire_prop.set_ignition_point(ignitionpoint2)
fire_prop.set_ignition_point(ignitionpoint3)
# fire_prop.set_ignition_point(TimedPoint(575000.0, 6200000.0, 0))
propagation_end_time = now + 60 * 60 * 4
export_front = (7000, 7100)
fire_prop.propagate(propagation_end_time)
fire_prop.ignitions().write_to_file("demo_propagation_highres_firemap.tif")
filt_fire = fire_prop.ignitions().clone()
filt_fire.data["ignition"][filt_fire.data["ignition"] < export_front[0]] = 0
filt_fire.data["ignition"][filt_fire.data["ignition"] > export_front[1]] = 0
filt_fire.data["ignition"][filt_fire.data["ignition"].nonzero()] = 65535
filt_fire.write_to_image_file("demo_propagation_highres_firemap_" +
str(export_front[0]) + "-" +
str(export_front[1]) + ".png",
layer_name="ignition")
fire_env.raster.write_to_file("demo_propagation_highres_dem.tif", "elevation")
# Figure terrain + ignition contour + ignition point
landcover_path=FIRERS_LANDCOVER_DATA)
seganosa_fire_env = Environment([[2799134.0, 2805134.0], [2296388.0, 2302388.0]],
3, 0, the_world)
# 6km by 6km around seganosa. Wind: 3m/s W->E
# FIRE PROPAGATION
## Fire started 6 hours ago. Predict until 2 hours in the future
fire_prop = FirePropagation(seganosa_fire_env)
four_hours_ago = (datetime.datetime.now() - datetime.timedelta(hours=6)).timestamp()
now = datetime.datetime.now().timestamp()
four_hours_from_now = (datetime.datetime.now() + datetime.timedelta(hours=2)).timestamp()
fire_start = TimedPoint(2802134.0 - 1500.0, 2299388.0, four_hours_ago)
fire_prop.set_ignition_point(fire_start)
fire_prop.propagate(until=four_hours_from_now)
## Figure terrain + ignition contour + ignition point
gdd = display.GeoDataDisplay.pyplot_figure(
seganosa_fire_env.raster.combine(fire_prop.ignitions().slice(["ignition"])),
frame=(0., 0.))
gdd.draw_elevation_shade(with_colorbar=False, cmap=matplotlib.cm.terrain)
gdd.draw_wind_quiver()
gdd.draw_ignition_contour(with_labels=True, cmap=matplotlib.cm.plasma)
gdd.draw_ignition_points(fire_start)
# gdd.figure.show()
gdd.figure.savefig(".".join(("demo_seganosa_propagation", "svg")), dpi=150, bbox_inches='tight')
# PLANNING
## Define the initial plan