GeoArray.read_pointData() returns wrong array shape in case of 2D GeoArray instance and offside coordinates

from geoarray import GeoArray
gA = GeoArray(np.array([[1, 2], [3, 4]], ..., ...)
mapXY_points = np.array([[..., ...]])  # offsite coordinate
gA.read_pointData(mapXY_points)  # returns 3D array instead of (n_coordinates x 1)