log_gamma 1106 modules/imgproc/src/lsd.cpp double log1term = (double(n) + 1) - log_gamma(double(k) + 1) log_gamma 1107 modules/imgproc/src/lsd.cpp - log_gamma(double(n-k) + 1)