First of all, I think I understand now, why you consider the frequency.
p = hbar k = hbar (2Pi / lambda)
And if p changes, does lambda also change? (Did you have the same thought?)
Seems strange to me, now that I think about it...I get an headache

.
Because it would be somehow odd if the photon changed its wavelength after passing the slit. This would violate energy conservation.
I also ask myself now if my explanation with the slit is correct
because the photon gets a momentum 'kick' due to the position
measurement. But where does this kick come from (hmm...momentum conservation). I think it's rather the DIRECTION of the momentum vector
that is changed but not its magnitude. Therefore the x and y component of p
are only distributed in a another way. (My second reply to your post
doesnt seem correct then, because px is not changed in my explanation)
2 things that I also don't understand:
Firstly: The change of distribution in x and y means, there's not only
a Delta_py for position measurement Delta_y, but also
a Delta_px. Quite strange for me, because it's not contained
in the uncertainty equation, that Delta_px occurs too.
Secondly: it also follows that Delta_y is indeed bounded, like you've said.
I also don't understand this because the uncertainty equation says
Delta_y * Delta_py > hbar/2. So if I made the slit smaller and smaller,
then Delta_py would go to infinity...
Hmm...I hope someone else can clear things up, lol.
Nevertheless, thanks for this question, it's very inspiring.
-Edgardo