A model-independent Dalitz plot analysis of B<sup>±</sup>→DK<sup>±</sup> with D→KS0h+h- (h=π, K) decays and constraints on the CKM angle γ