MIC-Modeling-Polymer-Composites/Statisics_on_samples_normalized.m at gh-pages · Materials-Informatics-Class-Fall2014/MIC-Modeling-Polymer-Composites · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101

normalize = @(A)( A-min(A(:)) ) ./ ( max(A(:)) - min(A(:)) )
%% Spacial Stats on Original sample 2

data.('name') = 'Spatial_Stats_Original2';
data.('local') = '_data/Spatial_Stats_Original2.mat';
data.('header') = '_data/Spatial_Stats_Original2.mat.json';

data.('source') = '_data/Ncropped.mat';

load( data.source);
cropped = double( Ncropped );
org2.phase_norm = normalize(cropped);
org2.stats_norm_xy150 = SpatialStatsFFT( org2.phase_norm(:,:,150));
saveas(gcf, './assets/Spacial_Stats_Original_sample2_xy_normalized.png' );
data.('images').('stats')(4) = struct( 'src', '/assets/Spacial_Stats_Original_sample2_xy_normalized.png', 'description', 'Spatial statistics on original image N in X-Y direction.' );
org2.stats_norm_xz150 = SpatialStatsFFT( squeeze(org2.phase_norm(:,150,:)));
saveas(gcf, './assets/Spacial_Stats_Original_sample2_xz_normalized.png' );
data.('images').('stats')(5) = struct( 'src', '/assets/Spacial_Stats_Original_sample2_xz_normalized.png', 'description', 'Spatial statistics on original image N in X-Z direction.' );
org2.stats_norm_yz150 = SpatialStatsFFT( squeeze(org2.phase_norm(150,:,:)));
saveas(gcf, './assets/Spacial_Stats_Original_sample2_yz_normalized.png' );
data.('images').('stats')(6) = struct( 'src', '/assets/Spacial_Stats_Original_sample2_yz_normalized.png', 'description', 'Spatial statistics on original image N in Y-Z direction.' );


s = savejson( [], data );
fo = fopen( data.header, 'w'); fwrite( fo, s ); fclose(fo);

%% Spacial Stats on Original sample 1

data.('name') = 'Spatial_Stats_Original';
data.('local') = '_data/Spatial_Stats_Original.mat';
data.('header') = '_data/Spatial_Stats_Original.mat.json';

data.('source') = '_data/Cropped_fiber_matrix_400_400.mat';

load( data.source);
cropped = double( Scropped );
org.phase_norm = normalize(cropped);
org.stats_norm_xy150 = SpatialStatsFFT( org.phase_norm(:,:,150));
saveas(gcf, './assets/Spacial_Stats_Original_sample1_xy_normalized.png' );
data.('images').('stats')(4) = struct( 'src', '/assets/Spacial_Stats_Original_sample1_xy_normalized.png', 'description', 'Spatial statistics on original image S in X-Y direction.' );
org.stats_norm_xz150 = SpatialStatsFFT( squeeze(org.phase_norm(:,150,:)));
saveas(gcf, './assets/Spacial_Stats_Original_sample1_xz_normalized.png' );
data.('images').('stats')(5) = struct( 'src', '/assets/Spacial_Stats_Original_sample1_xz_normalized.png', 'description', 'Spatial statistics on original image in X-Z direction.' );
org.stats_norm_yz150 = SpatialStatsFFT( squeeze(org.phase_norm(150,:,:)));
saveas(gcf, './assets/Spacial_Stats_Original_sample1_yz_normalized.png' );
data.('images').('stats')(6) = struct( 'src', '/assets/Spacial_Stats_Original_sample1_yz_normalized.png', 'description', 'Spatial statistics on original image in Y-Z direction.' );

s = savejson( [], data );
fo = fopen( data.header, 'w'); fwrite( fo, s ); fclose(fo);

%% Spacial Stats on segmented sample 2

data.('name') = 'Spatial_Stats_Segmented2';
data.('local') = '_data/Spatial_Stats_Segmented2.mat';
data.('header') = '_data/Spatial_Stats_Segmented2.mat.json';

data.('source') = '_data/SegSample2P2.mat';

load( data.source);
cropped = double( I);
seg2.phase_norm = normalize(cropped);
seg2.stats_norm_xy150 = SpatialStatsFFT( seg2.phase_norm(:,:,150));
saveas(gcf, './assets/Spacial_Stats_Segmented_sample2_xy_normalized.png' );
data.('images').('stats')(4) = struct( 'src', '/assets/Spacial_Stats_Segmented_sample2_xy_normalized.png', 'description', 'Spatial statistics on segmented image N in X-Y direction.' );
seg2.stats_norm_xz150 = SpatialStatsFFT( squeeze(seg2.phase_norm(:,150,:)));
saveas(gcf, './assets/Spacial_Stats_Segmented_sample2_xz_normalized.png' );
data.('images').('stats')(5) = struct( 'src', '/assets/Spacial_Stats_Segmented_sample2_xz_normalized.png', 'description', 'Spatial statistics on segmented image N in X-Z direction.' );
seg2.stats_norm_yz150 = SpatialStatsFFT( squeeze(seg2.phase_norm(150,:,:)));
saveas(gcf, './assets/Spacial_Stats_Segmented_sample2_yz_normalized.png' );
data.('images').('stats')(6) = struct( 'src', '/assets/Spacial_Stats_Segmented_sample2_yz_normalized.png', 'description', 'Spatial statistics on segmented image N in Y-Z direction.' );


s = savejson( [], data );
fo = fopen( data.header, 'w'); fwrite( fo, s ); fclose(fo);

%% Spacial Stats on segmented sample 1

data.('name') = 'Spatial_Stats_Segmented';
data.('local') = '_data/Spatial_Stats_Segmented.mat';
data.('header') = '_data/Spatial_Stats_Segmented.mat.json';

data.('source') = '_data/SegSample1P2.mat';

load( data.source);
cropped = double( I);
seg.phase_norm = normalize(cropped);
seg.stats_norm_xy150 = SpatialStatsFFT( seg.phase_norm(:,:,150));
saveas(gcf, './assets/Spacial_Stats_Segmented_sample1_xy_normalized.png' );
data.('images').('stats')(4) = struct( 'src', '/assets/Spacial_Stats_Segmented_sample1_xy_normalized.png', 'description', 'Spatial statistics on segmented image S in X-Y direction.' );
seg.stats_norm_xz150 = SpatialStatsFFT( squeeze(seg.phase_norm(:,150,:)));
saveas(gcf, './assets/Spacial_Stats_Segmented_sample1_xz_normalized.png' );
data.('images').('stats')(5) = struct( 'src', '/assets/Spacial_Stats_Segmented_sample1_xz_normalized.png', 'description', 'Spatial statistics on segmented image S in X-Z direction.' );
seg.stats_norm_yz150 = SpatialStatsFFT( squeeze(seg.phase_norm(150,:,:)));
saveas(gcf, './assets/Spacial_Stats_Segmented_sample1_yz_normalized.png' );
data.('images').('stats')(6) = struct( 'src', '/assets/Spacial_Stats_Segmented_sample1_yz_normalized.png', 'description', 'Spatial statistics on segmented image S in Y-Z direction.' );


s = savejson( [], data );
fo = fopen( data.header, 'w'); fwrite( fo, s ); fclose(fo);