Cokriging for large spatial datasets: mapping soil properties at a region scale from airborne hyperspectral imagery