A method for deriving Divergence-free wavelets is presented. The approach, in 2D, produces the same wavelets as mentioned by Battle. However, Battle claims that to produce the 3D equivalent construct is still an open problem. Our method solves this problem in a simple and intuitive manner. We discuss potential applications for 2D and 3D divergence-free vector wavelets in image processing problems: particularly using volumetric CT or MRI data.