I need to correct the reflectance values of LANDSAT-8 data for steep terrain and find it too tedious and clunky in ArcGIS. The data I have is: LANSDAT-8 DEM (ASTER) Slope and aspect (Derived using DEM in ArcGIS) Sun Elevation and Azimuth (Obtained from L8 MTL file)
All of which are clipped to a polygon of my study area.
Two methods I wanted to use were the cosine and minneart corrections. I found a package for R called "landsat" - http://cran.r-project.org/web/packages/landsat/index.html thatcontains a function named "topocorr" - http://www.inside-r.org/packages/cran/landsat/docs/topocorr which takes the uncorrected image, slope, aspect, azimuth and elevation then produces a corrected image.
I'm struggling to use it however as the input data type for the input image is "matrix, data frame, or SpatialGridDataFrame" formats.
How can I convert L8 data into this format?