Форум программистов, компьютерный форум, киберфорум
Java SE (J2SE)
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 4.80/5: Рейтинг темы: голосов - 5, средняя оценка - 4.80
0 / 0 / 0
Регистрация: 02.02.2017
Сообщений: 18

Детектор Канни

10.12.2017, 15:17. Показов 923. Ответов 4
Метки нет (Все метки)

Студворк — интернет-сервис помощи студентам
Доброго времени суток, дорогие друзья; буду очень вам признателен, если поможете. В программировании я новичок и классы еще толком не изучал. Мне необходимо протестировать готовые алгоритмы детектора Канни. Нашел этот вариант на java:

Кликните здесь для просмотра всего текста
Java
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
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
package image;
 
import java.awt.Point;
import java.awt.image.BufferedImage;
import java.awt.image.Raster;
import java.util.Stack;
 
 
public class ColorCannyEdgeDetector {
 
    private final static double GAUSSIAN_CUT_OFF = 0.005;
    private int W;
    private int H;
    private int[] r, g, b;
    private double yrConv[];
    private double xrConv[];
    private double ygConv[];
    private double xgConv[];
    private double ybConv[];
    private double xbConv[];
 
 
    private double xGradient[];
    private double yGradient[];
    private double vMagnitude[];
 
    private int[] data;
 
 
    private double[] magnitude;
 
    double gaussianKernelRadius = 1;
    int gaussianKernelWidth = 16;
    private static final double lowThreshold = 3;
    private static final double highThreshold = 7;
 
    public ColorCannyEdgeDetector(BufferedImage image){
 
        this.W = image.getWidth();
        this.H = image.getHeight();
        int p = 0;
        r = new int[W*H];
        g = new int[W*H];
        b = new int[W*H];
 
        for(int x=0; x<W; x++){
            for(int y=0; y<H; y++){
                int rgb = image.getRGB(x, y);
                r[p] = (rgb >> 16)&0xff;
                g[p] = (rgb >> 8)&0xff;
                b[p++] = rgb&0xff;
            }
        }
 
        yrConv = new double[r.length];
        xrConv = new double[r.length];
        ygConv = new double[r.length];
        xgConv = new double[r.length];
        ybConv = new double[r.length];
        xbConv = new double[r.length];
 
        xGradient = new double[r.length];
        yGradient = new double[r.length];
        vMagnitude = new double[r.length];
 
        data = new int[r.length];
        magnitude = new double[r.length];
 
    }
 
 
 
 
    public BufferedImage findEdges() {
        computeColorGradients(gaussianKernelRadius, gaussianKernelWidth);
        double low = lowThreshold;
        double high = highThreshold;
        hysterezis(low, high);
        threshold();
//        double max = 0;
//        for(int i=0; i<magnitude.length; i++){
//            max = Math.max(max, magnitude[i]);
//        }
//        for(int i=0; i<data.length; i++){
//            data[i] = (int) (255 * magnitude[i] / max);
//        }
        return image(data, data, data, W, H);
    }
 
    private static BufferedImage image(int[] r, int[] g, int[] b, int W, int H) {
        BufferedImage image = new BufferedImage(W, H, BufferedImage.TYPE_INT_RGB);
        for(int y=0; y<H; y++){
            for(int x=0; x<W; x++){
                image.setRGB(x, y, r[x*H+y]<<16 | g[x*H+y]<<8 | b[x*H+y]);
            }
        }
        return image;
    }
 
