Now we save the intermediate Hessians as a pickle, as converting them to numpy arrays breaks all hell. Theres a new option to pass an output name (otherwise, it gives a default filename)
ImageOp Hessian calculations appear to work but take forever :(