An automated method for assessment of cell density and biomarker expressions in an interface zone IZ of tissues of different classes wherein at least a first tissue and a second tissue interact, the method comprising steps: a) acquiring data from digital image analysis DIA of one or more microscopy images containing tissue sections pixel-wise classified to at least said first and second tissue classes, and bio-markers expressed by immune and/or other cells/particles wherein the DIA data contains types and coordinates of said bio-markers, b) subsampling the provided DIA data by a regular grid wherein all grid elements are equally-sized and equally-spaced covering the entire DIA image, c) detecting a tissue edge TE across the grid elements, wherein a TE class is assigned to grid elements in which an interface between the the first and second tissue classes is determined; and classifying grid elements not on the TE to a tissue class having the maximum area fraction in that grid element, d) augmenting each grid element h with a rank(h) based on a spatial ranking of all grid elements by their shortest distance to the detected TE, wherein the rank(h) assigned according to the tissue class of h is as follows: rank h = def if class h = tissue of first class, if class h = tissue of second class, if class h = TE, rank = shortest distance to TE rank = − shortest distance to TE; rank = 0 e) extracting the interface zone IZ between the first and second tissue classes by selecting an IZ width including ranks on both sides of the TE, f) calculating quantities of the bio-markers in grid elements of the extracted IZ ranks, and further, assessing one or more profiles of the calculated biomarker quantities across the IZ using any gradient-type indicators of: Centre of Mass CM q = ∑ r i r i q r i ∑ r i q r i, where i indexes all ranks r i across a selected IZ, and q ( r i) denotes the calculated bio-marker quantities in grid elements at the ranks r i of said IZ, and Immunodrop ID q = q r − 1 q r 1 which represents calculated bio-marker quantity q ( r) change at the TE between ranks r -1 and r 1, in a selected IZ. a) acquiring data from digital image analysis DIA of one or more microscopy images containing tissue sections pixel-wise classified to at least said first and second tissue classes, and bio-markers expressed by immune and/or other cells/particles wherein the DIA data contains types and coordinates of said bio-markers, b) subsampling the provided DIA data by a regular grid wherein all grid elements are equally-sized and equally-spaced covering the entire DIA image, c) detecting a tissue edge TE across the grid elements, wherein a TE class is assigned to grid elements in which an interface between the the first and second tissue classes is determined; and classifying grid elements not on the TE to a tissue class having the maximum area fraction in that grid element, d) augmenting each grid element h with a rank(h) based on a spatial ranking of all grid elements by their shortest distance to the detected TE, wherein the rank(h) assigned according to the tissue class of h is as follows: rank h = def if class h = tissue of first class, if class h = tissue of second class, if class h = TE, rank = shortest distance to TE rank = − shortest distance to TE; rank = 0 e) extracting the interface zone IZ between the first and second tissue classes by selecting an IZ width including ranks on both sides of the TE, f) calculating quantities of the bio-markers in grid elements of the extracted IZ ranks, and further, assessing one or more profiles of the calculated biomarker quantities across the IZ using any gradient-type indicators of: Centre of Mass CM q = ∑ r i r i q r i ∑ r i q r i, where i indexes all ranks r i across a selected IZ, and q ( r i) denotes the calculated bio-marker quantities in grid elements at the ranks r i of said IZ, and Immunodrop ID q = q r − 1 q r 1 which represents calculated bio-marker quantity q ( r) change at the TE between ranks r -1 and r 1, in a selected IZ. Centre of Mass CM q = ∑ r i r i q r i ∑ r i q r i, where i indexes all ranks r i across a selected IZ, and q ( r i) denotes the calculated bio-marker quantities in grid elements at the ranks r i of said IZ, and Immunodrop ID q = q r − 1 q r 1 which represents calculated bio-marker quantity q ( r) change at the TE between ranks r -1 and r 1, in a selected IZ.