docs: def modis1kmto250m(lons1km, lats1km, cores=1):
      """Getting 250m geolocation for modis from 1km tiepoints.
  
      http://www.icare.univ-lille1.fr/tutorials/MODIS_geolocation
      """
      if cores > 1:
          return _multi(modis1kmto250m, lons1km, lats1km, 10, cores)
  
      cols1km = np.arange(1354)
      cols250m = np.arange(1354 * 4) / 4.0
  
      along_track_order = 1
      cross_track_order = 3
  
      lines = lons1km.shape[0]
      rows1km = np.arange(lines)
      rows250m = (np.arange(lines * 4) - 1.5) / 4.0
  
      satint = SatelliteInterpolator((lons1km, lats1km),
                                     (rows1km, cols1km),
                                     (rows250m, cols250m),
                                     along_track_order,
                                     cross_track_order,
                                     chunk_size=40)
      satint.fill_borders("y", "x")
      lons250m, lats250m = satint.interpolate()
  
      return lons250m, lats250m