angular_cell_size_deg = 7;
rho_cell_size = 0.1;
rho_min = -1;
rho_max = 1;

base = h3d_create_struct(angular_cell_size_deg, rho_cell_size, rho_min, rho_max);


points = h3d_test_shape('pyramid');
points = [points rot_x(deg2rad(-30))*points];

R1 = rot_z(deg2rad(-30))*rot_y(deg2rad(-30));

s1 = h3d_compute_cube_ht(base, R1*points);


figure;

for face=1:6
	subplot(1,6,face);
	imagesc(squeeze(s1.cube_hs(face,:,:)));
end

figure;