matlab实现,数据曲线毛刺光滑

实现效果

c 复制代码
[Curve] = Model_optimization2(1:size(M,1),M1);

figure(1);
hold on; 
plot(M1,'k-');
plot(Curve,'g-.');
legend("处理前","处理后");
cpp 复制代码
function [Curve] = Model_optimization(HU_curve,Curve)
%循环遍历相邻数据点,直到所有相邻数据点之间的的差值的绝对值都小于指定的阈值,遍历结束。从第一个点开始遍历,
% 判断当前点A和下一个点B的差值的绝对值,如果大于指定阈值,则继续寻找下下个点C,如果点A和点C差值满足条件,则将点B位置的值替
% 换为 (A+B)/2,此时开始点跳到C点,继续遍历,如果C点不满足条件,也即是存在大于等于2个点是存在较大的偏差的,继续遍历,
% 直到找到满足条件的C点。(允许寻找周围的最多7个点,否则终止寻找,继续下一个数据点的判定)

%%%%%%输入 

% % % tempCurve_b = BONE_curve;
% % % tempCurve = WATER_curve;

Th = 0.005;
index = 1;
end_label = 0;

while (index <= length(HU_curve)-2)
    index
    if(abs(Curve(index) - Curve(index+1))> Th && end_label ==0)
        %%%%%%%%% 判断第3个点是否满足条件
        if(abs(Curve(index) - Curve(index+2))< Th)
            Curve(index+1) = (Curve(index) + Curve(index+2))/2;
            index = index +2 ;
% %             disp("in");
        else
            %%%%如果第3个点不满足,开始计数,5个点内还找不到则放弃,如果5个点内找到满足条件的点,
            %%% 进一步判断,第一个点和当前点的大小关系,判断是单调增还是单调减,如果是单调增,且找了3个点才找到符合条件的点,就按照每次递增1/3的第一个点和当前
            %%% 符合条件的点的绝对值之差
            counter =2;
%             for kkk =3:7
            for kkk =3:40
                if (index+kkk) > length(HU_curve)
                    end_label = 1;
                    break;
                end
                
                if(abs(Curve(index) - Curve(index+kkk))< Th)    
                    break;
                end
                counter = counter +1;
            end
            
            if end_label == 0
                %%%%%%%% 改变对应位置的数据点
                del = abs(Curve(index) - Curve(index+kkk))/counter;
                for ii =1:counter
                    %%% 判断递增还是递减
                    if (index+kkk) <= length(HU_curve)
                        if(Curve(index) > Curve(index+kkk))  %%%递减
                            Curve(index+ii) = Curve(index) - del*ii;
                        else  %%%递减
                            Curve(index+ii) = Curve(index) + del*ii;
                        end
                    end
                end
                index = index + kkk;
            end
        end
    elseif(abs(Curve(index) - Curve(index+1))<= Th && end_label ==0)
         index = index +1;
    else
        break;          
    end
end

end
相关推荐
The Future is mine11 分钟前
Python计算经纬度两点之间距离
开发语言·python
Enti7c12 分钟前
HTML5和CSS3的一些特性
开发语言·css3
斯汤雷13 分钟前
Matlab绘图案例,设置图片大小,坐标轴比例为黄金比
数据库·人工智能·算法·matlab·信息可视化
爱吃巧克力的程序媛19 分钟前
在 Qt 创建项目时,Qt Quick Application (Compat) 和 Qt Quick Application
开发语言·qt
独好紫罗兰1 小时前
洛谷题单3-P5719 【深基4.例3】分类平均-python-流程图重构
开发语言·python·算法
篝火悟者1 小时前
自学-C语言-基础-数组、函数、指针、结构体和共同体、文件
c语言·开发语言
genispan2 小时前
QT/C++ 多线程并发下载实践
开发语言·c++·qt
-代号95272 小时前
【JavaScript】十三、事件监听与事件类型
开发语言·javascript·ecmascript
写代码的小王吧2 小时前
【Java可执行命令】(十)JAR文件签名工具 jarsigner:通过数字签名及验证保证代码信任与安全,深入解析 Java的 jarsigner命令~
java·开发语言·网络·安全·web安全·网络安全·jar
小卡皮巴拉2 小时前
【力扣刷题实战】矩阵区域和
开发语言·c++·算法·leetcode·前缀和·矩阵