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') #item ID: SkySat

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]


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

Was this article helpful?
0 out of 0 found this helpful



Please sign in to leave a comment.