I'm trying to write a raw image converter, thanks to LibRaw and the ecosystem around it which made my journey easier.

I've read the image processing code from LibRaw, but I'm not sure I understand the procedure correctly. If someone can point me to the right direction, it would be very helpful.

Let's assume the raw image is RGB rather than using CFA. This is what I do when a pixel is acquired:

• 1. Black level subtraction. All three channels minus the corresponding black level.
• 2. Scale the values. The channels scale according to the white and black level range.
• 3. Clipping. Clip the pixel values in case they go beyond the target boundaries, e.g. 0 or max.
• 4. White balancing. Multiply the white balance acquired by the raw library, i.e., wb_coeffs.
• 5. Colour space conversion. Convert the values from camera colour space to the target colour space, e.g. sRGB. First, I query the adobe matrix table using the camera ID, I call this matrix "xyz_2_cam". Second, I take "rgb_2_xyz" from LibRaw and multiply it with "xyz_2_cam", then use "pseudo_inverse" to inverse it getting "cam_2_xyz".

// rgb_2_xyz is the matrix for sRGB to xyz conversion:
static const Matrix<double> rgb_2_xyz = {
                {0.4124564, 0.3575761, 0.1804375},
                {0.2126729, 0.7151522, 0.0721750},
                {0.0193339, 0.1191920, 0.9503041}

This process doesn't produce the colour I wanted (LibRaw colour).

I don't know a better place for seeking help, thank you!