This code finds the Landsat WRS-2 path/row which contains a given latitude/longitude co-ordinate.

Dependencies

• GDAL/ORG
• Shapely
• Landsat WRS-2 Path/Row Shapefiles - download from USGS site

How to use

>>> from get_wrs import ConvertToWRS
>>> conv = ConvertToWRS()
>>> conv.get_wrs(50.14, -1.7)
[{'path': 202, 'row': 25}]
>>> conv.get_wrs(50.14, -1.43)
[{'path': 201, 'row': 25}, {'path': 202, 'row': 25}]

If you can't put the shapefile in the same folder as the code, then you can pass the path to the shapefile as an argument - for example: ConvertToWRS(shapefile="\path\to\shapefile.shp").

More details are available in the docstrings (try ConvertToWRS? or conv.get_wrs? in IPython) and in my original blog post.

License

This code is released under the 3-clause BSD license (see the LICENSE file for more information).