Effect of rotation in the fits to symmetric ellipsoids

PENTACAM DATA

I generate datos_fig_pap_pentacam_rot by running fig_pap_correlpentacam_03, and changing in the code c_cuadros=5:5.

>> load datos_fig_pap_pentacam_rot
>> R_rot=R; Q_rot=Q;
>> load datos_fig_pap_pentacam

First subject:
>> close all
subplot(2,2,1)
hist(R{1,1})
title(‘Cara anterior, sin rotación’)
subplot(2,2,2)
hist(R_rot{1,1})
title(‘Cara anterior, con rotación’)
subplot(2,2,3)
hist(R{1,2})
title(‘Cara posterior, sin rotación’)
subplot(2,2,4)
hist(R_rot{1,2})
title(‘Cara posterior, con rotación’)

Second subject:

>> close all
subplot(2,2,1)
hist(R{2,1})
title(‘Cara anterior, sin rotación’)
subplot(2,2,2)
hist(R_rot{2,1})
title(‘Cara anterior, con rotación’)
subplot(2,2,3)
hist(R{2,2})
title(‘Cara posterior, sin rotación’)
subplot(2,2,4)
hist(R_rot{2,2})
title(‘Cara posterior, con rotación’)

Third subject:

close all
subplot(2,2,1)
hist(R{3,1})
title(‘Cara anterior, sin rotación’)
subplot(2,2,2)
hist(R_rot{3,1})
title(‘Cara anterior, con rotación’)
subplot(2,2,3)
hist(R{3,2})
title(‘Cara posterior, sin rotación’)
subplot(2,2,4)
hist(R_rot{3,2})
title(‘Cara posterior, con rotación’)

Asphericities, for the three subjects:

>> close all
subplot(2,2,1)
hist(Q{1,1})
title(‘Cara anterior, sin rotación’)
subplot(2,2,2)
hist(Q_rot{1,1})
title(‘Cara anterior, con rotación’)
subplot(2,2,3)
hist(Q{1,2})
title(‘Cara posterior, sin rotación’)
subplot(2,2,4)
hist(Q_rot{1,2})
title(‘Cara posterior, con rotación’)
figure
subplot(2,2,1)
hist(Q{2,1})
title(‘Cara anterior, sin rotación’)
subplot(2,2,2)
hist(Q_rot{2,1})
title(‘Cara anterior, con rotación’)
subplot(2,2,3)
hist(Q{2,2})
title(‘Cara posterior, sin rotación’)
subplot(2,2,4)
hist(Q_rot{2,2})
title(‘Cara posterior, con rotación’)
figure
subplot(2,2,1)
hist(Q{3,1})
title(‘Cara anterior, sin rotación’)
subplot(2,2,2)
hist(Q_rot{3,1})
title(‘Cara anterior, con rotación’)
subplot(2,2,3)
hist(Q{3,2})
title(‘Cara posterior, sin rotación’)
subplot(2,2,4)
hist(Q_rot{3,2})
title(‘Cara posterior, con rotación’)


ATLAS DATA

Only horizontal meridian.

First subject:

subplot(2,2,1)
hist(datos.quadric(:,1))
title(‘Radio, sin rotación’)
subplot(2,2,2)
hist(datos.quadric_rot(:,1))
title(‘Radio, con rotación’)
subplot(2,2,3)
hist(datos.quadric(:,3))
title(‘Asfericidad, sin rotación’)
subplot(2,2,4)
hist(datos.quadric_rot(:,3))
title(‘Asfericidad, con rotación’)

Second subject:

>> close all
>> subplot(2,2,1)
hist(datos.quadric_laura(:,1))
title(‘Radio, sin rotación’)
subplot(2,2,2)
hist(datos.quadric_laura_rot(:,1))
title(‘Radio, con rotación’)
subplot(2,2,3)
hist(datos.quadric_laura(:,3))
title(‘Asfericidad, sin rotación’)
subplot(2,2,4)
hist(datos.quadric_laura_rot(:,3))
title(‘Asfericidad, con rotación’)

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s

%d bloggers like this: