Sync working changes into develop

.gitignore

@@ -2,6 +2,7 @@
 .ccls-cache/
 
 build/
+compile_commands.json
 
 # Prerequisites
 *.d

CMakeLists.txt

@@ -0,0 +1,16 @@
+cmake_minimum_required(VERSION 3.10)
+project(libezw)
+
+set(CMAKE_C_FLAGS "${CMAKE_C_FLAGS} -Wall")
+set(CMAKE_C_FLAGS_DEBUG "${CMAKE_C_FLAGS_DEBUG} -Wall")
+set(CMAKE_C_FLAGS_RELEASE "${CMAKE_C_FLAGS_RELEASE} -Wall")
+set(SOURCE_DIR ${CMAKE_CURRENT_SOURCE_DIR}/src)
+
+file(GLOB SOURCE_FILES "${SOURCE_DIR}/*.c")
+
+include_directories(include wavelib/header)
+add_executable(main ${SOURCE_FILES})
+
+target_link_libraries(main ${CMAKE_CURRENT_SOURCE_DIR}/libwavelib.a)
+target_link_libraries(main m) # link math.h
+

include/utils.h

@@ -0,0 +1,9 @@
+#ifndef UTILS_H
+#define UTILS_H
+#include<stdio.h>
+
+void read_binary_file(FILE *in, unsigned char *dest, int rows, int cols);
+void write_binary_file(FILE *out, unsigned char *dest, int rows, int cols);
+void peek_binary_file();
+
+#endif

src/main.c

@@ -0,0 +1,55 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include "utils.h"
+#include "wavelib.h"
+
+#define ROWS 512
+#define COLS 512
+
+int main()
+{
+    FILE *fptr = fopen("/home/ramprakash/rpdata/nie/study-material/sem6/miniproject/libezw/utils/raw_lichtenstein_img_processing_test_bw.bin", "rb");
+    unsigned char *pix_arr = calloc(ROWS*COLS, 1); // array of 8 bit vals
+    if(!fptr) {
+        fprintf(stderr, "Unable to open file\n");
+    }
+    else {
+        read_binary_file(fptr, pix_arr, ROWS, COLS);
+        fclose(fptr);
+        wave_object obj;
+        wt2_object wt;
+        double *inp, *wavecoeffs, *oup, *cLL;
+        int N = ROWS*COLS;
+        int ir, ic;
+
+        char *name = "db2";
+        obj = wave_init(name);// Initialize the wavelet
+        inp = (double*) calloc(N, sizeof(double));
+        oup = (double*) calloc(N, sizeof(double));
+
+        int J = 1; // number of decompositions
+
+        wt = wt2_init(obj, "dwt", ROWS, COLS, J);
+
+        for (int i = 0; i < ROWS; ++i) {
+            for (int k = 0; k < COLS; ++k) {
+                inp[i*COLS + k] = (double) pix_arr[i*COLS + k];
+                oup[i*COLS + k] = 0.0;
+            }
+        }
+        wavecoeffs = dwt2(wt, inp);
+        cLL = getWT2Coeffs(wt, wavecoeffs, 1, "A", &ir, &ic);
+
+        dispWT2Coeffs(cLL, ir, ic);
+
+        wt2_summary(wt);
+
+        // clean up
+        wave_free(obj);
+        wt2_free(wt);
+        free(inp);
+        free(wavecoeffs);
+        free(oup);
+    }
+    return 0;
+}

src/utils.c

@@ -0,0 +1,6 @@
+#include "utils.h"
+
+void read_binary_file(FILE *in, unsigned char* dest, int rows, int cols)
+{
+    fread(dest, 1, rows*cols, in);
+}

utils/get_gscale_raw.py

@@ -5,9 +5,20 @@
 WIDTH, HEIGHT = img.size
 
 data = list(img.getdata())
-data = [data[offset:offset+WIDTH] for offset in range(0, WIDTH*HEIGHT, WIDTH)]
-for row in data:
-    for value in row:
-        #  print(' '.join('{:3}'.format(value)), end='')
-        print(str(value) + ',', end='')
-    print()
+data_bytes = bytes(img.getdata())
+#  data = [data[offset:offset+WIDTH] for offset in range(0, WIDTH*HEIGHT, WIDTH)]
+
+with open("raw_lichtenstein_img_processing_test_bw.bin", "wb") as f:
+    f.write(data_bytes)
+#  for row in data:
+#      for value in row:
+#          #  print(' '.join('{:3}'.format(value)), end='')
+#          print(str(value) + ',', end='')
+#      print()
+
+#  chars = '@%#*+=-:. '  # Change as desired.
+#  scale = (len(chars)-1)/255.
+#  for row in data:
+#      for value in row:
+#          print(chars[int(value*scale)], end='')
+#      print()