docs: def metop20kmto1km(lons20km, lats20km):
      """Getting 1km geolocation for metop avhrr from 20km tiepoints.
      """
      cols20km = np.array([0] + list(range(4, 2048, 20)) + [2047])
      cols1km = np.arange(2048)
      lines = lons20km.shape[0]
      rows20km = np.arange(lines)
      rows1km = np.arange(lines)
  
      along_track_order = 1
      cross_track_order = 3
  
      satint = SatelliteInterpolator((lons20km, lats20km),
                                     (rows20km, cols20km),
                                     (rows1km, cols1km),
                                     along_track_order,
                                     cross_track_order)
      return satint.interpolate()