import numpy as np
from PIL import Image
from scipy.ndimage import zoom

def ricostruttore_v3():
    mag = np.load('matrice_sasso_v3.npy')
    dims = np.load('proporzioni.npy')
    
    # Ricostruzione proporzionale
    h_ratio = dims[1] / mag.shape[0]
    w_ratio = dims[0] / mag.shape[1]
    
    # Sintesi non lineare: ricostruiamo i bordi
    # Usiamo un'interpolazione bicubica ad alta densità
    ricostruzione = zoom(mag, (h_ratio, w_ratio), order=3)
    
    # Applichiamo un effetto di "diffusione" per simulare il campo azzurro
    # (In questa versione ricostruiamo la struttura, il colore seguirà)
    final_img = Image.fromarray(ricostruzione.astype(np.uint8))
    final_img.save("volto_strutturale_v3.png")
    print("V3: Ricostruzione proporzionale ultimata.")

if __name__ == "__main__":
    ricostruttore_v3()