veloren_world/site/

mod.rs

pub mod economy;
mod generation;
pub mod genstat;
pub mod namegen;
pub mod plot;
mod tile;
pub mod util;

use self::tile::{HazardKind, KeepKind, RoofKind, TILE_SIZE, Tile, TileGrid};
pub use self::{
    economy::Economy,
    generation::{Fill, Painter, Primitive, PrimitiveRef, Structure, aabr_with_z},
    genstat::{GenStatPlotKind, GenStatSiteKind, SitesGenMeta},
    plot::{Plot, PlotKind, foreach_plot},
    tile::TileKind,
};
use crate::{
    Canvas, IndexRef, Land,
    config::CONFIG,
    sim::Path,
    util::{CARDINALS, DHashSet, Grid, SQUARE_4, SQUARE_9, attempt},
};
use common::{
    astar::Astar,
    calendar::Calendar,
    comp::Alignment,
    generation::{EntityInfo, EntitySpawn},
    lottery::Lottery,
    map::MarkerKind,
    spiral::Spiral2d,
    store::{Id, Store},
    terrain::{
        Block, BlockKind, SiteKindMeta, SpriteKind, TerrainChunkSize,
        site::{DungeonKindMeta, SettlementKindMeta},
    },
    util::Dir2,
    vol::RectVolSize,
};
use hashbrown::DefaultHashBuilder;
use namegen::NameGen;
use rand::{SeedableRng, prelude::*, seq::IndexedRandom};
use rand_chacha::{ChaCha8Rng, ChaChaRng};
use std::ops::Range;
use vek::*;

/// Seed a new RNG from an old RNG, thereby making the old RNG independent of
/// changing use of the new RNG. The practical effect of this is to reduce the
/// extent to which changes to child generation algorithm produce a 'butterfly
/// effect' on their parent generators, meaning that generators will be less
/// likely to produce entirely different outcomes if some detail of a generation
/// algorithm changes slightly. This is generally good and makes worldgen code
/// easier to maintain and less liable to breaking changes.
fn reseed<R: Rng>(rng: &mut R) -> impl Rng + use<R> {
    ChaChaRng::from_seed(rng.random::<[u8; 32]>())
}

pub struct SpawnRules {
    pub trees: bool,
    pub max_warp: f32,
    pub paths: bool,
    pub waypoints: bool,
    // (alt, total_weight, max_factor)
    // factor = 0.0: no effect
    // factor = 1.0: exact alt
    pub preferred_alt: (f32, f32, f32),
}

impl SpawnRules {
    /// Specify that the column these rules relates to should prefer to use the
    /// given `alt` as an altitude (causing the terrain to shift up or down).
    ///
    /// The `weight` parameter indicates how strong this pull should be: 0.0
    /// indicates no change, >= 1.0 indicates *exactly* the given altitude.
    ///
    /// Successive calls to this method will collect a weighted average. If you
    /// wish to strongly override the preferred altitude, a weight > 1.0 might
    /// be desirable.
    pub fn prefer_alt(&mut self, alt: f32, weight: f32) {
        self.preferred_alt.0 += alt * weight;
        self.preferred_alt.1 += weight;
        self.preferred_alt.2 = self.preferred_alt.2.max(weight);
    }

    // Returns (alt, factor)
    pub fn get_preferred_alt(&self) -> (f32, f32) {
        // + 0.1 is a hack, prevents precision issues
        (
            self.preferred_alt.0 / self.preferred_alt.1.max(0.0001) + 0.1,
            self.preferred_alt.2,
        )
    }
}

impl Default for SpawnRules {
    fn default() -> Self {
        Self {
            trees: true,
            max_warp: 1.0,
            paths: true,
            waypoints: true,
            preferred_alt: (0.0, 0.0, f32::NEG_INFINITY),
        }
    }
}

#[derive(Debug, Clone, Copy, PartialEq, Eq)]
pub enum SiteKind {
    Refactor,
    CliffTown,
    SavannahTown,
    DesertCity,
    ChapelSite,
    DwarvenMine,
    CoastalTown,
    Citadel,
    Terracotta,
    GiantTree,
    Gnarling,
    Bridge(Vec2<i32>, Vec2<i32>),
    Adlet,
    Haniwa,
    PirateHideout,
    JungleRuin,
    RockCircle,
    TrollCave,
    Camp,
    Cultist,
    Sahagin,
    VampireCastle,
    GliderCourse,
    Myrmidon,
}

impl SiteKind {
    pub fn meta(&self) -> Option<SiteKindMeta> {
        match self {
            SiteKind::Refactor => Some(SiteKindMeta::Settlement(SettlementKindMeta::Default)),
            SiteKind::CliffTown => Some(SiteKindMeta::Settlement(SettlementKindMeta::CliffTown)),
            SiteKind::SavannahTown => {
                Some(SiteKindMeta::Settlement(SettlementKindMeta::SavannahTown))
            },
            SiteKind::CoastalTown => {
                Some(SiteKindMeta::Settlement(SettlementKindMeta::CoastalTown))
            },
            SiteKind::DesertCity => Some(SiteKindMeta::Settlement(SettlementKindMeta::DesertCity)),
            SiteKind::Gnarling => Some(SiteKindMeta::Dungeon(DungeonKindMeta::Gnarling)),
            SiteKind::Adlet => Some(SiteKindMeta::Dungeon(DungeonKindMeta::Adlet)),
            SiteKind::Terracotta => Some(SiteKindMeta::Dungeon(DungeonKindMeta::Terracotta)),
            SiteKind::Haniwa => Some(SiteKindMeta::Dungeon(DungeonKindMeta::Haniwa)),
            SiteKind::Myrmidon => Some(SiteKindMeta::Dungeon(DungeonKindMeta::Myrmidon)),
            SiteKind::DwarvenMine => Some(SiteKindMeta::Dungeon(DungeonKindMeta::DwarvenMine)),
            SiteKind::ChapelSite => Some(SiteKindMeta::Dungeon(DungeonKindMeta::SeaChapel)),
            SiteKind::Cultist => Some(SiteKindMeta::Dungeon(DungeonKindMeta::Cultist)),
            SiteKind::Sahagin => Some(SiteKindMeta::Dungeon(DungeonKindMeta::Sahagin)),
            SiteKind::VampireCastle => Some(SiteKindMeta::Dungeon(DungeonKindMeta::VampireCastle)),

            _ => None,
        }
    }

    pub fn should_do_economic_simulation(&self) -> bool {
        matches!(
            self,
            SiteKind::Refactor
                | SiteKind::CliffTown
                | SiteKind::SavannahTown
                | SiteKind::CoastalTown
                | SiteKind::DesertCity
        )
    }

    pub fn marker(&self) -> Option<MarkerKind> {
        match self {
            SiteKind::Refactor
            | SiteKind::CliffTown
            | SiteKind::SavannahTown
            | SiteKind::CoastalTown
            | SiteKind::DesertCity => Some(MarkerKind::Town),
            SiteKind::Citadel => Some(MarkerKind::Castle),
            SiteKind::Bridge(_, _) => Some(MarkerKind::Bridge),
            SiteKind::GiantTree => Some(MarkerKind::Tree),
            SiteKind::Gnarling => Some(MarkerKind::Gnarling),
            SiteKind::DwarvenMine => Some(MarkerKind::DwarvenMine),
            SiteKind::ChapelSite => Some(MarkerKind::ChapelSite),
            SiteKind::Terracotta => Some(MarkerKind::Terracotta),
            SiteKind::GliderCourse => Some(MarkerKind::GliderCourse),
            SiteKind::Cultist => Some(MarkerKind::Cultist),
            SiteKind::Sahagin => Some(MarkerKind::Sahagin),
            SiteKind::Myrmidon => Some(MarkerKind::Myrmidon),
            SiteKind::Adlet => Some(MarkerKind::Adlet),
            SiteKind::Haniwa => Some(MarkerKind::Haniwa),
            SiteKind::VampireCastle => Some(MarkerKind::VampireCastle),

            SiteKind::PirateHideout
            | SiteKind::JungleRuin
            | SiteKind::RockCircle
            | SiteKind::TrollCave
            | SiteKind::Camp => None,
        }
    }
}

#[derive(Default)]
pub struct Site {
    pub origin: Vec2<i32>,
    name: Option<String>,
    // NOTE: Do we want these to be public?
    pub tiles: TileGrid,
    pub plots: Store<Plot>,
    pub plazas: Vec<Id<Plot>>,
    pub roads: Vec<Id<Plot>>,
    pub economy: Option<Box<Economy>>,
    pub kind: Option<SiteKind>,
}

