使用matlab对一组数据的处理以及研究思路路线（9）--output invision

%% 在最小容忍 
[P_X_1,P_Y_1,P_Z_1] = Get_Range_Max(1.5,Delta_k_1,precision,tolerance_Minimal);
[P_X_2,P_Y_2,P_Z_2] = Get_Range_Max(3  ,Delta_k_2,precision,tolerance_Minimal);
[P_X_3,P_Y_3,P_Z_3] = Get_Range_Max(5  ,Delta_k_3,precision,tolerance_Minimal);
[P_X_4,P_Y_4,P_Z_4] = Get_Range_Max(10 ,Delta_k_4,precision,tolerance_Minimal);
[P_X_5,P_Y_5,P_Z_5] = Get_Range_Max(15 ,Delta_k_5,precision,tolerance_Minimal);
%% 得到对应的 位置 
[Delta_select_1] = Get_Range_Max_Delta_select(tolerance_Minimal,Delta_k_1);
[Delta_select_2] = Get_Range_Max_Delta_select(tolerance_Minimal,Delta_k_2);
[Delta_select_3] = Get_Range_Max_Delta_select(tolerance_Minimal,Delta_k_3);
[Delta_select_4] = Get_Range_Max_Delta_select(tolerance_Minimal,Delta_k_4);
[Delta_select_5] = Get_Range_Max_Delta_select(tolerance_Minimal,Delta_k_5);
P_X_all=[P_X_1;P_X_2;P_X_3;P_X_4;P_X_5];
P_Y_all=[P_Y_1;P_Y_2;P_Y_3;P_Y_4;P_Y_5];
P_Z_all=[P_Z_1;P_Z_2;P_Z_3;P_Z_4;P_Z_5];
Delta_select_All=[Delta_select_1;Delta_select_2;Delta_select_3;Delta_select_4;Delta_select_5];
ALL=[P_X_all,P_Y_all,P_Z_all,Delta_select_All];
%% 容忍 插值点云图
figure(N);
N=N+1;
subplot(1,2,1);
plot3(P_X_all,P_Y_all,P_Z_all,'.','markersize',4);
xlabel('x/b'), ylabel('y/b'), zlabel('z/b');
grid on;
%%
%***************************************************************************************************************************************************
%   ************************************************************************************************************
%***************************************************************************************************************************************************
subplot(1,2,2);
F = scatteredInterpolant(ALL(:,1:3),ALL(:,4)); 
[xq,yq,zq] = meshgrid(min(ALL(:,1)):50*precision:max(ALL(:,1))...
                     ,min(ALL(:,2)):precision: max(ALL(:,2))  ...
                     ,min(ALL(:,3)):precision: max(ALL(:,3)));
F_Delta_select_All=F(xq,yq,zq);  %%%%%%%%%%%  由插值的集合 映射关系
xslice = (min(ALL(:,1)):50*precision:max(ALL(:,1)));
yslice=(min(ALL(:,2)):precision: max(ALL(:,2)));
zslice=(min(ALL(:,3)):precision: max(ALL(:,3)));
slice(xq,yq,zq,F_Delta_select_All,xslice,yslice,zslice)
xlabel('x/b'), ylabel('y/b'), zlabel('z/b')
%%
% 没意义
% subplot(1,2,2);
% [yi,zi] = meshgrid(-0.375:precision:0.375, -0.375:precision:0.375);
% xi = griddata(P_Y_all,P_Z_all,P_X_all,yi,zi,'cubic');
% surf(xi,yi,zi);
% xlabel('x/b'), ylabel('y/b'), zlabel('z/b')
% clear xi yi zi;