vec4 color = tex3D(volume, p);
color.b = tex3D(vessels, p);
color.a = (1.0 - color.b)*0.03;
return color;