Hi DataCubers,
A quick question about discovering what data are available in the datacube (Landsat).
First let me clarify that I am not using the Python API at the moment.
What I do is search by filename. For example if I want L8 for a given tile and a given year and month, I go to the corresponding folder and look all the filenames that match my requirements.
In most cases I get a right answer but in some cases there are files missing
Example: for L8 in the tile 145 -35 and January 2015, I go to /g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS8_OLI_TIRS/145_-035/2015/
And then search using the string _2015-01-
In that particular case I find 3 NBAR :
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS8_OLI_TIRS/145_-035/2015/LS8_OLI_TIRS_NBAR_145_-035_2015-01-08T00-08-18.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS8_OLI_TIRS/145_-035/2015/LS8_OLI_TIRS_NBAR_145_-035_2015-01-17T00-02-28.tif
- But there are no PQA files for them!! -
Another example, LS7, 2015, December, tile 145, -36
I find only two NBAR
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_NBAR_145_-036_2015-12-02T00-09-34.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_NBAR_145_-036_2015-12-02T00-09-58.tif
But 5 FC (which need NBAR to be done first)
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_FC_145_-036_2015-12-02T00-09-34.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_FC_145_-036_2015-12-02T00-09-58.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_FC_145_-036_2015-12-11T00-03-34.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_FC_145_-036_2015-12-11T00-03-58.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_FC_145_-036_2015-12-27T00-03-47.tif
And 8 PQA!!
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_PQA_145_-036_2015-12-02T00-09-34.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_PQA_145_-036_2015-12-02T00-09-58.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_PQA_145_-036_2015-12-11T00-03-34.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_PQA_145_-036_2015-12-11T00-03-58.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_PQA_145_-036_2015-12-18T00-09-51.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_PQA_145_-036_2015-12-18T00-10-15.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_PQA_145_-036_2015-12-27T00-03-47.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_PQA_145_-036_2015-12-27T00-04-11.tif
Any idea why and if I am doing anything wrong? Is there a python function that does this?
Or is there any other easier way to find what’s available?
Cheers
Juan
Hi DataCubers,
A quick question about discovering what data are available in the datacube (Landsat).
First let me clarify that I am not using the Python API at the moment.
What I do is search by filename. For example if I want L8 for a given tile and a given year and month, I go to the corresponding folder and look all the filenames that match my requirements.
In most cases I get a right answer but in some cases there are files missing
Example: for L8 in the tile 145 -35 and January 2015, I go to /g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS8_OLI_TIRS/145_-035/2015/
And then search using the string _2015-01-
In that particular case I find 3 NBAR :
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS8_OLI_TIRS/145_-035/2015/LS8_OLI_TIRS_NBAR_145_-035_2015-01-08T00-08-18.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS8_OLI_TIRS/145_-035/2015/LS8_OLI_TIRS_NBAR_145_-035_2015-01-17T00-02-28.tif
Another example, LS7, 2015, December, tile 145, -36
I find only two NBAR
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_NBAR_145_-036_2015-12-02T00-09-34.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_NBAR_145_-036_2015-12-02T00-09-58.tif
But 5 FC (which need NBAR to be done first)
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_FC_145_-036_2015-12-02T00-09-34.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_FC_145_-036_2015-12-02T00-09-58.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_FC_145_-036_2015-12-11T00-03-34.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_FC_145_-036_2015-12-11T00-03-58.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_FC_145_-036_2015-12-27T00-03-47.tif
And 8 PQA!!
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_PQA_145_-036_2015-12-02T00-09-34.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_PQA_145_-036_2015-12-02T00-09-58.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_PQA_145_-036_2015-12-11T00-03-34.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_PQA_145_-036_2015-12-11T00-03-58.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_PQA_145_-036_2015-12-18T00-09-51.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_PQA_145_-036_2015-12-18T00-10-15.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_PQA_145_-036_2015-12-27T00-03-47.tif
/g/data2/rs0/tiles/EPSG4326_1deg_0.00025pixel/LS7_ETM/145_-036/2015/LS7_ETM_PQA_145_-036_2015-12-27T00-04-11.tif
Any idea why and if I am doing anything wrong? Is there a python function that does this?
Or is there any other easier way to find what’s available?
Cheers
Juan