Форум программистов, компьютерный форум CyberForum.ru

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

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 10, средняя оценка - 4.70
Aldragon
2 / 2 / 0
Регистрация: 28.09.2012
Сообщений: 34
24.03.2013, 13:03     Прямое и обратное преобразование Хаара #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);
Буду благодарен, за любой совет.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
24.03.2013, 13:03     Прямое и обратное преобразование Хаара
Посмотрите здесь:

C++ прямое присвоение двух массивов в С++
C++ Прямое обращение к жесткому диску
C++ Каскады хаара
Прямое и обратное преобразование чисел в Код Грея C++
Товарищи помогите с вейвлетом Хаара C++
C++ Вейвлет Хаара 2.Ввод и вывод.
Обратное преобразование Лапласа C++
Прямое отображение слов C++

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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
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){
но восстановить все ровно не получается.
Никто ничего больше не видит?
Yandex
Объявления
25.03.2013, 13:40     Прямое и обратное преобразование Хаара
Ответ Создать тему
Опции темы

Текущее время: 01:49. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru