1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41

use geo::Coord;
use geo::CoordFloat;
use num_traits::FloatConst;

use super::builder::types::BuilderAntimeridianResampleNoClip;
use super::builder_conic::Builder;
use super::builder_conic::ParallelsSet;
use super::equal_area::EqualArea;
use super::CenterSet;
use super::RawBase;
use super::RotateSet;
use super::ScaleSet;
use super::TranslateSet;

type Output<DRAIN, T> = Builder<BuilderAntimeridianResampleNoClip<DRAIN, EqualArea<T>, T>, T>;

/// Albers - [``ConicEqualArea``](crate::projection::conic_equal_area::ConicEqualArea) centered on the U.S.
///
/// # Panics
/// unwrap() is used here but a panic will never happen as constants  will
/// always be converted into T.
#[must_use]
pub fn albers<DRAIN, T>() -> Output<DRAIN, T>
where
    DRAIN: Clone,
    T: 'static + CoordFloat + Default + FloatConst,
{
    let mut b = EqualArea::<T>::builder();
    b.parallels_set(T::from(29.5_f64).unwrap(), T::from(45.5_f64).unwrap())
        .scale_set(T::from(1070_f64).unwrap())
        .translate_set(&Coord {
            x: T::from(480_f64).unwrap(),
            y: T::from(250_f64).unwrap(),
        })
        .rotate2_set(&[T::from(96_f64).unwrap(), T::zero()])
        .center_set(&Coord {
            x: T::from(-0.6_f64).unwrap(),
            y: T::from(38.7_f64).unwrap(),
        });
    b
}