/
raspi_dng.c
190 lines (158 loc) · 6.55 KB
/
raspi_dng.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
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
/* Read in jpeg from Raspberry Pi camera captured using 'raspistill --raw'
and extract raw file with 10-bit values stored left-justified at 16 bpp
in Adobe DNG (TIFF-EP) format, convert with 'ufraw out.dng' for example
John Beale 26 May 2013
and others
Contains code written by Dave Coffin for Berkeley Engineering and Research.
Free for all uses.
Requires LibTIFF 3.8.0 plus a patch, see http://www.cybercom.net/~dcoffin/dcraw/
Compile with:
gcc -o raspi_dng raspi_dng.c -O4 -Wall -lm -ljpeg -ltiff
*/
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <sys/stat.h>
#include <time.h>
#include <math.h>
#include <tiffio.h>
#include <errno.h>
#define LINELEN 256 // how long a string we need to hold the filename
#define RAWBLOCKSIZE 6404096
#define HEADERSIZE 32768
#define ROWSIZE 3264 // number of bytes per row of pixels, including 24 'other' bytes at end
#define IDSIZE 4 // number of bytes in raw header ID string
#define HPIXELS 2592 // number of horizontal pixels on OV5647 sensor
#define VPIXELS 1944 // number of vertical pixels on OV5647 sensor
int main (int argc, char **argv)
{
static const short CFARepeatPatternDim[] = { 2,2 };
// this color matrix is definitely inaccurate, TODO: calibrate
static const float cam_xyz[] = {
// R G B
1.000, 0.000, 0.000, // R
0.000, 1.000, 0.000, // G
0.000, 0.000, 1.000 // B
};
static const float neutral[] = { 1.0, 1.0, 1.0 }; // TODO calibrate
long sub_offset=0, white=0xffff;
int status=1, i, j, row, col;
unsigned short curve[256];
struct stat st;
struct tm tm;
char datetime[64];
FILE *ifp;
TIFF *tif;
const char* fname = argv[1];
unsigned long fileLen; // number of bytes in file
unsigned long offset; // offset into file to start reading pixel data
unsigned char *buffer;
unsigned short pixel[HPIXELS]; // array holds 16 bits per pixel
unsigned char split; // single byte with 4 pairs of low-order bits
if (argc != 3) {
fprintf (stderr, "Usage: %s infile outfile\n"
"Example: %s rpi.jpg output.dng\n", argv[0], argv[0]);
return 1;
}
if (!(ifp = fopen (fname, "rb"))) {
perror (argv[2]);
return 1;
}
stat (fname, &st);
gmtime_r (&st.st_mtime, &tm);
sprintf (datetime, "%04d:%02d:%02d %02d:%02d:%02d",
tm.tm_year+1900,tm.tm_mon+1,tm.tm_mday,tm.tm_hour,tm.tm_min,tm.tm_sec);
//Get file length
fseek(ifp, 0, SEEK_END);
fileLen=ftell(ifp);
if (fileLen < RAWBLOCKSIZE) {
fprintf(stderr, "File %s too short to contain expected 6MB RAW data.\n", fname);
exit(1);
}
offset = (fileLen - RAWBLOCKSIZE) ; // location in file the raw header starts
fseek(ifp, offset, SEEK_SET);
//printf("File length = %d bytes.\n",fileLen);
//printf("offset = %d:",offset);
//Allocate memory for one line of pixel data
buffer=(unsigned char *)malloc(ROWSIZE+1);
if (!buffer)
{
fprintf(stderr, "Memory error!");
status = ENOMEM;
goto fail;
}
if (!(tif = TIFFOpen (argv[2], "w"))) goto fail;
//fprintf(stderr, "Writing TIFF header...\n");
TIFFSetField (tif, TIFFTAG_SUBFILETYPE, 1);
TIFFSetField (tif, TIFFTAG_IMAGEWIDTH, HPIXELS >> 4);
TIFFSetField (tif, TIFFTAG_IMAGELENGTH, VPIXELS >> 4);
TIFFSetField (tif, TIFFTAG_BITSPERSAMPLE, 8);
TIFFSetField (tif, TIFFTAG_COMPRESSION, COMPRESSION_NONE);
TIFFSetField (tif, TIFFTAG_PHOTOMETRIC, PHOTOMETRIC_RGB);
TIFFSetField (tif, TIFFTAG_MAKE, "Raspberry Pi");
TIFFSetField (tif, TIFFTAG_MODEL, "Model OV5647");
TIFFSetField (tif, TIFFTAG_ORIENTATION, ORIENTATION_TOPLEFT);
TIFFSetField (tif, TIFFTAG_SAMPLESPERPIXEL, 3);
TIFFSetField (tif, TIFFTAG_PLANARCONFIG, PLANARCONFIG_CONTIG);
TIFFSetField (tif, TIFFTAG_SOFTWARE, "raspi_dng");
TIFFSetField (tif, TIFFTAG_DATETIME, datetime);
TIFFSetField (tif, TIFFTAG_SUBIFD, 1, &sub_offset);
TIFFSetField (tif, TIFFTAG_DNGVERSION, "\001\001\0\0");
TIFFSetField (tif, TIFFTAG_DNGBACKWARDVERSION, "\001\0\0\0");
TIFFSetField (tif, TIFFTAG_UNIQUECAMERAMODEL, "Raspberry Pi - OV5647");
TIFFSetField (tif, TIFFTAG_COLORMATRIX1, 9, cam_xyz);
TIFFSetField (tif, TIFFTAG_ASSHOTNEUTRAL, 3, neutral);
TIFFSetField (tif, TIFFTAG_CALIBRATIONILLUMINANT1, 21);
TIFFSetField (tif, TIFFTAG_ORIGINALRAWFILENAME, fname);
// fprintf(stderr, "Writing TIFF thumbnail...\n");
memset (pixel, 0, HPIXELS); // all-black thumbnail
for (row=0; row < VPIXELS >> 4; row++)
TIFFWriteScanline (tif, pixel, row, 0);
TIFFWriteDirectory (tif);
// fprintf(stderr, "Writing TIFF header for main image...\n");
TIFFSetField (tif, TIFFTAG_SUBFILETYPE, 0);
TIFFSetField (tif, TIFFTAG_IMAGEWIDTH, HPIXELS);
TIFFSetField (tif, TIFFTAG_IMAGELENGTH, VPIXELS);
TIFFSetField (tif, TIFFTAG_BITSPERSAMPLE, 16);
TIFFSetField (tif, TIFFTAG_PHOTOMETRIC, PHOTOMETRIC_CFA);
TIFFSetField (tif, TIFFTAG_SAMPLESPERPIXEL, 1);
TIFFSetField (tif, TIFFTAG_PLANARCONFIG, PLANARCONFIG_CONTIG);
TIFFSetField (tif, TIFFTAG_CFAREPEATPATTERNDIM, CFARepeatPatternDim);
TIFFSetField (tif, TIFFTAG_CFAPATTERN, 4, "\002\001\001\0"); // 0 = Red, 1 = Green, 2 = Blue, 3 = Cyan, 4 = Magenta, 5 = Yellow, 6 = White
//TIFFSetField (tif, TIFFTAG_LINEARIZATIONTABLE, 256, curve);
TIFFSetField (tif, TIFFTAG_WHITELEVEL, 1, &white);
fprintf(stderr, "Processing RAW data...\n");
// for one file, (TotalFileLength:11112983 - RawBlockSize:6404096) + Header:32768 = 4741655
// The pixel data is arranged in the file in rows, with 3264 bytes per row.
// with 3264 bytes per row x 1944 rows we have 6345216 bytes, that is the full 2592x1944 image area.
//Read one line of pixel data into buffer
fread(buffer, IDSIZE, 1, ifp);
// now on to the pixel data
offset = (fileLen - RAWBLOCKSIZE) + HEADERSIZE; // location in file the raw pixel data starts
fseek(ifp, offset, SEEK_SET);
for (row=0; row < VPIXELS; row++) { // iterate over pixel rows
fread(buffer, ROWSIZE, 1, ifp); // read next line of pixel data
j = 0; // offset into buffer
for (col = 0; col < HPIXELS; col+= 4) { // iterate over pixel columns
pixel[col+0] = buffer[j++] << 8;
pixel[col+1] = buffer[j++] << 8;
pixel[col+2] = buffer[j++] << 8;
pixel[col+3] = buffer[j++] << 8;
split = buffer[j++]; // low-order packed bits from previous 4 pixels
pixel[col+0] += (split & 0b11000000); // unpack them bits, add to 16-bit values, left-justified
pixel[col+1] += (split & 0b00110000)<<2;
pixel[col+2] += (split & 0b00001100)<<4;
pixel[col+3] += (split & 0b00000011)<<6;
}
if (TIFFWriteScanline (tif, pixel, row, 0) != 1) {
fprintf(stderr, "Error writing TIFF scanline.");
exit(1);
}
} // end for(k..)
free(buffer); // free up that memory we allocated
TIFFClose (tif);
status = 0;
fail:
fclose (ifp);
return status;
}