pub fn morton_encode(x: u32, y: u32, z: u32) -> u32
Compute a 30-bit Morton code from 3D integer coordinates.
Each coordinate must be in [0, 1023].
[0, 1023]