native.landmass_steradians.argtypes = [ctypes.POINTER(ctypes.c_double), ctypes.c_int, ctypes.c_int, ctypes.c_int, ctypes.c_int]
native.landmass_steradians.restype = ctypes.c_double
landmass_steradians = native.landmass_steradians

const SEA_LEVEL: f64 = 0.1;

    const SEA_LEVEL: f64 = 0.1;

        }
    }
}
#[no_mangle]
pub extern "C" fn landmass_steradians(source: *const f64, width: u32, height: u32, x: u32, y: u32) -> f64 {
    use std::collections::HashMap;
    use std::collections::HashSet;
    use std::collections::VecDeque;
    use std::f64::consts::PI;
    let source = unsafe {
        std::slice::from_raw_parts(source, (width as usize) * (height as usize))
    };
    let mut patch = HashSet::new();
    let mut q = VecDeque::new();
    q.push_back((x,y));
    while let Some((x,y)) = q.pop_front() {
        if source[(y * width + x) as usize] >= SEA_LEVEL {
            if patch.insert((x,y)) {
                let y = y as i32;
                let x = x as i32;
                for (x, y) in [(x - 1, y - 1), (x, y - 1), (x + 1, y - 1), (x + 1, y), (x + 1, y + 1), (x, y + 1), (x - 1, y + 1), (x - 1, y)].iter() {
                    if *y > 0 && (*y as u32) < height {
                        let np = ((x % width as i32) as u32, *y as u32);
                        q.push_back(np);
                    }
                }
            }

    };
    std::mem::drop(q);
    let mut lines: HashMap<u32, u32> = HashMap::new();
    for (_,y) in patch.iter() {
        lines.entry(*y)
            .and_modify(|w| *w += 1)
            .or_insert(1);

    std::mem::drop(patch);
    let dy = PI / (height as f64);
    let dx = PI * 2. / (width as f64);
    lines.iter().map(|(y,w)| dy * ((*y as f64) * dy).sin() * (*w as f64) * dx).sum()

        outputs = tf.math.multiply(255, colourize_heightmap(tf.reshape(outputs, (512, 1024, outputs.shape[1]))))
        image = tensor_to_surface(outputs)

        image = tensor_to_surface(tf.math.multiply(255, colourize_heightmap(tf.reshape(outputs, (512, 1024, outputs.shape[1])))))

              solid_angle = landmass_steradians(tf.reshape(outputs, (512,1024,1)), event.button.x, event.button.y)
              print('Clicked land mass: %s square Km | %s square miles' % (solid_angle * 6372.0**2, solid_angle * 3959**2))

def landmass_steradians(source, x, y):
    shape = source.shape
    source = tf.cast(source, tf.float64).numpy()
    return helpers.landmass_steradians(ctypes.cast(ctypes.c_voidp(source.ctypes.data), ctypes.POINTER(ctypes.c_double)), ctypes.c_int(shape[1]), ctypes.c_int(shape[0]), ctypes.c_int(x), ctypes.c_int(y))