GeoArray of harmonized output does not contain GeoTransform, projection and other metadata.
When calling SpectralHomogenizer as documented, the metadata from the input dataset is not copied to the output data.
from spechomo import SpectralHomogenizer
print(image_l8.gt)
# get an instance of SpectralHomogenizer class:
SH = SpectralHomogenizer()
# run the spectral homogenization
image_s2, errors = SH.predict_by_machine_learner(
arrcube=image_l8[:,:,:7],
method='LR',
n_clusters=50,
src_satellite='Landsat-8',
src_sensor='OLI_TIRS',
src_LBA=['1', '2', '3', '4', '5', '6', '7'], # must be passed as list of strings and match the band numbers of the input
print(image_s2.gt)
)
(380100.0, 30.0, 0.0, 3696300.0, 0.0, -30.0)
[0, 1, 0, 0, 0, -1]