A computer implemented method of processing seismic image of the subsurface, said seismic image comprising pixels, each pixel having a color value, said seismic image comprising columns of pixels, wherein the method comprises: - determining a first surface and a second surface in the seismic image, the first surface and the second surface being determined based on geological information extracted from well data from at least one well (202, 203, 204) in the subsurface; - determining first dip components for pixels of the first surface and second dip components for pixels of the second surface; characterized in that the method further comprises: - for each pixel of the seismic image between the first surface and the second surface, determining interpolated dip components based on at least the first dip components and the second dip components. - for some pixels of the seismic image between the first surface and the second surface, determining an interpolated surface coming through the pixel and based on at least the interpolated dip components for pixels between the first surface and the second surface; - for each pixel of the seismic image between the first surface and the second surface, determining a density value based on a number of intersections of the determined surfaces with a column of the seismic image comprising said pixel, said intersections being below said pixel in said column; - determining a geological-time image comprising columns of pixels, the columns of said geological-time image are computed based: - on the density values of pixels of a corresponding column of the seismic image; - on the color value of pixels of the corresponding column of the seismic image. - determining a first surface and a second surface in the seismic image, the first surface and the second surface being determined based on geological information extracted from well data from at least one well (202, 203, 204) in the subsurface; - determining first dip components for pixels of the first surface and second dip components for pixels of the second surface; - for each pixel of the seismic image between the first surface and the second surface, determining interpolated dip components based on at least the first dip components and the second dip components. - for some pixels of the seismic image between the first surface and the second surface, determining an interpolated surface coming through the pixel and based on at least the interpolated dip components for pixels between the first surface and the second surface; - for each pixel of the seismic image between the first surface and the second surface, determining a density value based on a number of intersections of the determined surfaces with a column of the seismic image comprising said pixel, said intersections being below said pixel in said column; - determining a geological-time image comprising columns of pixels, the columns of said geological-time image are computed based: - on the density values of pixels of a corresponding column of the seismic image; - on the color value of pixels of the corresponding column of the seismic image. - on the density values of pixels of a corresponding column of the seismic image; - on the color value of pixels of the corresponding column of the seismic image.