Decision trees have become a widely adapted tool for mapping soil properties in geographic space. However, it is problematic to implement spatial relationships in the models. We present a new method which uses geographic coordinates along several axes tilted at oblique angles in the models. We test this method on four spatial datasets. The results show that the new method is at least as accurate as other proposed alternatives, has a computational advantage and is flexible and interpretable.