Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

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

Open
chlsl opened this issue Oct 29, 2019 · 2 comments
Open

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

chlsl opened this issue Oct 29, 2019 · 2 comments

Comments

@chlsl
Copy link
Collaborator

chlsl commented Oct 29, 2019

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"
}
@carlodef
Copy link
Member

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).

@chlsl
Copy link
Collaborator Author

chlsl commented Oct 29, 2019

Thank you for your answer @carlodef.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants