Segmenting white matter and grey matter is a nontrivial task. There is probably literature describing different methods to do this.
You should note as well that imread does not return the image in HU, even dicomread has to be scaled. SO the thresholding you're doing here is not actually connected to what you describe.
