We develop the noncommutative harmonic space (NHS) analysis to study the problem of solving the non-linear constraints eqs of noncommutative Yang-Mills selfduality in four-dimensions. We show that this space, denoted also as NHS(η, θ), has two SU(2) isovector deformations η (ij) and θ (ij) parametrising respectively two noncommutative harmonic subspaces NHS(η,0) and NHS(0,θ) used to study the selfdual and anti self-dual noncommutative Yang-Mills solutions. We reformulate the Yang-Mills self-dual constraint eqs on NHS(η,0) by extending the idea of harmonic analyticity to linearize them. Then we give a perturbative self dual solution recovering the ordinary one. Finally we present the explicit computation of an exact self-dual solution.