//stim::image<float> gaussian_derivative_filter_odd(stim::image<float> image, int r, unsigned int sigma_n, float theta);
//stim::image<float> func_mPb_theta(stim::image<float> img, float theta, int* r, float* alpha, int s);
//stim::image<float> func_mPb(stim::image<float> img, unsigned int theta_n, int* r, float* alpha, int s);

stim::image<float> Gd1(stim::image<float> image, int r, unsigned int sigma_n);
stim::image<float> Gd2(stim::image<float> image, int r, unsigned int sigma_n);
stim::image<float> Gd_odd(stim::image<float> image, int r, unsigned int sigma_n, float theta);
stim::image<float> Gd_even(stim::image<float> image, int r, unsigned int sigma_n, float theta);
stim::image<float> Gd_center(stim::image<float> image, int r, unsigned int sigma_n);

stim::image<float> textons(stim::image<float> image, unsigned int theta_n);
stim::image<float> kmeans(stim::image<float> textons, unsigned int K);
stim::image<float> Pb(stim::image<float> image, int r, unsigned int sigma_n);
stim::image<float> cPb(stim::image<float> img, int* r, float* alpha, int s);
stim::image<float> tPb(stim::image<float> img, int* r, float* alpha, unsigned int theta_n, unsigned int bin_n, int s, unsigned int K);