impl Site {
    pub fn filter_plots<'a, F: FnMut(&'a Plot) -> bool>(
        &'a self,
        mut f: F,
    ) -> std::iter::Filter<impl ExactSizeIterator<Item = &'a Plot>, impl FnMut(&&'a Plot) -> bool>
    {
        self.plots.values().filter(move |p| f(p))
    }

    pub fn any_plot<F: FnMut(&Plot) -> bool>(&self, f: F) -> bool { self.plots.values().any(f) }

    pub fn meta(&self) -> Option<SiteKindMeta> { self.kind.and_then(|s| s.meta()) }

    pub fn economy_mut(&mut self) -> &mut Economy { self.economy.get_or_insert_default() }

    pub fn do_economic_simulation(&self) -> bool {
        self.kind.is_some_and(|s| s.should_do_economic_simulation())
    }

    pub fn trade_information(&self, id: Id<Site>) -> Option<common::trade::SiteInformation> {
        self.economy
            .as_ref()
            .map(|econ| common::trade::SiteInformation {
                id: id.id(),
                unconsumed_stock: econ.get_available_stock(),
            })
    }

    pub fn radius(&self) -> f32 {
        ((self
            .tiles
            .bounds
            .min
            .map(|e| e.abs())
            .reduce_max()
            .max(self.tiles.bounds.max.map(|e| e.abs()).reduce_max())
            // Temporary solution for giving giant_tree's leaves enough space to be painted correctly
            // TODO: This will have to be replaced by a system as described on discord :
            // https://discord.com/channels/449602562165833758/450064928720814081/937044837461536808
            + if self
                .plots
                .values()
                .any(|p| matches!(&p.kind, PlotKind::GiantTree(_)))
            {
                // 25 Seems to be big enough for the current scale of 4.0
                25
            } else {
                5
            })
            * TILE_SIZE as i32) as f32
    }

    pub fn spawn_rules(&self, spawn_rules: &mut SpawnRules, land: &Land, wpos: Vec2<i32>) {
        let tile_pos = self.wpos_tile_pos(wpos);
        let max_warp = SQUARE_9
            .iter()
            .filter_map(|rpos| {
                let tile_pos = tile_pos + rpos;
                if self.tiles.get(tile_pos).is_natural() {
                    None
                } else {
                    let clamped =
                        wpos.clamped(self.tile_wpos(tile_pos), self.tile_wpos(tile_pos + 1) - 1);
                    Some(clamped.as_::<f32>().distance_squared(wpos.as_::<f32>()))
                }
            })
            .min_by_key(|d2| *d2 as i32)
            .map(|d2| d2.sqrt() / TILE_SIZE as f32)
            .unwrap_or(1.0);
        SQUARE_9.iter().for_each(|rpos| {
            let tile_pos = tile_pos + rpos;
            let tile = self.tiles.get(tile_pos);
            let clamped = wpos.clamped(self.tile_wpos(tile_pos), self.tile_wpos(tile_pos + 1) - 1);
            let dist = clamped.as_::<f32>().distance(wpos.as_::<f32>()) / TILE_SIZE as f32;
            let weight = (1.0 - dist).clamped(0.001, 1.0);
            // Special case: roads get their altitude from the tile, not the plot
            // TODO: Should this be true of all tiles?
            if let TileKind::Road { alt, .. } = &tile.kind {
                spawn_rules.prefer_alt(*alt, weight.powi(4) * 25.0);
            } else if let TileKind::Path { closest_pos, path } = &tile.kind
                && let Some(hard_alt) = tile.hard_alt
            {
                spawn_rules.prefer_alt(
                    hard_alt as f32,
                    (1.0 - (closest_pos.distance(wpos.as_()) - path.width * 1.25).max(0.0) * 0.15)
                        .clamped(0.001, 1.0)
                        .powi(2),
                );
            } else if let Some(plot) = tile.plot {
                self.plot(plot).spawn_rules(spawn_rules, land, wpos, weight);
            }
        });

        spawn_rules.trees &= max_warp == 1.0;
        spawn_rules.max_warp = spawn_rules.max_warp.min(max_warp);
    }

    pub fn bounds(&self) -> Aabr<i32> {
        let border = 1;
        Aabr {
            min: self.tile_wpos(self.tiles.bounds.min - border),
            max: self.tile_wpos(self.tiles.bounds.max + 1 + border),
        }
    }

    pub fn plot(&self, id: Id<Plot>) -> &Plot { &self.plots[id] }

    pub fn plots(&self) -> impl ExactSizeIterator<Item = &Plot> + '_ { self.plots.values() }

    pub fn plazas(&self) -> impl ExactSizeIterator<Item = Id<Plot>> + '_ {
        self.plazas.iter().copied()
    }

    pub fn create_plot(&mut self, plot: Plot) -> Id<Plot> { self.plots.insert(plot) }

    pub fn blit_aabr(&mut self, aabr: Aabr<i32>, tile: Tile) {
        for y in 0..aabr.size().h {
            for x in 0..aabr.size().w {
                self.tiles.set(aabr.min + Vec2::new(x, y), tile.clone());
            }
        }
    }

    pub fn create_road(
        &mut self,
        land: &Land,
        a: Vec2<i32>,
        b: Vec2<i32>,
        w: u16,
        kind: plot::RoadKind,
        (src, src_alt): (Id<Plot>, f32),
        (dst, dst_alt): (Id<Plot>, f32),
    ) {
        const MAX_ITERS: usize = 4096;
        let range = &(-(w as i32) / 2..w as i32 - (w as i32 + 1) / 2);
        // Manhattan distance.
        let heuristic =
            |(tile, _): &(Vec2<i32>, Vec2<i32>)| (tile - b).map(|e| e.abs()).sum() as f32;
        let Some((path, _cost)) =
            Astar::new(MAX_ITERS, (a, Vec2::zero()), DefaultHashBuilder::default())
                .poll(
                    MAX_ITERS,
                    &heuristic,
                    |(tile, prev_dir)| {
                        let tile = *tile;
                        let prev_dir = *prev_dir;
                        let this = &self;
                        CARDINALS.iter().map(move |dir| {
                            let neighbor = (tile + *dir, *dir);

                            // Transition cost
                            let alt_a = land.get_alt_approx(this.tile_center_wpos(tile));
                            let alt_b = land.get_alt_approx(this.tile_center_wpos(neighbor.0));
                            let mut cost = 1.0
                                + (alt_a - alt_b).abs() / TILE_SIZE as f32
                                + (prev_dir != *dir) as i32 as f32;

                            for i in range.clone() {
                                let orth = dir.yx() * i;
                                let tile = this.tiles.get(neighbor.0 + orth);
                                if tile.is_obstacle() {
                                    cost += 1000.0;
                                } else if !tile.is_empty() && !tile.is_road() {
                                    cost += 25.0;
                                }
                            }

                            (neighbor, cost)
                        })
                    },
                    |(tile, _)| self.tiles.get(*tile).plot == Some(dst),
                )
                .into_path()
        else {
            return;
        };

        // Skip tiles in the src plot
        let mut path = path.nodes();
        while let Some((tile, _)) = path.first()
            && self.tiles.get(*tile).plot == Some(src)
        {
            path = &path[1..];
        }
        let Some((a, _)) = path.first().copied() else {
            return;
        };

        // Start and end altitudes, to clamp road tile altitudes relative to
        // let [alt_a, alt_b] = [a, b].map(|tile_pos|
        //     match_some!(&self.tiles.get(tile_pos).kind, TileKind::Road { alt, .. } =>
        // *alt)     .unwrap_or_else(||
        // land.get_alt_approx(self.tile_center_wpos(tile_pos))));

        let plot = self.create_plot(Plot {
            kind: PlotKind::Road(plot::Road {
                path: path.iter().map(|(tile, _)| *tile).collect(),
                kind,
            }),
            root_tile: a,
            tiles: path.iter().map(|(tile, _)| *tile).collect(),
        });

        self.roads.push(plot);

        for (i, (tile, _)) in path.iter().enumerate() {
            for y in range.clone() {
                for x in range.clone() {
                    let tile = tile + Vec2::new(x, y);
                    let old_tile = self.tiles.get(tile);
                    if matches!(
                        old_tile.kind,
                        TileKind::Empty | TileKind::Path { .. } /* | TileKind::Road { .. } */
                    ) {
                        self.tiles.set(tile, Tile {
                            kind: TileKind::Road {
                                a: i.saturating_sub(1) as u16,
                                b: (i + 1).min(path.len() - 1) as u16,
                                w,
                                alt: land.get_alt_approx(self.tile_center_wpos(tile))
                                    // Clamp roads to an altitude that means they can always reach their destination with less than a 45 degree slope
                                    .clamp(
                                        src_alt - i as f32 * tile::TILE_SIZE as f32,
                                        src_alt + i as f32 * tile::TILE_SIZE as f32,
                                    )
                                    .clamp(
                                        dst_alt - path.len().saturating_sub(i + 1) as f32 * tile::TILE_SIZE as f32,
                                        dst_alt + path.len().saturating_sub(i + 1) as f32 * tile::TILE_SIZE as f32,
                                    ),//Lerp::lerp(alt_a, alt_b, i as f32 / path.len() as f32),
                            },
                            plot: old_tile.plot.or(Some(plot)),
                            hard_alt: Some(land.get_alt_approx(self.tile_center_wpos(tile)) as i32),
                        });
                    }
                }
            }
        }
    }

    pub fn find_aabr(
        &mut self,
        search_pos: Vec2<i32>,
        area_range: Range<u32>,
        min_dims: Extent2<u32>,
    ) -> Option<(Aabr<i32>, Vec2<i32>, Vec2<i32>, Option<i32>)> {
        let ((aabr, (door_dir, hard_alt)), door_pos) =
            self.tiles.find_near(search_pos, |center, _| {
                let dir = CARDINALS
                    .iter()
                    .find(|dir| self.tiles.get(center + *dir).is_road())?;
                let hard_alt = self.tiles.get(center + *dir).hard_alt.or(self
                    .tiles
                    .get(center + *dir)
                    .plot
                    .and_then(|plot| {
                        if let PlotKind::Plaza(p) = self.plots.get(plot).kind() {
                            Some(p.hard_alt.unwrap_or(p.alt))
                        } else {
                            None
                        }
                    }));
                self.tiles
                    .grow_aabr(center, area_range.clone(), min_dims)
                    .ok()
                    .zip(Some((*dir, hard_alt)))
            })?;
        Some((aabr, door_pos, door_dir, hard_alt))
    }

    pub fn find_roadside_aabr(
        &mut self,
        rng: &mut impl Rng,
        area_range: Range<u32>,
        min_dims: Extent2<u32>,
    ) -> Option<(Aabr<i32>, Vec2<i32>, Vec2<i32>, Option<i32>)> {
        let dir = Vec2::<f32>::zero()
            .map(|_| rng.random_range(-1.0..1.0))
            .normalized();
        let search_pos = if rng.random() {
            let plot = self.plot(*self.plazas.choose(rng)?);
            let sz = plot.find_bounds().size();
            plot.root_tile + dir.map(|e: f32| e.round() as i32) * (sz + 1)
        } else if let PlotKind::Road(plot::Road { path, .. }) =
            &self.plot(*self.roads.choose(rng)?).kind
        {
            *path.nodes().choose(rng)? + (dir * 1.0).map(|e: f32| e.round() as i32)
        } else {
            return None;
        };

        let (aabr, door_pos, door_dir, hard_alt) =
            self.find_aabr(search_pos, area_range, min_dims)?;

        let alt = if let TileKind::Road { alt, .. } = &self.tiles.get(door_pos + door_dir).kind {
            Some(*alt as i32 - 1)
        } else {
            None
        };

        Some((aabr, door_pos, door_dir, hard_alt.or(alt)))
    }

    pub fn find_rural_aabr(
        &mut self,
        rng: &mut impl Rng,
        area_range: Range<u32>,
        min_dims: Extent2<u32>,
    ) -> Option<(Aabr<i32>, Vec2<i32>, Vec2<i32>, Option<i32>)> {
        // Choose a random plaza as the center of our search
        let search_center = self
            .plazas
            .choose(rng)
            .map(|&p| self.plot(p).root_tile)
            .unwrap_or_default();

        // Search in a random direction from the plaza
        let dir = Vec2::<f32>::zero()
            .map(|_| rng.random_range(-1.0..1.0))
            .normalized();
        let search_offset = dir.map2(min_dims.into(), |e: f32, sz: u32| {
            (e * sz as f32 * 0.75 + 10.0).round() as i32
        });

        self.find_aabr(search_center + search_offset, area_range, min_dims)
    }

    pub fn make_plaza_at(
        &mut self,
        land: &Land,
        index: IndexRef,
        tile_aabr: Aabr<i32>,
        rng: &mut impl Rng,
        road_kind: plot::RoadKind,
    ) -> Option<Id<Plot>> {
        let tpos = tile_aabr.center();
        let plaza_alt = land.get_alt_approx(self.tile_center_wpos(tpos)) as i32;

        let plaza = self.create_plot(Plot {
            kind: PlotKind::Plaza(plot::Plaza::generate(
                tile_aabr, road_kind, self, land, index, rng,
            )),
            root_tile: tpos,
            tiles: aabr_tiles(tile_aabr).collect(),
        });
        self.plazas.push(plaza);
        self.blit_aabr(tile_aabr, Tile {
            kind: TileKind::Road {
                a: 0,
                b: 0,
                w: 0,
                alt: plaza_alt as f32,
            },
            plot: Some(plaza),
            hard_alt: Some(plaza_alt),
        });

        let mut already_pathed = vec![];
        // One major, one minor road
        for _ in (0..rng.random_range(1.25..2.25) as u16).rev() {
            if let Some(&p) = self
                .plazas
                .iter()
                .filter(|&&p| {
                    !already_pathed.contains(&p)
                        && p != plaza
                        && already_pathed.iter().all(|&ap| {
                            (self.plot(ap).root_tile - tpos)
                                .map(|e| e as f32)
                                .normalized()
                                .dot(
                                    (self.plot(p).root_tile - tpos)
                                        .map(|e| e as f32)
                                        .normalized(),
                                )
                                < 0.0
                        })
                })
                .min_by_key(|&&p| self.plot(p).root_tile.distance_squared(tpos))
                && let PlotKind::Plaza(src_plaza) = &self.plot(p).kind
                && let PlotKind::Plaza(dst_plaza) = &self.plot(plaza).kind
            {
                self.create_road(
                    land,
                    self.plot(p).root_tile,
                    tpos,
                    2, /* + i */
                    road_kind,
                    (p, src_plaza.alt as f32),
                    (plaza, dst_plaza.alt as f32),
                );
                already_pathed.push(p);
            } else {
                break;
            }
        }

        Some(plaza)
    }

    pub fn make_plaza(
        &mut self,
        land: &Land,
        index: IndexRef,
        rng: &mut impl Rng,
        generator_stats: &mut SitesGenMeta,
        site_name: &str,
        road_kind: plot::RoadKind,
    ) -> Option<Id<Plot>> {
        generator_stats.attempt(site_name, GenStatPlotKind::Plaza);
        let plaza_radius = rng.random_range(1..3);
        let plaza_dist = 6.5 + plaza_radius as f32 * 3.0;
        let aabr = attempt(32, || {
            self.plazas
                .choose(rng)
                .map(|&p| {
                    self.plot(p).root_tile
                        + (Vec2::new(rng.random_range(-1.0..1.0), rng.random_range(-1.0..1.0))
                            .normalized()
                            * plaza_dist)
                            .map(|e| e as i32)
                })
                .or_else(|| Some(Vec2::zero()))
                .map(|center_tile| Aabr {
                    min: center_tile + Vec2::broadcast(-plaza_radius),
                    max: center_tile + Vec2::broadcast(plaza_radius + 1),
                })
                .filter(|&aabr| {
                    rng.random_range(0..48) > aabr.center().map(|e| e.abs()).reduce_max()
                        && aabr_tiles(aabr).all(|tile| !self.tiles.get(tile).is_obstacle())
                })
                .filter(|&aabr| {
                    self.plazas.iter().all(|&p| {
                        let dist_sqr = if let PlotKind::Plaza(plaza) = &self.plot(p).kind {
                            let intersection = plaza.aabr.intersection(aabr);
                            // If the size of the intersection is greater than zero they intersect
                            // on that axis and the distance on that axis is 0.
                            intersection
                                .size()
                                .map(|e| e.min(0) as f32)
                                .magnitude_squared()
                        } else {
                            let r = self.plot(p).root_tile();
                            let closest_point = aabr.projected_point(r);
                            closest_point.as_::<f32>().distance_squared(r.as_::<f32>())
                        };
                        dist_sqr > (plaza_dist * 0.85).powi(2)
                    })
                })
        })?;
        generator_stats.success(site_name, GenStatPlotKind::Plaza);
        self.make_plaza_at(land, index, aabr, rng, road_kind)
    }

    pub fn demarcate_obstacles(&mut self, land: &Land) {
        const SEARCH_RADIUS: u32 = 96;

        Spiral2d::new()
            .take((SEARCH_RADIUS * 2 + 1).pow(2) as usize)
            .for_each(|tile| {
                let wpos = self.tile_center_wpos(tile);
                if let Some(kind) = Spiral2d::new()
                    .take(9)
                    .find_map(|rpos| wpos_is_hazard(land, wpos + rpos))
                {
                    for &rpos in &SQUARE_4 {
                        // `get_mut` doesn't increase generation bounds
                        self.tiles
                            .get_mut(tile - rpos - 1)
                            .filter(|tile| tile.is_natural())
                            .map(|tile| tile.kind = TileKind::Hazard(kind));
                    }
                }
                if let Some((_, path_wpos, path, _)) = land.get_nearest_path(wpos) {
                    let tile_aabr = Aabr {
                        min: self.tile_wpos(tile),
                        max: self.tile_wpos(tile + 1) - 1,
                    };

                    if tile_aabr
                        .projected_point(path_wpos.as_())
                        .as_()
                        .distance_squared(path_wpos)
                        < path.width.powi(2)
                    {
                        self.tiles
                            .get_mut(tile)
                            .filter(|tile| tile.is_natural())
                            .map(|tile| {
                                tile.kind = TileKind::Path {
                                    closest_pos: path_wpos,
                                    path,
                                };
                                tile.hard_alt = Some(land.get_alt_approx(path_wpos.as_()) as i32);
                            });
                    }
                }
            });
    }

    /// The find_roadside_aabr function wants to have an existing plaza or road.
    /// This function is used to find a suitable location for the first plaza in
    /// a town, which has the side-effect of creating at least one road.
    /// This function is more expensive than the make_plaza function but
    /// fails to find a plaza location only if there are no suitable
    /// locations within the entire search radius.
    ///
    /// It works by exhaustively finding all tiles within a ring pattern around
    /// the town center where the tile and all surrounding tiles to the
    /// plaza radius are not hazards or roads. It then chooses the tile with
    /// the minimum distance from the town center as the plaza location. See
    /// the comments in common/src/spiral.rs for more information on the spiral
    /// ring pattern.
    ///
    /// demarcate_obstacles() should be called before this function to mark the
    /// obstacles and roads. (Cliff Towns are an exception).
    pub fn make_initial_plaza(
        &mut self,
        land: &Land,
        index: IndexRef,
        rng: &mut impl Rng,
        plaza_radius: u32,
        search_inner_radius: u32,
        search_width: u32,
        generator_stats: &mut SitesGenMeta,
        site_name: &str,
        road_kind: plot::RoadKind,
    ) -> Option<Id<Plot>> {
        generator_stats.attempt(site_name, GenStatPlotKind::InitialPlaza);
        // Find all the suitable locations for a plaza.
        let mut plaza_locations = vec![];
        // Search over a spiral ring pattern
        Spiral2d::with_ring(search_inner_radius, search_width).for_each(|tpos| {
            let aabr = Aabr {
                min: tpos - Vec2::broadcast(plaza_radius as i32),
                max: tpos + Vec2::broadcast(plaza_radius as i32 + 1),
            };
            // if all the tiles in the proposed plaza location are also not hazards or roads
            // then add the tile as a candidate for a plaza location
            if aabr_tiles(aabr).all(|tpos| self.tiles.get(tpos).is_empty()) {
                plaza_locations.push(aabr);
            }
        });
        if plaza_locations.is_empty() {
            // No suitable plaza locations were found, it's unlikely that the town will be
            // able to be generated, but we can try to make a plaza anyway with
            // the original make_plaza function.
            self.make_plaza(land, index, rng, generator_stats, site_name, road_kind)
        } else {
            // Choose the minimum distance from the town center.
            plaza_locations.sort_by_key(|&aabr| {
                aabr.min
                    .map2(aabr.max, |a, b| a.abs().min(b.abs()))
                    .magnitude_squared()
            });
            // use the first plaza location as the plaza position
            let aabr = plaza_locations.first()?;
            generator_stats.success(site_name, GenStatPlotKind::InitialPlaza);
            self.make_plaza_at(land, index, *aabr, rng, road_kind)
        }
    }

    /// This is make_initial_plaza with default options/parameters. This calls
    /// make_initial_plaza with the default parameters for the plaza_radius
    /// and search_inner_radius. The plaza_radius will be in the range 1-3,
    /// and the search_inner_radius will be 7 + plaza_radius. The search_width
    /// will be PLAZA_MAX_SEARCH_RADIUS - search_inner_radius. The
    /// search_inner_radius is approximately the same distance
    /// from the center of town as for the original make_plaza function, so this
    /// function will place the initial plaza and roads near where the
    /// original make_plaza function would place them in the case where the site
    /// is clear of hazards.
    ///
    /// This default plaza generation function is used for generating cities,
    /// cliff towns, savannah towns, and coastal towns. The other town types
    /// (terracotta, myrmidon, desert city) have a central feature so they
    /// use specific plaza generation parameters and call the make_initial_plaza
    /// function directly.
    ///
    /// demarcate_obstacles() should be called before this function to mark the
    /// obstacles and roads.
    pub fn make_initial_plaza_default(
        &mut self,
        land: &Land,
        index: IndexRef,
        rng: &mut impl Rng,
        generator_stats: &mut SitesGenMeta,
        site_name: &str,
        road_kind: plot::RoadKind,
    ) -> Option<Id<Plot>> {
        // The plaza radius can be 1, 2, or 3.
        let plaza_radius = rng.random_range(1..3);
        // look for plaza locations within a ring with an outer dimension
        // of 24 tiles and an inner dimension that will offset the plaza from the town
        // center.
        let search_inner_radius = 7 + plaza_radius;
        const PLAZA_MAX_SEARCH_RADIUS: u32 = 24;
        self.make_initial_plaza(
            land,
            index,
            rng,
            plaza_radius,
            search_inner_radius,
            PLAZA_MAX_SEARCH_RADIUS - search_inner_radius,
            generator_stats,
            site_name,
            road_kind,
        )
    }

    pub fn name(&self) -> Option<&str> { self.name.as_deref() }

    pub fn generate_mine(land: &Land, rng: &mut impl Rng, origin: Vec2<i32>) -> Self {
        let mut rng = reseed(rng);
        let mut site = Site {
            origin,
            kind: Some(SiteKind::DwarvenMine),
            ..Site::default()
        };

        let size = 60.0;

        let aabr = Aabr {
            min: Vec2::broadcast(-size as i32),
            max: Vec2::broadcast(size as i32),
        };

        let wpos: Vec2<i32> = [1, 2].into();

        let dwarven_mine =
            plot::DwarvenMine::generate(land, &mut reseed(&mut rng), &site, wpos, aabr);
        site.name = Some(dwarven_mine.name().to_string());
        let plot = site.create_plot(Plot {
            kind: PlotKind::DwarvenMine(dwarven_mine),
            root_tile: aabr.center(),
            tiles: aabr_tiles(aabr).collect(),
        });

        site.blit_aabr(aabr, Tile {
            kind: TileKind::Empty,
            plot: Some(plot),
            hard_alt: Some(1_i32),
        });

        site
    }

    pub fn generate_citadel(land: &Land, rng: &mut impl Rng, origin: Vec2<i32>) -> Self {
        let mut rng = reseed(rng);
        let mut site = Site {
            origin,
            kind: Some(SiteKind::Citadel),
            ..Site::default()
        };
        site.demarcate_obstacles(land);
        let citadel = plot::Citadel::generate(origin, land, &mut rng);
        site.name = Some(citadel.name().to_string());
        let size = citadel.radius() / tile::TILE_SIZE as i32;
        let aabr = Aabr {
            min: Vec2::broadcast(-size),
            max: Vec2::broadcast(size),
        };
        let plot = site.create_plot(Plot {
            kind: PlotKind::Citadel(citadel),
            root_tile: aabr.center(),
            tiles: aabr_tiles(aabr).collect(),
        });
        site.blit_aabr(aabr, Tile {
            kind: TileKind::Building,
            plot: Some(plot),
            hard_alt: None,
        });
        site
    }

    pub fn generate_gnarling(land: &Land, rng: &mut impl Rng, origin: Vec2<i32>) -> Self {
        let mut rng = reseed(rng);
        let mut site = Site {
            origin,
            kind: Some(SiteKind::Gnarling),
            ..Site::default()
        };
        site.demarcate_obstacles(land);
        let gnarling_fortification = plot::GnarlingFortification::generate(origin, land, &mut rng);
        site.name = Some(gnarling_fortification.name().to_string());
        let size = gnarling_fortification.radius() / TILE_SIZE as i32;
        let aabr = Aabr {
            min: Vec2::broadcast(-size),
            max: Vec2::broadcast(size),
        };
        let plot = site.create_plot(Plot {
            kind: PlotKind::Gnarling(gnarling_fortification),
            root_tile: aabr.center(),
            tiles: aabr_tiles(aabr).collect(),
        });
        site.blit_aabr(aabr, Tile {
            kind: TileKind::GnarlingFortification,
            plot: Some(plot),
            hard_alt: None,
        });
        site
    }

    pub fn generate_adlet(
        land: &Land,
        rng: &mut impl Rng,
        origin: Vec2<i32>,
        index: IndexRef,
    ) -> Self {
        let mut rng = reseed(rng);
        let mut site = Site {
            origin,
            kind: Some(SiteKind::Adlet),
            ..Site::default()
        };
        site.demarcate_obstacles(land);
        let adlet_stronghold = plot::AdletStronghold::generate(origin, land, &mut rng, index);
        site.name = Some(adlet_stronghold.name().to_string());
        let (cavern_aabr, wall_aabr) = adlet_stronghold.plot_tiles(origin);
        let plot = site.create_plot(Plot {
            kind: PlotKind::Adlet(adlet_stronghold),
            root_tile: cavern_aabr.center(),
            tiles: aabr_tiles(cavern_aabr)
                .chain(aabr_tiles(wall_aabr))
                .collect(),
        });
        site.blit_aabr(cavern_aabr, Tile {
            kind: TileKind::AdletStronghold,
            plot: Some(plot),
            hard_alt: None,
        });
        site.blit_aabr(wall_aabr, Tile {
            kind: TileKind::AdletStronghold,
            plot: Some(plot),
            hard_alt: None,
        });
        site
    }

    pub fn generate_terracotta(
        land: &Land,
        index: IndexRef,
        rng: &mut impl Rng,
        origin: Vec2<i32>,
        generator_stats: &mut SitesGenMeta,
    ) -> Self {
        let mut rng = reseed(rng);
        let gen_name = NameGen::location(&mut rng).generate_terracotta();
        let suffix = [
            "Tombs",
            "Necropolis",
            "Ruins",
            "Mausoleum",
            "Cemetery",
            "Burial Grounds",
            "Remains",
            "Temples",
            "Gardens",
        ]
        .choose(&mut rng)
        .unwrap();
        let name = match rng.random_range(0..2) {
            0 => format!("{} {}", gen_name, suffix),
            _ => format!("{} of {}", suffix, gen_name),
        };
        let mut site = Site {
            origin,
            name: Some(name.clone()),
            kind: Some(SiteKind::Terracotta),
            ..Site::default()
        };

        // place the initial plaza
        site.demarcate_obstacles(land);
        // The terracotta_palace is 15 tiles in radius, so the plaza should be outside
        // the palace.
        const TERRACOTTA_PLAZA_RADIUS: u32 = 3;
        const TERRACOTTA_PLAZA_SEARCH_INNER: u32 = 17;
        const TERRACOTTA_PLAZA_SEARCH_WIDTH: u32 = 12;
        generator_stats.add(site.name(), GenStatSiteKind::Terracotta);
        let road_kind = plot::RoadKind {
            lights: plot::RoadLights::Terracotta,
            material: plot::RoadMaterial::Sandstone,
        };
        site.make_initial_plaza(
            land,
            index,
            &mut rng,
            TERRACOTTA_PLAZA_RADIUS,
            TERRACOTTA_PLAZA_SEARCH_INNER,
            TERRACOTTA_PLAZA_SEARCH_WIDTH,
            generator_stats,
            &name,
            road_kind,
        );

        let size = 15.0 as i32;
        let aabr = Aabr {
            min: Vec2::broadcast(-size),
            max: Vec2::broadcast(size),
        };
        {
            let terracotta_palace =
                plot::TerracottaPalace::generate(land, &mut reseed(&mut rng), &site, aabr);
            let terracotta_palace_alt = terracotta_palace.alt;
            let plot = site.create_plot(Plot {
                kind: PlotKind::TerracottaPalace(terracotta_palace),
                root_tile: aabr.center(),
                tiles: aabr_tiles(aabr).collect(),
            });

            site.blit_aabr(aabr, Tile {
                kind: TileKind::Building,
                plot: Some(plot),
                hard_alt: Some(terracotta_palace_alt),
            });
        }
        let build_chance = Lottery::from(vec![(12.0, 1), (4.0, 2)]);
        for _ in 0..16 {
            match *build_chance.choose_seeded(rng.random()) {
                1 => {
                    // TerracottaHouse
                    generator_stats.attempt(site.name(), GenStatPlotKind::House);
                    let size = (9.0 + rng.random::<f32>().powf(5.0) * 1.5).round() as u32;
                    if let Some((aabr, _, _, alt)) = attempt(32, || {
                        site.find_roadside_aabr(
                            &mut rng,
                            9..(size + 1).pow(2),
                            Extent2::broadcast(size),
                        )
                    }) {
                        let terracotta_house = plot::TerracottaHouse::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            aabr,
                            alt,
                        );
                        let terracotta_house_alt = terracotta_house.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::TerracottaHouse(terracotta_house),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(terracotta_house_alt),
                        });

                        generator_stats.success(site.name(), GenStatPlotKind::House);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },

                2 => {
                    // TerracottaYard
                    generator_stats.attempt(site.name(), GenStatPlotKind::Yard);
                    let size = (9.0 + rng.random::<f32>().powf(5.0) * 1.5).round() as u32;
                    if let Some((aabr, _, _, alt)) = attempt(32, || {
                        site.find_roadside_aabr(
                            &mut rng,
                            9..(size + 1).pow(2),
                            Extent2::broadcast(size),
                        )
                    }) {
                        let terracotta_yard = plot::TerracottaYard::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            aabr,
                            alt,
                        );
                        let terracotta_yard_alt = terracotta_yard.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::TerracottaYard(terracotta_yard),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(terracotta_yard_alt),
                        });

                        generator_stats.success(site.name(), GenStatPlotKind::Yard);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },
                _ => {},
            }
        }
        site
    }

    pub fn generate_myrmidon(
        land: &Land,
        index: IndexRef,
        rng: &mut impl Rng,
        origin: Vec2<i32>,
        generator_stats: &mut SitesGenMeta,
    ) -> Self {
        let mut rng = reseed(rng);
        let gen_name = NameGen::location(&mut rng).generate_danari();
        let suffix = ["City", "Metropolis"].choose(&mut rng).unwrap();
        let name = match rng.random_range(0..2) {
            0 => format!("{} {}", gen_name, suffix),
            _ => format!("{} of {}", suffix, gen_name),
        };
        let mut site = Site {
            origin,
            name: Some(name.clone()),
            kind: Some(SiteKind::Myrmidon),
            ..Site::default()
        };

        let road_kind = plot::RoadKind {
            lights: plot::RoadLights::Default,
            material: plot::RoadMaterial::Dirt,
        };

        // place the initial plaza
        site.demarcate_obstacles(land);
        // The myrmidon_arena is 16 tiles in radius, so the plaza should be outside the
        // palace.
        const MYRMIDON_PLAZA_RADIUS: u32 = 3;
        const MYRMIDON_PLAZA_SEARCH_INNER: u32 = 18;
        const MYRMIDON_PLAZA_SEARCH_WIDTH: u32 = 12;
        generator_stats.add(site.name(), GenStatSiteKind::Myrmidon);
        generator_stats.attempt(site.name(), GenStatPlotKind::InitialPlaza);
        site.make_initial_plaza(
            land,
            index,
            &mut rng,
            MYRMIDON_PLAZA_RADIUS,
            MYRMIDON_PLAZA_SEARCH_INNER,
            MYRMIDON_PLAZA_SEARCH_WIDTH,
            generator_stats,
            &name,
            road_kind,
        );

        let size = 16.0 as i32;
        let aabr = Aabr {
            min: Vec2::broadcast(-size),
            max: Vec2::broadcast(size),
        };
        {
            let myrmidon_arena =
                plot::MyrmidonArena::generate(land, &mut reseed(&mut rng), &site, aabr);
            let myrmidon_arena_alt = myrmidon_arena.alt;
            let plot = site.create_plot(Plot {
                kind: PlotKind::MyrmidonArena(myrmidon_arena),
                root_tile: aabr.center(),
                tiles: aabr_tiles(aabr).collect(),
            });

            site.blit_aabr(aabr, Tile {
                kind: TileKind::Building,
                plot: Some(plot),
                hard_alt: Some(myrmidon_arena_alt),
            });
        }
        for _ in 0..30 {
            // MyrmidonHouse
            generator_stats.attempt(site.name(), GenStatPlotKind::House);
            let size = (9.0 + rng.random::<f32>().powf(5.0) * 1.5).round() as u32;
            if let Some((aabr, _, _, alt)) = attempt(32, || {
                site.find_roadside_aabr(&mut rng, 9..(size + 1).pow(2), Extent2::broadcast(size))
            }) {
                let myrmidon_house =
                    plot::MyrmidonHouse::generate(land, &mut reseed(&mut rng), &site, aabr, alt);
                let myrmidon_house_alt = myrmidon_house.alt;
                let plot = site.create_plot(Plot {
                    kind: PlotKind::MyrmidonHouse(myrmidon_house),
                    root_tile: aabr.center(),
                    tiles: aabr_tiles(aabr).collect(),
                });

                site.blit_aabr(aabr, Tile {
                    kind: TileKind::Building,
                    plot: Some(plot),
                    hard_alt: Some(myrmidon_house_alt),
                });

                generator_stats.success(site.name(), GenStatPlotKind::House);
            } else {
                site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
            }
        }

        site
    }

    pub fn generate_giant_tree(land: &Land, rng: &mut impl Rng, origin: Vec2<i32>) -> Self {
        let mut rng = reseed(rng);
        let mut site = Site {
            origin,
            kind: Some(SiteKind::GiantTree),
            ..Site::default()
        };
        site.demarcate_obstacles(land);
        let giant_tree = plot::GiantTree::generate(&site, Vec2::zero(), land, &mut rng);
        site.name = Some(giant_tree.name().to_string());
        let size = (giant_tree.radius() / TILE_SIZE as f32).ceil() as i32;
        let aabr = Aabr {
            min: Vec2::broadcast(-size),
            max: Vec2::broadcast(size) + 1,
        };
        let plot = site.create_plot(Plot {
            kind: PlotKind::GiantTree(giant_tree),
            root_tile: aabr.center(),
            tiles: aabr_tiles(aabr).collect(),
        });
        site.blit_aabr(aabr, Tile {
            kind: TileKind::Building,
            plot: Some(plot),
            hard_alt: None,
        });
        site
    }

    // Size is 0..1
    pub fn generate_city(
        land: &Land,
        index: IndexRef,
        rng: &mut impl Rng,
        origin: Vec2<i32>,
        size: f32,
        calendar: Option<&Calendar>,
        generator_stats: &mut SitesGenMeta,
    ) -> Self {
        let mut rng = reseed(rng);
        let name = NameGen::location(&mut rng).generate_town();
        let mut site = Site {
            origin,
            name: Some(name.clone()),
            kind: Some(SiteKind::Refactor),
            ..Site::default()
        };
        let road_kind = plot::RoadKind {
            lights: plot::RoadLights::Default,
            material: plot::RoadMaterial::Cobblestone,
        };

        // place the initial plaza
        site.demarcate_obstacles(land);
        generator_stats.add(site.name(), GenStatSiteKind::City);
        site.make_initial_plaza_default(land, index, &mut rng, generator_stats, &name, road_kind);

        let build_chance = Lottery::from(vec![
            (64.0, 1), // house
            (5.0, 2),  // guard tower
            (25.0, 3), // field
            //(32.0, 4), // castle
            (5.0, 5),  // workshop
            (15.0, 6), // airship dock
            (15.0, 7), // tavern
            (5.0, 8),  // barn
        ]);

        // These plots have minimums or limits.
        let mut workshops = 0;
        let mut castles = 0;
        let mut taverns = 0;
        let mut airship_docks = 0;

        for _ in 0..(size * 200.0) as i32 {
            match *build_chance.choose_seeded(rng.random()) {
                // Workshop
                n if (n == 5 && workshops < (size * 5.0) as i32) || workshops == 0 => {
                    generator_stats.attempt(site.name(), GenStatPlotKind::Workshop);
                    let size = (3.0 + rng.random::<f32>().powf(5.0) * 1.5).round() as u32;
                    if let Some((aabr, door_tile, door_dir, alt)) = attempt(32, || {
                        site.find_roadside_aabr(
                            &mut rng,
                            4..(size + 1).pow(2),
                            Extent2::broadcast(size),
                        )
                    }) {
                        let workshop = plot::Workshop::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            door_tile,
                            door_dir,
                            aabr,
                            alt,
                        );
                        let workshop_alt = workshop.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::Workshop(workshop),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(workshop_alt),
                        });
                        workshops += 1;
                        generator_stats.success(site.name(), GenStatPlotKind::Workshop);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },
                // House
                1 => {
                    let size = (1.5 + rng.random::<f32>().powf(5.0) * 1.0).round() as u32;
                    generator_stats.attempt(site.name(), GenStatPlotKind::House);
                    if let Some((aabr, door_tile, door_dir, alt)) = attempt(32, || {
                        site.find_roadside_aabr(
                            &mut rng,
                            4..(size + 1).pow(2),
                            Extent2::broadcast(size),
                        )
                    }) {
                        let house = plot::House::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            door_tile,
                            door_dir,
                            aabr,
                            calendar,
                            alt,
                        );
                        let house_alt = house.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::House(house),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(house_alt),
                        });
                        generator_stats.success(site.name(), GenStatPlotKind::House);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },
                // Guard tower
                2 => {
                    generator_stats.attempt(site.name(), GenStatPlotKind::GuardTower);
                    if let Some((_aabr, _, _door_dir, _)) = attempt(10, || {
                        site.find_roadside_aabr(&mut rng, 4..4, Extent2::new(2, 2))
                    }) {
                        // let plot = site.create_plot(Plot {
                        //     kind: PlotKind::Castle(plot::Castle::generate(
                        //         land,
                        //         &mut rng,
                        //         &site,
                        //         aabr,
                        //     )),
                        //     root_tile: aabr.center(),
                        //     tiles: aabr_tiles(aabr).collect(),
                        // });

                        // site.blit_aabr(aabr, Tile {
                        //     kind: TileKind::Castle,
                        //     plot: Some(plot),
                        //     hard_alt: None,
                        // });
                    }
                },
                // Field
                3 => {
                    Self::generate_farm(false, &mut rng, &mut site, land);
                },
                // Castle
                4 if size > 0.2 && castles < 1 => {
                    generator_stats.attempt(site.name(), GenStatPlotKind::Castle);
                    if let Some((aabr, _entrance_tile, _door_dir, _alt)) = attempt(64, || {
                        site.find_rural_aabr(&mut rng, 16 * 16..18 * 18, Extent2::new(16, 16))
                    }) {
                        let offset = rng.random_range(5..(aabr.size().w.min(aabr.size().h) - 4));
                        let gate_aabr = Aabr {
                            min: Vec2::new(aabr.min.x + offset - 1, aabr.min.y),
                            max: Vec2::new(aabr.min.x + offset + 2, aabr.min.y + 1),
                        };
                        let castle = plot::Castle::generate(land, &mut rng, &site, aabr, gate_aabr);
                        let castle_alt = castle.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::Castle(castle),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        let wall_north = Tile {
                            kind: TileKind::Wall(Dir2::Y),
                            plot: Some(plot),
                            hard_alt: Some(castle_alt),
                        };

                        let wall_east = Tile {
                            kind: TileKind::Wall(Dir2::X),
                            plot: Some(plot),
                            hard_alt: Some(castle_alt),
                        };
                        for x in 0..aabr.size().w {
                            site.tiles
                                .set(aabr.min + Vec2::new(x, 0), wall_east.clone());
                            site.tiles.set(
                                aabr.min + Vec2::new(x, aabr.size().h - 1),
                                wall_east.clone(),
                            );
                        }
                        for y in 0..aabr.size().h {
                            site.tiles
                                .set(aabr.min + Vec2::new(0, y), wall_north.clone());
                            site.tiles.set(
                                aabr.min + Vec2::new(aabr.size().w - 1, y),
                                wall_north.clone(),
                            );
                        }

                        let gate = Tile {
                            kind: TileKind::Gate,
                            plot: Some(plot),
                            hard_alt: Some(castle_alt),
                        };
                        let tower_parapet = Tile {
                            kind: TileKind::Tower(RoofKind::Parapet),
                            plot: Some(plot),
                            hard_alt: Some(castle_alt),
                        };
                        let tower_pyramid = Tile {
                            kind: TileKind::Tower(RoofKind::Pyramid),
                            plot: Some(plot),
                            hard_alt: Some(castle_alt),
                        };

                        site.tiles.set(
                            Vec2::new(aabr.min.x + offset - 2, aabr.min.y),
                            tower_parapet.clone(),
                        );
                        site.tiles
                            .set(Vec2::new(aabr.min.x + offset - 1, aabr.min.y), gate.clone());
                        site.tiles
                            .set(Vec2::new(aabr.min.x + offset, aabr.min.y), gate.clone());
                        site.tiles
                            .set(Vec2::new(aabr.min.x + offset + 1, aabr.min.y), gate.clone());
                        site.tiles.set(
                            Vec2::new(aabr.min.x + offset + 2, aabr.min.y),
                            tower_parapet.clone(),
                        );

                        site.tiles
                            .set(Vec2::new(aabr.min.x, aabr.min.y), tower_parapet.clone());
                        site.tiles
                            .set(Vec2::new(aabr.max.x - 1, aabr.min.y), tower_parapet.clone());
                        site.tiles
                            .set(Vec2::new(aabr.min.x, aabr.max.y - 1), tower_parapet.clone());
                        site.tiles.set(
                            Vec2::new(aabr.max.x - 1, aabr.max.y - 1),
                            tower_parapet.clone(),
                        );

                        // Courtyard
                        site.blit_aabr(
                            Aabr {
                                min: aabr.min + 1,
                                max: aabr.max - 1,
                            },
                            Tile {
                                kind: TileKind::Road {
                                    a: 0,
                                    b: 0,
                                    w: 0,
                                    alt: castle_alt as f32,
                                },
                                plot: Some(plot),
                                hard_alt: Some(castle_alt),
                            },
                        );

                        // Keep
                        site.blit_aabr(
                            Aabr {
                                min: aabr.center() - 3,
                                max: aabr.center() + 3,
                            },
                            Tile {
                                kind: TileKind::Wall(Dir2::Y),
                                plot: Some(plot),
                                hard_alt: Some(castle_alt),
                            },
                        );
                        site.tiles.set(
                            Vec2::new(aabr.center().x + 2, aabr.center().y + 2),
                            tower_pyramid.clone(),
                        );
                        site.tiles.set(
                            Vec2::new(aabr.center().x + 2, aabr.center().y - 3),
                            tower_pyramid.clone(),
                        );
                        site.tiles.set(
                            Vec2::new(aabr.center().x - 3, aabr.center().y + 2),
                            tower_pyramid.clone(),
                        );
                        site.tiles.set(
                            Vec2::new(aabr.center().x - 3, aabr.center().y - 3),
                            tower_pyramid.clone(),
                        );

                        site.blit_aabr(
                            Aabr {
                                min: aabr.center() - 2,
                                max: aabr.center() + 2,
                            },
                            Tile {
                                kind: TileKind::Keep(KeepKind::Middle),
                                plot: Some(plot),
                                hard_alt: Some(castle_alt),
                            },
                        );

                        castles += 1;
                        generator_stats.success(site.name(), GenStatPlotKind::Castle);
                    }
                },
                //airship dock
                6 if (size > 0.125 && airship_docks == 0) => {
                    generator_stats.attempt(site.name(), GenStatPlotKind::AirshipDock);
                    // The airship dock rendered size is 10x10.
                    // The rendering code allows for a margin of 1 tile
                    // and we want to keep the dock footprint as small as possible.
                    // The area range for the aabr is fixed at size squared (81) since the
                    // dock structure is square.
                    let size = 9u32;
                    if let Some((aabr, door_tile, door_dir, _)) = attempt(32, || {
                        site.find_roadside_aabr(&mut rng, 81..82, Extent2::broadcast(size))
                    }) {
                        let airship_dock = plot::AirshipDock::generate(
                            land,
                            index,
                            &mut reseed(&mut rng),
                            &site,
                            door_tile,
                            door_dir,
                            aabr,
                        );
                        let airship_dock_alt = airship_dock.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::AirshipDock(airship_dock),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(airship_dock_alt),
                        });
                        airship_docks += 1;
                        generator_stats.success(site.name(), GenStatPlotKind::AirshipDock);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },
                7 if (size > 0.125 && taverns < 2) => {
                    generator_stats.attempt(site.name(), GenStatPlotKind::Tavern);
                    let size = (4.5 + rng.random::<f32>().powf(5.0) * 2.0).round() as u32;
                    if let Some((aabr, door_tile, door_dir, alt)) = attempt(32, || {
                        site.find_roadside_aabr(
                            &mut rng,
                            8..(size + 1).pow(2),
                            Extent2::broadcast(size),
                        )
                    }) {
                        let tavern = plot::Tavern::generate(
                            land,
                            index,
                            &mut reseed(&mut rng),
                            &site,
                            door_tile,
                            Dir2::from_vec2(door_dir),
                            aabr,
                            alt,
                        );
                        let tavern_alt = tavern.door_wpos.z;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::Tavern(tavern),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(tavern_alt),
                        });

                        taverns += 1;
                        generator_stats.success(site.name(), GenStatPlotKind::Tavern);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },
                8 => {
                    Self::generate_barn(false, &mut rng, &mut site, land, index);
                },
                _ => {},
            }
        }

        site
    }

    pub fn generate_glider_course(
        land: &Land,
        _index: IndexRef,
        rng: &mut impl Rng,
        origin: Vec2<i32>,
    ) -> Self {
        let mut rng = reseed(rng);
        let mut site = Site {
            origin,
            kind: Some(SiteKind::GliderCourse),
            ..Site::default()
        };

        // TODO use the nearest peak name. Unfortunately this requires `Civs` but we
        // only have access to `WorldSim`
        site.name = Some(NameGen::location(&mut rng).generate_town() + " Glider Course");

        // Pick the starting downhill direction based on the average drop over
        // two chunks in the four cardinal directions
        let origin_alt = land.get_alt_approx(origin);
        let alt_drops: Vec<f32> = CARDINALS
            .iter()
            .map(|c| {
                origin_alt
                    - 0.5
                        * (land.get_alt_approx(origin + *c * TerrainChunkSize::RECT_SIZE.x as i32)
                            + land.get_alt_approx(
                                origin + 2 * *c * TerrainChunkSize::RECT_SIZE.x as i32,
                            ))
            })
            .collect();
        let mut cardinal = 0;
        let mut max_drop = 0.0;
        for (i, drop) in alt_drops.iter().enumerate() {
            if *drop > max_drop {
                max_drop = *drop;
                cardinal = i;
            }
        }
        let dir = match cardinal {
            0 => Dir2::X,
            1 => Dir2::Y,
            2 => Dir2::NegX,
            3 => Dir2::NegY,
            _ => Dir2::X,
        };
        let size = 2.0;

        let mut valid_course = true;
        let mut positions = Vec::new();

        // Platform
        let pos = origin;
        let tile_pos: Vec2<i32> = Vec2::zero();
        positions.push((pos, tile_pos));

        // This defines the distance between rings
        // An offset of 5 results in courses that are about 1 minute long
        // An offset of 6+ results in not all plots being in range of the site
        const CHUNK_OFFSET: usize = 5;
        // WARNING: This assumes x and y lengths of a chunk are the same!!!
        let offset = CHUNK_OFFSET as i32 * TerrainChunkSize::RECT_SIZE.x as i32;
        // Always convert to tiles then back to wpos to remove any integer division
        // artifacts
        let tile_offset = offset / TILE_SIZE as i32;
        let pos_offset = tile_offset * TILE_SIZE as i32;

        // Loop 1 is always straight forward from the launch platform
        let pos = origin + pos_offset * dir.to_vec2();
        let tile_pos = tile_offset * dir.to_vec2();
        positions.push((pos, tile_pos));

        // Loops 2-9 follow the downhill path of terrain chunks
        // In the future it may be desirable to follow ridges and the like but that
        // would be a future MR
        let mut last_pos = pos;
        let mut last_tile_pos = tile_pos;
        for j in 1..(CHUNK_OFFSET * 9 + 1) {
            let c_downhill = land.get_chunk_wpos(last_pos).and_then(|c| c.downhill);
            if let Some(downhill) = c_downhill {
                let downhill_chunk =
                    downhill.map2(TerrainChunkSize::RECT_SIZE, |e, sz: u32| e / (sz as i32));
                let downhill_chunk_pos = TerrainChunkSize::center_wpos(downhill_chunk);
                let downhill_vec = downhill_chunk_pos - last_pos;
                // Convert to tiles first, then back to wpos to ensure coordinates align, as
                // chunks are not tile aligned
                let tile_offset = downhill_vec / (TILE_SIZE as i32);
                let pos_offset = tile_offset * TILE_SIZE as i32;
                let pos = last_pos + pos_offset;
                let tile_pos = last_tile_pos + tile_offset;
                last_pos = pos;
                last_tile_pos = tile_pos;
                // Only want to save positions with large enough chunk offsets, not every chunk
                // position
                if j % CHUNK_OFFSET == 0 {
                    positions.push((pos, tile_pos));
                }
            } else {
                valid_course = false;
            }
        }
        // Currently there is no check to ensure the delta z between rings is
        // sufficient to successfully fly through the course. This should cause
        // no glider course site to be created. Right now it just doesn't spawn
        // one in the world (similar to towns when placed near/on bodies of water).
        // In the future maybe the generate functions should return an `Option`
        // instead of a `Site`
        if valid_course && positions.len() > 1 {
            for (i, window) in positions.windows(2).enumerate() {
                if !window.is_empty() {
                    let [(pos, tile_pos), (next_pos, next_tile_pos)] = window else {
                        panic!(
                            "previous condition required positions Vec to have at least two \
                             elements"
                        );
                    };
                    if i == 0 {
                        // Launch platform
                        let aabr = Aabr {
                            min: Vec2::broadcast(-size as i32),
                            max: Vec2::broadcast(size as i32),
                        };
                        let glider_platform = plot::GliderPlatform::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            *pos,
                            dir,
                        );
                        let alt = glider_platform.alt - 5;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::GliderPlatform(glider_platform),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });
                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(alt),
                        });
                    } else if i < 9 {
                        // Point each ring after 1 towards the next ring
                        // This provides a subtle guide through the course
                        let dir = if i > 1 {
                            Dir2::from_vec2(next_pos - pos)
                        } else {
                            dir
                        };
                        let aabr = Aabr {
                            min: Vec2::broadcast(-size as i32) + tile_pos,
                            max: Vec2::broadcast(size as i32) + tile_pos,
                        };
                        let glider_ring = plot::GliderRing::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            pos,
                            i,
                            dir,
                        );
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::GliderRing(glider_ring),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });
                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: None,
                        });
                    } else if i == 9 {
                        // last ring (ring 9) and finish platform
                        // Separate condition due to window iterator to ensure
                        // the finish platform is generated
                        let dir = Dir2::from_vec2(next_pos - pos);
                        let aabr = Aabr {
                            min: Vec2::broadcast(-size as i32) + tile_pos,
                            max: Vec2::broadcast(size as i32) + tile_pos,
                        };
                        let glider_ring = plot::GliderRing::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            pos,
                            i,
                            dir,
                        );
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::GliderRing(glider_ring),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });
                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: None,
                        });
                        // Finish
                        let size = 10.0;
                        let aabr = Aabr {
                            min: Vec2::broadcast(-size as i32) + next_tile_pos,
                            max: Vec2::broadcast(size as i32) + next_tile_pos,
                        };
                        let glider_finish = plot::GliderFinish::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            *next_pos,
                        );
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::GliderFinish(glider_finish),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });
                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: None,
                        });
                    }
                }
            }
        }

        site
    }

    pub fn generate_cliff_town(
        land: &Land,
        index: IndexRef,
        rng: &mut impl Rng,
        origin: Vec2<i32>,
        generator_stats: &mut SitesGenMeta,
    ) -> Self {
        let mut rng = reseed(rng);
        let name = NameGen::location(&mut rng).generate_arabic();
        let mut site = Site {
            origin,
            name: Some(name.clone()),
            kind: Some(SiteKind::CliffTown),
            ..Site::default()
        };
        let mut campfires = 0;
        let road_kind = plot::RoadKind {
            lights: plot::RoadLights::Default,
            material: plot::RoadMaterial::Sandstone,
        };

        // place the initial plaza
        generator_stats.add(site.name(), GenStatSiteKind::CliffTown);
        site.make_initial_plaza_default(land, index, &mut rng, generator_stats, &name, road_kind);

        let build_chance = Lottery::from(vec![(30.0, 1), (50.0, 2)]);
        let mut airship_docks = 0;
        for _ in 0..80 {
            match *build_chance.choose_seeded(rng.random()) {
                1 => {
                    // CliffTower
                    let size = (9.0 + rng.random::<f32>().powf(5.0) * 1.0).round() as u32;
                    generator_stats.attempt(site.name(), GenStatPlotKind::House);
                    let campfire = campfires < 4;
                    if let Some((aabr, door_tile, door_dir, alt)) = attempt(32, || {
                        site.find_roadside_aabr(
                            &mut rng,
                            8..(size + 1).pow(2),
                            Extent2::broadcast(size),
                        )
                    }) {
                        let cliff_tower = plot::CliffTower::generate(
                            land,
                            index,
                            &mut reseed(&mut rng),
                            &site,
                            door_tile,
                            door_dir,
                            aabr,
                            campfire,
                            alt,
                        );
                        let cliff_tower_alt = cliff_tower.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::CliffTower(cliff_tower),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });
                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(cliff_tower_alt),
                        });
                        campfires += 1;
                        generator_stats.success(site.name(), GenStatPlotKind::House);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },
                2 if airship_docks < 1 => {
                    // CliffTownAirshipDock
                    let size = 25u32;
                    generator_stats.attempt(site.name(), GenStatPlotKind::AirshipDock);
                    if let Some((aabr, door_tile, door_dir, _)) = attempt(32, || {
                        site.find_roadside_aabr(&mut rng, 625..626, Extent2::broadcast(size))
                    }) {
                        let cliff_town_airship_dock = plot::CliffTownAirshipDock::generate(
                            land,
                            index,
                            &mut reseed(&mut rng),
                            &site,
                            door_tile,
                            door_dir,
                            aabr,
                        );
                        let cliff_town_airship_dock_alt = cliff_town_airship_dock.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::CliffTownAirshipDock(cliff_town_airship_dock),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(cliff_town_airship_dock_alt),
                        });
                        airship_docks += 1;
                        generator_stats.success(site.name(), GenStatPlotKind::AirshipDock);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },
                _ => {},
            }
        }

        site.demarcate_obstacles(land);
        site
    }

    pub fn generate_savannah_town(
        land: &Land,
        index: IndexRef,
        rng: &mut impl Rng,
        origin: Vec2<i32>,
        generator_stats: &mut SitesGenMeta,
    ) -> Self {
        let mut rng = reseed(rng);
        let name = NameGen::location(&mut rng).generate_savannah_custom();
        let mut site = Site {
            origin,
            name: Some(name.clone()),
            kind: Some(SiteKind::SavannahTown),
            ..Site::default()
        };
        let road_kind = plot::RoadKind {
            lights: plot::RoadLights::Default,
            material: plot::RoadMaterial::Dirt,
        };

        // place the initial plaza
        site.demarcate_obstacles(land);
        generator_stats.add(site.name(), GenStatSiteKind::SavannahTown);
        site.make_initial_plaza_default(land, index, &mut rng, generator_stats, &name, road_kind);

        let mut workshops = 0;
        let mut airship_dock = 0;
        let build_chance = Lottery::from(vec![
            (25.0, 1),
            (5.0, 2),
            (5.0, 3),
            (15.0, 4),
            (5.0, 5),
            (5.0, 6),
        ]);

        for _ in 0..50 {
            match *build_chance.choose_seeded(rng.random()) {
                n if (n == 2 && workshops < 3) || workshops == 0 => {
                    // SavannahWorkshop
                    let size = (4.0 + rng.random::<f32>().powf(5.0) * 1.5).round() as u32;
                    generator_stats.attempt(site.name(), GenStatPlotKind::Workshop);
                    if let Some((aabr, door_tile, door_dir, alt)) = attempt(32, || {
                        site.find_roadside_aabr(
                            &mut rng,
                            4..(size + 1).pow(2),
                            Extent2::broadcast(size),
                        )
                    }) {
                        let savannah_workshop = plot::SavannahWorkshop::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            door_tile,
                            door_dir,
                            aabr,
                            alt,
                        );
                        let savannah_workshop_alt = savannah_workshop.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::SavannahWorkshop(savannah_workshop),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(savannah_workshop_alt),
                        });
                        workshops += 1;
                        generator_stats.success(site.name(), GenStatPlotKind::Workshop);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },
                1 => {
                    // SavannahHut
                    let size = (4.0 + rng.random::<f32>().powf(5.0) * 1.5).round() as u32;
                    generator_stats.attempt(site.name(), GenStatPlotKind::House);
                    if let Some((aabr, door_tile, door_dir, alt)) = attempt(32, || {
                        site.find_roadside_aabr(
                            &mut rng,
                            4..(size + 1).pow(2),
                            Extent2::broadcast(size),
                        )
                    }) {
                        let savannah_hut = plot::SavannahHut::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            door_tile,
                            door_dir,
                            aabr,
                            alt,
                        );
                        let savannah_hut_alt = savannah_hut.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::SavannahHut(savannah_hut),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(savannah_hut_alt),
                        });
                        generator_stats.success(site.name(), GenStatPlotKind::House);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },
                3 if airship_dock < 1 => {
                    // SavannahAirshipDock
                    let size = 9u32;
                    generator_stats.attempt(site.name(), GenStatPlotKind::AirshipDock);
                    if let Some((aabr, door_tile, _, _)) = attempt(48, || {
                        site.find_roadside_aabr(&mut rng, 81..82, Extent2::broadcast(size))
                    }) {
                        let savannah_airship_dock = plot::SavannahAirshipDock::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            door_tile,
                            aabr,
                        );
                        let savannah_airship_dock_alt = savannah_airship_dock.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::SavannahAirshipDock(savannah_airship_dock),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(savannah_airship_dock_alt),
                        });
                        airship_dock += 1;
                        generator_stats.success(site.name(), GenStatPlotKind::AirshipDock);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },
                // Field
                4 => {
                    Self::generate_farm(false, &mut rng, &mut site, land);
                },
                5 => {
                    Self::generate_barn(false, &mut rng, &mut site, land, index);
                },
                6 => {
                    // SavannahGuardHut
                    let size = (4.0 + rng.random::<f32>().powf(5.0) * 1.5).round() as u32;
                    generator_stats.attempt(site.name(), GenStatPlotKind::House);
                    if let Some((aabr, door_tile, door_dir, alt)) = attempt(32, || {
                        site.find_roadside_aabr(
                            &mut rng,
                            2..(size + 1).pow(2),
                            Extent2::broadcast(size),
                        )
                    }) {
                        let savannah_guard_hut = plot::SavannahGuardHut::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            door_tile,
                            door_dir,
                            aabr,
                            alt,
                        );
                        let savannah_guard_hut_alt = savannah_guard_hut.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::SavannahGuardHut(savannah_guard_hut),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(savannah_guard_hut_alt),
                        });
                        generator_stats.success(site.name(), GenStatPlotKind::House);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },
                _ => {},
            }
        }
        site
    }

    pub fn generate_coastal_town(
        land: &Land,
        index: IndexRef,
        rng: &mut impl Rng,
        origin: Vec2<i32>,
        generator_stats: &mut SitesGenMeta,
    ) -> Self {
        let mut rng = reseed(rng);
        let name = NameGen::location(&mut rng).generate_danari();
        let mut site = Site {
            origin,
            name: Some(name.clone()),
            kind: Some(SiteKind::CoastalTown),
            ..Site::default()
        };
        let road_kind = plot::RoadKind {
            lights: plot::RoadLights::Default,
            material: plot::RoadMaterial::Marble,
        };

        // place the initial plaza
        site.demarcate_obstacles(land);
        generator_stats.add(site.name(), GenStatSiteKind::CoastalTown);
        site.make_initial_plaza_default(land, index, &mut rng, generator_stats, &name, road_kind);

        let mut workshops = 0;
        let build_chance = Lottery::from(vec![(38.0, 1), (5.0, 2), (15.0, 3), (15.0, 4), (5.0, 5)]);
        let mut airship_docks = 0;
        for _ in 0..55 {
            match *build_chance.choose_seeded(rng.random()) {
                n if (n == 2 && workshops < 3) || workshops == 0 => {
                    // CoastalWorkshop
                    let size = (7.0 + rng.random::<f32>().powf(5.0) * 1.5).round() as u32;
                    generator_stats.attempt(site.name(), GenStatPlotKind::Workshop);
                    if let Some((aabr, door_tile, door_dir, alt)) = attempt(32, || {
                        site.find_roadside_aabr(
                            &mut rng,
                            7..(size + 1).pow(2),
                            Extent2::broadcast(size),
                        )
                    }) {
                        let coastal_workshop = plot::CoastalWorkshop::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            door_tile,
                            door_dir,
                            aabr,
                            alt,
                        );
                        let coastal_workshop_alt = coastal_workshop.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::CoastalWorkshop(coastal_workshop),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(coastal_workshop_alt),
                        });
                        workshops += 1;
                        generator_stats.success(site.name(), GenStatPlotKind::Workshop);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },
                1 => {
                    // CoastalHouse
                    let size = (7.0 + rng.random::<f32>().powf(5.0) * 1.5).round() as u32;
                    generator_stats.attempt(site.name(), GenStatPlotKind::House);
                    if let Some((aabr, door_tile, door_dir, alt)) = attempt(32, || {
                        site.find_roadside_aabr(
                            &mut rng,
                            7..(size + 1).pow(2),
                            Extent2::broadcast(size),
                        )
                    }) {
                        let coastal_house = plot::CoastalHouse::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            door_tile,
                            door_dir,
                            aabr,
                            alt,
                        );
                        let coastal_house_alt = coastal_house.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::CoastalHouse(coastal_house),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(coastal_house_alt),
                        });

                        generator_stats.success(site.name(), GenStatPlotKind::House);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },
                3 if airship_docks < 1 => {
                    // CoastalAirshipDock
                    // The airship dock rendered size is 10x10.
                    // The rendering code allows for a margin of 1 tile
                    // and we want to keep the dock footprint as small as possible.
                    // The area range for the aabr is fixed at size squared (81) since the
                    // dock structure is square.
                    let size = 9u32;
                    generator_stats.attempt(site.name(), GenStatPlotKind::AirshipDock);
                    if let Some((aabr, door_tile, _, _)) = attempt(32, || {
                        site.find_roadside_aabr(&mut rng, 81..82, Extent2::broadcast(size))
                    }) {
                        let coastal_airship_dock = plot::CoastalAirshipDock::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            door_tile,
                            aabr,
                        );
                        let coastal_airship_dock_alt = coastal_airship_dock.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::CoastalAirshipDock(coastal_airship_dock),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(coastal_airship_dock_alt),
                        });
                        airship_docks += 1;
                        generator_stats.success(site.name(), GenStatPlotKind::AirshipDock);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },
                // Field
                4 => {
                    Self::generate_farm(false, &mut rng, &mut site, land);
                },
                5 => {
                    Self::generate_barn(false, &mut rng, &mut site, land, index);
                },
                _ => {},
            }
        }
        site
    }

    pub fn generate_desert_city(
        land: &Land,
        index: IndexRef,
        rng: &mut impl Rng,
        origin: Vec2<i32>,
        generator_stats: &mut SitesGenMeta,
    ) -> Self {
        let mut rng = reseed(rng);

        let name = NameGen::location(&mut rng).generate_arabic();
        let mut site = Site {
            origin,
            name: Some(name.clone()),
            kind: Some(SiteKind::DesertCity),
            ..Site::default()
        };
        let road_kind = plot::RoadKind {
            lights: plot::RoadLights::Default,
            material: plot::RoadMaterial::Sandstone,
        };

        // place the initial plaza
        site.demarcate_obstacles(land);
        // The desert_city_arena is 17 tiles in radius, so the plaza should be outside
        // the palace.
        const DESERT_CITY_PLAZA_RADIUS: u32 = 3;
        const DESERT_CITY_PLAZA_SEARCH_INNER: u32 = 19;
        const DESERT_CITY_PLAZA_SEARCH_WIDTH: u32 = 12;
        generator_stats.add(site.name(), GenStatSiteKind::DesertCity);
        site.make_initial_plaza(
            land,
            index,
            &mut rng,
            DESERT_CITY_PLAZA_RADIUS,
            DESERT_CITY_PLAZA_SEARCH_INNER,
            DESERT_CITY_PLAZA_SEARCH_WIDTH,
            generator_stats,
            &name,
            road_kind,
        );

        let size = 17.0 as i32;
        let aabr = Aabr {
            min: Vec2::broadcast(-size),
            max: Vec2::broadcast(size),
        };

        let desert_city_arena =
            plot::DesertCityArena::generate(land, &mut reseed(&mut rng), &site, aabr);

        let desert_city_arena_alt = desert_city_arena.alt;
        let plot = site.create_plot(Plot {
            kind: PlotKind::DesertCityArena(desert_city_arena),
            root_tile: aabr.center(),
            tiles: aabr_tiles(aabr).collect(),
        });

        site.blit_aabr(aabr, Tile {
            kind: TileKind::Building,
            plot: Some(plot),
            hard_alt: Some(desert_city_arena_alt),
        });

        let build_chance =
            Lottery::from(vec![(20.0, 1), (10.0, 2), (15.0, 3), (10.0, 4), (0.0, 5)]);

        let mut temples = 0;
        let mut airship_docks = 0;
        let mut campfires = 0;

        for _ in 0..35 {
            match *build_chance.choose_seeded(rng.random()) {
                // DesertCityMultiplot
                1 => {
                    let size = (9.0 + rng.random::<f32>().powf(5.0) * 1.5).round() as u32;
                    generator_stats.attempt(site.name(), GenStatPlotKind::MultiPlot);
                    let campfire = campfires < 4;
                    if let Some((aabr, door_tile, door_dir, alt)) = attempt(32, || {
                        site.find_roadside_aabr(
                            &mut rng,
                            8..(size + 1).pow(2),
                            Extent2::broadcast(size),
                        )
                    }) {
                        let desert_city_multi_plot = plot::DesertCityMultiPlot::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            door_tile,
                            door_dir,
                            aabr,
                            campfire,
                            alt,
                        );
                        let desert_city_multi_plot_alt = desert_city_multi_plot.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::DesertCityMultiPlot(desert_city_multi_plot),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(desert_city_multi_plot_alt),
                        });
                        campfires += 1;
                        generator_stats.success(site.name(), GenStatPlotKind::MultiPlot);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },
                // DesertCityTemple
                2 if temples < 1 => {
                    let size = (9.0 + rng.random::<f32>().powf(5.0) * 1.5).round() as u32;
                    generator_stats.attempt(site.name(), GenStatPlotKind::Temple);
                    if let Some((aabr, door_tile, door_dir, alt)) = attempt(32, || {
                        site.find_roadside_aabr(
                            &mut rng,
                            8..(size + 1).pow(2),
                            Extent2::broadcast(size),
                        )
                    }) {
                        let desert_city_temple = plot::DesertCityTemple::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            door_tile,
                            door_dir,
                            aabr,
                            alt,
                        );
                        let desert_city_temple_alt = desert_city_temple.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::DesertCityTemple(desert_city_temple),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(desert_city_temple_alt),
                        });
                        temples += 1;
                        generator_stats.success(site.name(), GenStatPlotKind::Temple);
                    }
                },
                3 if airship_docks < 1 => {
                    // DesertCityAirshipDock
                    // The airship dock rendered size is 10x10.
                    // The rendering code allows for a margin of 1 tile
                    // and we want to keep the dock footprint as small as possible.
                    // The area range for the aabr is fixed at size squared (81) since the
                    // dock structure is square.
                    let size = 9u32;
                    generator_stats.attempt(site.name(), GenStatPlotKind::AirshipDock);
                    if let Some((aabr, door_tile, door_dir, alt)) = attempt(100, || {
                        site.find_roadside_aabr(&mut rng, 81..82, Extent2::broadcast(size))
                    }) {
                        let desert_city_airship_dock = plot::DesertCityAirshipDock::generate(
                            land,
                            &mut reseed(&mut rng),
                            &site,
                            door_tile,
                            door_dir,
                            aabr,
                            alt,
                        );
                        let desert_city_airship_dock_alt = desert_city_airship_dock.alt;
                        let plot = site.create_plot(Plot {
                            kind: PlotKind::DesertCityAirshipDock(desert_city_airship_dock),
                            root_tile: aabr.center(),
                            tiles: aabr_tiles(aabr).collect(),
                        });

                        site.blit_aabr(aabr, Tile {
                            kind: TileKind::Building,
                            plot: Some(plot),
                            hard_alt: Some(desert_city_airship_dock_alt),
                        });
                        airship_docks += 1;
                        generator_stats.success(site.name(), GenStatPlotKind::AirshipDock);
                    } else {
                        site.make_plaza(land, index, &mut rng, generator_stats, &name, road_kind);
                    }
                },
                // cactus farm
                4 => {
                    Self::generate_farm(true, &mut rng, &mut site, land);
                },
                // desert barn - disabled for now (0.0 spawn chance)
                // need desert-variant sprite
                5 => {
                    Self::generate_barn(true, &mut rng, &mut site, land, index);
                },
                _ => {},
            }
        }
        site
    }

    pub fn generate_farm(
        is_desert: bool,
        mut rng: &mut impl Rng,
        site: &mut Site,
        land: &Land,
    ) -> bool {
        let size = (3.0 + rng.random::<f32>().powf(5.0) * 6.0).round() as u32;
        if let Some((aabr, door_tile, door_dir, _alt)) = attempt(32, || {
            site.find_rural_aabr(&mut rng, 6..(size + 1).pow(2), Extent2::broadcast(size))
        }) {
            let field = plot::FarmField::generate(
                land,
                &mut reseed(&mut rng),
                site,
                door_tile,
                door_dir,
                aabr,
                is_desert,
            );

            let field_alt = field.alt;
            let plot = site.create_plot(Plot {
                kind: PlotKind::FarmField(field),
                root_tile: aabr.center(),
                tiles: aabr_tiles(aabr).collect(),
            });

            site.blit_aabr(aabr, Tile {
                kind: TileKind::Field,
                plot: Some(plot),
                hard_alt: Some(field_alt),
            });
            true
        } else {
            false
        }
    }

    pub fn generate_barn(
        is_desert: bool,
        mut rng: &mut impl Rng,
        site: &mut Site,
        land: &Land,
        index: IndexRef,
    ) -> bool {
        let size = (7.0 + rng.random::<f32>().powf(5.0) * 1.5).round() as u32;
        if let Some((aabr, door_tile, door_dir, _alt)) = attempt(32, || {
            site.find_rural_aabr(&mut rng, 7..(size + 1).pow(2), Extent2::broadcast(size))
        }) {
            let bounds = Aabr {
                min: site.tile_wpos(aabr.min),
                max: site.tile_wpos(aabr.max),
            };

            // Barns don't place very well in hilly sections of land due to
            // their size. They need a relatively flat area.
            let gradient_avg = get_gradient_average(bounds, land);

            if gradient_avg > 0.5 {
                false
            } else {
                let barn = plot::Barn::generate(
                    land,
                    index,
                    &mut reseed(&mut rng),
                    site,
                    door_tile,
                    door_dir,
                    aabr,
                    is_desert,
                );
                let barn_alt = barn.alt;
                let plot = site.create_plot(Plot {
                    kind: PlotKind::Barn(barn),
                    root_tile: aabr.center(),
                    tiles: aabr_tiles(aabr).collect(),
                });

                site.blit_aabr(aabr, Tile {
                    kind: TileKind::Building,
                    plot: Some(plot),
                    hard_alt: Some(barn_alt),
                });

                true
            }
        } else {
            false
        }
    }

    pub fn generate_haniwa(land: &Land, rng: &mut impl Rng, origin: Vec2<i32>) -> Self {
        let mut rng = reseed(rng);
        let mut site = Site {
            origin,
            name: Some(format!(
                "{} {}",
                NameGen::location(&mut rng).generate_haniwa(),
                [
                    "Catacombs",
                    "Crypt",
                    "Tomb",
                    "Gravemound",
                    "Tunnels",
                    "Vault",
                    "Chambers",
                    "Halls",
                    "Tumulus",
                    "Barrow",
                ]
                .choose(&mut rng)
                .unwrap()
            )),
            kind: Some(SiteKind::Haniwa),
            ..Site::default()
        };
        let size = 24.0 as i32;
        let aabr = Aabr {
            min: Vec2::broadcast(-size),
            max: Vec2::broadcast(size),
        };
        {
            let haniwa = plot::Haniwa::generate(land, &mut reseed(&mut rng), &site, aabr);
            let haniwa_alt = haniwa.alt;
            let plot = site.create_plot(Plot {
                kind: PlotKind::Haniwa(haniwa),
                root_tile: aabr.center(),
                tiles: aabr_tiles(aabr).collect(),
            });

            site.blit_aabr(aabr, Tile {
                kind: TileKind::Building,
                plot: Some(plot),
                hard_alt: Some(haniwa_alt),
            });
        }
        site
    }

    pub fn generate_chapel_site(land: &Land, rng: &mut impl Rng, origin: Vec2<i32>) -> Self {
        let mut rng = reseed(rng);
        let mut site = Site {
            origin,
            name: Some(NameGen::location(&mut rng).generate_danari()),
            kind: Some(SiteKind::ChapelSite),
            ..Site::default()
        };

        // SeaChapel
        let size = 10.0 as i32;
        let aabr = Aabr {
            min: Vec2::broadcast(-size),
            max: Vec2::broadcast(size),
        };
        {
            let sea_chapel = plot::SeaChapel::generate(land, &mut reseed(&mut rng), &site, aabr);
            let sea_chapel_alt = sea_chapel.alt;
            let plot = site.create_plot(Plot {
                kind: PlotKind::SeaChapel(sea_chapel),
                root_tile: aabr.center(),
                tiles: aabr_tiles(aabr).collect(),
            });

            site.blit_aabr(aabr, Tile {
                kind: TileKind::Building,
                plot: Some(plot),
                hard_alt: Some(sea_chapel_alt),
            });
        }
        site
    }

    pub fn generate_pirate_hideout(land: &Land, rng: &mut impl Rng, origin: Vec2<i32>) -> Self {
        let mut rng = reseed(rng);
        let mut site = Site {
            origin,
            name: None,
            kind: Some(SiteKind::PirateHideout),
            ..Site::default()
        };

        let size = 8.0 as i32;
        let aabr = Aabr {
            min: Vec2::broadcast(-size),
            max: Vec2::broadcast(size),
        };
        {
            let pirate_hideout =
                plot::PirateHideout::generate(land, &mut reseed(&mut rng), &site, aabr);
            let pirate_hideout_alt = pirate_hideout.alt;
            let plot = site.create_plot(Plot {
                kind: PlotKind::PirateHideout(pirate_hideout),
                root_tile: aabr.center(),
                tiles: aabr_tiles(aabr).collect(),
            });

            site.blit_aabr(aabr, Tile {
                kind: TileKind::Building,
                plot: Some(plot),
                hard_alt: Some(pirate_hideout_alt),
            });
        }
        site
    }

    pub fn generate_jungle_ruin(land: &Land, rng: &mut impl Rng, origin: Vec2<i32>) -> Self {
        let mut rng = reseed(rng);
        let mut site = Site {
            origin,
            name: None,
            kind: Some(SiteKind::JungleRuin),
            ..Site::default()
        };
        let size = 8.0 as i32;
        let aabr = Aabr {
            min: Vec2::broadcast(-size),
            max: Vec2::broadcast(size),
        };
        {
            let jungle_ruin = plot::JungleRuin::generate(land, &mut reseed(&mut rng), &site, aabr);
            let jungle_ruin_alt = jungle_ruin.alt;
            let plot = site.create_plot(Plot {
                kind: PlotKind::JungleRuin(jungle_ruin),
                root_tile: aabr.center(),
                tiles: aabr_tiles(aabr).collect(),
            });

            site.blit_aabr(aabr, Tile {
                kind: TileKind::Building,
                plot: Some(plot),
                hard_alt: Some(jungle_ruin_alt),
            });
        }
        site
    }

    pub fn generate_rock_circle(land: &Land, rng: &mut impl Rng, origin: Vec2<i32>) -> Self {
        let mut rng = reseed(rng);
        let mut site = Site {
            origin,
            kind: Some(SiteKind::RockCircle),
            ..Site::default()
        };
        let size = 8.0 as i32;
        let aabr = Aabr {
            min: Vec2::broadcast(-size),
            max: Vec2::broadcast(size),
        };
        {
            let rock_circle = plot::RockCircle::generate(land, &mut reseed(&mut rng), &site, aabr);
            let rock_circle_alt = rock_circle.alt;
            let plot = site.create_plot(Plot {
                kind: PlotKind::RockCircle(rock_circle),
                root_tile: aabr.center(),
                tiles: aabr_tiles(aabr).collect(),
            });

            site.blit_aabr(aabr, Tile {
                kind: TileKind::Building,
                plot: Some(plot),
                hard_alt: Some(rock_circle_alt),
            });
        }
        site
    }

    pub fn generate_troll_cave(land: &Land, rng: &mut impl Rng, origin: Vec2<i32>) -> Self {
        let mut rng = reseed(rng);
        let mut site = Site {
            origin,
            name: None,
            kind: Some(SiteKind::TrollCave),
            ..Site::default()
        };
        let size = 2.0 as i32;
        let aabr = Aabr {
            min: Vec2::broadcast(-size),
            max: Vec2::broadcast(size),
        };
        let site_temp = temp_at_wpos(land, origin);
        {
            let troll_cave =
                plot::TrollCave::generate(land, &mut reseed(&mut rng), &site, aabr, site_temp);
            let troll_cave_alt = troll_cave.alt;
            let plot = site.create_plot(Plot {
                kind: PlotKind::TrollCave(troll_cave),
                root_tile: aabr.center(),
                tiles: aabr_tiles(aabr).collect(),
            });

            site.blit_aabr(aabr, Tile {
                kind: TileKind::Building,
                plot: Some(plot),
                hard_alt: Some(troll_cave_alt),
            });
        }
        site
    }

    pub fn generate_camp(land: &Land, rng: &mut impl Rng, origin: Vec2<i32>) -> Self {
        let mut rng = reseed(rng);
        let mut site = Site {
            origin,
            name: None,
            kind: Some(SiteKind::Camp),
            ..Site::default()
        };
        let size = 2.0 as i32;
        let aabr = Aabr {
            min: Vec2::broadcast(-size),
            max: Vec2::broadcast(size),
        };
        let site_temp = temp_at_wpos(land, origin);
        {
            let camp = plot::Camp::generate(land, &mut reseed(&mut rng), &site, aabr, site_temp);
            let camp_alt = camp.alt;
            let plot = site.create_plot(Plot {
                kind: PlotKind::Camp(camp),
                root_tile: aabr.center(),
                tiles: aabr_tiles(aabr).collect(),
            });

            site.blit_aabr(aabr, Tile {
                kind: TileKind::Building,
                plot: Some(plot),
                hard_alt: Some(camp_alt),
            });
        }
        site
    }

    pub fn generate_cultist(land: &Land, rng: &mut impl Rng, origin: Vec2<i32>) -> Self {
        let mut rng = reseed(rng);
        let mut site = Site {
            origin,
            name: Some({
                let name = NameGen::location(&mut rng).generate();
                match rng.random_range(0..5) {
                    0 => format!("{} Dungeon", name),
                    1 => format!("{} Lair", name),
                    2 => format!("{} Crib", name),
                    3 => format!("{} Catacombs", name),
                    _ => format!("{} Pit", name),
                }
            }),
            kind: Some(SiteKind::Cultist),
            ..Site::default()
        };
        let size = 22.0 as i32;
        let aabr = Aabr {
            min: Vec2::broadcast(-size),
            max: Vec2::broadcast(size),
        };
        {
            let cultist = plot::Cultist::generate(land, &mut reseed(&mut rng), &site, aabr);
            let cultist_alt = cultist.alt;
            let plot = site.create_plot(Plot {
                kind: PlotKind::Cultist(cultist),
                root_tile: aabr.center(),
                tiles: aabr_tiles(aabr).collect(),
            });

            site.blit_aabr(aabr, Tile {
                kind: TileKind::Building,
                plot: Some(plot),
                hard_alt: Some(cultist_alt),
            });
        }
        site
    }

    pub fn generate_sahagin(
        land: &Land,
        index: IndexRef,
        rng: &mut impl Rng,
        origin: Vec2<i32>,
    ) -> Self {
        let mut rng = reseed(rng);
        let mut site = Site {
            origin,
            name: Some({
                let name = NameGen::location(&mut rng).generate();
                match rng.random_range(0..5) {
                    0 => format!("{} Isle", name),
                    1 => format!("{} Islet", name),
                    2 => format!("{} Key", name),
                    3 => format!("{} Cay", name),
                    _ => format!("{} Rock", name),
                }
            }),
            kind: Some(SiteKind::Sahagin),
            ..Site::default()
        };
        let size = 16.0 as i32;
        let aabr = Aabr {
            min: Vec2::broadcast(-size),
            max: Vec2::broadcast(size),
        };
        {
            let sahagin = plot::Sahagin::generate(land, index, &mut reseed(&mut rng), &site, aabr);
            let sahagin_alt = sahagin.alt;
            let plot = site.create_plot(Plot {
                kind: PlotKind::Sahagin(sahagin),
                root_tile: aabr.center(),
                tiles: aabr_tiles(aabr).collect(),
            });

            site.blit_aabr(aabr, Tile {
                kind: TileKind::Building,
                plot: Some(plot),
                hard_alt: Some(sahagin_alt),
            });
        }
        site
    }

    pub fn generate_vampire_castle(land: &Land, rng: &mut impl Rng, origin: Vec2<i32>) -> Self {
        let mut rng = reseed(rng);
        let mut site = Site {
            origin,
            name: Some({
                let name = NameGen::location(&mut rng).generate_vampire();
                match rng.random_range(0..4) {
                    0 => format!("{} Keep", name),
                    1 => format!("{} Chateau", name),
                    2 => format!("{} Manor", name),
                    _ => format!("{} Palace", name),
                }
            }),
            kind: Some(SiteKind::VampireCastle),
            ..Site::default()
        };
        let size = 22.0 as i32;
        let aabr = Aabr {
            min: Vec2::broadcast(-size),
            max: Vec2::broadcast(size),
        };
        {
            let vampire_castle =
                plot::VampireCastle::generate(land, &mut reseed(&mut rng), &site, aabr);
            let vampire_castle_alt = vampire_castle.alt;
            let plot = site.create_plot(Plot {
                kind: PlotKind::VampireCastle(vampire_castle),
                root_tile: aabr.center(),
                tiles: aabr_tiles(aabr).collect(),
            });

            site.blit_aabr(aabr, Tile {
                kind: TileKind::Building,
                plot: Some(plot),
                hard_alt: Some(vampire_castle_alt),
            });
        }
        site
    }

    pub fn generate_bridge(
        land: &Land,
        index: IndexRef,
        rng: &mut impl Rng,
        start_chunk: Vec2<i32>,
        end_chunk: Vec2<i32>,
    ) -> Self {
        let mut rng = reseed(rng);
        let start = TerrainChunkSize::center_wpos(start_chunk);
        let end = TerrainChunkSize::center_wpos(end_chunk);
        let origin = (start + end) / 2;

        let mut site = Site {
            origin,
            name: Some(format!(
                "Bridge of {}",
                NameGen::location(&mut rng).generate_town()
            )),
            kind: Some(SiteKind::Bridge(start_chunk, end_chunk)),
            ..Site::default()
        };

        let start_tile = site.wpos_tile_pos(start);
        let end_tile = site.wpos_tile_pos(end);

        let width = 1;

        let orth = (start_tile - end_tile).yx().map(|dir| dir.signum().abs());

        let start_aabr = Aabr {
            min: start_tile.map2(end_tile, |a, b| a.min(b)) - orth * width,
            max: start_tile.map2(end_tile, |a, b| a.max(b)) + 1 + orth * width,
        };

        let bridge = plot::Bridge::generate(land, index, &mut rng, &site, start_tile, end_tile);

        let start_tile = site.wpos_tile_pos(bridge.start.xy());
        let end_tile = site.wpos_tile_pos(bridge.end.xy());

        let width = (bridge.width() + TILE_SIZE as i32 / 2) / TILE_SIZE as i32;
        let aabr = Aabr {
            min: start_tile.map2(end_tile, |a, b| a.min(b)) - orth * width,
            max: start_tile.map2(end_tile, |a, b| a.max(b)) + 1 + orth * width,
        };

        let line = LineSegment2 {
            start: site.tile_wpos(bridge.dir.select_aabr_with(start_aabr, start_aabr.center())),
            end: site.tile_wpos(
                bridge
                    .dir
                    .opposite()
                    .select_aabr_with(start_aabr, start_aabr.center()),
            ),
        }
        .as_();

        for y in start_aabr.min.y..start_aabr.max.y {
            for x in start_aabr.min.x..start_aabr.max.x {
                let tpos = Vec2::new(x, y);
                let tile_aabr = Aabr {
                    min: site.tile_wpos(tpos),
                    max: site.tile_wpos(tpos + 1) - 1,
                };
                if let Some(tile) = site.tiles.get_mut(tpos) {
                    let closest_point = line.projected_point(tile_aabr.center().as_());
                    let w = TILE_SIZE as f32;
                    if tile_aabr
                        .as_()
                        .projected_point(closest_point)
                        .distance_squared(closest_point)
                        < w.powi(2)
                    {
                        tile.kind = TileKind::Path {
                            closest_pos: closest_point,
                            path: Path { width: w },
                        };
                    }
                }
            }
        }

        let plot = site.create_plot(Plot {
            kind: PlotKind::Bridge(bridge),
            root_tile: start_tile,
            tiles: aabr_tiles(aabr).collect(),
        });

        site.blit_aabr(aabr, Tile {
            kind: TileKind::Bridge,
            plot: Some(plot),
            hard_alt: None,
        });

        site
    }

    pub fn wpos_tile_pos(&self, wpos2d: Vec2<i32>) -> Vec2<i32> {
        (wpos2d - self.origin).map(|e| e.div_euclid(TILE_SIZE as i32))
    }

    pub fn wpos_tile(&self, wpos2d: Vec2<i32>) -> &Tile {
        self.tiles.get(self.wpos_tile_pos(wpos2d))
    }

    pub fn tile_wpos(&self, tile: Vec2<i32>) -> Vec2<i32> { self.origin + tile * TILE_SIZE as i32 }

    pub fn tile_center_wpos(&self, tile: Vec2<i32>) -> Vec2<i32> {
        self.origin + tile * TILE_SIZE as i32 + TILE_SIZE as i32 / 2
    }

    pub fn render_tile(&self, canvas: &mut Canvas, tpos: Vec2<i32>) {
        let tile = self.tiles.get(tpos);
        let twpos = self.tile_wpos(tpos);
        let border = TILE_SIZE as i32;
        let cols = (-border..TILE_SIZE as i32 + border).flat_map(|y| {
            (-border..TILE_SIZE as i32 + border)
                .map(move |x| (twpos + Vec2::new(x, y), Vec2::new(x, y)))
        });
        if let TileKind::Path { closest_pos, path } = &tile.kind {
            let near_connections = CARDINALS.iter().filter_map(|rpos| {
                let npos = tpos + rpos;
                let tile = self.tiles.get(npos);
                let tile_aabr = Aabr {
                    min: self.tile_wpos(tpos).map(|e| e as f32),
                    max: self.tile_wpos(tpos + 1).map(|e| e as f32) - 1.0,
                };
                match tile.kind {
                    TileKind::Road { a, b, w, .. } => {
                        if let Some(PlotKind::Road(road)) = tile.plot.map(|p| &self.plot(p).kind) {
                            let start = road.path.nodes[a as usize];
                            let end = road.path.nodes[b as usize];
                            let dir = Dir2::from_vec2(end - start);
                            let orth = dir.orthogonal();
                            let aabr = Aabr {
                                min: self.tile_center_wpos(start)
                                    - w as i32 * 2 * orth.to_vec2()
                                    - dir.to_vec2() * TILE_SIZE as i32 / 2,
                                max: self.tile_center_wpos(end)
                                    + w as i32 * 2 * orth.to_vec2()
                                    + dir.to_vec2() * TILE_SIZE as i32 / 2,
                            }
                            .made_valid()
                            .as_();
                            Some(aabr)
                        } else {
                            Some(tile_aabr)
                        }
                    },
                    TileKind::Bridge | TileKind::Plaza => Some(tile_aabr),
                    _ => tile
                        .plot
                        .and_then(|plot| self.plot(plot).door_tile())
                        .is_some_and(|door_tile| door_tile == npos)
                        .then_some(tile_aabr),
                }
            });
            cols.for_each(|(wpos2d, _offs)| {
                let wpos2df = wpos2d.map(|e| e as f32);

                if closest_pos.distance_squared(wpos2d.as_()) < path.width.powi(2)
                    || near_connections
                        .clone()
                        .map(|aabr| aabr.distance_to_point(wpos2df))
                        .min_by_key(|d| (*d * 100.0) as i32)
                        .is_some_and(|d| d <= 1.5)
                {
                    let alt = canvas.col(wpos2d).map_or(0, |col| col.alt as i32);
                    let sub_surface_color = canvas
                        .col(wpos2d)
                        .map_or(Rgb::zero(), |col| col.sub_surface_color);
                    for z in -8..6 {
                        let wpos = Vec3::new(wpos2d.x, wpos2d.y, alt + z);
                        canvas.map(wpos, |b| {
                            if b.kind() == BlockKind::Snow {
                                b.into_vacant()
                            } else if b.is_filled() {
                                if b.is_terrain() {
                                    Block::new(
                                        BlockKind::Earth,
                                        path.surface_color((sub_surface_color * 255.0).as_(), wpos),
                                    )
                                } else {
                                    b
                                }
                            } else {
                                b.into_vacant()
                            }
                        })
                    }
                }
            });
        }
    }

    pub fn render(&self, canvas: &mut Canvas, dynamic_rng: &mut ChaCha8Rng) {
        let tile_aabr = Aabr {
            min: self.wpos_tile_pos(canvas.wpos()) - 1,
            max: self
                .wpos_tile_pos(canvas.wpos() + TerrainChunkSize::RECT_SIZE.map(|e| e as i32) + 2)
                + 3, // Round up, uninclusive, border
        };

        // Don't double-generate the same plot per chunk!
        let mut plots = DHashSet::default();

        for y in tile_aabr.min.y..tile_aabr.max.y {
            for x in tile_aabr.min.x..tile_aabr.max.x {
                self.render_tile(canvas, Vec2::new(x, y));

                if let Some(plot) = self.tiles.get(Vec2::new(x, y)).plot {
                    plots.insert(plot);
                }
            }
        }

        canvas.foreach_col(|canvas, wpos2d, col| {
            let tile = self.wpos_tile(wpos2d);
            for z_off in (-2..4).rev() {
                if let Some(plot) = tile.plot.map(|p| self.plot(p)) {
                    canvas.map_resource(
                        wpos2d.with_z(plot.rel_terrain_offset(col) + z_off),
                        |block| {
                            plot.terrain_surface_at(wpos2d, block, dynamic_rng, col, z_off, self)
                                .unwrap_or(block)
                        },
                    );
                }
            }
        });

        // TODO: Solve the 'trees are too big' problem and remove this
        for (id, plot) in self.plots.iter() {
            if matches!(&plot.kind, PlotKind::GiantTree(_)) {
                plots.insert(id);
            }
        }

        let mut plots_to_render = plots.into_iter().collect::<Vec<_>>();
        // First sort by priority, then id.
        plots_to_render.sort_unstable_by_key(|plot| (self.plot(*plot).render_ordering(), *plot));

        let wpos2d = canvas.info().wpos();
        let chunk_aabr = Aabr {
            min: wpos2d,
            max: wpos2d + TerrainChunkSize::RECT_SIZE.as_::<i32>(),
        };

        let info = canvas.info();

        for plot in plots_to_render {
            let (prim_tree, fills, mut entities) = self.plot(plot).render_collect(self, canvas);

            let mut spawn = |pos, last_block| {
                if let Some(entity) = match &self.plots[plot].kind {
                    PlotKind::GiantTree(tree) => tree.entity_at(pos, &last_block, dynamic_rng),
                    _ => None,
                } {
                    entities.push(entity);
                }
            };

            let mut entities_from_structure_blocks = Vec::<EntityInfo>::new();

            for (prim, fill) in fills {
                for mut aabb in Fill::get_bounds_disjoint(&prim_tree, prim) {
                    aabb.min = Vec2::max(aabb.min.xy(), chunk_aabr.min).with_z(aabb.min.z);
                    aabb.max = Vec2::min(aabb.max.xy(), chunk_aabr.max).with_z(aabb.max.z);

                    for x in aabb.min.x..aabb.max.x {
                        for y in aabb.min.y..aabb.max.y {
                            let wpos = Vec2::new(x, y);
                            let col_tile = self.wpos_tile(wpos);
                            if
                            /* col_tile.is_building() && */
                            col_tile
                                .plot
                                .and_then(|p| self.plots[p].z_range())
                                .zip(self.plots[plot].z_range())
                                .is_some_and(|(a, b)| a.end > b.end)
                            {
                                continue;
                            }
                            let mut last_block = None;

                            let col = canvas
                                .col(wpos)
                                .map(|col| col.get_info())
                                .unwrap_or_default();

                            for z in aabb.min.z..aabb.max.z {
                                let pos = Vec3::new(x, y, z);

                                let mut sprite_cfg = None;

                                let map = |block| {
                                    let (current_block, _sb, entity_path) = fill.sample_at(
                                        &prim_tree,
                                        prim,
                                        pos,
                                        &info,
                                        block,
                                        &mut sprite_cfg,
                                        &col,
                                    );

                                    if let Some(spec) = entity_path {
                                        let entity = EntityInfo::at(pos.as_());
                                        let mut loadout_rng = rand::rng();
                                        entities_from_structure_blocks.push(
                                            entity.with_asset_expect(&spec, &mut loadout_rng, None),
                                        );
                                    };

                                    if let (Some(last_block), None) = (last_block, current_block) {
                                        spawn(pos, last_block);
                                    }
                                    last_block = current_block;
                                    current_block.unwrap_or(block)
                                };

                                match fill {
                                    Fill::ResourceSprite { .. } | Fill::Prefab(..) => {
                                        canvas.map_resource(pos, map)
                                    },
                                    _ => canvas.map(pos, map),
                                };

                                if let Some(sprite_cfg) = sprite_cfg {
                                    canvas.set_sprite_cfg(pos, sprite_cfg);
                                }
                            }
                            if let Some(block) = last_block {
                                spawn(Vec3::new(x, y, aabb.max.z), block);
                            }
                        }
                    }
                }
            }

            for entity in entities {
                canvas.spawn(EntitySpawn::Entity(Box::new(entity)));
            }

            for entity in entities_from_structure_blocks {
                canvas.spawn(EntitySpawn::Entity(Box::new(entity)));
            }
        }
    }

    pub fn apply_supplement(
        &self,
        dynamic_rng: &mut impl Rng,
        wpos2d: Vec2<i32>,
        supplement: &mut crate::ChunkSupplement,
    ) {
        for (_, plot) in self.plots.iter() {
            match &plot.kind {
                PlotKind::Gnarling(g) => g.apply_supplement(dynamic_rng, wpos2d, supplement),
                PlotKind::Adlet(a) => a.apply_supplement(dynamic_rng, wpos2d, supplement),
                _ => {},
            }
        }
    }
}

