A method for performing super-resolution processing of a low resolution input data structure (So) of digital 1 D, 2D or 3D data, comprising steps of - filtering the input data structure (S 0) by a first low-pass filter (F l,0), wherein a low-frequency input data structure (L 0) is obtained; - calculating a difference between the input data structure (S 0) and the low-frequency input data structure (L 0), whereby a high-frequency input data structure (H 0) is generated; - upscaling the input data structure (So), and filtering the upscaled input data structure by a second low-pass filter (F l,1), wherein a low-frequency upscaled data structure (L 1) is obtained; - determining in the low-frequency upscaled data structure (L 1) a first patch (P n,L1) at a first position; - searching (152,154) in the low-frequency input data structure (L 0) a first block (B n,L0) that matches the first patch (P n,L1) best, and determining the position of said first block (B n,L0) within the low-frequency input data structure (L 0); - selecting a second block (B n,H0) in the high-frequency input data structure (H 0) at the position that was determined for said first block (B n,L0) within the low-frequency input data structure (L 0); - accumulating, by pixelwise addition, pixel data of the selected second block (B n,H0) to a second patch (P n,H1), the second patch being a patch in a high-frequency upscaled data structure (H 1,acc) at the first position, wherein the high-frequency upscaled data structure (H 1,acc) is initially empty; - repeating the steps of determining a new patch (P n,L1) in the low-frequency upscaled data structure (L 1), searching (152,154) in the low-frequency input data structure (L 0) a block (B n,L0) that matches the selected patch (P n,L1) best, selecting a corresponding block (B n,H0) in the high-frequency input data structure (H 0) and accumulating pixel data of the selected corresponding block (B n,H0) to a patch (P n,H1) in the high-frequency upscaled data structure (H 1,acc) at the position of said new patch (P n,L1); - normalizing the accumulated pixel values in the high-frequency upscaled data structure (H 1,acc), whereby a normalized high-frequency upscaled data structure (H 1) is obtained and wherein the step of normalizing comprises, for each pixel, dividing the accumulated value of said pixel by the number of contributions of pixel values to the accumulated value of said pixel; and - adding the normalized high-frequency upscaled data structure (H 1) to the low-frequency upscaled data structure (L 1), whereby a super-resolved data structure (S 1) is obtained. - filtering the input data structure (S 0) by a first low-pass filter (F l,0), wherein a low-frequency input data structure (L 0) is obtained; - calculating a difference between the input data structure (S 0) and the low-frequency input data structure (L 0), whereby a high-frequency input data structure (H 0) is generated; - upscaling the input data structure (So), and filtering the upscaled input data structure by a second low-pass filter (F l,1), wherein a low-frequency upscaled data structure (L 1) is obtained; - determining in the low-frequency upscaled data structure (L 1) a first patch (P n,L1) at a first position; - searching (152,154) in the low-frequency input data structure (L 0) a first block (B n,L0) that matches the first patch (P n,L1) best, and determining the position of said first block (B n,L0) within the low-frequency input data structure (L 0); - selecting a second block (B n,H0) in the high-frequency input data structure (H 0) at the position that was determined for said first block (B n,L0) within the low-frequency input data structure (L 0); - accumulating, by pixelwise addition, pixel data of the selected second block (B n,H0) to a second patch (P n,H1), the second patch being a patch in a high-frequency upscaled data structure (H 1,acc) at the first position, wherein the high-frequency upscaled data structure (H 1,acc) is initially empty; - repeating the steps of determining a new patch (P n,L1) in the low-frequency upscaled data structure (L 1), searching (152,154) in the low-frequency input data structure (L 0) a block (B n,L0) that matches the selected patch (P n,L1) best, selecting a corresponding block (B n,H0) in the high-frequency input data structure (H 0) and accumulating pixel data of the selected corresponding block (B n,H0) to a patch (P n,H1) in the high-frequency upscaled data structure (H 1,acc) at the position of said new patch (P n,L1); - normalizing the accumulated pixel values in the high-frequency upscaled data structure (H 1,acc), whereby a normalized high-frequency upscaled data structure (H 1) is obtained and wherein the step of normalizing comprises, for each pixel, dividing the accumulated value of said pixel by the number of contributions of pixel values to the accumulated value of said pixel; and - adding the normalized high-frequency upscaled data structure (H 1) to the low-frequency upscaled data structure (L 1), whereby a super-resolved data structure (S 1) is obtained.