Bidirectional DEM relief shading method for extraction of gully shoulder line in loess tableland area. (4th July 2018)