function [X,Y,Z]=func_restruction(I1,I2,I3,is_real);
II1 = func_getbk(I1,is_real);
II2 = func_getbk(I2,is_real);
II3 = func_getbk(I3,is_real);
[Sa1,Sb1,H1,W1] = func_featrure(II1);
[Sa2,Sb2,H2,W2] = func_featrure(II2);
[Sa3,Sb3,H3,W3] = func_featrure(II3);
if H3/W3 > 0.9 & abs(Sb3/Sa3 - 0.5) <0.05
t = H3/W3*2 :-0.1: 0;
[x,y,z] = cylinder(t);
X = W3*x/4;
Y = W3*y/4;
Z = H3*(z);
end
if H3/W3 < 1 & H3/W3 > 0.85 & abs(Sb3/Sa3 - 0.75) <0.05
RR = round((H3+W3)/4);
[x,y,z] = sphere(40);
X = RR*x;
Y = RR*y;
Z = RR*z;
end
if H3/W3 < 1 & H3/W3 > 0.85 & abs(Sb3/Sa3 - 0.6) <0.07
RR = round((H3+W3)/4);
[x,y,z]=yuanpan(RR,RR/2);
X = RR*x;
Y = RR*y;
Z = RR*z;
end
if H3/W3 < 0.6 & H3/W3 > 0.5 & abs(Sb3/Sa3 - 0.8) <0.05
XR = W3;
YR = W3/2;
ZR = H3/2;
[x,y,z]=ellipsoid(0,0,0,XR,YR,ZR,40);
X = x;
Y = y;
Z = z;
end