Hi Micha,

Thank you for reporting this.

I tried to get the same error with a large area, but I could not reproduce it. My process worked fine and downloaded a bigger area than your 75km limit.

My download took many hours, and SITS successfully renewed the MPC access tokens when needed.

This means the problem is most likely related to the version of SITS you are using. A newer version of SITS probably fixed the limit you found.

I suggest the following action plan for you:

1. Update SITS
   The first thing you should do is update your SITS package to the newest version from CRAN. This is the simplest way to fix the problem. Use this command in R to update SITS:

install.packages("sits")

After you update, please try your original code again.

1. Alternative strategy if the issue continues:
   Try to reduce the total size of the request by asking for data year by year. This avoids any single large request that the MPC API might still reject. Make sure you copy all the data into the same local directory. Once all the annual data is copied to the local folder, SITS can open the full 10-year data cube from your local disk using the

sits_cube(source = "MPC", ..., local_dir = "your_local_folder")

command. See more details on local cubes in SITS documentation: ?sita_cube

Let us know if the update fixes the issue!

I second partitioning the search to be annual and then merging the results client-side. I'm not familiar with this sits library, but I know that enumerating large numbers of STAC items through the STAC API is prone to failures if the service experiences extreme load during your enumeration.

If you need to enumerate all of the STAC items for Landsat over a period of time, I suggest you instead use the stac-geoparquet exported version of the collection to bypass the API altogether: https://planetarycomputer.microsoft.com/docs/quickstarts/stac-geoparquet/

This may not be easy to do in R, I understand if that's not an option.