cpsd_disc function

optimal spatial data discretization based on SPADE q-statistics