public DNGFile(string aFileName)
: base(aFileName)
{
byte a = mFileReader.ReadByte();
byte b = mFileReader.ReadByte();
bool fileIsLittleEndian;
if (a == b && b == 'I')
{
fileIsLittleEndian = true;
}
else
if (a == b && b == 'M')
{
fileIsLittleEndian = false;
}
else
{
throw new FileLoadException("Could not determine file endianess. Is this a proper TIFF/PEF file?", aFileName);
}
mEndianSwap = fileIsLittleEndian != BitConverter.IsLittleEndian;
ushort magicNumber = ReadUI2();
if (magicNumber != 42)
{
throw new FileLoadException("This is not a valid TIFF/PEF file: Magic number is not 42.", aFileName);
}
uint offsetToFirstIFD = ReadUI4();
mFile.Seek(offsetToFirstIFD, SeekOrigin.Begin);
List <ImageFileDirectory> ifds = new List <ImageFileDirectory>();
while (true)
{
ImageFileDirectory ifd = new ImageFileDirectory(this);
ifds.Add(ifd);
uint offsetToNext = ReadUI4();
if (offsetToNext == 0)
{
break;
}
Seek(offsetToNext, System.IO.SeekOrigin.Begin);
}
IFDDNGVersion version = ifds[0].GetEntry <IFDDNGVersion>();
isDNGVersionLarge = version == null;
if (version != null)
{
isDNGVersionLarge = version.Value[1] > 1;
}
IFDSubIFDs sub = ifds[0].GetEntry <IFDSubIFDs>();
mWidth = (int)sub.Value[0].GetEntry <IFDImageWidth>().Value;
mHeight = (int)sub.Value[0].GetEntry <IFDImageLength>().Value;
IFDExif exif = ifds[0].GetEntry <IFDExif>();
mISO = exif.GetEntry <ExifEntry.ExifISOSpeedRatings>().Value;
mBitDepth = sub.Value[0].GetEntry <IFDBitsPerSample>().Value;
mBayerPattern = sub.Value[0].GetEntry <IFDCFAPattern>().Value;
mBayerWidth = sub.Value[0].GetEntry <IFDCFARepeatPatternDim>().Value[0];
mBayerHeight = sub.Value[0].GetEntry <IFDCFARepeatPatternDim>().Value[1];
mWhiteLevel = sub.Value[0].GetEntry <IFDDNGWhiteLevel>();
mBlackLevel = sub.Value[0].GetEntry <IFDDNGBlackLevel>();
mAsShotNeutral = sub.Value[0].GetEntry <IFDDNGAsShotNeutral>();
mLinearizationTable = sub.Value[0].GetEntry <IFDDNGLinearizationTable>();
mMinVal = new float[3];
mMinVal[0] = (float)mBlackLevel.Value[0].Value;
mMinVal[1] = (float)mBlackLevel.Value[mBlackLevel.Value.Length > 1 ? 1 : 0].Value;
mMinVal[2] = (float)mBlackLevel.Value[mBlackLevel.Value.Length > 3 ? 3 : 0].Value;
mMaxVal = new float[3];
mMaxVal[0] = (float)mWhiteLevel.Value[0];
mMaxVal[1] = (float)mWhiteLevel.Value[mWhiteLevel.Value.Length > 1 ? 1 : 0];
mMaxVal[2] = (float)mWhiteLevel.Value[mWhiteLevel.Value.Length > 3 ? 3 : 0];
int tileWidth = (int)sub.Value[0].GetEntry <IFDTileWidth>().Value;
int tileHeight = (int)sub.Value[0].GetEntry <IFDTileLength>().Value;
uint[] offsets = sub.Value[0].GetEntry <IFDTileOffsets>().Value;
uint[] byteCounts = sub.Value[0].GetEntry <IFDTileByteCounts>().Value;
mRawImage = new ushort[mHeight, mWidth];
int row = 0;
int col = 0;
for (int tile = 0; tile < offsets.Length; tile++)
{
byte[] data;
Seek(offsets[tile], SeekOrigin.Begin);
data = mFileReader.ReadBytes((int)byteCounts[tile]);
MemoryStream ms = new MemoryStream(data);
ms.Seek(0, SeekOrigin.Begin);
jhead jh = new jhead();
int ret = ljpeg_start(jh, 0, ms);
int jrow, jcol;
if (ret > 0 && jh != null)
{
int jwide = jh.wide;
jwide *= jh.clrs;
unsafe
{
if (jh.algo == 0xc3) //lossless JPEG
{
for (jrow = 0; jrow < jh.high; jrow++)
{
fixed(ushort *ptr = jh.row)
{
jh.rowPtr = ptr;
ushort *rp = ljpeg_row(jrow, jh, ms);
for (jcol = 0; jcol < jwide; jcol++)
{
if (jcol + col < mWidth && jrow + row < mHeight)
{
if (mLinearizationTable != null)
{
mRawImage[row + jrow, col + jcol] = mLinearizationTable.Value[rp[jcol] < mLinearizationTable.Value.Length ? rp[jcol] : mLinearizationTable.Value.Length - 1];
}
else
{
mRawImage[row + jrow, col + jcol] = rp[jcol];
}
}
}
jh.rowPtr = null;
}
}
}
}
}
col += tileWidth;
if (col > mWidth)
{
col = 0;
row += tileHeight;
}
}
}
public DNGFile(string aFileName, bool HeaderOnly)
: base(aFileName)
{
byte a = mFileReader.ReadByte();
byte b = mFileReader.ReadByte();
bool fileIsLittleEndian;
if (a == b && b == 'I')
{
fileIsLittleEndian = true;
}
else
if (a == b && b == 'M')
{
fileIsLittleEndian = false;
}
else
{
throw new FileLoadException("Could not determine file endianess. Is this a proper TIFF/PEF file?", aFileName);
}
mEndianSwap = fileIsLittleEndian != BitConverter.IsLittleEndian;
ushort magicNumber = ReadUI2();
if (magicNumber != 42)
{
throw new FileLoadException("This is not a valid TIFF/PEF file: Magic number is not 42.", aFileName);
}
uint offsetToFirstIFD = ReadUI4();
mFile.Seek(offsetToFirstIFD, SeekOrigin.Begin);
List <ImageFileDirectory> ifds = new List <ImageFileDirectory>();
while (true)
{
ImageFileDirectory ifd = new ImageFileDirectory(this);
ifds.Add(ifd);
uint offsetToNext = ReadUI4();
if (offsetToNext == 0)
{
break;
}
Seek(offsetToNext, System.IO.SeekOrigin.Begin);
}
IFDSubIFDs sub = ifds[0].GetEntry <IFDSubIFDs>();
mWidth = (int)sub.Value[0].GetEntry <IFDImageWidth>().Value;
mHeight = (int)sub.Value[0].GetEntry <IFDImageLength>().Value;
IFDExif exif = ifds[0].GetEntry <IFDExif>();
mISO = exif.GetEntry <ExifEntry.ExifISOSpeedRatings>().Value;
mBitDepth = sub.Value[0].GetEntry <IFDBitsPerSample>().Value;
mBayerPattern = sub.Value[0].GetEntry <IFDCFAPattern>().Value;
mBayerWidth = sub.Value[0].GetEntry <IFDCFARepeatPatternDim>().Value[0];
mBayerHeight = sub.Value[0].GetEntry <IFDCFARepeatPatternDim>().Value[1];
mWhiteLevel = sub.Value[0].GetEntry <IFDDNGWhiteLevel>();
mBlackLevel = sub.Value[0].GetEntry <IFDDNGBlackLevel>();
mLinearizationTable = sub.Value[0].GetEntry <IFDDNGLinearizationTable>();
mMinVal = new float[3];
mMinVal[0] = (float)mBlackLevel.Value[0].Value;
mMinVal[1] = (float)mBlackLevel.Value[mBlackLevel.Value.Length > 1 ? 1 : 0].Value;
mMinVal[2] = (float)mBlackLevel.Value[mBlackLevel.Value.Length > 3 ? 3 : 0].Value;
mMaxVal = new float[3];
mMaxVal[0] = (float)mWhiteLevel.Value[0];
mMaxVal[1] = (float)mWhiteLevel.Value[mWhiteLevel.Value.Length > 1 ? 1 : 0];
mMaxVal[2] = (float)mWhiteLevel.Value[mWhiteLevel.Value.Length > 3 ? 3 : 0];
mAsShotNeutral = ifds[0].GetEntry <IFDDNGAsShotNeutral>();
}
private ImageFileDirectory SetMembers()
{
ImageFileDirectory rawIFD = null;
//find the IFD with the RAW Bayer image
for (int i = 0; i < mIfds.Count; i++)
{
//well, it should actually be somewhere in IFD0...
if (mIfds[i].GetEntry <IFDPhotometricInterpretation>() != null)
{
if (mIfds[i].GetEntry <IFDPhotometricInterpretation>().Value ==
IFDPhotometricInterpretation.PhotometricInterpretation.CFA)
{
rawIFD = mIfds[i];
break;
}
}
}
//no root IFD seems to contain RAW bayer, search for Sub-IFDs
if (rawIFD == null)
{
//find the IFD with the RAW Bayer image
for (int i = 0; i < mIfds.Count; i++)
{
IFDSubIFDs subIFD = mIfds[i].GetEntry <IFDSubIFDs>();
if (subIFD == null)
{
continue;
}
for (int j = 0; j < subIFD.Value.Count; j++)
{
if (subIFD.Value[j].GetEntry <IFDPhotometricInterpretation>().Value ==
IFDPhotometricInterpretation.PhotometricInterpretation.CFA)
{
rawIFD = subIFD.Value[j];
break;
}
}
}
}
if (rawIFD == null)
{
throw new ArgumentException("Can't find IFD with Bayer RAW image.");
}
mLinearizationTable = rawIFD.GetEntry <IFDDNGLinearizationTable>();
mWidth = (int)rawIFD.GetEntry <IFDImageWidth>().Value;
mHeight = (int)rawIFD.GetEntry <IFDImageLength>().Value;
if (mIfds[0].GetEntry <IFDDateTime>() != null)
{
mRecordingDate = mIfds[0].GetEntry <IFDDateTime>().Value;
}
//in case of Pentax this will have succes:
IFDDNGPrivateData privateData = mIfds[0].GetEntry <IFDDNGPrivateData>();
if (privateData != null)
{
MNLevelInfo levelInfo = privateData.PentaxMakerNotes.GetEntry <MNLevelInfo>();
if (levelInfo != null)
{
mRollAngle = levelInfo.Value.RollAngle;
mRollAnglePresent = true;
}
}
IFDExif exif = mIfds[0].GetEntry <IFDExif>();
if (exif != null)
{
mISO = exif.GetEntry <ExifISOSpeedRatings>().Value;
mExposureTime = exif.GetEntry <ExifExposureTime>().Value;
mRecordingDate = (exif.GetEntry <ExifDateTimeDigitized>()?.Value).GetValueOrDefault(mRecordingDate);
}
else if (rawIFD.GetEntry <IFDISOSpeedRatings>() != null)
{
mISO = rawIFD.GetEntry <IFDISOSpeedRatings>().Value;
//very likely that exposure time is also present
mExposureTime = rawIFD.GetEntry <IFDExposureTime>().Value;
}
else if (mIfds[0].GetEntry <IFDISOSpeedRatings>() != null)
{
mISO = mIfds[0].GetEntry <IFDISOSpeedRatings>().Value;
//very likely that exposure time is also present
mExposureTime = mIfds[0].GetEntry <IFDExposureTime>().Value;
}
mBitDepth = rawIFD.GetEntry <IFDBitsPerSample>().Value[0];
int bayerWidth = rawIFD.GetEntry <IFDCFARepeatPatternDim>().Value[0];
int bayerHeight = rawIFD.GetEntry <IFDCFARepeatPatternDim>().Value[1];
if (bayerHeight != 2 || bayerWidth != 2)
{
throw new ArgumentException("This file has a bayer pattern size different than 2x2. Can't decode that.");
}
ExifCFAPattern.BayerColor[] bayer = rawIFD.GetEntry <IFDCFAPattern>().Value;
mBayerPattern = new BayerColor[bayer.Length];
for (int i = 0; i < bayer.Length; i++)
{
mBayerPattern[i] = (BayerColor)(int)bayer[i];
}
IFDDNGCFAPlaneColor planeColor = rawIFD.GetEntry <IFDDNGCFAPlaneColor>();
int[] planeColorHelper = new int[] { 0, 1, 2 };
if (planeColor != null) //did it ever differ from 0,1,2?
{
//0 = red, 1 = gree, 2 = blue.
//The debayer algo creates images with plane order red/green/blue.
//If this order differs, we need to re-order the planes in order to
//have the color matrices correct.
//reset colorTwist matrix:
mColorTwist = new float[3, 4];
if (planeColor.Value.Length > 3 || planeColor.Value.Length < 3)
{
throw new ArgumentException("This image doesn't contain three color planes.");
}
for (int i = 0; i < planeColor.Value.Length; i++)
{
int color = planeColor.Value[i];
planeColorHelper[i] = color;
if (color > 2)
{
throw new ArgumentException("This image contains colors different than red/green/blue.");
}
mColorTwist[color, i] = 1;
if (color != i)
{
mColorTwistIsIdentity = false;
}
}
}
mColorSpec = new DNGColorSpec(3, mIfds[0], rawIFD);
IFDDNGWhiteLevel whiteLevel = rawIFD.GetEntry <IFDDNGWhiteLevel>();
IFDDNGBlackLevel blackLevel = rawIFD.GetEntry <IFDDNGBlackLevel>();
mBlackLevel = new float[3];
if (blackLevel != null)
{
//only one value for all colors
if (blackLevel.Value.Length == 1)
{
mBlackLevel[0] = (float)blackLevel.Value[0].Value;
mBlackLevel[1] = (float)blackLevel.Value[0].Value;
mBlackLevel[2] = (float)blackLevel.Value[0].Value;
}
//values per color channel
if (blackLevel.Value.Length == 3)
{
mBlackLevel[planeColorHelper[0]] = (float)blackLevel.Value[0].Value;
mBlackLevel[planeColorHelper[1]] = (float)blackLevel.Value[1].Value;
mBlackLevel[planeColorHelper[2]] = (float)blackLevel.Value[2].Value;
}
//values per color bayer pattern
if (blackLevel.Value.Length == 4)
{
//red
int indexR = -1;
for (int i = 0; i < mBayerPattern.Length; i++)
{
if (mBayerPattern[i] == BayerColor.Red)
{
indexR = i;
break;
}
}
mBlackLevel[0] = (float)blackLevel.Value[indexR].Value;
//blue
int indexB = -1;
for (int i = 0; i < mBayerPattern.Length; i++)
{
if (mBayerPattern[i] == BayerColor.Blue)
{
indexB = i;
break;
}
}
mBlackLevel[2] = (float)blackLevel.Value[indexB].Value;
//green, the two remaining indices
int indexG1 = -1, indexG2 = -1;
for (int i = 0; i < mBayerPattern.Length; i++)
{
if (mBayerPattern[i] == BayerColor.Green && indexG1 == -1)
{
indexG1 = i;
}
if (mBayerPattern[i] == BayerColor.Green && indexG1 != -1)
{
indexG2 = i;
}
}
float g1 = (float)blackLevel.Value[indexG1].Value;
float g2 = (float)blackLevel.Value[indexG2].Value;
mBlackLevel[1] = Math.Max(g1, g2); //well, one could distinguish the two greens, but what for?
}
}
mWhiteLevel = new float[3];
mWhiteLevel[0] = (float)whiteLevel.Value[0];
mWhiteLevel[1] = (float)whiteLevel.Value[0];
mWhiteLevel[2] = (float)whiteLevel.Value[0];
//subtract black level from white level
mWhiteLevel[0] -= mBlackLevel[0];
mWhiteLevel[1] -= mBlackLevel[1];
mWhiteLevel[2] -= mBlackLevel[2];
//get white balance from color spec
DNGVector wb = mColorSpec.CameraWhite;
mWhiteBalance = new float[3];
mWhiteBalance[planeColorHelper[0]] = 1.0f / (float)wb[0];
mWhiteBalance[planeColorHelper[1]] = 1.0f / (float)wb[1];
mWhiteBalance[planeColorHelper[2]] = 1.0f / (float)wb[2];
//look for orientation tag. If RAW ifd has the tag, choose that one
if (rawIFD.GetEntry <IFDOrientation>() != null)
{
mOrientation = new DNGOrientation(rawIFD.GetEntry <IFDOrientation>().Value);
}
else if (mIfds[0].GetEntry <IFDOrientation>() != null)
{
mOrientation = new DNGOrientation(mIfds[0].GetEntry <IFDOrientation>().Value);
}
else
{
//no tag found, use default
mOrientation = new DNGOrientation(DNGOrientation.Orientation.Normal);
}
//default Values:
int cropLeft = 0;
int cropTop = 0;
int croppedWidth = mWidth;
int croppedHeight = mHeight;
IFDDNGActiveArea activeArea = rawIFD.GetEntry <IFDDNGActiveArea>();
//if active area is defined:
if (activeArea != null)
{
int top, left, bottom, right;
top = (int)activeArea.Value[0];
left = (int)activeArea.Value[1];
bottom = (int)activeArea.Value[2];
right = (int)activeArea.Value[3];
cropLeft += left;
cropTop += top;
croppedWidth = right - left;
croppedHeight = bottom - top;
//CFA pattern is defined on active area. If top/left is uneven we need to shift the CFA pattern accordingly
if (top % 2 != 0)
{
BayerColor bayer0 = BayerPattern[0];
BayerColor bayer1 = BayerPattern[1];
BayerPattern[0] = BayerPattern[2];
BayerPattern[1] = BayerPattern[3];
BayerPattern[2] = bayer0;
BayerPattern[3] = bayer1;
}
if (left % 2 != 0)
{
BayerColor bayer0 = BayerPattern[0];
BayerColor bayer2 = BayerPattern[2];
BayerPattern[0] = BayerPattern[1];
BayerPattern[2] = BayerPattern[3];
BayerPattern[1] = bayer0;
BayerPattern[3] = bayer2;
}
}
IFDDNGDefaultCropOrigin cropOrigin = rawIFD.GetEntry <IFDDNGDefaultCropOrigin>();
IFDDNGDefaultCropSize cropSize = rawIFD.GetEntry <IFDDNGDefaultCropSize>();
if (cropOrigin != null && cropSize != null)
{
int top, left, width, height;
left = (int)(cropOrigin.Value[0].Value);
top = (int)(cropOrigin.Value[0].Value);
width = (int)(cropSize.Value[0].Value);
height = (int)(cropSize.Value[0].Value);
cropLeft += left;
cropTop += top;
croppedWidth = width;
croppedHeight = height;
}
//we always crop at least two pixels because of our algos...
mCropLeft = Math.Max(2, cropLeft);
mCropTop = Math.Max(2, cropTop);
mCroppedWidth = croppedWidth - Math.Max(0, (cropLeft + croppedWidth) - (mWidth - 2));
mCroppedHeight = croppedHeight - Math.Max(0, (cropTop + croppedHeight) - (mHeight - 2));
IFDDNGNoiseProfile noise = rawIFD.GetEntry <IFDDNGNoiseProfile>();
if (noise == null)
{
noise = mIfds[0].GetEntry <IFDDNGNoiseProfile>();
}
if (noise != null)
{
//if noise level is given for all channels,
//take the green one as it is usually scalled to one
if (noise.Value.Length > 2)
{
mNoiseModelAlpha = (float)noise.Value[planeColorHelper[1] * 2];
mNoiseModelBeta = (float)noise.Value[planeColorHelper[1] * 2 + 1];
}
else
{
mNoiseModelAlpha = (float)noise.Value[0];
mNoiseModelBeta = (float)noise.Value[1];
}
}
mMake = mIfds[0].GetEntry <IFDMake>().Value;
mUniqueModelName = mIfds[0].GetEntry <IFDDNGUniqueCameraModel>().Value;
return(rawIFD);
}
