tsd icon indicating copy to clipboard operation
tsd copied to clipboard

Published 20 hours ago • verified publisher icon centreborelli

--satellite Landsat-7 option actually downloads Landsat-8 images

Open chlsl opened this issue 5 years ago • 2 comments

I can't manage to download Landsat-7 images with tsd. Using the option --satellite Landsat-7 doesn't seem to work.

With the command python3 tsd/tsd/get_landsat.py --satellite Landsat-7 --geom rio_bounds.geojson --start-date 2000-06-01 --end-date 2000-10-01 -o test, tsd can't find any images, although on earthexplorer I found 9 datasets for the same place and dates.

Changing the dates: python3 tsd/tsd/get_landsat.py --satellite Landsat-7 --geom rio_bounds.geojson --start-date 2013-06-01 --end-date 2013-10-01 -o test, tsd finds 8 images and downloads them:

$ ls -1 test
2013-06-01_L8_OLITIRS_LC81990262013152LGN01_band_B8.tif
2013-06-01_L8_OLITIRS_LC81990262013152LGN01_band_BQA.tif
2013-06-24_L8_OLITIRS_LC82000262013175LGN01_band_B8.tif
2013-06-24_L8_OLITIRS_LC82000262013175LGN01_band_BQA.tif
2013-07-19_L8_OLITIRS_LC81990262013200LGN01_band_B8.tif
2013-07-19_L8_OLITIRS_LC81990262013200LGN01_band_BQA.tif
2013-07-26_L8_OLITIRS_LC82000262013207LGN01_band_B8.tif
2013-07-26_L8_OLITIRS_LC82000262013207LGN01_band_BQA.tif
2013-08-04_L8_OLITIRS_LC81990262013216LGN01_band_B8.tif
2013-08-04_L8_OLITIRS_LC81990262013216LGN01_band_BQA.tif
2013-08-20_L8_OLITIRS_LC81990262013232LGN01_band_B8.tif
2013-08-20_L8_OLITIRS_LC81990262013232LGN01_band_BQA.tif
2013-09-05_L8_OLITIRS_LC81990262013248LGN02_band_B8.tif
2013-09-05_L8_OLITIRS_LC81990262013248LGN02_band_BQA.tif

These images are from Landsat-8, as confirmed by gdalinfo:

$ gdalinfo test/2013-06-01_L8_OLITIRS_LC81990262013152LGN01_band_B8.tif
Driver: GTiff/GeoTIFF
Files: test/2013-06-01_L8_OLITIRS_LC81990262013152LGN01_band_B8.tif
Size is 340, 344
Coordinate System is:
PROJCS["WGS 84 / UTM zone 31N",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",3],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH],
    AUTHORITY["EPSG","32631"]]
