simple mask test

Tianshu Cheng
1 parent 95272d90
Showing 11 changed files with 9 additions and 6 deletions Show diff stats
fun_mPb_theta.cpp
gauss_derivative_odd.cpp
test result/0829_01_vs_gradient_theta=0.PNG
test result/0829_02_vs_gradient_theta=0.PNG
test result/0829_03_vs_gradient_theta=0_r=3,5,10.PNG
test result/0829_03_vs_gradient_theta=90_r=3,5,10.PNG
test result/0829_04_compare_colora_theta=90_r=3,5,10.PNG
test result/0829_04_compare_colorb_theta=90_r=3,5,10.PNG
test result/0829_04_compare_vs_gradient_theta=90_r=3,5,10.PNG
test result/0831_01_compare_vs_mPb.PNG
test result/0831_01_compare_vs_mPb_time.PNG
@@ -28,8 +28,8 @@ stim::image&lt;float&gt; func_mPb_theta(stim::image&lt;float&gt; lab, float theta, unsigned 
  
 	stim::image<float> l1,l2,l3,a1,a2,a3,b1,b2,b3;
  
-	l1 = gaussian_derivative_filter_odd(pic_light, sigma, sigma_n, r1 * 2, theta, w, h);
-	stim::cpu2image(l1.data(), "data_output/l1_tex2.bmp", w, h, stim::cmBrewer);
+	l1 = gaussian_derivative_filter_odd(pic_light, sigma, sigma_n, r3 * 2, theta, w, h);
+	stim::cpu2image(l1.data(), "data_output/testl_tex5.bmp", w, h, stim::cmBrewer);
 	l2 = gaussian_derivative_filter_odd(pic_light, sigma, sigma_n, r2 * 2, theta, w, h);
 	stim::cpu2image(l2.data(), "data_output/l2_tex.bmp", w, h, stim::cmBrewer);
 	l3 = gaussian_derivative_filter_odd(pic_light, sigma, sigma_n, r3 * 2, theta, w, h);
@@ -12,12 +12,14 @@ void array_multiply(float* lhs, float rhs, unsigned int N);
 // winsize = 2 * r, side of mask = winsize + 1
 stim::image<float> gaussian_derivative_filter_odd(stim::image<float> image, float sigma, unsigned int sigma_n, unsigned int winsize, float theta, unsigned int w, unsigned int h){
  
-	stim::image<float> mask_x(winsize+1, winsize+1), mask_y(winsize+1, winsize+1), mask_theta(winsize+1, winsize+1), derivative_x, derivative_y, derivative_theta(w, h);
+	stim::image<float> mask_x(winsize+1, winsize+1), mask_y(winsize+1, winsize+1), mask_theta(winsize+1, winsize+1), mask_test(winsize+1, winsize+1, 1), derivative_x, derivative_y, derivative_theta(w, h);
 	//float* ptr = mask_x.data();
  
-	//mask_x.load("101087.bmp");
-	//float s[169];
-	//float *ptr = s;
+
+	memset ( mask_test.data(), 0,  mask_test.size() * sizeof(float));
+	mask_test.data()[winsize*(winsize+2)/2] = 1;
+	stim::cpu2image(mask_test.data(), "data_output/mask_test.bmp", winsize+1, winsize+1, stim::cmBrewer);
+
  
 	// set parameters
 	unsigned N = w * h;
@@ -51,6 +53,7 @@ stim::image&lt;float&gt; gaussian_derivative_filter_odd(stim::image&lt;float&gt; image, floa
 	// 2D convolution
 	//derivative_theta = image.convolve2(mask_theta);
 	//stim::cpu2image(derivative_theta.data(), "data_output/derivative_theta1.bmp", w, h, stim::cmBrewer);
+	//conv2(image.data(), mask_test.data(), derivative_theta.data(), w, h, winsize+1);
 	conv2(image.data(), mask_theta.data(), derivative_theta.data(), w, h, winsize+1);
 	//stim::cpu2image(derivative_theta.data(), "data_output/derivative_theta_tex1.bmp", w, h, stim::cmBrewer);