Function features 

pub fn features(
    p: &Pyramid,
    bi: usize,
    pbi: Option<usize>,
    q: &[i32],
    z: usize,
    y: usize,
    x: usize,
) -> Feat

Gather the features of position (z, y, x) in band bi.