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
相关推荐
264玫瑰资源库16 分钟前
问道数码兽 怀旧剧情回合手游源码搭建教程(反查重优化版)
java·开发语言·前端·游戏
普if加的帕41 分钟前
java Springboot使用扣子Coze实现实时音频对话智能客服
java·开发语言·人工智能·spring boot·实时音视频·智能客服
安冬的码畜日常1 小时前
【AI 加持下的 Python 编程实战 2_10】DIY 拓展:从扫雷小游戏开发再探问题分解与 AI 代码调试能力(中)
开发语言·前端·人工智能·ai·扫雷游戏·ai辅助编程·辅助编程
朝阳5812 小时前
Rust项目GPG签名配置指南
开发语言·后端·rust
朝阳5812 小时前
Rust实现高性能目录扫描工具ll的技术解析
开发语言·后端·rust
程高兴2 小时前
基于Matlab的车牌识别系统
开发语言·matlab
牛马baby2 小时前
Java高频面试之并发编程-07
java·开发语言·面试
CodeWithMe2 小时前
【C++】STL之deque
开发语言·c++
炯哈哈3 小时前
【上位机——MFC】运行时类信息机制
开发语言·c++·mfc·上位机
T.Ree.5 小时前
【数据结构】_树和二叉树
c语言·开发语言·数据结构