… | |
… | |
68 | |
68 | |
69 | *h &= 255; |
69 | *h &= 255; |
70 | } |
70 | } |
71 | } |
71 | } |
72 | |
72 | |
|
|
73 | struct feature { |
|
|
74 | float v1, v2, v3; // mean, square, cube |
|
|
75 | int n; |
|
|
76 | }; |
|
|
77 | |
|
|
78 | static void |
|
|
79 | feature_init (struct feature *f) |
|
|
80 | { |
|
|
81 | f->v1 = 0.; |
|
|
82 | f->v2 = 0.; |
|
|
83 | f->v3 = 0.; |
|
|
84 | f->n = 0; |
|
|
85 | } |
|
|
86 | |
|
|
87 | // didn't find an algorithm to neatly do mean, variance and skew in one pass. |
|
|
88 | // elmex ist schuld. |
|
|
89 | static void |
|
|
90 | feature_update_pass_1 (struct feature *f, unsigned int v) |
|
|
91 | { |
|
|
92 | f->v1 += v; |
|
|
93 | f->n += 1; |
|
|
94 | } |
|
|
95 | |
|
|
96 | static void |
|
|
97 | feature_finish_pass_1 (struct feature *f) |
|
|
98 | { |
|
|
99 | if (f->n < 1) |
|
|
100 | return; |
|
|
101 | |
|
|
102 | f->v1 /= f->n; |
|
|
103 | } |
|
|
104 | |
|
|
105 | static void |
|
|
106 | feature_update_pass_2 (struct feature *f, unsigned int v) |
|
|
107 | { |
|
|
108 | float d = v - f->v1; |
|
|
109 | |
|
|
110 | f->v2 += d * d; |
|
|
111 | f->v3 += d * d * d; |
|
|
112 | } |
|
|
113 | |
|
|
114 | static void |
|
|
115 | feature_finish_pass_2 (struct feature *f) |
|
|
116 | { |
|
|
117 | if (f->n < 1) |
|
|
118 | return; |
|
|
119 | |
|
|
120 | f->v2 /= f->n; |
|
|
121 | f->v3 /= f->n; |
|
|
122 | |
|
|
123 | f->v1 /= 255.; |
|
|
124 | f->v2 /= 255. * 255.; f->v2 = sqrtf (f->v2); |
|
|
125 | f->v3 /= 255. * 255. * 255.; f->v3 = powf (fabsf (f->v3), 1./3.); |
|
|
126 | } |
|
|
127 | |
73 | static guint32 a85_val; |
128 | static guint32 a85_val; |
74 | static guint a85_cnt; |
129 | static guint a85_cnt; |
75 | static guchar a85_buf[LINELENGTH], *a85_ptr; |
130 | static guchar a85_buf[LINELENGTH], *a85_ptr; |
76 | |
131 | |
77 | static void |
132 | static void |
… | |
… | |
586 | ############################################################################# |
641 | ############################################################################# |
587 | |
642 | |
588 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster |
643 | MODULE = Gtk2::CV PACKAGE = Gtk2::CV::Plugin::RCluster |
589 | |
644 | |
590 | SV * |
645 | SV * |
|
|
646 | extract_features (SV *ar) |
|
|
647 | CODE: |
|
|
648 | { |
|
|
649 | int i; |
|
|
650 | AV *av, *result; |
|
|
651 | |
|
|
652 | if (!SvROK (ar) || SvTYPE (SvRV (ar)) != SVt_PVAV) |
|
|
653 | croak ("Not an array ref as first argument to make_histogram"); |
|
|
654 | |
|
|
655 | av = (AV *) SvRV (ar); |
|
|
656 | result = newAV (); |
|
|
657 | |
|
|
658 | for (i = 0; i <= av_len (av); ++i) |
|
|
659 | { |
|
|
660 | SV *sv = *av_fetch (av, i, 1); |
|
|
661 | SV *histsv = newSV (9 * sizeof (float) + 1); |
|
|
662 | |
|
|
663 | SvPOK_on (histsv); |
|
|
664 | SvCUR_set (histsv, 9 * sizeof (float)); |
|
|
665 | float *hist = (float *)SvPVX (histsv); |
|
|
666 | |
|
|
667 | struct feature f_h, f_s, f_v; |
|
|
668 | feature_init (&f_h); |
|
|
669 | feature_init (&f_s); |
|
|
670 | feature_init (&f_v); |
|
|
671 | |
|
|
672 | { |
|
|
673 | STRLEN len; |
|
|
674 | unsigned char *buf = (unsigned char *)SvPVbyte (sv, len); |
|
|
675 | while (len >= 3) |
|
|
676 | { |
|
|
677 | unsigned int r, g, b, h, s, v; |
|
|
678 | r = *buf++; g = *buf++; b = *buf++; |
|
|
679 | rgb_to_hsv (r, g, b, &h, &s, &v); |
|
|
680 | |
|
|
681 | feature_update_pass_1 (&f_h, h); |
|
|
682 | feature_update_pass_1 (&f_s, s); |
|
|
683 | feature_update_pass_1 (&f_v, v); |
|
|
684 | |
|
|
685 | len -= 3; |
|
|
686 | } |
|
|
687 | |
|
|
688 | feature_finish_pass_1 (&f_h); |
|
|
689 | feature_finish_pass_1 (&f_s); |
|
|
690 | feature_finish_pass_1 (&f_v); |
|
|
691 | } |
|
|
692 | |
|
|
693 | { |
|
|
694 | STRLEN len; |
|
|
695 | unsigned char *buf = (unsigned char *)SvPVbyte (sv, len); |
|
|
696 | while (len >= 3) |
|
|
697 | { |
|
|
698 | unsigned int r, g, b, h, s, v; |
|
|
699 | r = *buf++; g = *buf++; b = *buf++; |
|
|
700 | rgb_to_hsv (r, g, b, &h, &s, &v); |
|
|
701 | |
|
|
702 | feature_update_pass_2 (&f_h, h); |
|
|
703 | feature_update_pass_2 (&f_s, s); |
|
|
704 | feature_update_pass_2 (&f_v, v); |
|
|
705 | |
|
|
706 | len -= 3; |
|
|
707 | } |
|
|
708 | |
|
|
709 | feature_finish_pass_2 (&f_h); |
|
|
710 | feature_finish_pass_2 (&f_s); |
|
|
711 | feature_finish_pass_2 (&f_v); |
|
|
712 | } |
|
|
713 | |
|
|
714 | hist [0] = f_h.v1 * 2.; hist [1] = f_h.v2 * 2.; hist [2] = f_h.v3 * 2.; |
|
|
715 | hist [3] = f_s.v1 ; hist [4] = f_s.v2 ; hist [5] = f_s.v3 ; |
|
|
716 | hist [6] = f_v.v1 ; hist [7] = f_v.v2 ; hist [8] = f_v.v3 ; |
|
|
717 | |
|
|
718 | av_push (result, histsv); |
|
|
719 | } |
|
|
720 | |
|
|
721 | RETVAL = newRV_noinc ((SV *)result); |
|
|
722 | } |
|
|
723 | OUTPUT: |
|
|
724 | RETVAL |
|
|
725 | |
|
|
726 | SV * |
591 | make_histograms (SV *ar) |
727 | make_histograms (SV *ar) |
592 | CODE: |
728 | CODE: |
593 | { |
729 | { |
594 | int i; |
730 | int i; |
595 | AV *av, *result; |
731 | AV *av, *result; |