| 47 | } |
| 48 | |
| 49 | std::vector<double> GaussianBlur::boxesForGauss(double sigma, int n) { |
| 50 | double wIdeal = sqrt((12 * sigma * sigma / n) + 1); |
| 51 | int wl = floor(wIdeal); |
| 52 | |
| 53 | if(wl % 2 == 0) |
| 54 | wl--; |
| 55 | |
| 56 | const int wu = wl + 2; |
| 57 | |
| 58 | const double mIdeal = (12 * sigma * sigma - n * wl * wl - 4 * n * wl - 3 * n) / (-4 * wl - 4); |
| 59 | const int m = round(mIdeal); |
| 60 | |
| 61 | std::vector<double> sizes; |
| 62 | |
| 63 | for(int i = 0; i < n; i++) { |
| 64 | sizes.push_back(i < m ? wl : wu); |
| 65 | } |
| 66 | |
| 67 | return sizes; |
| 68 | } |
| 69 | |
| 70 | void GaussianBlur::boxBlur(IntensityMap &input, IntensityMap &result, double radius, bool tileable) { |
| 71 | for(int i = 0; i < input.getWidth() * input.getHeight(); i++) { |
nothing calls this directly
no outgoing calls
no test coverage detected