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

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

Восстановить пароль Регистрация
Другие темы раздела
C++ SDL & OpenGL http://www.cyberforum.ru/cpp-beginners/thread817003.html
Можно при использовании связки SDL и OpenGL также использовать сторонние библиотеки для создания пользовательского интерфейса (например эту)? Или лучше использовать Qt?
C++ Создание программы на начальном этапе Подскажите, пожалуйста, какого типа проекта в С++ надо создать и куда именно вставить следующий программный код, что бы программа запустилась и можно было посмотреть на результат? #include <stdio.h> int main() { double b = 213.492222232; printf("%i",(int)(b+.5)); } http://www.cyberforum.ru/cpp-beginners/thread816977.html
C++ Распаковка PAK Файлов
Есть одна игра, новая Dance - MMO, я все ресурсы игры - модели,картинки тексты и т.д находятся в этом файле. Я пытался открыть этот файл с помощью всех программ которые открывает файлы .PAK но никакая программа не открыла этот файл. В игре есть такие 4 файла, у одного размер 2GB, у второго размер 750MB, у третьего 350MB, и у четвертого 1GB, так что кинуть вам файлы не смогу, так как много...
C++ В двумерном массиве определить первый отрицательный элемент на интервале [a,b]
Подскажите пожалуйста что мне нужно про писать в 5 строке чтобы определить первый отрицательный элемент массива не используя break for (i=1; i<3; i++) for (j=1; j<3; j++) if(a<= M && M<=b) cout <<"Первый отрицательный элемент массива в интервале ="<<z<<endl;
C++ Подключение flex/bison или antlr к с++ http://www.cyberforum.ru/cpp-beginners/thread816955.html
Подключение flex/bison или antlr к с++. Кто-нибудь с этим сталкивался? ничего толкового не нашёл, слышал ещё через CMake можно подключать, но у меня не получилось(
C++ Определить номер строки и номер столбца матрицы, на пересечении которых находятся требуемые элементы В двумерном массиве определить номер строки и номер столбца, на пересечении которых находятся требуемые элементы массива: 1) Последний отрицательный элемент; 3) Первый отрицательный попавший в интервал Добавлено через 5 часов 14 минут Не правильно работает строка 40 . помогите исправить ошибку #include <iostream> #include <conio.h> #include <clocale> using namespace std; подробнее

Показать сообщение отдельно
Aldragon
2 / 2 / 0
Регистрация: 28.09.2012
Сообщений: 34
24.03.2013, 13:03     Прямое и обратное преобразование Хаара
Здравствуйте. Написал вейвлет-сжатие Хаара, чисто визуально все работало правильно. Для надежности решил написать обратное преобразование и восстановить изображение.
Прямое преобразование
Кликните здесь для просмотра всего текста

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);
Буду благодарен, за любой совет.
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
Текущее время: 05:21. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru