Форум программистов, компьютерный форум, киберфорум
Наши страницы

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Рейтинг: Рейтинг темы: голосов - 10, средняя оценка - 4.70
Aldragon
2 / 2 / 0
Регистрация: 28.09.2012
Сообщений: 34
#1

Прямое и обратное преобразование Хаара - C++

24.03.2013, 13:03. Просмотров 1548. Ответов 1
Метки нет (Все метки)

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

C++
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
#include "stdafx.h"
#include "libopencv.h"
#include <fstream>
#include <math.h>
 
#define N 32
 
using namespace std;
using namespace cv;
 
 
int main(int argc, char* argv[]){
    if(argc != 2){
        cout << "Too few arguments" << endl;
        system("pause");
        return(-1);
    }
    Mat Image(N, N, CV_32FC1, Scalar(0, 0, 0, 0));
    ifstream in(argv[1]);
    if(!(in.is_open())){
        cout << "Error of opening file" << endl;
        system("pause");
        return(-1);
    }
    for(int i=0; i<N; i++){
        for(int j=0; j<N; j++){
            float tmp;
            in >> tmp;
            Image.at<float>(i,j) = tmp;
        }
    }
    normalize(Image, Image, 0, 1, NORM_MINMAX);
    namedWindow("SourceImage", CV_WINDOW_FREERATIO);
    imshow("SourceImage", Image);
    int Row = Image.rows, Col = Image.cols;
    Mat Zero(Row, Col, CV_32FC1, Scalar(0, 0, 0));
    Mat BufMat(Row, Col, CV_32FC1, Scalar(0, 0, 0));
    for (int k=Row; k>=1; k/=2){
        for (int i=0; i<Row; i++){
            for(int j=0; j<Col/2; j++){
                BufMat.at<float>(i,j) = (Image.at<float>(i,2*j) + Image.at<float>(i,2*j+1))/2;
                BufMat.at<float>(i, Col/2+j) = fabs((Image.at<float>(i,2*j) - Image.at<float>(i,2*j+1))/2);
            }
        }
        Image = BufMat;
        BufMat = Zero;
        for(int i=0; i<Col; i++){
            for(int j=0; j<Row/2; j++){
                BufMat.at<float>(j,i) = (Image.at<float>(2*j,i) + Image.at<float>(2*j+1,i))/2;
                BufMat.at<float>(Row/2+j, i) = fabs((Image.at<float>(2*j,i) - Image.at<float>(2*j+1,i))/2);
            }
        }
        Image = BufMat;
        Mat TMP;
        normalize(BufMat, TMP, 0, 1, NORM_MINMAX);
        namedWindow("Step", CV_WINDOW_FREERATIO);
        imshow("Step", TMP);
        cout << k << endl;
        waitKey();
    }
    ofstream Res("ResOfExpansion.dat");
    for(int i=0; i<Row; i++){
        for(int j=0; j<Col; j++){
            Res << Image.at<float>(i,j) << " ";
        }
        Res << endl;
    }
    Res.close();
    in.close();
    return(0);
}

Обратное преобразование.
Кликните здесь для просмотра всего текста
C++
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
#include "stdafx.h"
#include "libopencv.h"
#include <fstream>
#include <math.h>
 
#define N 32
 
