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:
Comments
Please sign in to leave a comment.