Preserve shape for simple region reference selections

Previously, region reference selection always resulted in a 1D output shape.  This PR deduces the shape of simple rectangular selections (including strided selections) and provides a correctly-shaped output array.