hsd2geotiff_py2.py

import himawari.HimawariScene as hsd
from sys import argv
import rasterio as rio
import numpy as np

scaleup = 64

scene = hsd.HimawariScene(argv[1])

rad = scene.radiances()
rad[rad.mask] = 0
rad = np.clip(rad.data, 0, 65535/scaleup)

rad = (rad * scaleup).astype(np.uint16)

meta = {
    'driver': 'GTiff',
    'count': 1,
    'height': rad.shape[0],
    'width': rad.shape[1],
    'dtype': np.uint16,
    'compress': 'LZW'
}

with rio.open(argv[2], 'w', **meta) as dst:
    dst.write(rad, 1)