a python tool to process landsat images
TRANSCRIPT
![Page 1: A Python Tool to Process Landsat Images](https://reader030.vdocument.in/reader030/viewer/2022012423/61779fcb1ce0ca3e5a354687/html5/thumbnails/1.jpg)
A Python Tool to Process Landsat Images
William J. Kramber Idaho Department of Water Resources
Esri International Users Conference, San Diego, CA July 22, 2015
Landsat 8 NASA/USGS
![Page 2: A Python Tool to Process Landsat Images](https://reader030.vdocument.in/reader030/viewer/2022012423/61779fcb1ce0ca3e5a354687/html5/thumbnails/2.jpg)
Landsat USGS/NASA mission L1 launched 1972
Current satellites
L7 launched 1999 L8 launched 2013
15/30/100-120 meter (Pan/Multispectral/Thermal)
16 day cycle 106 by 115 miles Landsat 9 launch 2023
NASA
NASA
Landsat 7
![Page 3: A Python Tool to Process Landsat Images](https://reader030.vdocument.in/reader030/viewer/2022012423/61779fcb1ce0ca3e5a354687/html5/thumbnails/3.jpg)
Landsat at IDWR Since 1975
Land use and land cover mapping Water rights evaluation Conservation Reserve Enhancement Program Mapping Evapotranspiration (ET)
Hydrologic modeling Endangered Species Planning Water Delivery Call Curtailment Order
CREP
2006 Evapotranspiration
Eastern Snake Plain Aquifer Model
2010 Irrigated Land Use
![Page 4: A Python Tool to Process Landsat Images](https://reader030.vdocument.in/reader030/viewer/2022012423/61779fcb1ce0ca3e5a354687/html5/thumbnails/4.jpg)
Landsat Archive
USGS EROS Center, Sioux Falls, SD Over 4 million scenes July 1972 to present (thermal since 1984) Free Search and download sites:
Earth Explorer (earthexplorer.usgs.gov) GloVis (glovis.usgs.gov)
![Page 5: A Python Tool to Process Landsat Images](https://reader030.vdocument.in/reader030/viewer/2022012423/61779fcb1ce0ca3e5a354687/html5/thumbnails/5.jpg)
Download Landsat data (path 40 rows 30 and 31, 2014)
![Page 6: A Python Tool to Process Landsat Images](https://reader030.vdocument.in/reader030/viewer/2022012423/61779fcb1ce0ca3e5a354687/html5/thumbnails/6.jpg)
Process Landsat 8 Tool
![Page 7: A Python Tool to Process Landsat Images](https://reader030.vdocument.in/reader030/viewer/2022012423/61779fcb1ce0ca3e5a354687/html5/thumbnails/7.jpg)
Starting 1 of 26
![Page 8: A Python Tool to Process Landsat Images](https://reader030.vdocument.in/reader030/viewer/2022012423/61779fcb1ce0ca3e5a354687/html5/thumbnails/8.jpg)
Finished 26 of 26 (11.5 hours)
Dell OptiPlex 790 3.4 GHz Core i7 processor 16 GB memory Two 1 TB hard drives AMD 1GB Radeon 6450 dual video card
![Page 9: A Python Tool to Process Landsat Images](https://reader030.vdocument.in/reader030/viewer/2022012423/61779fcb1ce0ca3e5a354687/html5/thumbnails/9.jpg)
Output Folders
![Page 10: A Python Tool to Process Landsat Images](https://reader030.vdocument.in/reader030/viewer/2022012423/61779fcb1ce0ca3e5a354687/html5/thumbnails/10.jpg)
p40r30_l8_20140610
![Page 11: A Python Tool to Process Landsat Images](https://reader030.vdocument.in/reader030/viewer/2022012423/61779fcb1ce0ca3e5a354687/html5/thumbnails/11.jpg)
p40r30_l8_20140610
![Page 12: A Python Tool to Process Landsat Images](https://reader030.vdocument.in/reader030/viewer/2022012423/61779fcb1ce0ca3e5a354687/html5/thumbnails/12.jpg)
Tool Summary Create list of gz files for gz in gzList: Check Directory and create Subdirectory Copy gz to Subdirectory Uncompress gz and tar files Project band 8 and BQA to IDTM83 Project other bands to IDTM83 Trim bands to 0 where any of the bands has 0 Create mask of 0 and 1 Create 10-band composite Create 3-band (564 and 543) composites Apply 2.2 SD Contrast Stretch to 3-band composites Compute TOA Reflectance Compute NDVI End of Tool ArcGIS Tools used
ProjectRaster CellStatistics Con SetNull CompositeBands ExtractByMask CalculateStatistics CopyRaster
![Page 13: A Python Tool to Process Landsat Images](https://reader030.vdocument.in/reader030/viewer/2022012423/61779fcb1ce0ca3e5a354687/html5/thumbnails/13.jpg)
Set Null: change gaps to nodata Zonal Statistics: compute mean of polygons Is Null & Con: replace gaps with mean value
L7 gaps filled using ArcGIS Tools
Gaps in Landsat 7 Gaps filled
![Page 14: A Python Tool to Process Landsat Images](https://reader030.vdocument.in/reader030/viewer/2022012423/61779fcb1ce0ca3e5a354687/html5/thumbnails/14.jpg)
For more information Earth Explorer earthexplorer.usgs.gov IDWR idwr.idaho.gov Mapping Evapotranspiration idwr.idaho.gov/GeographicInfo/METRIC/et.htm
Seasonal Evapotranspiration for Cropland in 2000
![Page 15: A Python Tool to Process Landsat Images](https://reader030.vdocument.in/reader030/viewer/2022012423/61779fcb1ce0ca3e5a354687/html5/thumbnails/15.jpg)
Landsat 9
2023 or better yet 2021
![Page 16: A Python Tool to Process Landsat Images](https://reader030.vdocument.in/reader030/viewer/2022012423/61779fcb1ce0ca3e5a354687/html5/thumbnails/16.jpg)
Future work
ProcessLandsat8 Tool Pan-sharpen 3-band composites
Smart Eliminate Tool for ArcGIS
MRLC NLCD Mapping Tool www.mrlc.gov/resources.php