    private void computeColorGradients(double kernelRadius, int kernelWidth) {
        double kernel[] = new double[kernelWidth];
        double diffKernel[] = new double[kernelWidth];
        int kWidth;
        for (kWidth = 0; kWidth < kernelWidth; kWidth++) {
            double g1 = gaussian(kWidth, kernelRadius);
            if (g1 <= GAUSSIAN_CUT_OFF && kWidth >= 2) break;
            double g2 = gaussian(kWidth - 0.5, kernelRadius);
            double g3 = gaussian(kWidth + 0.5, kernelRadius);
            kernel[kWidth] = (g1 + g2 + g3) / 3. / (2 * Math.PI * kernelRadius * kernelRadius);
            diffKernel[kWidth] = g3 - g2;
        }
 
        int initX = kWidth - 1;
        int maxX = W - (kWidth - 1);
        int initY = (kWidth - 1);
        int maxY = (H - (kWidth - 1));
 
        for (int x = initX; x < maxX; x++) {
            for (int y = initY; y < maxY; y++) {
                double sumrX = r[x * H + y] * kernel[0];
                double sumrY = r[x * H + y] * kernel[0];
                double sumgX = g[x * H + y] * kernel[0];
                double sumgY = g[x * H + y] * kernel[0];
                double sumbX = b[x * H + y] * kernel[0];
                double sumbY = b[x * H + y] * kernel[0];
 
                for (int ri = 1; ri < kWidth; ri++) {
                    sumrY += kernel[ri] * (r[x * H + y - ri] + r[x * H + y + ri]);
                    sumrX += kernel[ri] * (r[(x - ri) * H + y] + r[(x + ri) * H + y]);
                    sumgY += kernel[ri] * (g[x * H + y - ri] + g[x * H + y + ri]);
                    sumgX += kernel[ri] * (g[(x - ri) * H + y] + g[(x + ri) * H + y]);
                    sumbY += kernel[ri] * (b[x * H + y - ri] + b[x * H + y + ri]);
                    sumbX += kernel[ri] * (b[(x - ri) * H + y] + b[(x + ri) * H + y]);
 
                }
 
                yrConv[x * H + y] = sumrY;
                xrConv[x * H + y] = sumrX;
                ygConv[x * H + y] = sumgY;
                xgConv[x * H + y] = sumgX;
                ybConv[x * H + y] = sumbY;
                xbConv[x * H + y] = sumbX;
            }
 
        }
 
 
        initX = kWidth;
        maxX = W - kWidth;
        initY = kWidth;
        maxY = (H - kWidth);
        for (int x = initX; x < maxX; x++) {
            for (int y = initY; y < maxY; y++) {
                double rx = (xrConv[(x + 1) * H + y] - xrConv[(x - 1) * H + y]) / 2.;
                double gx = (xgConv[(x + 1) * H + y] - xgConv[(x - 1) * H + y]) / 2.;
                double bx = (xbConv[(x + 1) * H + y] - xbConv[(x - 1) * H + y]) / 2.;
 
                double ry = (yrConv[x * H + y - 1] - yrConv[x * H + y + 1]) / 2.;
                double gy = (ygConv[x * H + y - 1] - ygConv[x * H + y + 1]) / 2.;
                double by = (ybConv[x * H + y - 1] - ybConv[x * H + y + 1]) / 2.;
 
 
                double q1 = rx * rx + gx * gx + bx * bx;
                double q2 = rx * ry + gx * gy + bx * by;
                double q4 = ry * ry + gy * gy + by * by;
 
                double[] eigen = eigen(q1, q2, q2, q4);
                double fe1 = eigen[0];
                double fe2 = eigen[1];
                double fv = eigen[2];
 
                xGradient[x * H + y] = fe1;
                yGradient[x * H + y] = fe2;
                vMagnitude[x * H + y] = Math.sqrt(fv);
            }
 
        }
 
 
        for (int x = initX+1; x < maxX-1; x++) {
            for (int y = initY+1; y < maxY-1; y++) {
 
                double a1 = vMagnitude[(x - 1) * H + y - 1];
                double a2 = vMagnitude[x * H + y - 1];
                double a3 = vMagnitude[(x + 1) * H + y - 1];
 
                double a4 = vMagnitude[(x - 1) * H + y];
                double a5 = vMagnitude[x * H + y];
                double a6 = vMagnitude[(x + 1) * H + y];
 
                double a7 = vMagnitude[(x - 1) * H + y + 1];
                double a8 = vMagnitude[x * H + y + 1];
                double a9 = vMagnitude[(x + 1) * H + y + 1];
 
                double gx = xGradient[x * H + y];
                double gy = yGradient[x * H + y];
 
                double angle = Math.atan2(gy, gx);
                double h, v;
 
                if (ishorizonatal(angle)) {
                    h = a4;
                    v = a6;
                } else if (isvertical(angle)) {
                    h = a2;
                    v = a8;
                } else if (isdiag1(angle)) {
                    h = a3;
                    v = a7;
                } else if (isdiag2(angle)) {
                    h = a1;
                    v = a9;
                } else {
                   throw new RuntimeException("Incorrect angle");
                }
 
                double val = 0;
                if (a5 > h && a5 >= v)
                    val = a5;
 
                magnitude[x * H + y] = val;
            }
        }
    }
 
