A method, comprising: obtaining input data as 2D seismic data or 3D cubes of seismic data from data sampling locations from a seismic survey, along with corresponding dip and azimuth data of a geological layer; selecting a data sampling location as a central computation location, defining an operator size for an adaptive structure-oriented operator to be created, respectively in a dip direction (Rdip), strike direction (Rstrike) and normal direction (Rnorm); converting the dip and azimuth data of the data sampling locations to displacement vector coordinates projected on the crossline, inline, and Z axis, and translating the dip and azimuth data into discrete steps for each input sample; determining a normalized dip vector, representing a projected distance covered on the crossline, inline and Z axis for each step in the dip direction; determining normalized strike vectors (S1, S2), representing a projected distance covered on the crossline, inline and Z axis for each step in the strike direction; determining normal vector, representing a projected distance covered on the crossline, inline and Z axis for each step in the normal direction; creating, with a computer, an adaptive structure-oriented operator for the central computational location; wherein creating the adaptive structure-oriented operator, comprises: - starting from the central computational location, selecting neighboring samples (n) in the normal, strike and dip direction by using the determined normalized dip vector, normalized strike vectors (S1, S2) and normal vector; - starting from one or more of the selected neighboring samples (n), selecting new neighboring samples (n) in the normal, strike and dip direction by using the determined normalized dip vector, normalized strike vectors (S1, S2) and normal vector; - stopping to select neighboring samples when the samples distance from the central computational location corresponds to the defined size for the operator respectively in the dip direction (Rdip), strike direction (Rstrike), and normal direction (Rnorm); populating the operator by associating to each of the selected samples (n) a corresponding amplitude from the seismic input data; construct an image of a subsurface using the amplitude values from the seismic input data associated to the operator; using the constructed image of a subsurface in exploration and recovery of hydrocarbons. obtaining input data as 2D seismic data or 3D cubes of seismic data from data sampling locations from a seismic survey, along with corresponding dip and azimuth data of a geological layer; selecting a data sampling location as a central computation location, defining an operator size for an adaptive structure-oriented operator to be created, respectively in a dip direction (Rdip), strike direction (Rstrike) and normal direction (Rnorm); converting the dip and azimuth data of the data sampling locations to displacement vector coordinates projected on the crossline, inline, and Z axis, and translating the dip and azimuth data into discrete steps for each input sample; determining a normalized dip vector, representing a projected distance covered on the crossline, inline and Z axis for each step in the dip direction; determining normalized strike vectors (S1, S2), representing a projected distance covered on the crossline, inline and Z axis for each step in the strike direction; determining normal vector, representing a projected distance covered on the crossline, inline and Z axis for each step in the normal direction; creating, with a computer, an adaptive structure-oriented operator for the central computational location; wherein creating the adaptive structure-oriented operator, comprises: - starting from the central computational location, selecting neighboring samples (n) in the normal, strike and dip direction by using the determined normalized dip vector, normalized strike vectors (S1, S2) and normal vector; - starting from one or more of the selected neighboring samples (n), selecting new neighboring samples (n) in the normal, strike and dip direction by using the determined normalized dip vector, normalized strike vectors (S1, S2) and normal vector; - stopping to select neighboring samples when the samples distance from the central computational location corresponds to the defined size for the operator respectively in the dip direction (Rdip), strike direction (Rstrike), and normal direction (Rnorm); - starting from the central computational location, selecting neighboring samples (n) in the normal, strike and dip direction by using the determined normalized dip vector, normalized strike vectors (S1, S2) and normal vector; - starting from one or more of the selected neighboring samples (n), selecting new neighboring samples (n) in the normal, strike and dip direction by using the determined normalized dip vector, normalized strike vectors (S1, S2) and normal vector; - stopping to select neighboring samples when the samples distance from the central computational location corresponds to the defined size for the operator respectively in the dip direction (Rdip), strike direction (Rstrike), and normal direction (Rnorm); populating the operator by associating to each of the selected samples (n) a corresponding amplitude from the seismic input data; construct an image of a subsurface using the amplitude values from the seismic input data associated to the operator; using the constructed image of a subsurface in exploration and recovery of hydrocarbons.