find Landsat Product ID from old Scene ID

The whole Landsat archive has been reprocessed and renamed, this file finds the new product names from a list of old sceneIDs

The USGS hosts the largest Landsat archive and has reprocessed these scenes. Prior to this reprocessing, the "scene-id" was used. However, now a "product-id" is used. Searching through the archive is done through USGS EarthExplorer (https://earthexplorer.usgs.gov/), but this portal is only capable of searching for the newer "product-id"s. Because the new "product-id" has an additional processing date, it is not possible to transfer codes prior to the processing, to the new format.

Thus this script looks at a text file with old "scene-id"s and downloads the metadata for the associated Landsat satellites. Hereafter, it mates the list with the database and produces a resulting text file with "product-id"s.

Beware! Because the metadata files of Landsat are very big, it might take some time on your computer to process.

A "scene-id" looks like: LXSSPPPRRRYYYYMMDD
A "product-id" looks like: LXSS_LLLL_PPPRRR_YYYYMMDD_yyyymmdd_CC_TX

Where:
L = Landsat
X = Sensor (“C”=OLI/TIRS combined, “O”=OLI-only, “T”=TIRS-only, “E”=ETM+, “T”=“TM, “M”=MSS)
SS = Satellite (”07”=Landsat 7, “08”=Landsat 8)
LLL = Processing correction level (L1TP/L1GT/L1GS)
PPP = WRS path
RRR = WRS row
YYYYMMDD = Acquisition year, month, day
yyyymmdd - Processing year, month, day
CC = Collection number (01, 02, …)
TX = Collection category (“RT”=Real-Time, “T1”=Tier 1, “T2”=Tier 2)