Origin = (431572.500000000000000,5409247.500000000000000)
Pixel Size = (15.000000000000000,-15.000000000000000)
Metadata:
  AREA_OR_POINT=Point
  aws_base_url=https://s3-us-west-2.amazonaws.com/landsat-pds/c1/L8/199/026/LC08_L1TP_199026_20130601_20170504_01_T1
  cloud_cover=30
  date=2013-06-01 10:42:53.077784
  downloaded_by=TSD on 2019-10-29T14:19:23.146459
  filename=2013-06-01_L8_OLITIRS_LC81990262013152LGN01
  metadata_source=devseed
  path=199
  product_id=LC08_L1TP_199026_20130601_20170504_01_T1
  row=026
  satellite=L8
  scene_id=LC81990262013152LGN01
  sensor=OLITIRS
  thumbnail=https://s3-us-west-2.amazonaws.com/landsat-pds/c1/L8/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_thumb_large.jpg
  urls={'aws': {'B1': 'https://s3-us-west-2.amazonaws.com/landsat-pds/c1/L8/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B1.TIF', 'B2': 'https://s3-us-west-2.amazonaws.com/landsat-pds/c1/L8/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B2.TIF', 'B3': 'https://s3-us-west-2.amazonaws.com/landsat-pds/c1/L8/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B3.TIF', 'B4': 'https://s3-us-west-2.amazonaws.com/landsat-pds/c1/L8/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B4.TIF', 'B5': 'https://s3-us-west-2.amazonaws.com/landsat-pds/c1/L8/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B5.TIF', 'B6': 'https://s3-us-west-2.amazonaws.com/landsat-pds/c1/L8/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B6.TIF', 'B7': 'https://s3-us-west-2.amazonaws.com/landsat-pds/c1/L8/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B7.TIF', 'B8': 'https://s3-us-west-2.amazonaws.com/landsat-pds/c1/L8/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B8.TIF', 'B9': 'https://s3-us-west-2.amazonaws.com/landsat-pds/c1/L8/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B9.TIF', 'B10': 'https://s3-us-west-2.amazonaws.com/landsat-pds/c1/L8/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B10.TIF', 'B11': 'https://s3-us-west-2.amazonaws.com/landsat-pds/c1/L8/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B11.TIF', 'BQA': 'https://s3-us-west-2.amazonaws.com/landsat-pds/c1/L8/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_BQA.TIF'}, 'gcloud': {'B1': 'https://storage.googleapis.com/gcp-public-data-landsat/LC08/01/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B1.TIF', 'B2': 'https://storage.googleapis.com/gcp-public-data-landsat/LC08/01/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B2.TIF', 'B3': 'https://storage.googleapis.com/gcp-public-data-landsat/LC08/01/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B3.TIF', 'B4': 'https://storage.googleapis.com/gcp-public-data-landsat/LC08/01/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B4.TIF', 'B5': 'https://storage.googleapis.com/gcp-public-data-landsat/LC08/01/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B5.TIF', 'B6': 'https://storage.googleapis.com/gcp-public-data-landsat/LC08/01/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B6.TIF', 'B7': 'https://storage.googleapis.com/gcp-public-data-landsat/LC08/01/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B7.TIF', 'B8': 'https://storage.googleapis.com/gcp-public-data-landsat/LC08/01/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B8.TIF', 'B9': 'https://storage.googleapis.com/gcp-public-data-landsat/LC08/01/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B9.TIF', 'B10': 'https://storage.googleapis.com/gcp-public-data-landsat/LC08/01/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B10.TIF', 'B11': 'https://storage.googleapis.com/gcp-public-data-landsat/LC08/01/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_B11.TIF', 'BQA': 'https://storage.googleapis.com/gcp-public-data-landsat/LC08/01/199/026/LC08_L1TP_199026_20130601_20170504_01_T1/LC08_L1TP_199026_20130601_20170504_01_T1_BQA.TIF'}}
Image Structure Metadata:
  COMPRESSION=DEFLATE
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  431572.500, 5409247.500) (  2d 4' 3.30"E, 48d49'56.76"N)
Lower Left  (  431572.500, 5404087.500) (  2d 4' 6.39"E, 48d47' 9.67"N)
Upper Right (  436672.500, 5409247.500) (  2d 8'13.44"E, 48d49'58.70"N)
Lower Right (  436672.500, 5404087.500) (  2d 8'16.31"E, 48d47'11.61"N)
Center      (  434122.500, 5406667.500) (  2d 6' 9.86"E, 48d48'34.21"N)
Band 1 Block=340x12 Type=UInt16, ColorInterp=Gray

Note, the rio_bounds.geojson file used in the commands above:

$ cat rio_bounds.geojson | jq                                                                                                                                     ✘ 2
{
  "bbox": [
    2.0685128788396763,
    48.78649911256507,
    2.137141255915403,
    48.832360818215854
  ],
  "features": [
    {
      "bbox": [
        2.0685128788396763,
        48.78649911256507,
        2.137141255915403,
        48.832360818215854
      ],
      "geometry": {
        "coordinates": [
          [
            [
              2.0685128788396763,
              48.78649911256507
            ],
            [
              2.137141255915403,
              48.78649911256507
            ],
            [
              2.137141255915403,
              48.832360818215854
            ],
            [
              2.0685128788396763,
              48.832360818215854
            ],
            [
              2.0685128788396763,
              48.78649911256507
            ]
          ]
        ],
        "type": "Polygon"
      },
      "properties": {
        "filename": "2018-01-01_S2A_orbit_051_tile_31UDQ_L1C_band_B02.tif",
        "id": "0",
        "title": "./L1C/2018-01-01_S2A_orbit_051_tile_31UDQ_L1C_band_B02.tif"
      },
      "type": "Feature"
    }
  ],
  "type": "FeatureCollection"
}

chlsl avatar Oct 29 '19 14:10 chlsl

Thank you @chlsl for reporting this issue. For Landsat-7 and older Landsats the only available catalog is GCloud, which you can use by adding the option --api gcloud to your command (but it requires credentials).

carlodef avatar Oct 29 '19 15:10 carlodef

Thank you for your answer @carlodef.

chlsl avatar Oct 29 '19 15:10 chlsl