| [45] | 1 | % make rotational invariant features |
|---|
| 2 | % |
|---|
| 3 | % The whole procedure goes like: |
|---|
| 4 | % 1. readcells.m -> cells.mat |
|---|
| 5 | % 1a. scalecellimage -> cellscaled.mat |
|---|
| 6 | % 2. gencellgabor.m -> cells_gabor_sc.mat |
|---|
| 7 | % 3. makerotinv |
|---|
| 8 | % 3a. gencellsimple.m (max) -> cellssimplemax.mat |
|---|
| 9 | % 3b. gencellsimple.m (mean) -> cellssimplemean.mat |
|---|
| 10 | % 3c. gencellsimple.m (var) -> cellssimplevar.mat |
|---|
| 11 | |
|---|
| 12 | [a,b]=gencellsimple(x,z,'max'); |
|---|
| 13 | xnew = sumgabororientations(a); |
|---|
| 14 | znew = sumgabororientations(b); |
|---|
| 15 | [a,b]=gencellsimple(x,z,'mean'); |
|---|
| 16 | xnew = [xnew sumgabororientations(a)]; |
|---|
| 17 | znew = [znew sumgabororientations(b)]; |
|---|
| 18 | [a,b]=gencellsimple(x,z,'var'); |
|---|
| 19 | xnew = [xnew sumgabororientations(a)]; |
|---|
| 20 | znew = [znew sumgabororientations(b)]; |
|---|
| 21 | |
|---|
| 22 | % fix a wrong labeling in the training set: |
|---|
| 23 | lab = getlabels(xnew); |
|---|
| 24 | I = strmatch('centromeric',lab); |
|---|
| 25 | lab(I,:) = repmat('centromere ',length(I),1); |
|---|
| 26 | xnew = remclass(setlabels(xnew,lab)); |
|---|
| 27 | |
|---|
| 28 | a = xnew; b = znew; |
|---|
| 29 | save cellssimple_rotinv a b; |
|---|
| 30 | |
|---|
| 31 | |
|---|
| 32 | |
|---|
