The method is used for creating an all-in-focus image of a portion of an iridocorneal angle of an eye starting from a plurality of images (610), these images being shallow field-depth images of that portion taken with different focus planes but with a same image frame; the method comprises the steps of: A) pre-processing (620) said plurality of images so to obtain a plurality of single-channel images if the images of said plurality are color images, B) calculating (630) an initial focus map starting from said plurality of single-channel images, C) creating (640) an initial estimated all-in-focus image starting from said initial focus map and said plurality of single-channel images, D) calculating (650) an initial edge map starting from said initial focus map and said plurality of single-channel images, E) iterating (660) a regularization process until a stop regularization condition is met, wherein the regularization process comprises removing one or more edge points of the initial edge map by considering a reference model of an iridocorneal angle of an eye, F) creating (670) said all-in-focus image starting from said plurality of images and a focus map last calculated in step E; the reference model expresses at least one geometrical parameter of a surface of an iridocorneal angle of an eye.