# -*- 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()
