Based on our newly developed methods and the XMM-Newton large program of SN1006, we extract and analyze the spectra from 3596 tessellated regions of this SNR each with 0.3-8 keV counts $>10^4$. For the first time, we map out multiple physical parameters, such as the temperature ($kT$), electron density ($n_e$), ionization parameter ($n_et$), ionization age ($t_{ion}$), metal abundances, as well as the radio-to-X-ray slope ($alpha$) and cutoff frequency ($ u_{cutoff}$) of the synchrotron emission. We construct probability distribution functions of $kT$ and $n_et$, and model them with several Gaussians, in order to characterize the average thermal and ionization states of such an extended source. We construct equivalent width (EW) maps based on continuum interpolation with the spectral model of each regions. We then compare the EW maps of OVII, OVIII, OVII K$delta-zeta$, Ne, Mg, SiXIII, SiXIV, and S lines constructed with this method to those constructed with linear interpolation. We further extract spectra from larger regions to confirm the features revealed by parameter and EW maps, which are often not directly detectable on X-ray intensity images. For example, O abundance is consistent with solar across the SNR, except for a low-abundance hole in the center. This O Hole has enhanced OVII K$delta-zeta$ and Fe emissions, indicating recently reverse shocked ejecta, but also has the highest $n_et$, indicating forward shocked ISM. Therefore, a multi-temperature model is needed to decompose these components. The asymmetric metal distributions suggest there is either an asymmetric explosion of the SN or an asymmetric distribution of the ISM.