    private static int conv(byte col) {
        return col >= 0 ? col : (int) col + 256;
    }
 
    private void hysterezis(double low, double high) {
        for (int y = 0; y < H; y++) {
            for (int x = 0; x < W; x++) {
                if (data[x * H + y] == 0 && magnitude[x * H + y] >= high) {
                    track(x, y, low);
                }
            }
        }
    }
 
    private void threshold() {
        for (int y = 0; y < H; y++) {
            for (int x = 0; x < W; x++) {
                data[x * H + y] = data[x * H + y] > 0 ? 255 : 0;
            }
        }
    }
 
    private void track(int x_, int y_, double threshold) {
        Stack<Point> stack = new Stack<Point>();
        stack.push(new Point(x_, y_));
        data[x_*H+y_] = (int) magnitude[x_*H+y_];
        while(!stack.empty()) {
            Point n = stack.pop();
            int x1 = n.x;
            int y1 = n.y;
            int x0 = x1 == 0 ? 0 : x1 - 1;
            int x2 = x1 == W-1 ? W-1 : x1 + 1;
            int y0 = y1 == 0 ? 0 : y1 - 1;
            int y2 = y1 == H-1 ? H-1 : y1 + 1;
 
            for (int x = x0; x <= x2; x++) {
                for (int y = y0; y <= y2; y++) {
                    if ((y != y1 || x != x1) && data[x * H + y] == 0 && magnitude[x * H + y] >= threshold) {
                        data[x*H+y] = (int) magnitude[x*H+y] + 1;
                        stack.push(new Point(x, y));
                    }
                }
            }
        }
    }
 
    private double[] eigen(double q1, double q2, double q3, double q4) {
 
        double a = 1;
        double b = -(q1 + q4);
        double c = q1 * q4 - q2 * q3;
        double v = ((-b + Math.sqrt(b * b - 4 * a * c)) / (2 * a));
        double e1 = -q2;
        double e2 = q1 - v;
 
        if (v < 0) {
            v = -v;
            e1 = -e1;
            e2 = -e2;
        }
 
        return new double[]{e1, e2, v};
    }
 
    private boolean isdiag2(double angle) {
        return ((angle >= -3 * Math.PI / 8) && (angle < -(Math.PI / 8)))
                || ((angle >= 5 * Math.PI / 8) && (angle < 7 * Math.PI / 8));
    }
 
    private boolean isdiag1(double angle) {
        return ((angle >= -7 * Math.PI / 8) && (angle < -5 * Math.PI / 8))
                || ((angle >= Math.PI / 8) && (angle < 3 * Math.PI / 8));
    }
 
    private boolean isvertical(double angle) {
        return ((angle >= -5 * Math.PI / 8) && (angle < -3 * Math.PI / 8))
                || ((angle >= 3 * Math.PI / 8) && (angle < 5 * Math.PI / 8));
    }
 
    private boolean ishorizonatal(double angle) {
        return (angle < -7 * Math.PI / 8) || ((angle >= -(Math.PI / 8))
                && (angle < Math.PI / 8)) || (angle >= 7 * Math.PI / 8);
    }
 
 
    private double gaussian(double x, double sigma) {
        return Math.exp(-(x * x) / (2f * sigma * sigma));
    }
}


Я так понимаю, здесь отсутствует обращение к классу? Или нет.. Кто может подсказать, как запустить этот код? Компилятор, пока, ругается на отсутствие каких-то файлов. Использую InteliJ IDEA
0
cpp_developer
Эксперт
20123 / 5690 / 1417
Регистрация: 09.04.2010
Сообщений: 22,546
Блог
10.12.2017, 15:17
Ответы с готовыми решениями:

Градиент в детекторе Канни
Собственно, судя по информации в интернете, угол направления в градиенте принимает значения от 0 до 180. Но, при вычислении этого...

Детектор напряжений
Здравствуйте. прошу подсказки. есть вот такая задачка: есть 3 входящие линии с напряжением 5 вольт и необходимо вывести одну исходящую но...

Пороговый детектор
Добрый день. Столкнулся с задачей детектирования сигнала на фоне шума. Решающей статистикой в моем случае является выход согласованного...

4
Эксперт Java
3639 / 2971 / 918
Регистрация: 05.07.2013
Сообщений: 14,220
10.12.2017, 23:16
напиши hello world, например.
0
 Аватар для Aviz__
