Extracting reflectance coefficients from a SkySat analytic image by using Python

Before you can convert to reflectance values, you must extract the conversion coefficients included in the GeoTIFF Header, under ‘TIFFTAG_IMAGEDESCRIPTION’ in the SkySat analytic image.

import rasterio
import json

filename = rasterio.open('20210913_211131_ssc15_u0001_analytic.tif')

tags = filename.tags()['TIFFTAG_IMAGEDESCRIPTION']
jtags = json.loads(tags)
coeffs = jtags['properties']['reflectance_coefficients']
# print ("reflectance coefficients:", coeffs)

for coeff in [0, 1, 2, 3]:
print('reflectance coefficient band', coeff + 1, ':' , coeffs[coeff]

Screen_Shot_2021-09-24_at_9.30.00_PM.png

 

For more information about Convert Radiance values to Reflectance please check this document:

 

War dieser Beitrag hilfreich?
0 von 0 fanden dies hilfreich

Kommentare

0 Kommentare

Bitte melden Sie sich an, um einen Kommentar zu hinterlassen.