Finding Optimal Integral Sampling Lattices for a given Frequency Support in Multidimensions