# -*- coding: utf-8 -*- # script01.py # 2014.11.20 K. Kuwata from image import landsat # Band1 limage = landsat("/Users/ken/Data/Landsat/LE71070352000329EDC00_B3.TIF", "/Users/ken/Data/Landsat/LE71070352000329EDC00_MTL.txt", 7, 3) array = limage.Band2Array() Parameters = limage.ReturnParameters() Radiance = limage.DN2Radiance(array, Parameters) Reflectance3 = limage.Radiance2Reflectance(Radiance, Parameters) limage.WriteArrayAsImage("Reflectance_B1.tif", Reflectance3) # Band3 limage = landsat("/Users/ken/Data/Landsat/LE71070352000329EDC00_B4.TIF", "/Users/ken/Data/Landsat/LE71070352000329EDC00_MTL.txt", 7, 4) array = limage.Band2Array() Parameters = limage.ReturnParameters() Radiance = limage.DN2Radiance(array, Parameters) Reflectance4 = limage.Radiance2Reflectance(Radiance, Parameters)
# -*- coding: utf-8 -*- # script01.py # 2014.11.20 K. Kuwata from image import landsat name= "LE70240322008153EDC00" # Calculating reflectance of Band1 limage = landsat("/home/faizan/USA_data/landsat/4/LE70240322008153EDC00_B1.TIF", "/home/faizan/USA_data/landsat/4/LE70240322008153EDC00_MTL.txt", 7, 1) # Get pixel data as array array = limage.Band2Array() # Get information from a parameter file Parameters = limage.ReturnParameters() # Calculate Radiance from DN Radiance = limage.DN2Radiance(array, Parameters) # Calculate Reflectance from Radiance Reflectance1 = limage.Radiance2Reflectance(Radiance, Parameters) # Make a Geotiff file limage.WriteArrayAsImage("4_Reflectance_B1_B3.tif", Reflectance1) # Calculating reflectance of Band3 limage = landsat("/home/faizan/USA_data/landsat/4/LE70240322008153EDC00_B3.TIF", "/home/faizan/USA_data/landsat/4/LE70240322008153EDC00_MTL.txt", 7, 3) # Get pixel data as array array = limage.Band2Array() # Get information from a parameter file Parameters = limage.ReturnParameters()