using namespace std;
using namespace cv;
 
 
int main(int argc, char* argv[]){
    //if(argc != 2){
    //  cout << "Too few arguments" << endl;
    //  system("pause");
    //  return(-1);
    //}
    Mat Image(N, N, CV_32FC1, Scalar(0, 0, 0, 0));
    Mat TMP;
    ifstream in("img-0.dat");
    if(!(in.is_open())){
        cout << "Error of opening file" << endl;
        system("pause");
        return(-1);
    }
    for(int i=0; i<N; i++){
        for(int j=0; j<N; j++){
            float tmp;
            in >> tmp;
            Image.at<float>(i,j) = tmp;
        }
    }
    normalize(Image, TMP, 0, 1, NORM_MINMAX);
    namedWindow("SourceImage", CV_WINDOW_FREERATIO);
    imshow("SourceImage", TMP);
    int Row = Image.rows, Col = Image.cols;
    Mat Zero(Row, Col, CV_32FC1, Scalar(0, 0, 0));
    Mat BufMat(Row, Col, CV_32FC1, Scalar(0, 0, 0));
    for (int k=1; k<=Row; k*=2){
        for (int i=0; i<Row; i++){
            for(int j=0; j<Col/2; j++){
                BufMat.at<float>(i,2*j) = Image.at<float>(i,j) + Image.at<float>(i,j+Col/2);
                BufMat.at<float>(i,2*j+1) = Image.at<float>(i,j) - Image.at<float>(i,j+Col/2);
            }
        }
        Image = BufMat;
        BufMat = Zero;
        for(int i=0; i<Col; i++){
            for(int j=0; j<Row/2; j++){
                BufMat.at<float>(2*j,i) = Image.at<float>(j,i) + Image.at<float>(j+Row/2,i);
                BufMat.at<float>(2*j+1,i) = Image.at<float>(j,i) - Image.at<float>(j+Row/2,i);
            }
        }
        Image = BufMat;
        normalize(BufMat, TMP, 0, 1, NORM_MINMAX);
        namedWindow("Step", CV_WINDOW_FREERATIO);
        imshow("Step", TMP);
        cout << k << endl;
        waitKey();
    }
    in.close();
    return(0);
}

Для визуализации использую OpenCV.

Если выполнить всего один шаг сжатия(можно сказать, изображение делится на 4 части), все отлично восстанавливается. Но уже для двух шагов после восстановления получается какая то билиберда. Не могу понять в какой программе ошибка.
Сильно сомневаюсь в необходимости функции fabs(), в первой программе, хотя без нее получается еще большая билеберда.
C++
1
BufMat.at<float>(Row/2+j, i) = fabs((Image.at<float>(2*j,i) - Image.at<float>(2*j+1,i))/2);
Буду благодарен, за любой совет.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
24.03.2013, 13:03
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Прямое и обратное преобразование Хаара (C++):

Прямое и обратное преобразование функции - C++
Всем привет! Нужно написать программу прямого и обратного преобразования функции (скорее всего дискретное преобразование Фурье). У кого...

Прямое и обратное преобразование чисел в Код Грея - C++
Здравствуйте. Есть вот такая интересная задача, надо прямое и обратное преобразование чисел в Код Грея. Допустим числа от 1 до 10. ...

Прямое и обратное отображение набранных символов в консоли - C++
#include &lt;iostream&gt; using namespace std; int main() { string s; int i; cin&gt;&gt;s; while (i=length(s)&lt;=1) cout&lt;&lt;s; ...

Обратное преобразование Лапласа - C++
Доброго времени суток. Подскажите пожалуйста как запрограммировать обратное преобразование Лапласа. К примеру у меня есть число 0,7 и к...

Обратное быстрое преобразование Фурье - C++
Прямое преобразование работает корректно, а вот обратное возвращает странную дугу в независимсти от исходного сигнала. &quot;Не быстрый&quot;...

Обратное вейвлет-преобразование - C++
Здравствуйте. Пишу код двумерного вейвлет-преобразования изображений. Использую вейвлет Хаара. Коэффициенты матриц разложения считаются...

1
Aldragon
2 / 2 / 0
Регистрация: 28.09.2012
Сообщений: 34
25.03.2013, 13:40  [ТС] #2
нашел одну ошибку в разложении. в 38 строке, была лишняя итерация.
C++
1
for (int k=Row; k>1; k/=2){
но восстановить все ровно не получается.
Никто ничего больше не видит?
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
25.03.2013, 13:40
Привет! Вот еще темы с ответами:

Каскады хаара - C++
Всем привет. Кто нибудь пожалуйста дайте ссылку для скачивания cvhaartraining.lib, очень нужна эта библа. Заранее спасибо!

Вейвлет Хаара 2.Ввод и вывод. - C++
Есть две программы. 1-ая это прямое преобразование Хаара. 2-ая это собственно обратное. Что нужно бы сделать. Нужно считать файл input.dat...

Товарищи помогите с вейвлетом Хаара - C++
#include &lt;iostream&gt; #include &lt;time.h&gt; using namespace std; int main () { int n; std::cout &lt;&lt; &quot;Enter n:\n&gt;&quot;; ...

Прямое отображение слов - C++
Доброго времени суток. Дано слово pass. А отведено 13 элементов. Нужно сделать,чтобы получилось passpasspassp есть ровно 13...


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

Или воспользуйтесь поиском по форуму:
2
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.