My guess is that your 3dNwarpApply is effectively just removing the obliquity info from the header and keeping everything on the closest cardinal equivalent usually. You can do the same with "3drefit -deoblique". Note, the data in this case is not transformed; only the header info is removed.
