A data processing device (100, 200) for extracting a desired vital signal, which contains a physiological information component pertaining to a subject of interest, from sensor data that includes time-dependent first sensor data in the form of photoplethysmography data, which is indicative of an amount of electromagnetic radiation reflected from or transmitted through a sensed region of a subject of interest as a function of time in at least one first spectral channel that is sensitive to blood volume variations in the sensed region, the time-dependent first sensor data comprising the physiological information component and at least one motion artifact component, and that includes time-dependent second sensor data that is indicative of a position, a velocity or an acceleration of the sensed region as a function of time in one or more spatial dimensions, wherein the second sensor data forms a sequence of samples each allocated to a respective one of a sequence of time-base values; the data processing device comprising: - a decomposition unit (104, 204), which is configured to receive the second sensor data, to decompose the second sensor data into at least two components of decomposed sensor data and to provide, based on the decomposed second sensor data, at least two different sets of motion reference data (MR11,..., MR1m, MR21,..., MR2m) in at least two different motion reference data channels, and - an artifact removal unit (106, 206), which is configured to receive the first sensor data and the at least two different motion reference data channels and to determine and provide at its output the vital signal formed from a combination of the first sensor data and the motion reference data of at least two of the motion reference data channels, characterized in that the decomposition unit is configured to decompose the second sensor data as a function of time by generating a plurality of time-shifted components (MR11,..., MR1m) from the second sensor data, each time-shifted component being determined from the second sensor data by shifting the allocation to the time-base values by a respective predetermined number of samples, so as to provide the at least two components of decomposed second sensor data (MR11,MR12,..., MR1m); and wherein the artifact removal unit is configured to determine individual weights (W 0, W1,... W L) of the respective spectral and motion reference data channels subject to the combination using a boundary condition, which requires that a weight vector that has as its vector components the individual weights of the spectral channels and of the motion reference data channels to be selected for the combination forms an optimum of a correspondence measure indicative of a correspondence of the vital signal to a prestored normalized correlation vector (b). - a decomposition unit (104, 204), which is configured to receive the second sensor data, to decompose the second sensor data into at least two components of decomposed sensor data and to provide, based on the decomposed second sensor data, at least two different sets of motion reference data (MR11,..., MR1m, MR21,..., MR2m) in at least two different motion reference data channels, and - an artifact removal unit (106, 206), which is configured to receive the first sensor data and the at least two different motion reference data channels and to determine and provide at its output the vital signal formed from a combination of the first sensor data and the motion reference data of at least two of the motion reference data channels, characterized in that the decomposition unit is configured to decompose the second sensor data as a function of time by generating a plurality of time-shifted components (MR11,..., MR1m) from the second sensor data, each time-shifted component being determined from the second sensor data by shifting the allocation to the time-base values by a respective predetermined number of samples, so as to provide the at least two components of decomposed second sensor data (MR11,MR12,..., MR1m); and wherein the artifact removal unit is configured to determine individual weights (W 0, W1,... W L) of the respective spectral and motion reference data channels subject to the combination using a boundary condition, which requires that a weight vector that has as its vector components the individual weights of the spectral channels and of the motion reference data channels to be selected for the combination forms an optimum of a correspondence measure indicative of a correspondence of the vital signal to a prestored normalized correlation vector (b).