code: def get_xy_1D(ds, stride=1, getval=False):
      
      gt = ds.GetGeoTransform()
      #stride = stride_m/gt[1]
      pX = np.arange(0, ds.RasterXSize, stride)
      pY = np.arange(0, ds.RasterYSize, stride)
      mX, dummy = pixelToMap(pX, pY[0], gt)
      dummy, mY = pixelToMap(pX[0], pY, gt)
      return mX, mY
docs: Return 1D arrays of x and y map coordinates for input GDAL Dataset