This function splits a certain phase into two subphases based on a size threshold, which makes sense in the reconstruction framework if the two subphases behave different structurally.
A 2D array.
The phase to modify.
The threshold for splitting (in pixels).
The phase to change pixels to.
The modified 2D array.