public static GrayscaleRenderOptions FromMinMax(DicomDataset dataset) { var bits = BitDepth.FromDataset(dataset); var options = new GrayscaleRenderOptions(bits); options.RescaleSlope = dataset.Get<double>(DicomTag.RescaleSlope, 1.0); options.RescaleIntercept = dataset.Get<double>(DicomTag.RescaleIntercept, 0.0); int padding = dataset.Get<int>(DicomTag.PixelPaddingValue, 0, Int32.MinValue); var transcoder = new DicomTranscoder(dataset.InternalTransferSyntax, DicomTransferSyntax.ExplicitVRLittleEndian); var pixels = transcoder.DecodePixelData(dataset, 0); var range = pixels.GetMinMax(padding); if (range.Minimum < bits.MinimumValue || range.Minimum == Double.MaxValue) range.Minimum = bits.MinimumValue; if (range.Maximum > bits.MaximumValue || range.Maximum == Double.MinValue) range.Maximum = bits.MaximumValue; options.WindowWidth = Math.Abs(range.Maximum - range.Minimum); options.WindowCenter = range.Minimum + (options.WindowWidth / 2.0); options.VOILUTFunction = dataset.Get<string>(DicomTag.VOILUTFunction, "LINEAR"); options.Monochrome1 = dataset.Get<PhotometricInterpretation>(DicomTag.PhotometricInterpretation) == PhotometricInterpretation.Monochrome1; return options; }
/// <summary> /// Create grayscale render options based on pixel data histogram. /// </summary> /// <param name="dataset">DICOM dataset from which render options should be obtained.</param> /// <param name="percent">Percentage of histogram window to include.</param> /// <returns>Grayscale render options based on pixel data histogram.</returns> public static GrayscaleRenderOptions FromHistogram(DicomDataset dataset, int percent = 90) { var bits = BitDepth.FromDataset(dataset); var options = new GrayscaleRenderOptions(bits); options.RescaleSlope = dataset.Get<double>(DicomTag.RescaleSlope, 1.0); options.RescaleIntercept = dataset.Get<double>(DicomTag.RescaleIntercept, 0.0); var transcoder = new DicomTranscoder( dataset.InternalTransferSyntax, DicomTransferSyntax.ExplicitVRLittleEndian); var pixels = transcoder.DecodePixelData(dataset, 0); var histogram = pixels.GetHistogram(0); int padding = dataset.Get<int>(DicomTag.PixelPaddingValue, 0, Int32.MinValue); if (padding != Int32.MinValue) histogram.Clear(padding); histogram.ApplyWindow(percent); var min = histogram.WindowStart * options.RescaleSlope + options.RescaleIntercept; var max = histogram.WindowEnd * options.RescaleSlope + options.RescaleIntercept; options.WindowWidth = Math.Abs(max - min); options.WindowCenter = (max + min) / 2.0; options.VOILUTFunction = dataset.Get<string>(DicomTag.VOILUTFunction, "LINEAR"); options.ColorMap = GetColorMap(dataset); return options; }
public static GrayscaleRenderOptions FromHistogram(DicomDataset dataset, int percent = 90) { var bits = BitDepth.FromDataset(dataset); var options = new GrayscaleRenderOptions(bits); options.RescaleSlope = dataset.Get<double>(DicomTag.RescaleSlope, 1.0); options.RescaleIntercept = dataset.Get<double>(DicomTag.RescaleIntercept, 0.0); var transcoder = new DicomTranscoder(dataset.InternalTransferSyntax, DicomTransferSyntax.ExplicitVRLittleEndian); var pixels = transcoder.DecodePixelData(dataset, 0); var histogram = pixels.GetHistogram(0); int padding = dataset.Get<int>(DicomTag.PixelPaddingValue, 0, Int32.MinValue); if (padding != Int32.MinValue) histogram.Clear(padding); histogram.ApplyWindow(percent); options.WindowWidth = histogram.WindowEnd - histogram.WindowStart; options.WindowCenter = histogram.WindowStart + (options.WindowWidth / 2.0); options.VOILUTFunction = dataset.Get<string>(DicomTag.VOILUTFunction, "LINEAR"); options.Monochrome1 = dataset.Get<PhotometricInterpretation>(DicomTag.PhotometricInterpretation) == PhotometricInterpretation.Monochrome1; return options; }
/// <summary> /// Create grayscale render options based on identified minimum and maximum pixel values. /// </summary> /// <param name="dataset">DICOM dataset from which render options should be obtained.</param> /// <returns>Grayscale render options based on identified minimum and maximum pixel values.</returns> public static GrayscaleRenderOptions FromMinMax(DicomDataset dataset) { var bits = BitDepth.FromDataset(dataset); var options = new GrayscaleRenderOptions(bits); options.RescaleSlope = dataset.Get<double>(DicomTag.RescaleSlope, 1.0); options.RescaleIntercept = dataset.Get<double>(DicomTag.RescaleIntercept, 0.0); int padding = dataset.Get<int>(DicomTag.PixelPaddingValue, 0, Int32.MinValue); var transcoder = new DicomTranscoder( dataset.InternalTransferSyntax, DicomTransferSyntax.ExplicitVRLittleEndian); var pixels = transcoder.DecodePixelData(dataset, 0); var range = pixels.GetMinMax(padding); if (range.Minimum < bits.MinimumValue || range.Minimum == Double.MaxValue) range.Minimum = bits.MinimumValue; if (range.Maximum > bits.MaximumValue || range.Maximum == Double.MinValue) range.Maximum = bits.MaximumValue; var min = range.Minimum * options.RescaleSlope + options.RescaleIntercept; var max = range.Maximum * options.RescaleSlope + options.RescaleIntercept; options.WindowWidth = Math.Abs(max - min); options.WindowCenter = (max + min) / 2.0; options.VOILUTFunction = dataset.Get<string>(DicomTag.VOILUTFunction, "LINEAR"); options.ColorMap = GetColorMap(dataset); return options; }