mirror of
https://github.com/qgis/QGIS.git
synced 2025-02-23 00:02:38 -05:00
54 lines
1.9 KiB
Common Lisp
54 lines
1.9 KiB
Common Lisp
|
#include "calcfirstder.cl"
|
||
|
|
|
||
|
|
// Aspect renderer for QGIS
|
||
|
|
|
||
|
|
__kernel void processNineCellWindow( __global float *scanLine1,
|
||
|
|
__global float *scanLine2,
|
||
|
|
__global float *scanLine3,
|
||
|
|
__global uchar4 *resultLine,
|
||
|
|
__global float *rasterParams
|
||
|
|
) {
|
||
|
|
|
||
|
|
// Get the index of the current element
|
||
|
|
const int i = get_global_id(0);
|
||
|
|
|
||
|
|
// Do the operation
|
||
|
|
//return (( (x31 - x11) + 2 * (x32 - x12) + (x33 - x13) ) / (8 * mCellSizeX))
|
||
|
|
float derX = calcFirstDer( scanLine1[i], scanLine2[i], scanLine3[i],
|
||
|
|
scanLine1[i+1], scanLine2[i+1], scanLine3[i+1],
|
||
|
|
scanLine1[i+2], scanLine2[i+2], scanLine3[i+2],
|
||
|
|
rasterParams[0], rasterParams[1], rasterParams[2], rasterParams[3]
|
||
|
|
);
|
||
|
|
//return (((x11 - x13) + 2 * (x21 - x23) + (x31 - x33)) / ( 8 * mCellSizeY));
|
||
|
|
float derY = calcFirstDer( scanLine1[i+2], scanLine1[i+1], scanLine1[i],
|
||
|
|
scanLine2[i+2], scanLine2[i+1], scanLine2[i],
|
||
|
|
scanLine3[i+2], scanLine3[i+1], scanLine3[i],
|
||
|
|
rasterParams[0], rasterParams[1], rasterParams[2], rasterParams[4]
|
||
|
|
);
|
||
|
|
|
||
|
|
|
||
|
|
float res;
|
||
|
|
if ( derX == rasterParams[1] || derY == rasterParams[1] ||
|
||
|
|
( derX == 0.0f && derY == 0.0f) )
|
||
|
|
{
|
||
|
|
res = rasterParams[1];
|
||
|
|
}
|
||
|
|
else
|
||
|
|
{
|
||
|
|
// 180.0 / M_PI = 57.29577951308232
|
||
|
|
float aspect = atan2( derX, derY ) * 57.29577951308232;
|
||
|
|
if ( aspect < 0 )
|
||
|
|
res = 90.0f - aspect;
|
||
|
|
else if (aspect > 90.0f)
|
||
|
|
// 360 + 90 = 450
|
||
|
|
res = 450.0f - aspect;
|
||
|
|
else
|
||
|
|
res = 90.0 - aspect;
|
||
|
|
}
|
||
|
|
|
||
|
|
res = res / 360 * 255;
|
||
|
|
|
||
|
|
resultLine[i] = (uchar4)(res, res, res, 255);
|
||
|
|
}
|
||
|
|
|