Image Edge Enhancing Coherence Filter Toolbox: [Jxx, Jxy, Jxz, Jyy, Jyz, Jzz]=StructureTensor3D(ux,uy,uz, rho) - File Exchange

Code covered by the BSD License

Highlights from
Image Edge Enhancing Coherence Filter Toolbox

CoherenceFilter(u,Options) This function COHERENCEFILTER will perform Anisotropic Diffusion of a
D=derivatives(I,option) Sobel like derivatives with Scharr rotation invariance stencil notations.
I=imgaussian(I,sigma,siz) IMGAUSSIAN filters an 1D, 2D color/greyscale or 3D image with an
[Dxx,Dxy,Dxz,Dyy,Dyz,Dzz]=Str... From Structure Tensor to Diffusion Tensor
[Dxx,Dxy,Dyy]=ConstructDiffus... Construct the edge preserving diffusion tensor D = [a,b;b,c] such as
[Iout]=ut_gauss(varargin) ut_gauss - 2-D filtering using Gaussian masks
[Jxx, Jxy, Jxz, Jyy, Jyz, Jzz... This function calculates the 3D Structure Tensor
[Jxx, Jxy, Jyy]=StructureTens... This function calculates the 2D
[Lambda1,Lambda2,I2x,I2y,I1x,... This function computes the eigenvectors and eigen values of the 2D image
[mu3,mu2,mu1,v3x,v3y,v3z,v2x,... This function calculates the eigen values and vectors, of the 3D image hessian.
[x11,x12,x22]=hessian(x)
a=clmp(a,mina,maxa) Function clmp to boundaries
diffusion_scheme_3D_standard(... Standard Discretization of 3D image diffusion.
err=error_diffusion_scheme_2D... Parameters must be positive. Otherwise the kernels will go to a more
err=error_diffusion_scheme_3D... Parameters must be positive. Otherwise the kernels will go to a more
err=fitgaussian(I,a,t)
fminlbfgs.m
optimize_scheme_2D.m This script finds optimized kernel values for in the novel discretization
optimize_scheme_3D.m Initial values for unknowns in optimal derivative kernels
showcs3(varargin) SHOWCS3 M-file for showcs3.fig
u=diffusion_scheme_2D_implici... Diffusion scheme as introduced by Weickert "Anisotropic Diffusion
u=diffusion_scheme_2D_non_neg... The Basic non_negativity diffusion equation. (Can be found in "A Scheme for
u=diffusion_scheme_2D_novel(u... par=[0.007520981141059, 0.049564810649554, 0.031509665995882, ...
u=diffusion_scheme_2D_rotatio... Most diffusion discretizations are not rotation-invariant, probably
u=diffusion_scheme_2D_standar... The standard diffusion equation. (Can be found in "A Scheme for
u=diffusion_scheme_3D_implici... This 3D Diffusion scheme is based on introduced by Weickert
u=diffusion_scheme_3D_non_neg... This is a basic non_negativity discretization of the 3D diffusion
u=diffusion_scheme_3D_novel(u...
u=diffusion_scheme_3D_novel_g...
u=diffusion_scheme_3D_rotatio... Most diffusion discretizations are not rotation-invariant. This is an
CoherenceFilterStep2D.m
CoherenceFilterStep3D.m
compile_c_files.m
View all files

from Image Edge Enhancing Coherence Filter Toolbox by Dirk-Jan Kroon
Advanced 2D/3D noise removal and edge enhancing with anisotropic diffusion filtering ( Weickert )

[Jxx, Jxy, Jxz, Jyy, Jyz, Jzz]=StructureTensor3D(ux,uy,uz, rho)

function [Jxx, Jxy, Jxz, Jyy, Jyz, Jzz]=StructureTensor3D(ux,uy,uz, rho)
% This function calculates the 3D Structure Tensor
% Jp ( grad(u) ) = conv ( Kp  , grad(u) * grad(u)^T )
% 
% Function is written by D.Kroon University of Twente (September 2009)

% J(grad u_sigma)
Jxx = ux.^2;
Jxy = ux.*uy;
Jxz = ux.*uz;
Jyy = uy.^2;
Jyz = uy.*uz;
Jzz = uz.^2;

% Do the gaussian smoothing
Jxx = imgaussian(Jxx,rho,4*rho);
Jxy = imgaussian(Jxy,rho,4*rho);
Jxz = imgaussian(Jxz,rho,4*rho);
Jyy = imgaussian(Jyy,rho,4*rho);
Jyz = imgaussian(Jyz,rho,4*rho);
Jzz = imgaussian(Jzz,rho,4*rho);

Highlights from Image Edge Enhancing Coherence Filter Toolbox

Highlights from
Image Edge Enhancing Coherence Filter Toolbox