forked from cruxopen/openISP
-
Notifications
You must be signed in to change notification settings - Fork 0
/
test_bnf.py
336 lines (313 loc) · 14.8 KB
/
test_bnf.py
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
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
from matplotlib import pyplot as plt
import numpy as np
import csv
# from model.dpc import DPC
# from model.blc import BLC
# from model.aaf import AAF
# from model.awb import WBGC
# from model.cnf import CNF
# from model.cfa import CFA
# from model.gac import GC
# from model.ccm import CCM
from model.csc import CSC
from model.bnf import BNF
# from model.eeh import EE
# from model.fcs import FCS
# from model.bcc import BCC
# from model.hsc import HSC
from model.nlm import NLM
raw_path = './raw/test.RAW'
# decomment for testing
# config_path = './config/config.csv'
config_path = './config/config_test.csv' # shrink the size of raw img for accelerating test result.
f = open(config_path, 'r', encoding='utf-8-sig')
with f:
reader = csv.reader(f, delimiter=',')
raw_h = 128 #for testing, original:1280
raw_w = 72 #for testing, original:720
dpc_thres = 30
dpc_mode = 'gradient'
dpc_clip = 1023
bl_r = 0
bl_gr = 0
bl_gb = 0
bl_b = 0
alpha = 0
beta = 0
blc_clip = 1023
bayer_pattern = 'rggb'
r_gain = 1.5
gr_gain = 1.0
gb_gain = 1.0
b_gain = 1.1
awb_clip = 1023
cfa_mode = 'malvar'
cfa_clip = 1023
ccm = np.zeros((3, 4))
csc = np.zeros((3, 4))
bnf_dw = np.zeros((5,5))
bnf_rw = [1, 1, 1, 1]
bnf_rthres = [32, 64, 128]
bnf_clip = 255
edge_filter = np.zeros((3, 5))
ee_gain = [32, 128]
ee_thres = [32, 64]
ee_emclip = [-64, 64]
fcs_edge = [32, 64]
fcs_gain = 32
fcs_intercept = 2
fcs_slope = 3
hue = 128
saturation = 256
hsc_clip = 255
brightness = 10 # [-255, 255]
contrast = 10 / pow(2, 5) # [-32,128]
bcc_clip = 255
nlm_h = 10
nlm_clip = 255
for row in reader:
parameter = row[0]
value = row[1]
description = row[2]
print(parameter, value, description)
if 'raw' in str(parameter):
raw_w = int(value) if '_w' in str(parameter) else raw_w
raw_h = int(value) if '_h' in str(parameter) else raw_h
# elif 'dpc' in str(parameter):
# dpc_thres = int(value) if '_thres' in str(parameter) else dpc_thres
# dpc_mode = str(value) if '_mode' in str(parameter) else dpc_mode
# dpc_clip = int(value) if '_clip' in str(parameter) else dpc_clip
# elif 'bl' in str(parameter):
# bl_r = int(value) if '_r' in str(parameter) else bl_r
# bl_gr = int(value) if '_gr' in str(parameter) else bl_gr
# bl_gb = int(value) if '_gb' in str(parameter) else bl_gb
# bl_b = int(value) if '_b' in str(parameter) else bl_b
# alpha = int(value) if '_alpha' in str(parameter) else alpha
# beta = int(value) if '_beta' in str(parameter) else beta
# blc_clip = int(value) if '_clip' in str(parameter) else beta
# elif 'bayer_pattern' in str(parameter):
# bayer_pattern = str(value)
# elif 'awb' in str(parameter):
# r_gain = int(value) if '_rgain' in str(parameter) else r_gain
# gr_gain = int(value) if '_grgain' in str(parameter) else gr_gain
# gb_gain = int(value) if '_gbgain' in str(parameter) else gb_gain
# b_gain = int(value) if '_bgain' in str(parameter) else b_gain
# awb_clip = int(value) if '_clip' in str(parameter) else awb_clip
# elif 'cfa' in str(parameter):
# cfa_mode = str(value) if '_mode' in str(parameter) else cfa_mode
# cfa_clip = int(value) if '_clip' in str(parameter) else cfa_clip
# elif 'ccm' in str(parameter):
# ccm[0][0] = int(value) if '_00' in str(parameter) else ccm[0][0]
# ccm[0][1] = int(value) if '_01' in str(parameter) else ccm[0][1]
# ccm[0][2] = int(value) if '_02' in str(parameter) else ccm[0][2]
# ccm[0][3] = int(value) if '_03' in str(parameter) else ccm[0][3]
# ccm[1][0] = int(value) if '_10' in str(parameter) else ccm[1][0]
# ccm[1][1] = int(value) if '_11' in str(parameter) else ccm[1][1]
# ccm[1][2] = int(value) if '_12' in str(parameter) else ccm[1][2]
# ccm[1][3] = int(value) if '_13' in str(parameter) else ccm[1][3]
# ccm[2][0] = int(value) if '_20' in str(parameter) else ccm[2][0]
# ccm[2][1] = int(value) if '_21' in str(parameter) else ccm[2][1]
# ccm[2][2] = int(value) if '_22' in str(parameter) else ccm[2][2]
# ccm[2][3] = int(value) if '_23' in str(parameter) else ccm[2][3]
# elif 'csc' in str(parameter):
# csc[0][0] = 1024 * float(value) if '_00' in str(parameter) else csc[0][0]
# csc[0][1] = 1024 * float(value) if '_01' in str(parameter) else csc[0][1]
# csc[0][2] = 1024 * float(value) if '_02' in str(parameter) else csc[0][2]
# csc[0][3] = 1024 * float(value) if '_03' in str(parameter) else csc[0][3]
# csc[1][0] = 1024 * float(value) if '_10' in str(parameter) else csc[1][0]
# csc[1][1] = 1024 * float(value) if '_11' in str(parameter) else csc[1][1]
# csc[1][2] = 1024 * float(value) if '_12' in str(parameter) else csc[1][2]
# csc[1][3] = 1024 * float(value) if '_13' in str(parameter) else csc[1][3]
# csc[2][0] = 1024 * float(value) if '_20' in str(parameter) else csc[2][0]
# csc[2][1] = 1024 * float(value) if '_21' in str(parameter) else csc[2][1]
# csc[2][2] = 1024 * float(value) if '_22' in str(parameter) else csc[2][2]
# csc[2][3] = 1024 * float(value) if '_23' in str(parameter) else csc[2][3]
elif 'bnf' in str(parameter):
bnf_dw[0][0] = int(value) if '_dw_00' in str(parameter) else bnf_dw[0][0]
bnf_dw[0][1] = int(value) if '_dw_01' in str(parameter) else bnf_dw[0][1]
bnf_dw[0][2] = int(value) if '_dw_02' in str(parameter) else bnf_dw[0][2]
bnf_dw[0][3] = int(value) if '_dw_03' in str(parameter) else bnf_dw[0][3]
bnf_dw[0][4] = int(value) if '_dw_04' in str(parameter) else bnf_dw[0][4]
bnf_dw[1][0] = int(value) if '_dw_10' in str(parameter) else bnf_dw[1][0]
bnf_dw[1][1] = int(value) if '_dw_11' in str(parameter) else bnf_dw[1][1]
bnf_dw[1][2] = int(value) if '_dw_12' in str(parameter) else bnf_dw[1][2]
bnf_dw[1][3] = int(value) if '_dw_13' in str(parameter) else bnf_dw[1][3]
bnf_dw[1][4] = int(value) if '_dw_14' in str(parameter) else bnf_dw[1][4]
bnf_dw[2][0] = int(value) if '_dw_20' in str(parameter) else bnf_dw[2][0]
bnf_dw[2][1] = int(value) if '_dw_21' in str(parameter) else bnf_dw[2][1]
bnf_dw[2][2] = int(value) if '_dw_22' in str(parameter) else bnf_dw[2][2]
bnf_dw[2][3] = int(value) if '_dw_23' in str(parameter) else bnf_dw[2][3]
bnf_dw[2][4] = int(value) if '_dw_24' in str(parameter) else bnf_dw[2][4]
bnf_rw[0] = int(value) if '_rw_0' in str(parameter) else bnf_rw[0]
bnf_rw[1] = int(value) if '_rw_1' in str(parameter) else bnf_rw[1]
bnf_rw[2] = int(value) if '_rw_2' in str(parameter) else bnf_rw[2]
bnf_rw[3] = int(value) if '_rw_3' in str(parameter) else bnf_rw[3]
bnf_rthres[0] = int(value) if '_rthres_0' in str(parameter) else bnf_rthres[0]
bnf_rthres[1] = int(value) if '_rthres_1' in str(parameter) else bnf_rthres[1]
bnf_rthres[2] = int(value) if '_rthres_2' in str(parameter) else bnf_rthres[2]
bnf_clip = int(value) if '_clip' in str(parameter) else bnf_clip
# elif 'edge_filter' in str(parameter):
# edge_filter[0][0] = int(value) if '_00' in str(parameter) else edge_filter[0][0]
# edge_filter[0][1] = int(value) if '_01' in str(parameter) else edge_filter[0][1]
# edge_filter[0][2] = int(value) if '_02' in str(parameter) else edge_filter[0][2]
# edge_filter[0][3] = int(value) if '_03' in str(parameter) else edge_filter[0][3]
# edge_filter[0][4] = int(value) if '_04' in str(parameter) else edge_filter[0][4]
# edge_filter[1][0] = int(value) if '_10' in str(parameter) else edge_filter[1][0]
# edge_filter[1][1] = int(value) if '_11' in str(parameter) else edge_filter[1][1]
# edge_filter[1][2] = int(value) if '_12' in str(parameter) else edge_filter[1][2]
# edge_filter[1][3] = int(value) if '_13' in str(parameter) else edge_filter[1][3]
# edge_filter[1][4] = int(value) if '_14' in str(parameter) else edge_filter[1][4]
# edge_filter[2][0] = int(value) if '_20' in str(parameter) else edge_filter[2][0]
# edge_filter[2][1] = int(value) if '_21' in str(parameter) else edge_filter[2][1]
# edge_filter[2][2] = int(value) if '_22' in str(parameter) else edge_filter[2][2]
# edge_filter[2][3] = int(value) if '_23' in str(parameter) else edge_filter[2][3]
# edge_filter[2][4] = int(value) if '_24' in str(parameter) else edge_filter[2][4]
# elif 'ee' in str(parameter):
# ee_gain[0] = int(value) if 'gain_min' in str(parameter) else ee_gain[0]
# ee_gain[1] = int(value) if 'gain_max' in str(parameter) else ee_gain[1]
# ee_thres[0] = int(value) if 'thres_min' in str(parameter) else ee_thres[0]
# ee_thres[1] = int(value) if 'thres_max' in str(parameter) else ee_thres[1]
# ee_emclip[0] = int(value) if 'emclip_min' in str(parameter) else ee_emclip[0]
# ee_emclip[1] = int(value) if 'emclip_max' in str(parameter) else ee_emclip[1]
# elif 'fcs' in str(parameter):
# fcs_edge[0] = int(value) if 'edge_min' in str(parameter) else fcs_edge[0]
# fcs_edge[1] = int(value) if 'edge_min' in str(parameter) else fcs_edge[1]
# fcs_gain = int(value) if '_gain' in str(parameter) else fcs_gain
# fcs_intercept = int(value) if '_intercept' in str(parameter) else fcs_intercept
# fcs_slope = int(value) if '_slope' in str(parameter) else fcs_slope
# elif 'nlm' in str(parameter):
# nlm_h = int(value) if '_h' in str(parameter) else nlm_h
# nlm_clip = int(value) if '_clip' in str(parameter) else nlm_clip
# else:
# hue = int(value) if 'hue' in str(parameter) else hue
# saturation = int(value) if 'saturation' in str(parameter) else saturation
# hsc_clip = int(value) if 'hsc_clip' in str(parameter) else hsc_clip
# brightness = int(value) if 'brightness' in str(parameter) else brightness
# contrast = int(value) if 'contrast' in str(parameter) else contrast
# bcc_clip = int(value) if 'bcc_clip' in str(parameter) else bcc_clip
rawimg = np.fromfile(raw_path, dtype='uint16', sep='')
print(rawimg.shape)
# rawimg = rawimg.reshape([raw_h, raw_w])
# for testing
rawimg = rawimg.reshape([raw_h*10, raw_w*10])
rawimg = rawimg[:raw_h,:raw_w]
print(50*'-' + '\nLoading RAW Image Done......')
#plt.imshow(rawimg, cmap='gray')
#plt.show()
# # dead pixel correction
# dpc = DPC(rawimg, dpc_thres, dpc_mode, dpc_clip)
# rawimg_dpc = dpc.execute()
# print(50*'-' + '\nDead Pixel Correction Done......')
# #plt.imshow(rawimg_dpc, cmap='gray')
# #plt.show()
# # black level compensation
# parameter = [bl_r, bl_gr, bl_gb, bl_b, alpha, beta]
# blc = BLC(rawimg_dpc, parameter, bayer_pattern, blc_clip)
# rawimg_blc = blc.execute()
# print(50*'-' + '\nBlack Level Compensation Done......')
# #plt.imshow(rawimg_blc, cmap='gray')
# #plt.show()
# # lens shading correction
# # anti-aliasing filter
# aaf = AAF(rawimg_blc)
# rawimg_aaf = aaf.execute()
# print(50*'-' + '\nAnti-aliasing Filtering Done......')
# #plt.imshow(rawimg_aaf, cmap='gray')
# #plt.show()
# #rawimg_diff = rawimg_blc - rawimg_aaf
# #plt.imshow(rawimg_diff, cmap='gray')
# #plt.show()
# # white balance gain control
# parameter = [r_gain, gr_gain, gb_gain, b_gain]
# awb = WBGC(rawimg_aaf, parameter, bayer_pattern, awb_clip)
# rawimg_awb = awb.execute()
# print(50*'-' + '\nWhite Balance Gain Done......')
# #plt.imshow(rawimg_awb, cmap='gray')
# #plt.show()
# # chroma noise filtering
# cnf = CNF(rawimg_awb, bayer_pattern, 0, parameter, 1023)
# rawimg_cnf = cnf.execute()
# print(50*'-' + '\nChroma Noise Filtering Done......')
# #plt.imshow(rawimg_cnf/4, cmap='gray')
# #plt.show()
# # color filter array interpolation
# cfa = CFA(rawimg_cnf, cfa_mode, bayer_pattern, cfa_clip)
# rgbimg_cfa = cfa.execute()
# print(50*'-' + '\nDemosaicing Done......')
# #plt.imshow(rgbimg_cfa/4)
# #plt.show()
# # color correction matrix
# ccm = CCM(rgbimg_cfa, ccm)
# rgbimg_ccm = ccm.execute()
# print(50*'-' + '\nColor Correction Done......')
# #plt.imshow(rgbimg_ccm)
# #plt.show()
# # gamma correction
# # look up table
# bw = 10
# gamma = 0.5
# mode = 'rgb'
# maxval = pow(2,bw)
# ind = range(0, maxval)
# val = [round(pow(float(i)/maxval, gamma) * maxval) for i in ind]
# lut = dict(zip(ind, val))
# #print(ind, val, lut)
# gc = GC(rgbimg_ccm, lut, mode)
# rgbimg_gc = gc.execute()
# print(50*'-' + '\nGamma Correction Done......')
# #plt.imshow(rgbimg_gc)
# #plt.show()
# # color space conversion
# csc = CSC(rgbimg_ccm, csc)
# yuvimg_csc = csc.execute()
# print(50*'-' + '\nColor Space Conversion Done......')
# plt.imshow(yuvimg_csc[:,:,0], cmap='gray')
# #plt.show()
# non-local means denoising
# nlm = NLM(yuvimg_csc[:,:,0], 1, 4, nlm_h, nlm_clip)
# yuvimg_nlm = nlm.execute()
# print(50*'-' + '\nNon Local Means Denoising Done......')
#plt.imshow(yuvimg_nlm, cmap='gray')
#plt.show()
# bilateral filter
# bnf = BNF(yuvimg_nlm, bnf_dw, bnf_rw, bnf_rthres, bnf_clip)
# change yuvimg_nlm to rawimg
bnf = BNF(rawimg, bnf_dw, bnf_rw, bnf_rthres, bnf_clip)
# yuvimg_bnf = bnf.execute()
# change yuvimg to rawimg
rawimg_bnf = bnf.execute()
print(50*'-' + '\nBilateral Filtering Done......')
#plt.imshow(yuvimg_bnf, cmap='gray')
#plt.show()
# # edge enhancement
# ee = EE(yuvimg_bnf[:,:], edge_filter, ee_gain, ee_thres, ee_emclip)
# yuvimg_ee, yuvimg_edgemap = ee.execute()
# print(50*'-' + '\nEdge Enhancement Done......')
# #plt.imshow(yuvimg_ee)
# #plt.show()
# #plt.imshow(yuvimg_edgemap)
# #plt.show()
# # false color suppresion
# fcs = FCS(yuvimg_csc[:,:,1:3], yuvimg_edgemap, fcs_edge, fcs_gain, fcs_intercept, fcs_slope)
# yuvimg_fcs = fcs.execute()
# print(50*'-' + '\nFalse Color Suppresion Done......')
# #plt.imshow(yuvimg_fcs)
# #plt.show()
# # hue/saturation control
# hsc = HSC(yuvimg_fcs, hue, saturation, hsc_clip)
# yuvimg_hsc = hsc.execute()
# print(50*'-' + '\nHue/Saturation Adjustment Done......')
# #plt.imshow(yuvimg_hsc)
# #plt.show()
# # brighyness/contrast control
# contrast = contrast / pow(2,5) #[-32,128]
# bcc = BCC(yuvimg_ee, brightness, contrast, bcc_clip)
# yuvimg_bcc = bcc.execute()
# print(50*'-' + '\nBrightness/Contrast Adjustment Done......')
# #plt.imshow(yuvimg_bcc)
# #plt.show()
# yuvimg_out = np.empty((raw_h, raw_w, 3), dtype=np.uint8)
# yuvimg_out[:,:,0] = yuvimg_bcc
# yuvimg_out[:,:,1:3] = yuvimg_hsc
# #plt.imshow(yuvimg_out)
# #plt.show()