GeoArray.__getitem__() returns wrong shape in case of a 3D array with only one line indexed with a 2D array
Example:
gA = GeoArray(np.random.randint(1, 10, (1, 5, 3))
mask = np.random.randint(0, 2, (1, 5)).astype(bool)
out = gA[mask] # returns array of shape (5, 1) but should return (5, 3)
This only occurrs in case gA holds a 3D array with only one line.
Edited by Daniel Scheffler