How to intepret slope on Temperatures rasters

Alejandro Hinojosa Corona alhinc at cicese.mx
Mon Feb 10 19:17:46 EST 1997


Hello grass colleagues:

I am applying the slope formula proposed by Shapiro and Westervelt in the
"R.MAPCALC An Algebra for GIS and Image Processing" document to
a raster map of Sea surface temperatures.
The raster is in lat-long with a cell size of 5 minutes in both latitude
and longitud. The raster values are temperatures coded as 10ths of a degree
Celsius.

 The formula applied is:

slope = eval( \
              x= (DTM[-1,-1] + 2 * DTM[0,-1] + DTM[1,-1]          \
                 -DTM[-1,1]  - 2 * DTM[0,1]  - DTM[1,1])/8.0,    \                       
              y= (DTM[-1,-1] + 2 * DTM[-1,0] + DTM[-1,1]    \
                 -DTM[1,-1] -  2 * DTM[1,0]  - DTM[1,1])/ (8.0 ), \
             atan(sqrt(x * x + y *y)))    
 
The output  is a slope raster with values in degrees
ranging from  0 to 90.
My question is:
     -how to interpret raster output as change in temperature per cell.
     -how to transform the slope value in degrees to an increment in temperature
      per unit cell.

 Best Regards

Alejandro Hinojosa    
alhinc at cicese.mx    http://www.cicese.mx/~alhinc
Research Associate
Earth Sciences Division
CICESE 
Ensenada, Baja California, Mexico. 



More information about the grass-user mailing list