2755 / 2062 / 509
Регистрация: 17.02.2014
Сообщений: 9,491
11.12.2017, 17:47
Arthuric, вы привели код, который является составной частью чего-то большего((.
Как фара от автомобиля, например...
0
4 / 4 / 2
Регистрация: 23.12.2016
Сообщений: 124
11.12.2017, 20:24
Что бы работала программа нужен метод с именем main (главный метод).Сдесь нету такого метода, исходя из этого я предположу что вы вырвали что-то из контекста и сами не знаете что. Советую начинать заниматься с более простых заданий. Найти вы из сможете например на сайте e-olymp или других подобных.
0
 Аватар для Aviz__
2755 / 2062 / 509
Регистрация: 17.02.2014
Сообщений: 9,491
12.12.2017, 09:51
Цитата Сообщение от WWH Посмотреть сообщение
и сами не знаете что
слишком сильное утверждение, коллега! Учитывая, что мы не знаем человека совсем...

Добавлено через 2 минуты
Цитата Сообщение от Arthuric Посмотреть сообщение
отсутствие каких-то файлов
каких конкретно?
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
raxper
Эксперт
30234 / 6612 / 1498
Регистрация: 28.12.2010
Сообщений: 21,154
Блог
12.12.2017, 09:51
Помогаю со студенческими работами здесь

Частотный детектор
помогите обосновать зачем в данной схеме сопротивление R3. т.е. что произойдет если сопротивление не будет.

Детектор нуля
В общем возникли непредвиденные проблемы. На мой взгляд простейший компаратор и будет им. Помоделировал, получил, что и ожидалось. ...

Детектор положения
Добрый день, уважаемые участники форума. Преамбула: у многих охотников есть собака aka такса, которая забегает в нору и зажимает дичь в...

Детектор темноты
Ломаю голову как сделать детектор темноты выдающий 1 сек импульс. Есть фотодиод bpw34 и транзисторы. Напряжение питания 2.5В. ...

Детектор огибающей
Здравствуйте. Имеется сигнал - АКФ радиоимпульса. Необходимо получить его огибающую (как я понимаю, должен быть треугольник). Сам сигнал...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
5
Ответ Создать тему
Новые блоги и статьи
1С: Контроль уникальности заводского номера
Maks 23.03.2026
Алгоритм контроля уникальности заводского (или серийного) номера на примере документа выдачи шин для спецтехники с табличной частью. Данные берутся из регистра сведений, по которому настроено. . .
Хочу заставить корпорации вкладываться в здоровье сотрудников: делаю мат модель здравосохранения
anaschu 22.03.2026
e7EYtONaj8Y Z4Tv2zpXVVo https:/ / github. com/ shumilovas/ med2. git
1С: Программный отбор элементов справочника по группе
Maks 22.03.2026
Установка программного отбора элементов справочника "Номенклатура" из модуля формы документа. В качестве фильтра для отбора справочника служит группа номенклатуры. Отбор по наименованию группы. . .
Как я обхитрил таблицу Word
Alexander-7 21.03.2026
Когда мигает курсор у внешнего края таблицы, и нам надо перейти на новую строку, а при нажатии Enter создается новый ряд таблицы с ячейками, то мы вместо нервных нажатий Энтеров мы пишем любые буквы. . .
Krabik - рыболовный бот для WoW 3.3.5a
AmbA 21.03.2026
без регистрации и смс. Это не торговля, приложение не содержит рекламы. Выполняет свою непосредственную задачу - автоматизацию рыбалки в WoW - и ничего более. Однако если админы будут против -. . .
1С: Программный отбор элементов справочника по значению перечисления
Maks 21.03.2026
Установка программного отбора элементов справочника "Сотрудники" из модуля формы документа. В качестве фильтра для отбора служит значение перечислений. / / Событие "НачалоВыбора" реквизита на форме. . .
Переходник USB-CAN-GPIO
Eddy_Em 20.03.2026
Достаточно давно на работе возникла необходимость в переходнике CAN-USB с гальваноразвязкой, оный и был разработан. Однако, все меня терзала совесть, что аж 48-ногий МК используется так тупо: просто. . .
Оттенки серого
Argus19 18.03.2026
Оттенки серого Нашёл в интернете 3 прекрасных модуля: Модуль класса открытия диалога открытия/ сохранения файла на Win32 API; Модуль класса быстрого перекодирования цветного изображения в оттенки. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru