I'm working on project where I'm using Independent Components Analysis (ICA) on NxN patches from images to separate components in an image signal. Each column of the resulting MxN^2 unmixing matrix is reshaped back to an NxN filter containing the independent components. For my project, this filter is then convolved with new images for feature extraction purposes.
My question is, how do I remove the DC offset of these NxN filters that result from ICA? In order words how can I ensure that (in MATLAB parlance)
sum(sum(filter)) == 0 returns
1? Is it fair to just subtract the filter's sum from each component, or will that affect the filter's response?
If the filter sums to some non-zero value X, is it appropriate to just subtract
X/number_of_coefficients from each filter coefficient?