pub fn test_site() -> Site {
    let index = crate::index::Index::new(0);
    let index_ref = IndexRef {
        colors: &index.colors(),
        features: &index.features(),
        index: &index,
    };
    let mut gen_meta = SitesGenMeta::new(0);
    Site::generate_city(
        &Land::empty(),
        index_ref,
        &mut rand::rng(),
        Vec2::zero(),
        0.5,
        None,
        &mut gen_meta,
    )
}

fn wpos_is_hazard(land: &Land, wpos: Vec2<i32>) -> Option<HazardKind> {
    if land
        .get_chunk_wpos(wpos)
        .is_none_or(|c| c.river.near_water())
    {
        Some(HazardKind::Water)
    } else {
        Some(land.get_gradient_approx(wpos))
            .filter(|g| *g > 0.8)
            .map(|gradient| HazardKind::Hill { gradient })
    }
}

fn temp_at_wpos(land: &Land, wpos: Vec2<i32>) -> f32 {
    land.get_chunk_wpos(wpos)
        .map(|c| c.temp)
        .unwrap_or(CONFIG.temperate_temp)
}

pub fn aabr_tiles(aabr: Aabr<i32>) -> impl Iterator<Item = Vec2<i32>> {
    (0..aabr.size().h)
        .flat_map(move |y| (0..aabr.size().w).map(move |x| aabr.min + Vec2::new(x, y)))
}

/// Returns the average gradient of the chunks neighboring the center of the
/// given `aabr`, as well as the gradient of the chunk at the center of the
/// `aabr` itself.
fn get_gradient_average(aabr: Aabr<i32>, land: &Land) -> f32 {
    let chunk_size = TerrainChunkSize::RECT_SIZE.reduce_max() as i32;

    let mut gradient_sum = 0.0;
    let mut gradient_sample_count = 0;

    let aabr_center = aabr.center();
    let range_x_min = aabr_center.x - chunk_size;
    let range_x_max = aabr_center.x + chunk_size;
    let range_y_min = aabr_center.y - chunk_size;
    let range_y_max = aabr_center.y + chunk_size;

    for x_pos in (range_x_min..=range_x_max).step_by(chunk_size as usize) {
        for y_pos in (range_y_min..=range_y_max).step_by(chunk_size as usize) {
            let gradient_at_pos = land.get_gradient_approx(Vec2::new(x_pos, y_pos));
            gradient_sum += gradient_at_pos;
            gradient_sample_count += 1;
        }
    }

    gradient_sum / (gradient_sample_count as f32)
}