forked from jbresearch/cr2_scripts
-
Notifications
You must be signed in to change notification settings - Fork 0
/
cr2_replace.py
executable file
·96 lines (85 loc) · 3.11 KB
/
cr2_replace.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
#!/usr/bin/env python
# -*- coding: utf-8 -*-
#
# Copyright © 2015-2018 Johann A. Briffa
#
# This file is part of CR2_Scripts.
#
# CR2_Scripts is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# CR2_Scripts is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with CR2_Scripts. If not, see <http:https://www.gnu.org/licenses/>.
import sys
import os
import tempfile
import argparse
sys.path.append(os.path.join(os.path.dirname(os.path.realpath(__file__)),'pyshared'))
import jbtiff
import jbcr2
import jbimage
## main program
def main():
# interpret user options
parser = argparse.ArgumentParser()
parser.add_argument("-i", "--input", required=True,
help="input raw file to use as basis")
parser.add_argument("-d", "--decode",
help="decode and keep previous sensor image file")
parser.add_argument("-s", "--sensor", required=True,
help="sensor image file to replace input")
parser.add_argument("-o", "--output", required=True,
help="output CR2 file")
args = parser.parse_args()
# read input raw file
tiff = jbtiff.tiff_file(open(args.input, 'rb'))
# obtain required parameters from RAW file
width,height = tiff.get_sensor_size()
slices = tiff.get_slices()
# extract existing sensor image
IFD, ifd_offset, strips = tiff.data[3]
# save into a temporary file
fid, tmpfile = tempfile.mkstemp()
for strip in strips:
os.write(fid, strip)
os.close(fid)
# decode lossless JPEG encoded file to determine parameters
a, components, precision = jbcr2.decode_lossless_jpeg(tmpfile)
# delete temporary file
os.remove(tmpfile)
# store decoded sensor image file as needed
if args.decode:
# unslice image
I = jbcr2.unslice_image(a, width, height, slices)
# save result
jbimage.pnm_file.write(I.astype('>H'), open(args.decode,'wb'))
# read sensor image file
sensor = jbimage.image_file.read(args.sensor).squeeze()
# check input image parameters
assert len(sensor.shape) == 2 # must be a one-channel image
assert sensor.shape == (height,width) # image size must be exact
# slice image
a = jbcr2.slice_image(sensor, width, height, slices)
# encode to a temporary lossless JPEG output file
fid, tmpfile = tempfile.mkstemp()
os.close(fid)
parts = jbcr2.encode_lossless_jpeg(a, components, precision, tmpfile)
# read and delete temporary file
with open(tmpfile, 'rb') as fid:
data = fid.read()
os.remove(tmpfile)
# replace data strips for main sensor image (IFD#3)
jbcr2.replace_ifd(tiff, 3, data)
# save updated CR2 file
tiff.write(open(args.output,'wb'))
return
# main entry point
if __name__ == '__main__':
main()