1.基于Weibull分布的混合高阶变分超声图像去噪方法，其特征在于包括以下步骤：步骤1：对超声图像进行噪声分布拟合，采用Weibull概率密度函数，对选定的图像局部均匀区域进行灰度直方图拟合，并对拟合分布参数采用最大似然法进行估计；

步骤2：根据拟合的噪声分布，推导对数真实无噪图像的似然函数，将似然函数进行负的对数变换，从而得到最小化能量函数的数据保真项；

步骤3：将对数真实无噪图像一阶正则化表达式和二阶正则化表达式进行加权融合，构建最小化能量函数的混合高阶正则化项，且权值能根据图像边缘特征自适应调整；

步骤4：将步骤2的数据保真项和步骤3的混合高阶正则化项形成最小化能量函数模型，采用Split‑Bregman迭代方法进行快速求解，迭代收敛后，得到对数去噪超声图像，最后通过指数变换得到去噪超声图像；

所述步骤1中，真实无噪图像v，受到乘性噪声n的干扰，生成观察到的超声图像f，且满足超声图像成像模型f＝vn，其中，乘性噪声n在图像的局部灰度均匀区域近似为Weibull分布，即噪声的概率密度函数满足：其中，exp()为指数函数；λ和k分别为Weibull分布的尺度参数和形状参数，选取超声图像的局部灰度均匀区域，获得直方图拟合的样本数据，利用最大似然法对尺度参数λ和形状参数k进行估计；

所述步骤2中，对超声图像成像模型f＝vn，进行对数变换可得logf＝logv+logn；令g＝logf，u＝logv和w＝logn，则有g＝u+w；根据乘性噪声n的概率密度函数pN(n)推导对数变换后的噪声w满足：进一步，求得对数变换后的对数真实无噪图像u的似然函数为：将似然函数进行负的对数变换，且忽略常数部分，可得：

‑k k(g‑u)

‑logpG|U(g|u)＝‑k(g‑u)+λ e ；

假定图像区域Ω内的所有像素g(x,y)∈Ω满足独立同分布，则能构建最小化能量函数的数据保真项如下：所述步骤3中，将对数真实无噪图像u的一阶正则化表达式 和二阶正则化表达式 进行加权组合，得到混合高阶正则化项为：其中， 分别代表x和y方向的偏导数，||·

||代表向量或矩阵的Euclid范数，加权系数γ1、γ2根据对数真实无噪图像u的边缘特征自适应变化，具体表达式为：γ2＝min(γ1)+max(γ1)‑γ1；

其中，η为控制函数平滑程度的阈值，目标轮廓越明显取值越大；

所述步骤4中，将步骤2的数据保真项和步骤3的混合高阶正则化项组合，引入数据保真项权系数α，得到最小化能量函数模型为：