First, I get the image of the bands of the desired area from the site https://earthexplorer.usgs.gov/, then I import the images into Mathmatica by addressing them as below:

Select dataset directory:

imgdir = 
    "Sacramento Area Data 07-10-2018"}];

Import bands without Atmospheric Correction in DN (with Radiometric correction):

{bandIR, bandR, bandG, bandB} = 
       "LC08_L1TP_044033_20181007_20200830_02_T1_B" <> ToString[#] <> 
        ".TIF"}], "Data"] & /@ {5, 4, 3, 2};

Show green band raw image:

5 ImageAdjust[Image[bandG, ImageSize -> Medium]]

enter image description here

This operation is not automatic and must be downloaded from the Usgs and then entered into Mathematica by addressing. My question is, is it possible to go directly from Mathematica to https://earthexplorer.usgs.gov/ and get a satellite image?

  • $\begingroup$ There was this API that has LandsatImages as a request option, but I don't know anything about and have never used it mathematica.stackexchange.com/a/163639/72953 $\endgroup$
    – ydd
    Commented Aug 15, 2023 at 15:44
  • 1
    $\begingroup$ It seems that USGS provide an API to access their data inventories. There is also a python library which you may use to call from Mathematica ... Read the documentation of the API, and then use HTTPRequest etc. to access it. $\endgroup$
    – Domen
    Commented Aug 16, 2023 at 7:51


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.