Apologies for the flurry, one more change:
float GDALSlopeHydroAlg (float* afWin, float fDstNoDataValue, void* pData)
{
// Hydrologic Slope is the max
// local slope btw center cell and adjacent cells
const double radiansToDegrees = 180.0 / M_PI;
GDALSlopeAlgData* psData = (GDALSlopeAlgData*)pData;
double key;
float pafLocalMax = afWin[0];
for ( int k = 1; k < 9; k++)
{
if (afWin[k] > pafLocalMax)
{
pafLocalMax=afWin[k];
}
}
key = pafLocalMax - afWin[4];
if (psData->slopeFormat == 1)
return (float) (atan(sqrt(key) / (2*psData->scale)) *
radiansToDegrees);
else
return (float) (100*(sqrt(key) / (2*psData->scale)));
}
Best,
Steve
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
