diff --git a/.editorconfig b/.editorconfig
index c28089d720..af1e5b44c1 100644
--- a/.editorconfig
+++ b/.editorconfig
@@ -104,8 +104,8 @@ dotnet_style_parentheses_in_relational_binary_operators = always_for_clarity:war
dotnet_style_parentheses_in_other_binary_operators = always_for_clarity:warning
dotnet_style_parentheses_in_other_operators = always_for_clarity:suggestion
# Expression-level preferences
-dotnet_style_object_initializer = true:warning
-dotnet_style_collection_initializer = true:warning
+dotnet_style_object_initializer = true:error
+dotnet_style_collection_initializer = true:error
dotnet_style_explicit_tuple_names = true:warning
dotnet_style_prefer_inferred_tuple_names = true:warning
dotnet_style_prefer_inferred_anonymous_type_member_names = true:warning
@@ -135,9 +135,9 @@ csharp_style_prefer_null_check_over_type_check = true:warning
# https://docs.microsoft.com/dotnet/fundamentals/code-analysis/style-rules/language-rules#c-style-rules
[*.{cs,csx,cake}]
# 'var' preferences
-csharp_style_var_for_built_in_types = false:warning
-csharp_style_var_when_type_is_apparent = false:warning
-csharp_style_var_elsewhere = false:warning
+csharp_style_var_for_built_in_types = false:error
+csharp_style_var_when_type_is_apparent = false:error
+csharp_style_var_elsewhere = false:error
# Expression-bodied members
csharp_style_expression_bodied_methods = true:warning
csharp_style_expression_bodied_constructors = true:warning
@@ -160,7 +160,7 @@ csharp_style_pattern_local_over_anonymous_function = true:warning
csharp_style_deconstructed_variable_declaration = true:warning
csharp_style_prefer_index_operator = true:warning
csharp_style_prefer_range_operator = true:warning
-csharp_style_implicit_object_creation_when_type_is_apparent = true:warning
+csharp_style_implicit_object_creation_when_type_is_apparent = true:error
# "Null" checking preferences
csharp_style_throw_expression = true:warning
csharp_style_conditional_delegate_call = true:warning
diff --git a/.gitattributes b/.gitattributes
index b5f742ab47..f7bd4d061e 100644
--- a/.gitattributes
+++ b/.gitattributes
@@ -136,3 +136,10 @@
*.ico filter=lfs diff=lfs merge=lfs -text
*.cur filter=lfs diff=lfs merge=lfs -text
*.ani filter=lfs diff=lfs merge=lfs -text
+*.heic filter=lfs diff=lfs merge=lfs -text
+*.hif filter=lfs diff=lfs merge=lfs -text
+*.avif filter=lfs diff=lfs merge=lfs -text
+###############################################################################
+# Handle ICC files by git lfs
+###############################################################################
+*.icc filter=lfs diff=lfs merge=lfs -text
diff --git a/shared-infrastructure b/shared-infrastructure
index 1dbfb576c8..5e13cde851 160000
--- a/shared-infrastructure
+++ b/shared-infrastructure
@@ -1 +1 @@
-Subproject commit 1dbfb576c83507645265c79e03369b66cdc0379f
+Subproject commit 5e13cde851a3d6e95d0dfdde2a57071f1efda9c3
diff --git a/src/ImageSharp/ColorProfiles/CieLab.cs b/src/ImageSharp/ColorProfiles/CieLab.cs
index 377cc20a99..ca72dd745a 100644
--- a/src/ImageSharp/ColorProfiles/CieLab.cs
+++ b/src/ImageSharp/ColorProfiles/CieLab.cs
@@ -35,7 +35,6 @@ public CieLab(float l, float a, float b)
/// The vector representing the l, a, b components.
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public CieLab(Vector3 vector)
- : this()
{
this.L = vector.X;
this.A = vector.Y;
@@ -82,6 +81,49 @@ public CieLab(Vector3 vector)
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static bool operator !=(CieLab left, CieLab right) => !left.Equals(right);
+ ///
+ public Vector4 ToScaledVector4()
+ {
+ Vector3 v3 = default;
+ v3 += this.AsVector3Unsafe();
+ v3 += new Vector3(0, 128F, 128F);
+ v3 /= new Vector3(100F, 255F, 255F);
+ return new Vector4(v3, 1F);
+ }
+
+ ///
+ public static CieLab FromScaledVector4(Vector4 source)
+ {
+ Vector3 v3 = source.AsVector3();
+ v3 *= new Vector3(100F, 255, 255);
+ v3 -= new Vector3(0, 128F, 128F);
+ return new CieLab(v3);
+ }
+
+ ///
+ public static void ToScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = source[i].ToScaledVector4();
+ }
+ }
+
+ ///
+ public static void FromScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = FromScaledVector4(source[i]);
+ }
+ }
+
///
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static CieLab FromProfileConnectingSpace(ColorConversionOptions options, in CieXyz source)
@@ -136,7 +178,7 @@ public CieXyz ToProfileConnectingSpace(ColorConversionOptions options)
float yr = l > CieConstants.Kappa * CieConstants.Epsilon ? Numerics.Pow3((l + 16F) / 116F) : l / CieConstants.Kappa;
float zr = fz3 > CieConstants.Epsilon ? fz3 : ((116F * fz) - 16F) / CieConstants.Kappa;
- CieXyz whitePoint = options.WhitePoint;
+ CieXyz whitePoint = options.SourceWhitePoint;
Vector3 wxyz = new(whitePoint.X, whitePoint.Y, whitePoint.Z);
Vector3 xyzr = new(xr, yr, zr);
diff --git a/src/ImageSharp/ColorProfiles/CieLch.cs b/src/ImageSharp/ColorProfiles/CieLch.cs
index 1319783406..e62aa2ba23 100644
--- a/src/ImageSharp/ColorProfiles/CieLch.cs
+++ b/src/ImageSharp/ColorProfiles/CieLch.cs
@@ -42,6 +42,17 @@ public CieLch(Vector3 vector)
this.H = vector.Z;
}
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+#pragma warning disable SA1313 // Parameter names should begin with lower-case letter
+ private CieLch(Vector3 vector, bool _)
+#pragma warning restore SA1313 // Parameter names should begin with lower-case letter
+ {
+ vector = Vector3.Clamp(vector, Min, Max);
+ this.L = vector.X;
+ this.C = vector.Y;
+ this.H = vector.Z;
+ }
+
///
/// Gets the lightness dimension.
/// A value ranging between 0 (black), 100 (diffuse white) or higher (specular white).
@@ -50,7 +61,7 @@ public CieLch(Vector3 vector)
///
/// Gets the a chroma component.
- /// A value ranging from 0 to 200.
+ /// A value ranging from -200 to 200.
///
public float C { get; }
@@ -82,6 +93,49 @@ public CieLch(Vector3 vector)
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static bool operator !=(CieLch left, CieLch right) => !left.Equals(right);
+ ///
+ public Vector4 ToScaledVector4()
+ {
+ Vector3 v3 = default;
+ v3 += this.AsVector3Unsafe();
+ v3 += new Vector3(0, 200, 0);
+ v3 /= new Vector3(100, 400, 360);
+ return new Vector4(v3, 1F);
+ }
+
+ ///
+ public static CieLch FromScaledVector4(Vector4 source)
+ {
+ Vector3 v3 = source.AsVector3();
+ v3 *= new Vector3(100, 400, 360);
+ v3 -= new Vector3(0, 200, 0);
+ return new CieLch(v3, true);
+ }
+
+ ///
+ public static void ToScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = source[i].ToScaledVector4();
+ }
+ }
+
+ ///
+ public static void FromScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = FromScaledVector4(source[i]);
+ }
+ }
+
///
public static CieLch FromProfileConnectingSpace(ColorConversionOptions options, in CieLab source)
{
diff --git a/src/ImageSharp/ColorProfiles/CieLchuv.cs b/src/ImageSharp/ColorProfiles/CieLchuv.cs
index 7fd95feb19..5478752ddc 100644
--- a/src/ImageSharp/ColorProfiles/CieLchuv.cs
+++ b/src/ImageSharp/ColorProfiles/CieLchuv.cs
@@ -35,7 +35,6 @@ public CieLchuv(float l, float c, float h)
/// The vector representing the l, c, h components.
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public CieLchuv(Vector3 vector)
- : this()
{
vector = Vector3.Clamp(vector, Min, Max);
this.L = vector.X;
@@ -43,6 +42,16 @@ public CieLchuv(Vector3 vector)
this.H = vector.Z;
}
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+#pragma warning disable SA1313 // Parameter names should begin with lower-case letter
+ private CieLchuv(Vector3 vector, bool _)
+#pragma warning restore SA1313 // Parameter names should begin with lower-case letter
+ {
+ this.L = vector.X;
+ this.C = vector.Y;
+ this.H = vector.Z;
+ }
+
///
/// Gets the lightness dimension.
/// A value ranging between 0 (black), 100 (diffuse white) or higher (specular white).
@@ -51,7 +60,7 @@ public CieLchuv(Vector3 vector)
///
/// Gets the a chroma component.
- /// A value ranging from 0 to 200.
+ /// A value ranging from -200 to 200.
///
public float C { get; }
@@ -81,6 +90,49 @@ public CieLchuv(Vector3 vector)
///
public static bool operator !=(CieLchuv left, CieLchuv right) => !left.Equals(right);
+ ///
+ public Vector4 ToScaledVector4()
+ {
+ Vector3 v3 = default;
+ v3 += this.AsVector3Unsafe();
+ v3 += new Vector3(0, 200, 0);
+ v3 /= new Vector3(100, 400, 360);
+ return new Vector4(v3, 1F);
+ }
+
+ ///
+ public static CieLchuv FromScaledVector4(Vector4 source)
+ {
+ Vector3 v3 = source.AsVector3();
+ v3 *= new Vector3(100, 400, 360);
+ v3 -= new Vector3(0, 200, 0);
+ return new CieLchuv(v3, true);
+ }
+
+ ///
+ public static void ToScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = source[i].ToScaledVector4();
+ }
+ }
+
+ ///
+ public static void FromScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = FromScaledVector4(source[i]);
+ }
+ }
+
///
public static CieLchuv FromProfileConnectingSpace(ColorConversionOptions options, in CieXyz source)
{
diff --git a/src/ImageSharp/ColorProfiles/CieLuv.cs b/src/ImageSharp/ColorProfiles/CieLuv.cs
index 97e2826f74..b17c433313 100644
--- a/src/ImageSharp/ColorProfiles/CieLuv.cs
+++ b/src/ImageSharp/ColorProfiles/CieLuv.cs
@@ -37,7 +37,6 @@ public CieLuv(float l, float u, float v)
/// The vector representing the l, u, v components.
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public CieLuv(Vector3 vector)
- : this()
{
this.L = vector.X;
this.U = vector.Y;
@@ -84,6 +83,18 @@ public CieLuv(Vector3 vector)
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static bool operator !=(CieLuv left, CieLuv right) => !left.Equals(right);
+ ///
+ public Vector4 ToScaledVector4() => throw new NotImplementedException();
+
+ ///
+ public static CieLuv FromScaledVector4(Vector4 source) => throw new NotImplementedException();
+
+ ///
+ public static void ToScaledVector4(ReadOnlySpan source, Span destination) => throw new NotImplementedException();
+
+ ///
+ public static void FromScaledVector4(ReadOnlySpan source, Span destination) => throw new NotImplementedException();
+
///
public static CieLuv FromProfileConnectingSpace(ColorConversionOptions options, in CieXyz source)
{
@@ -143,7 +154,7 @@ public CieXyz ToProfileConnectingSpace(ColorConversionOptions options)
// Use doubles here for accuracy.
// Conversion algorithm described here:
// http://www.brucelindbloom.com/index.html?Eqn_Luv_to_XYZ.html
- CieXyz whitePoint = options.WhitePoint;
+ CieXyz whitePoint = options.SourceWhitePoint;
double l = this.L, u = this.U, v = this.V;
diff --git a/src/ImageSharp/ColorProfiles/CieXyy.cs b/src/ImageSharp/ColorProfiles/CieXyy.cs
index 62873df147..744b6195e9 100644
--- a/src/ImageSharp/ColorProfiles/CieXyy.cs
+++ b/src/ImageSharp/ColorProfiles/CieXyy.cs
@@ -35,7 +35,6 @@ public CieXyy(float x, float y, float yl)
/// The vector representing the x, y, Y components.
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public CieXyy(Vector3 vector)
- : this()
{
// Not clamping as documentation about this space only indicates "usual" ranges
this.X = vector.X;
@@ -83,6 +82,38 @@ public CieXyy(Vector3 vector)
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static bool operator !=(CieXyy left, CieXyy right) => !left.Equals(right);
+ ///
+ public Vector4 ToScaledVector4()
+ => new(this.AsVector3Unsafe(), 1F);
+
+ ///
+ public static CieXyy FromScaledVector4(Vector4 source)
+ => new(source.AsVector3());
+
+ ///
+ public static void ToScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = source[i].ToScaledVector4();
+ }
+ }
+
+ ///
+ public static void FromScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = FromScaledVector4(source[i]);
+ }
+ }
+
///
public static CieXyy FromProfileConnectingSpace(ColorConversionOptions options, in CieXyz source)
{
diff --git a/src/ImageSharp/ColorProfiles/CieXyz.cs b/src/ImageSharp/ColorProfiles/CieXyz.cs
index 07f9b47f9b..94fcfb21bb 100644
--- a/src/ImageSharp/ColorProfiles/CieXyz.cs
+++ b/src/ImageSharp/ColorProfiles/CieXyz.cs
@@ -34,7 +34,6 @@ public CieXyz(float x, float y, float z)
///
/// The vector representing the x, y, z components.
public CieXyz(Vector3 vector)
- : this()
{
this.X = vector.X;
this.Y = vector.Y;
@@ -81,12 +80,85 @@ public CieXyz(Vector3 vector)
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static bool operator !=(CieXyz left, CieXyz right) => !left.Equals(right);
- ///
- /// Returns a new representing this instance.
- ///
- /// The .
[MethodImpl(MethodImplOptions.AggressiveInlining)]
- public Vector3 ToVector3() => new(this.X, this.Y, this.Z);
+ internal Vector3 ToVector3() => new(this.X, this.Y, this.Z);
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ internal Vector4 ToVector4()
+ {
+ Vector3 v3 = default;
+ v3 += this.AsVector3Unsafe();
+ return new Vector4(v3, 1F);
+ }
+
+ ///
+ public Vector4 ToScaledVector4()
+ {
+ Vector3 v3 = default;
+ v3 += this.AsVector3Unsafe();
+ v3 *= 32768F / 65535;
+ return new Vector4(v3, 1F);
+ }
+
+ internal static CieXyz FromVector4(Vector4 source)
+ {
+ Vector3 v3 = source.AsVector3();
+ return new CieXyz(v3);
+ }
+
+ ///
+ public static CieXyz FromScaledVector4(Vector4 source)
+ {
+ Vector3 v3 = source.AsVector3();
+ v3 *= 65535 / 32768F;
+ return new CieXyz(v3);
+ }
+
+ ///
+ public static void ToScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = source[i].ToScaledVector4();
+ }
+ }
+
+ ///
+ public static void FromScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = FromScaledVector4(source[i]);
+ }
+ }
+
+ internal static void FromVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = FromVector4(source[i]);
+ }
+ }
+
+ internal static void ToVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = source[i].ToVector4();
+ }
+ }
///
public static CieXyz FromProfileConnectingSpace(ColorConversionOptions options, in CieXyz source)
@@ -127,5 +199,5 @@ public static void ToProfileConnectionSpace(ColorConversionOptions options, Read
public bool Equals(CieXyz other)
=> this.AsVector3Unsafe() == other.AsVector3Unsafe();
- private Vector3 AsVector3Unsafe() => Unsafe.As(ref Unsafe.AsRef(in this));
+ internal Vector3 AsVector3Unsafe() => Unsafe.As(ref Unsafe.AsRef(in this));
}
diff --git a/src/ImageSharp/ColorProfiles/Cmyk.cs b/src/ImageSharp/ColorProfiles/Cmyk.cs
index e924904497..ee81ff9f7e 100644
--- a/src/ImageSharp/ColorProfiles/Cmyk.cs
+++ b/src/ImageSharp/ColorProfiles/Cmyk.cs
@@ -9,6 +9,7 @@ namespace SixLabors.ImageSharp.ColorProfiles;
///
/// Represents an CMYK (cyan, magenta, yellow, keyline) color.
+///
///
[StructLayout(LayoutKind.Sequential)]
public readonly struct Cmyk : IColorProfile
@@ -36,7 +37,18 @@ public Cmyk(float c, float m, float y, float k)
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public Cmyk(Vector4 vector)
{
- vector = Numerics.Clamp(vector, Min, Max);
+ vector = Vector4.Clamp(vector, Min, Max);
+ this.C = vector.X;
+ this.M = vector.Y;
+ this.Y = vector.Z;
+ this.K = vector.W;
+ }
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+#pragma warning disable SA1313 // Parameter names should begin with lower-case letter
+ private Cmyk(Vector4 vector, bool _)
+#pragma warning restore SA1313 // Parameter names should begin with lower-case letter
+ {
this.C = vector.X;
this.M = vector.Y;
this.Y = vector.Z;
@@ -89,16 +101,42 @@ public Cmyk(Vector4 vector)
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static bool operator !=(Cmyk left, Cmyk right) => !left.Equals(right);
+ ///
+ public Vector4 ToScaledVector4()
+ {
+ Vector4 v4 = default;
+ v4 += this.AsVector4Unsafe();
+ return v4;
+ }
+
+ ///
+ public static Cmyk FromScaledVector4(Vector4 source)
+ => new(source, true);
+
+ ///
+ public static void ToScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+ MemoryMarshal.Cast(source).CopyTo(destination);
+ }
+
+ ///
+ public static void FromScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+ MemoryMarshal.Cast(source).CopyTo(destination);
+ }
+
///
public static Cmyk FromProfileConnectingSpace(ColorConversionOptions options, in Rgb source)
{
// To CMY
- Vector3 cmy = Vector3.One - source.ToScaledVector3();
+ Vector3 cmy = Vector3.One - source.AsVector3Unsafe();
// To CMYK
Vector3 k = new(MathF.Min(cmy.X, MathF.Min(cmy.Y, cmy.Z)));
- if (MathF.Abs(k.X - 1F) < Constants.Epsilon)
+ if (k.X >= 1F - Constants.Epsilon)
{
return new Cmyk(0, 0, 0, 1F);
}
@@ -124,7 +162,7 @@ public static void FromProfileConnectionSpace(ColorConversionOptions options, Re
///
public Rgb ToProfileConnectingSpace(ColorConversionOptions options)
{
- Vector3 rgb = (Vector3.One - new Vector3(this.C, this.M, this.Y)) * (Vector3.One - new Vector3(this.K));
+ Vector3 rgb = (Vector3.One - new Vector3(this.C, this.M, this.Y)) * (1F - this.K);
return Rgb.FromScaledVector3(rgb);
}
@@ -134,8 +172,7 @@ public static void ToProfileConnectionSpace(ColorConversionOptions options, Read
// TODO: We can possibly optimize this by using SIMD
for (int i = 0; i < source.Length; i++)
{
- Cmyk cmyk = source[i];
- destination[i] = cmyk.ToProfileConnectingSpace(options);
+ destination[i] = source[i].ToProfileConnectingSpace(options);
}
}
diff --git a/src/ImageSharp/ColorProfiles/ColorConversionOptions.cs b/src/ImageSharp/ColorProfiles/ColorConversionOptions.cs
index 1eb118834a..44c7d2ac77 100644
--- a/src/ImageSharp/ColorProfiles/ColorConversionOptions.cs
+++ b/src/ImageSharp/ColorProfiles/ColorConversionOptions.cs
@@ -4,6 +4,7 @@
using System.Numerics;
using SixLabors.ImageSharp.ColorProfiles.WorkingSpaces;
using SixLabors.ImageSharp.Memory;
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
namespace SixLabors.ImageSharp.ColorProfiles;
@@ -13,11 +14,16 @@ namespace SixLabors.ImageSharp.ColorProfiles;
public class ColorConversionOptions
{
private Matrix4x4 adaptationMatrix;
+ private YCbCrMatrix yCbCrMatrix;
///
/// Initializes a new instance of the class.
///
- public ColorConversionOptions() => this.AdaptationMatrix = KnownChromaticAdaptationMatrices.Bradford;
+ public ColorConversionOptions()
+ {
+ this.AdaptationMatrix = KnownChromaticAdaptationMatrices.Bradford;
+ this.YCbCrMatrix = KnownYCbCrMatrices.BT601;
+ }
///
/// Gets the memory allocator.
@@ -27,7 +33,7 @@ public class ColorConversionOptions
///
/// Gets the source white point used for chromatic adaptation in conversions from/to XYZ color space.
///
- public CieXyz WhitePoint { get; init; } = KnownIlluminants.D50;
+ public CieXyz SourceWhitePoint { get; init; } = KnownIlluminants.D50;
///
/// Gets the destination white point used for chromatic adaptation in conversions from/to XYZ color space.
@@ -37,13 +43,36 @@ public class ColorConversionOptions
///
/// Gets the source working space used for companding in conversions from/to XYZ color space.
///
- public RgbWorkingSpace RgbWorkingSpace { get; init; } = KnownRgbWorkingSpaces.SRgb;
+ public RgbWorkingSpace SourceRgbWorkingSpace { get; init; } = KnownRgbWorkingSpaces.SRgb;
///
/// Gets the destination working space used for companding in conversions from/to XYZ color space.
///
public RgbWorkingSpace TargetRgbWorkingSpace { get; init; } = KnownRgbWorkingSpaces.SRgb;
+ ///
+ /// Gets the YCbCr matrix to used to perform conversions from/to RGB.
+ ///
+ public YCbCrMatrix YCbCrMatrix
+ {
+ get => this.yCbCrMatrix;
+ init
+ {
+ this.yCbCrMatrix = value;
+ this.TransposedYCbCrMatrix = value.Transpose();
+ }
+ }
+
+ ///
+ /// Gets the source ICC profile.
+ ///
+ public IccProfile? SourceIccProfile { get; init; }
+
+ ///
+ /// Gets the target ICC profile.
+ ///
+ public IccProfile? TargetIccProfile { get; init; }
+
///
/// Gets the transformation matrix used in conversion to perform chromatic adaptation.
/// for further information. Default is Bradford.
@@ -54,10 +83,12 @@ public Matrix4x4 AdaptationMatrix
init
{
this.adaptationMatrix = value;
- Matrix4x4.Invert(value, out Matrix4x4 inverted);
+ _ = Matrix4x4.Invert(value, out Matrix4x4 inverted);
this.InverseAdaptationMatrix = inverted;
}
}
+ internal YCbCrMatrix TransposedYCbCrMatrix { get; private set; }
+
internal Matrix4x4 InverseAdaptationMatrix { get; private set; }
}
diff --git a/src/ImageSharp/ColorProfiles/ColorProfileConverter.cs b/src/ImageSharp/ColorProfiles/ColorProfileConverter.cs
index 18b90a622a..d0afec4ab2 100644
--- a/src/ImageSharp/ColorProfiles/ColorProfileConverter.cs
+++ b/src/ImageSharp/ColorProfiles/ColorProfileConverter.cs
@@ -12,7 +12,7 @@ public class ColorProfileConverter
/// Initializes a new instance of the class.
///
public ColorProfileConverter()
- : this(new())
+ : this(new())
{
}
@@ -33,8 +33,8 @@ public ColorProfileConverter(ColorConversionOptions options)
where TTo : struct, IColorProfile
{
CieXyz sourceWhitePoint = TFrom.GetChromaticAdaptionWhitePointSource() == ChromaticAdaptionWhitePointSource.WhitePoint
- ? this.Options.WhitePoint
- : this.Options.RgbWorkingSpace.WhitePoint;
+ ? this.Options.SourceWhitePoint
+ : this.Options.SourceRgbWorkingSpace.WhitePoint;
CieXyz targetWhitePoint = TTo.GetChromaticAdaptionWhitePointSource() == ChromaticAdaptionWhitePointSource.WhitePoint
? this.Options.TargetWhitePoint
@@ -42,4 +42,7 @@ public ColorProfileConverter(ColorConversionOptions options)
return (sourceWhitePoint, targetWhitePoint);
}
+
+ internal bool ShouldUseIccProfiles()
+ => this.Options.SourceIccProfile != null && this.Options.TargetIccProfile != null;
}
diff --git a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieLabCieLab.cs b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieLabCieLab.cs
index 41ae4b08fa..a2dd5d9ced 100644
--- a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieLabCieLab.cs
+++ b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieLabCieLab.cs
@@ -12,6 +12,11 @@ public static TTo Convert(this ColorProfileConverter converter, in T
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ return converter.ConvertUsingIccProfile(source);
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS
@@ -33,6 +38,12 @@ public static void Convert(this ColorProfileConverter converter, Rea
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ converter.ConvertUsingIccProfile(source, destination);
+ return;
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS.
diff --git a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieLabCieXyz.cs b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieLabCieXyz.cs
index 04937e927e..096622564c 100644
--- a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieLabCieXyz.cs
+++ b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieLabCieXyz.cs
@@ -12,6 +12,11 @@ public static TTo Convert(this ColorProfileConverter converter, in T
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ return converter.ConvertUsingIccProfile(source);
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS
@@ -32,6 +37,12 @@ public static void Convert(this ColorProfileConverter converter, Rea
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ converter.ConvertUsingIccProfile(source, destination);
+ return;
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS.
diff --git a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieLabRgb.cs b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieLabRgb.cs
index 47e4d2a80a..51be13799c 100644
--- a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieLabRgb.cs
+++ b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieLabRgb.cs
@@ -12,6 +12,11 @@ public static TTo Convert(this ColorProfileConverter converter, in T
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ return converter.ConvertUsingIccProfile(source);
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS
@@ -33,6 +38,12 @@ public static void Convert(this ColorProfileConverter converter, Rea
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ converter.ConvertUsingIccProfile(source, destination);
+ return;
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS.
diff --git a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieXyzCieLab.cs b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieXyzCieLab.cs
index 6b1575d04c..3bab4e7b16 100644
--- a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieXyzCieLab.cs
+++ b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieXyzCieLab.cs
@@ -12,6 +12,11 @@ public static TTo Convert(this ColorProfileConverter converter, in T
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ return converter.ConvertUsingIccProfile(source);
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS
@@ -32,6 +37,12 @@ public static void Convert(this ColorProfileConverter converter, Rea
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ converter.ConvertUsingIccProfile(source, destination);
+ return;
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS.
diff --git a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieXyzCieXyz.cs b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieXyzCieXyz.cs
index 8f56a5a663..5188511476 100644
--- a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieXyzCieXyz.cs
+++ b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieXyzCieXyz.cs
@@ -12,6 +12,11 @@ public static TTo Convert(this ColorProfileConverter converter, in T
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ return converter.ConvertUsingIccProfile(source);
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS
@@ -29,6 +34,12 @@ public static void Convert(this ColorProfileConverter converter, Rea
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ converter.ConvertUsingIccProfile(source, destination);
+ return;
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS.
diff --git a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieXyzRgb.cs b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieXyzRgb.cs
index 9cc0bd9436..c56bf214b9 100644
--- a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieXyzRgb.cs
+++ b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsCieXyzRgb.cs
@@ -12,6 +12,11 @@ public static TTo Convert(this ColorProfileConverter converter, in T
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ return converter.ConvertUsingIccProfile(source);
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS
@@ -32,6 +37,12 @@ public static void Convert(this ColorProfileConverter converter, Rea
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ converter.ConvertUsingIccProfile(source, destination);
+ return;
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS.
diff --git a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsIcc.cs b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsIcc.cs
new file mode 100644
index 0000000000..c33f40001a
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsIcc.cs
@@ -0,0 +1,731 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+using System.Buffers;
+using System.Diagnostics.CodeAnalysis;
+using System.Numerics;
+using System.Runtime.CompilerServices;
+using System.Runtime.Intrinsics;
+using SixLabors.ImageSharp.ColorProfiles.Conversion.Icc;
+using SixLabors.ImageSharp.ColorProfiles.Icc;
+using SixLabors.ImageSharp.Memory;
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles;
+
+internal static class ColorProfileConverterExtensionsIcc
+{
+ private static readonly float[] PcsV2FromBlackPointScale =
+ [0.9965153F, 0.9965269F, 0.9965208F, 1F,
+ 0.9965153F, 0.9965269F, 0.9965208F, 1F,
+ 0.9965153F, 0.9965269F, 0.9965208F, 1F,
+ 0.9965153F, 0.9965269F, 0.9965208F, 1F];
+
+ private static readonly float[] PcsV2FromBlackPointOffset =
+ [0.00336F, 0.0034731F, 0.00287F, 0F,
+ 0.00336F, 0.0034731F, 0.00287F, 0F,
+ 0.00336F, 0.0034731F, 0.00287F, 0F,
+ 0.00336F, 0.0034731F, 0.00287F, 0F];
+
+ private static readonly float[] PcsV2ToBlackPointScale =
+ [1.0034969F, 1.0034852F, 1.0034913F, 1F,
+ 1.0034969F, 1.0034852F, 1.0034913F, 1F,
+ 1.0034969F, 1.0034852F, 1.0034913F, 1F,
+ 1.0034969F, 1.0034852F, 1.0034913F, 1F];
+
+ private static readonly float[] PcsV2ToBlackPointOffset =
+ [0.0033717495F, 0.0034852044F, 0.0028800198F, 0F,
+ 0.0033717495F, 0.0034852044F, 0.0028800198F, 0F,
+ 0.0033717495F, 0.0034852044F, 0.0028800198F, 0F,
+ 0.0033717495F, 0.0034852044F, 0.0028800198F, 0F];
+
+ internal static TTo ConvertUsingIccProfile(this ColorProfileConverter converter, in TFrom source)
+ where TFrom : struct, IColorProfile
+ where TTo : struct, IColorProfile
+ {
+ // TODO: Validation of ICC Profiles against color profile. Is this possible?
+ if (converter.Options.SourceIccProfile is null)
+ {
+ throw new InvalidOperationException("Source ICC profile is missing.");
+ }
+
+ if (converter.Options.TargetIccProfile is null)
+ {
+ throw new InvalidOperationException("Target ICC profile is missing.");
+ }
+
+ ConversionParams sourceParams = new(converter.Options.SourceIccProfile, toPcs: true);
+ ConversionParams targetParams = new(converter.Options.TargetIccProfile, toPcs: false);
+
+ ColorProfileConverter pcsConverter = new(new ColorConversionOptions
+ {
+ MemoryAllocator = converter.Options.MemoryAllocator,
+ SourceWhitePoint = new CieXyz(converter.Options.SourceIccProfile.Header.PcsIlluminant),
+ TargetWhitePoint = new CieXyz(converter.Options.TargetIccProfile.Header.PcsIlluminant),
+ });
+
+ // Normalize the source, then convert to the PCS space.
+ Vector4 sourcePcs = sourceParams.Converter.Calculate(source.ToScaledVector4());
+
+ // If both profiles need PCS adjustment, they both share the same unadjusted PCS space
+ // cancelling out the need to make the adjustment
+ // except if using TRC transforms, which always requires perceptual handling
+ // TODO: this does not include adjustment for absolute intent, which would double existing complexity, suggest throwing exception and addressing in future update
+ bool anyProfileNeedsPerceptualAdjustment = sourceParams.HasNoPerceptualHandling || targetParams.HasNoPerceptualHandling;
+ bool oneProfileHasV2PerceptualAdjustment = sourceParams.HasV2PerceptualHandling ^ targetParams.HasV2PerceptualHandling;
+
+ Vector4 targetPcs = anyProfileNeedsPerceptualAdjustment || oneProfileHasV2PerceptualAdjustment
+ ? GetTargetPcsWithPerceptualAdjustment(sourcePcs, sourceParams, targetParams, pcsConverter)
+ : GetTargetPcsWithoutAdjustment(sourcePcs, sourceParams, targetParams, pcsConverter);
+
+ return TTo.FromScaledVector4(targetParams.Converter.Calculate(targetPcs));
+ }
+
+ internal static void ConvertUsingIccProfile(this ColorProfileConverter converter, ReadOnlySpan source, Span destination)
+ where TFrom : struct, IColorProfile
+ where TTo : struct, IColorProfile
+ {
+ // TODO: Validation of ICC Profiles against color profile. Is this possible?
+ if (converter.Options.SourceIccProfile is null)
+ {
+ throw new InvalidOperationException("Source ICC profile is missing.");
+ }
+
+ if (converter.Options.TargetIccProfile is null)
+ {
+ throw new InvalidOperationException("Target ICC profile is missing.");
+ }
+
+ Guard.MustBeGreaterThanOrEqualTo(source.Length, destination.Length, nameof(destination));
+
+ ConversionParams sourceParams = new(converter.Options.SourceIccProfile, toPcs: true);
+ ConversionParams targetParams = new(converter.Options.TargetIccProfile, toPcs: false);
+
+ ColorProfileConverter pcsConverter = new(new ColorConversionOptions
+ {
+ MemoryAllocator = converter.Options.MemoryAllocator,
+ SourceWhitePoint = new CieXyz(converter.Options.SourceIccProfile.Header.PcsIlluminant),
+ TargetWhitePoint = new CieXyz(converter.Options.TargetIccProfile.Header.PcsIlluminant),
+ });
+
+ using IMemoryOwner pcsBuffer = converter.Options.MemoryAllocator.Allocate(source.Length);
+ Span pcs = pcsBuffer.GetSpan();
+
+ // Normalize the source, then convert to the PCS space.
+ TFrom.ToScaledVector4(source, pcs);
+ sourceParams.Converter.Calculate(pcs, pcs);
+
+ // If both profiles need PCS adjustment, they both share the same unadjusted PCS space
+ // cancelling out the need to make the adjustment
+ // except if using TRC transforms, which always requires perceptual handling
+ // TODO: this does not include adjustment for absolute intent, which would double existing complexity, suggest throwing exception and addressing in future update
+ bool anyProfileNeedsPerceptualAdjustment = sourceParams.HasNoPerceptualHandling || targetParams.HasNoPerceptualHandling;
+ bool oneProfileHasV2PerceptualAdjustment = sourceParams.HasV2PerceptualHandling ^ targetParams.HasV2PerceptualHandling;
+
+ if (anyProfileNeedsPerceptualAdjustment || oneProfileHasV2PerceptualAdjustment)
+ {
+ GetTargetPcsWithPerceptualAdjustment(pcs, sourceParams, targetParams, pcsConverter);
+ }
+ else
+ {
+ GetTargetPcsWithoutAdjustment(pcs, sourceParams, targetParams, pcsConverter);
+ }
+
+ // Convert to the target space.
+ targetParams.Converter.Calculate(pcs, pcs);
+ TTo.FromScaledVector4(pcs, destination);
+ }
+
+ private static Vector4 GetTargetPcsWithoutAdjustment(
+ Vector4 sourcePcs,
+ ConversionParams sourceParams,
+ ConversionParams targetParams,
+ ColorProfileConverter pcsConverter)
+ {
+ // Profile connecting spaces can only be Lab, XYZ.
+ // 16-bit Lab encodings changed from v2 to v4, but 16-bit LUTs always use the legacy encoding regardless of version
+ // so ensure that Lab is using the correct encoding when a 16-bit LUT is used
+ switch (sourceParams.PcsType)
+ {
+ // Convert from Lab to XYZ.
+ case IccColorSpaceType.CieLab when targetParams.PcsType is IccColorSpaceType.CieXyz:
+ {
+ sourcePcs = sourceParams.Is16BitLutEntry ? LabV2ToLab(sourcePcs) : sourcePcs;
+ CieLab lab = CieLab.FromScaledVector4(sourcePcs);
+ CieXyz xyz = pcsConverter.Convert(in lab);
+ return xyz.ToScaledVector4();
+ }
+
+ // Convert from XYZ to Lab.
+ case IccColorSpaceType.CieXyz when targetParams.PcsType is IccColorSpaceType.CieLab:
+ {
+ CieXyz xyz = CieXyz.FromScaledVector4(sourcePcs);
+ CieLab lab = pcsConverter.Convert(in xyz);
+ Vector4 targetPcs = lab.ToScaledVector4();
+ return targetParams.Is16BitLutEntry ? LabToLabV2(targetPcs) : targetPcs;
+ }
+
+ // Convert from XYZ to XYZ.
+ case IccColorSpaceType.CieXyz when targetParams.PcsType is IccColorSpaceType.CieXyz:
+ {
+ CieXyz xyz = CieXyz.FromScaledVector4(sourcePcs);
+ CieXyz targetXyz = pcsConverter.Convert(in xyz);
+ return targetXyz.ToScaledVector4();
+ }
+
+ // Convert from Lab to Lab.
+ case IccColorSpaceType.CieLab when targetParams.PcsType is IccColorSpaceType.CieLab:
+ {
+ // if both source and target LUT use same v2 LAB encoding, no need to correct them
+ if (sourceParams.Is16BitLutEntry && targetParams.Is16BitLutEntry)
+ {
+ CieLab sourceLab = CieLab.FromScaledVector4(sourcePcs);
+ CieLab targetLab = pcsConverter.Convert(in sourceLab);
+ return targetLab.ToScaledVector4();
+ }
+ else
+ {
+ sourcePcs = sourceParams.Is16BitLutEntry ? LabV2ToLab(sourcePcs) : sourcePcs;
+ CieLab sourceLab = CieLab.FromScaledVector4(sourcePcs);
+ CieLab targetLab = pcsConverter.Convert(in sourceLab);
+ Vector4 targetPcs = targetLab.ToScaledVector4();
+ return targetParams.Is16BitLutEntry ? LabToLabV2(targetPcs) : targetPcs;
+ }
+ }
+
+ default:
+ throw new ArgumentOutOfRangeException($"Source PCS {sourceParams.PcsType} to target PCS {targetParams.PcsType} is not supported");
+ }
+ }
+
+ private static void GetTargetPcsWithoutAdjustment(
+ Span pcs,
+ ConversionParams sourceParams,
+ ConversionParams targetParams,
+ ColorProfileConverter pcsConverter)
+ {
+ // Profile connecting spaces can only be Lab, XYZ.
+ // 16-bit Lab encodings changed from v2 to v4, but 16-bit LUTs always use the legacy encoding regardless of version
+ // so ensure that Lab is using the correct encoding when a 16-bit LUT is used
+ switch (sourceParams.PcsType)
+ {
+ // Convert from Lab to XYZ.
+ case IccColorSpaceType.CieLab when targetParams.PcsType is IccColorSpaceType.CieXyz:
+ {
+ if (sourceParams.Is16BitLutEntry)
+ {
+ LabV2ToLab(pcs, pcs);
+ }
+
+ using IMemoryOwner pcsFromBuffer = pcsConverter.Options.MemoryAllocator.Allocate(pcs.Length);
+ Span pcsFrom = pcsFromBuffer.GetSpan();
+
+ using IMemoryOwner pcsToBuffer = pcsConverter.Options.MemoryAllocator.Allocate(pcs.Length);
+ Span pcsTo = pcsToBuffer.GetSpan();
+
+ CieLab.FromScaledVector4(pcs, pcsFrom);
+ pcsConverter.Convert(pcsFrom, pcsTo);
+
+ CieXyz.ToScaledVector4(pcsTo, pcs);
+ break;
+ }
+
+ // Convert from XYZ to Lab.
+ case IccColorSpaceType.CieXyz when targetParams.PcsType is IccColorSpaceType.CieLab:
+ {
+ using IMemoryOwner pcsFromBuffer = pcsConverter.Options.MemoryAllocator.Allocate(pcs.Length);
+ Span pcsFrom = pcsFromBuffer.GetSpan();
+
+ using IMemoryOwner pcsToBuffer = pcsConverter.Options.MemoryAllocator.Allocate(pcs.Length);
+ Span pcsTo = pcsToBuffer.GetSpan();
+
+ CieXyz.FromScaledVector4(pcs, pcsFrom);
+ pcsConverter.Convert(pcsFrom, pcsTo);
+
+ CieLab.ToScaledVector4(pcsTo, pcs);
+
+ if (targetParams.Is16BitLutEntry)
+ {
+ LabToLabV2(pcs, pcs);
+ }
+
+ break;
+ }
+
+ // Convert from XYZ to XYZ.
+ case IccColorSpaceType.CieXyz when targetParams.PcsType is IccColorSpaceType.CieXyz:
+ {
+ using IMemoryOwner pcsFromToBuffer = pcsConverter.Options.MemoryAllocator.Allocate(pcs.Length);
+ Span pcsFromTo = pcsFromToBuffer.GetSpan();
+
+ CieXyz.FromScaledVector4(pcs, pcsFromTo);
+ pcsConverter.Convert(pcsFromTo, pcsFromTo);
+
+ CieXyz.ToScaledVector4(pcsFromTo, pcs);
+ break;
+ }
+
+ // Convert from Lab to Lab.
+ case IccColorSpaceType.CieLab when targetParams.PcsType is IccColorSpaceType.CieLab:
+ {
+ using IMemoryOwner pcsFromToBuffer = pcsConverter.Options.MemoryAllocator.Allocate(pcs.Length);
+ Span pcsFromTo = pcsFromToBuffer.GetSpan();
+
+ // if both source and target LUT use same v2 LAB encoding, no need to correct them
+ if (sourceParams.Is16BitLutEntry && targetParams.Is16BitLutEntry)
+ {
+ CieLab.FromScaledVector4(pcs, pcsFromTo);
+ pcsConverter.Convert(pcsFromTo, pcsFromTo);
+ CieLab.ToScaledVector4(pcsFromTo, pcs);
+ }
+ else
+ {
+ if (sourceParams.Is16BitLutEntry)
+ {
+ LabV2ToLab(pcs, pcs);
+ }
+
+ CieLab.FromScaledVector4(pcs, pcsFromTo);
+ pcsConverter.Convert(pcsFromTo, pcsFromTo);
+ CieLab.ToScaledVector4(pcsFromTo, pcs);
+
+ if (targetParams.Is16BitLutEntry)
+ {
+ LabToLabV2(pcs, pcs);
+ }
+ }
+
+ break;
+ }
+
+ default:
+ throw new ArgumentOutOfRangeException($"Source PCS {sourceParams.PcsType} to target PCS {targetParams.PcsType} is not supported");
+ }
+ }
+
+ ///
+ /// Effectively this is with an extra step in the middle.
+ /// It adjusts PCS by compensating for the black point used for perceptual intent in v2 profiles.
+ /// The adjustment needs to be performed in XYZ space, potentially an overhead of 2 more conversions.
+ /// Not required if both spaces need V2 correction, since they both have the same understanding of the PCS.
+ /// Not compatible with PCS adjustment for absolute intent.
+ ///
+ /// The source PCS values.
+ /// The source profile parameters.
+ /// The target profile parameters.
+ /// The converter to use for the PCS adjustments.
+ /// Thrown when the source or target PCS is not supported.
+ private static Vector4 GetTargetPcsWithPerceptualAdjustment(
+ Vector4 sourcePcs,
+ ConversionParams sourceParams,
+ ConversionParams targetParams,
+ ColorProfileConverter pcsConverter)
+ {
+ // all conversions are funneled through XYZ in case PCS adjustments need to be made
+ CieXyz xyz;
+
+ switch (sourceParams.PcsType)
+ {
+ // 16-bit Lab encodings changed from v2 to v4, but 16-bit LUTs always use the legacy encoding regardless of version
+ // so convert Lab to modern v4 encoding when returned from a 16-bit LUT
+ case IccColorSpaceType.CieLab:
+ sourcePcs = sourceParams.Is16BitLutEntry ? LabV2ToLab(sourcePcs) : sourcePcs;
+ CieLab lab = CieLab.FromScaledVector4(sourcePcs);
+ xyz = pcsConverter.Convert(in lab);
+ break;
+ case IccColorSpaceType.CieXyz:
+ xyz = CieXyz.FromScaledVector4(sourcePcs);
+ break;
+ default:
+ throw new ArgumentOutOfRangeException($"Source PCS {sourceParams.PcsType} is not supported");
+ }
+
+ bool oneProfileHasV2PerceptualAdjustment = sourceParams.HasV2PerceptualHandling ^ targetParams.HasV2PerceptualHandling;
+
+ // when converting from device to PCS with v2 perceptual intent
+ // the black point needs to be adjusted to v4 after converting the PCS values
+ if (sourceParams.HasNoPerceptualHandling ||
+ (oneProfileHasV2PerceptualAdjustment && sourceParams.HasV2PerceptualHandling))
+ {
+ Vector3 vector = xyz.ToVector3();
+
+ // when using LAB PCS, negative values are clipped before PCS adjustment (in DemoIccMAX)
+ if (sourceParams.PcsType == IccColorSpaceType.CieLab)
+ {
+ vector = Vector3.Max(vector, Vector3.Zero);
+ }
+
+ xyz = new CieXyz(AdjustPcsFromV2BlackPoint(vector));
+ }
+
+ // when converting from PCS to device with v2 perceptual intent
+ // the black point needs to be adjusted to v2 before converting the PCS values
+ if (targetParams.HasNoPerceptualHandling ||
+ (oneProfileHasV2PerceptualAdjustment && targetParams.HasV2PerceptualHandling))
+ {
+ Vector3 vector = AdjustPcsToV2BlackPoint(xyz.AsVector3Unsafe());
+
+ // when using XYZ PCS, negative values are clipped after PCS adjustment (in DemoIccMAX)
+ if (targetParams.PcsType == IccColorSpaceType.CieXyz)
+ {
+ vector = Vector3.Max(vector, Vector3.Zero);
+ }
+
+ xyz = new CieXyz(vector);
+ }
+
+ switch (targetParams.PcsType)
+ {
+ // 16-bit Lab encodings changed from v2 to v4, but 16-bit LUTs always use the legacy encoding regardless of version
+ // so convert Lab back to legacy encoding before using in a 16-bit LUT
+ case IccColorSpaceType.CieLab:
+ CieLab lab = pcsConverter.Convert(in xyz);
+ Vector4 targetPcs = lab.ToScaledVector4();
+ return targetParams.Is16BitLutEntry ? LabToLabV2(targetPcs) : targetPcs;
+ case IccColorSpaceType.CieXyz:
+ return xyz.ToScaledVector4();
+ default:
+ throw new ArgumentOutOfRangeException($"Target PCS {targetParams.PcsType} is not supported");
+ }
+ }
+
+ ///
+ /// Effectively this is with an extra step in the middle.
+ /// It adjusts PCS by compensating for the black point used for perceptual intent in v2 profiles.
+ /// The adjustment needs to be performed in XYZ space, potentially an overhead of 2 more conversions.
+ /// Not required if both spaces need V2 correction, since they both have the same understanding of the PCS.
+ /// Not compatible with PCS adjustment for absolute intent.
+ ///
+ /// The PCS values from the source.
+ /// The source profile parameters.
+ /// The target profile parameters.
+ /// The converter to use for the PCS adjustments.
+ /// Thrown when the source or target PCS is not supported.
+ private static void GetTargetPcsWithPerceptualAdjustment(
+ Span pcs,
+ ConversionParams sourceParams,
+ ConversionParams targetParams,
+ ColorProfileConverter pcsConverter)
+ {
+ // All conversions are funneled through XYZ in case PCS adjustments need to be made
+ using IMemoryOwner xyzBuffer = pcsConverter.Options.MemoryAllocator.Allocate(pcs.Length);
+ Span xyz = xyzBuffer.GetSpan();
+
+ switch (sourceParams.PcsType)
+ {
+ // 16-bit Lab encodings changed from v2 to v4, but 16-bit LUTs always use the legacy encoding regardless of version
+ // so convert Lab to modern v4 encoding when returned from a 16-bit LUT
+ case IccColorSpaceType.CieLab:
+ {
+ if (sourceParams.Is16BitLutEntry)
+ {
+ LabV2ToLab(pcs, pcs);
+ }
+
+ using IMemoryOwner pcsFromBuffer = pcsConverter.Options.MemoryAllocator.Allocate(pcs.Length);
+ Span pcsFrom = pcsFromBuffer.GetSpan();
+ CieLab.FromScaledVector4(pcs, pcsFrom);
+ pcsConverter.Convert(pcsFrom, xyz);
+ break;
+ }
+
+ case IccColorSpaceType.CieXyz:
+ CieXyz.FromScaledVector4(pcs, xyz);
+ break;
+ default:
+ throw new ArgumentOutOfRangeException($"Source PCS {sourceParams.PcsType} is not supported");
+ }
+
+ bool oneProfileHasV2PerceptualAdjustment = sourceParams.HasV2PerceptualHandling ^ targetParams.HasV2PerceptualHandling;
+
+ using IMemoryOwner vectorBuffer = pcsConverter.Options.MemoryAllocator.Allocate(pcs.Length);
+ Span vector = vectorBuffer.GetSpan();
+
+ // When converting from device to PCS with v2 perceptual intent
+ // the black point needs to be adjusted to v4 after converting the PCS values
+ if (sourceParams.HasNoPerceptualHandling ||
+ (oneProfileHasV2PerceptualAdjustment && sourceParams.HasV2PerceptualHandling))
+ {
+ CieXyz.ToVector4(xyz, vector);
+
+ // When using LAB PCS, negative values are clipped before PCS adjustment (in DemoIccMAX)
+ if (sourceParams.PcsType == IccColorSpaceType.CieLab)
+ {
+ ClipNegative(vector);
+ }
+
+ AdjustPcsFromV2BlackPoint(vector, vector);
+ CieXyz.FromVector4(vector, xyz);
+ }
+
+ // When converting from PCS to device with v2 perceptual intent
+ // the black point needs to be adjusted to v2 before converting the PCS values
+ if (targetParams.HasNoPerceptualHandling ||
+ (oneProfileHasV2PerceptualAdjustment && targetParams.HasV2PerceptualHandling))
+ {
+ CieXyz.ToVector4(xyz, vector);
+ AdjustPcsToV2BlackPoint(vector, vector);
+
+ // When using XYZ PCS, negative values are clipped after PCS adjustment (in DemoIccMAX)
+ if (targetParams.PcsType == IccColorSpaceType.CieXyz)
+ {
+ ClipNegative(vector);
+ }
+
+ CieXyz.FromVector4(vector, xyz);
+ }
+
+ switch (targetParams.PcsType)
+ {
+ // 16-bit Lab encodings changed from v2 to v4, but 16-bit LUTs always use the legacy encoding regardless of version
+ // so convert Lab back to legacy encoding before using in a 16-bit LUT
+ case IccColorSpaceType.CieLab:
+ {
+ using IMemoryOwner pcsToBuffer = pcsConverter.Options.MemoryAllocator.Allocate(pcs.Length);
+ Span pcsTo = pcsToBuffer.GetSpan();
+ pcsConverter.Convert(xyz, pcsTo);
+
+ CieLab.ToScaledVector4(pcsTo, pcs);
+
+ if (targetParams.Is16BitLutEntry)
+ {
+ LabToLabV2(pcs, pcs);
+ }
+
+ break;
+ }
+
+ case IccColorSpaceType.CieXyz:
+ CieXyz.ToScaledVector4(xyz, pcs);
+ break;
+ default:
+ throw new ArgumentOutOfRangeException($"Target PCS {targetParams.PcsType} is not supported");
+ }
+ }
+
+ // as per DemoIccMAX icPerceptual values in IccCmm.h
+ // refBlack = 0.00336F, 0.0034731F, 0.00287F
+ // refWhite = 0.9642F, 1.0000F, 0.8249F
+ // scale = 1 - (refBlack / refWhite)
+ // offset = refBlack
+ private static Vector3 AdjustPcsFromV2BlackPoint(Vector3 xyz)
+ => (xyz * new Vector3(0.9965153F, 0.9965269F, 0.9965208F)) + new Vector3(0.00336F, 0.0034731F, 0.00287F);
+
+ // as per DemoIccMAX icPerceptual values in IccCmm.h
+ // refBlack = 0.00336F, 0.0034731F, 0.00287F
+ // refWhite = 0.9642F, 1.0000F, 0.8249F
+ // scale = 1 / (1 - (refBlack / refWhite))
+ // offset = -refBlack * scale
+ private static Vector3 AdjustPcsToV2BlackPoint(Vector3 xyz)
+ => (xyz * new Vector3(1.0034969F, 1.0034852F, 1.0034913F)) - new Vector3(0.0033717495F, 0.0034852044F, 0.0028800198F);
+
+ private static void AdjustPcsFromV2BlackPoint(Span source, Span destination)
+ {
+ if (Vector.IsHardwareAccelerated && Vector.IsSupported &&
+ Vector.Count <= Vector512.Count &&
+ source.Length * 4 >= Vector.Count)
+ {
+ // TODO: Check our constants. They may require scaling.
+ Vector vScale = new(PcsV2FromBlackPointScale.AsSpan()[..Vector.Count]);
+ Vector vOffset = new(PcsV2FromBlackPointOffset.AsSpan()[..Vector.Count]);
+
+ // SIMD loop
+ int i = 0;
+ int simdBatchSize = Vector.Count / 4; // Number of Vector4 elements per SIMD batch
+ for (; i <= source.Length - simdBatchSize; i += simdBatchSize)
+ {
+ // Load the vector from source span
+ Vector v = Unsafe.ReadUnaligned>(ref Unsafe.As(ref source[i]));
+
+ // Scale and offset the vector
+ v *= vScale;
+ v += vOffset;
+
+ // Write the vector to the destination span
+ Unsafe.WriteUnaligned(ref Unsafe.As(ref destination[i]), v);
+ }
+
+ // Scalar fallback for remaining elements
+ for (; i < source.Length; i++)
+ {
+ Vector4 s = source[i];
+ s *= new Vector4(0.9965153F, 0.9965269F, 0.9965208F, 1F);
+ s += new Vector4(0.00336F, 0.0034731F, 0.00287F, 0F);
+ destination[i] = s;
+ }
+ }
+ else
+ {
+ // Scalar fallback if SIMD is not supported
+ for (int i = 0; i < source.Length; i++)
+ {
+ Vector4 s = source[i];
+ s *= new Vector4(0.9965153F, 0.9965269F, 0.9965208F, 1F);
+ s += new Vector4(0.00336F, 0.0034731F, 0.00287F, 0F);
+ destination[i] = s;
+ }
+ }
+ }
+
+ private static void AdjustPcsToV2BlackPoint(Span source, Span destination)
+ {
+ if (Vector.IsHardwareAccelerated && Vector.IsSupported &&
+ Vector.Count <= Vector512.Count &&
+ source.Length * 4 >= Vector.Count)
+ {
+ // TODO: Check our constants. They may require scaling.
+ Vector vScale = new(PcsV2ToBlackPointScale.AsSpan()[..Vector.Count]);
+ Vector vOffset = new(PcsV2ToBlackPointOffset.AsSpan()[..Vector.Count]);
+
+ // SIMD loop
+ int i = 0;
+ int simdBatchSize = Vector.Count / 4; // Number of Vector4 elements per SIMD batch
+ for (; i <= source.Length - simdBatchSize; i += simdBatchSize)
+ {
+ // Load the vector from source span
+ Vector v = Unsafe.ReadUnaligned>(ref Unsafe.As(ref source[i]));
+
+ // Scale and offset the vector
+ v *= vScale;
+ v -= vOffset;
+
+ // Write the vector to the destination span
+ Unsafe.WriteUnaligned(ref Unsafe.As(ref destination[i]), v);
+ }
+
+ // Scalar fallback for remaining elements
+ for (; i < source.Length; i++)
+ {
+ Vector4 s = source[i];
+ s *= new Vector4(1.0034969F, 1.0034852F, 1.0034913F, 1F);
+ s -= new Vector4(0.0033717495F, 0.0034852044F, 0.0028800198F, 0F);
+ destination[i] = s;
+ }
+ }
+ else
+ {
+ // Scalar fallback if SIMD is not supported
+ for (int i = 0; i < source.Length; i++)
+ {
+ Vector4 s = source[i];
+ s *= new Vector4(1.0034969F, 1.0034852F, 1.0034913F, 1F);
+ s -= new Vector4(0.0033717495F, 0.0034852044F, 0.0028800198F, 0F);
+ destination[i] = s;
+ }
+ }
+ }
+
+ private static void ClipNegative(Span source)
+ {
+ if (Vector.IsHardwareAccelerated && Vector.IsSupported && Vector.Count >= source.Length * 4)
+ {
+ // SIMD loop
+ int i = 0;
+ int simdBatchSize = Vector.Count / 4; // Number of Vector4 elements per SIMD batch
+ for (; i <= source.Length - simdBatchSize; i += simdBatchSize)
+ {
+ // Load the vector from source span
+ Vector v = Unsafe.ReadUnaligned>(ref Unsafe.As(ref source[i]));
+
+ v = Vector.Max(v, Vector.Zero);
+
+ // Write the vector to the destination span
+ Unsafe.WriteUnaligned(ref Unsafe.As(ref source[i]), v);
+ }
+
+ // Scalar fallback for remaining elements
+ for (; i < source.Length; i++)
+ {
+ ref Vector4 s = ref source[i];
+ s = Vector4.Max(s, Vector4.Zero);
+ }
+ }
+ else
+ {
+ // Scalar fallback if SIMD is not supported
+ for (int i = 0; i < source.Length; i++)
+ {
+ ref Vector4 s = ref source[i];
+ s = Vector4.Max(s, Vector4.Zero);
+ }
+ }
+ }
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ private static Vector4 LabToLabV2(Vector4 input)
+ => input * 65280F / 65535F;
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ private static Vector4 LabV2ToLab(Vector4 input)
+ => input * 65535F / 65280F;
+
+ private static void LabToLabV2(Span source, Span destination)
+ => LabToLab(source, destination, 65280F / 65535F);
+
+ private static void LabV2ToLab(Span source, Span destination)
+ => LabToLab(source, destination, 65535F / 65280F);
+
+ private static void LabToLab(Span source, Span destination, [ConstantExpected] float scale)
+ {
+ if (Vector.IsHardwareAccelerated && Vector.IsSupported)
+ {
+ Vector vScale = new(scale);
+ int i = 0;
+
+ // SIMD loop
+ int simdBatchSize = Vector.Count / 4; // Number of Vector4 elements per SIMD batch
+ for (; i <= source.Length - simdBatchSize; i += simdBatchSize)
+ {
+ // Load the vector from source span
+ Vector v = Unsafe.ReadUnaligned>(ref Unsafe.As(ref source[i]));
+
+ // Scale the vector
+ v *= vScale;
+
+ // Write the scaled vector to the destination span
+ Unsafe.WriteUnaligned(ref Unsafe.As(ref destination[i]), v);
+ }
+
+ // Scalar fallback for remaining elements
+ for (; i < source.Length; i++)
+ {
+ destination[i] = source[i] * scale;
+ }
+ }
+ else
+ {
+ // Scalar fallback if SIMD is not supported
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = source[i] * scale;
+ }
+ }
+ }
+
+ private class ConversionParams
+ {
+ private readonly IccProfile profile;
+
+ internal ConversionParams(IccProfile profile, bool toPcs)
+ {
+ this.profile = profile;
+ this.Converter = toPcs ? new IccDataToPcsConverter(profile) : new IccPcsToDataConverter(profile);
+ }
+
+ internal IccConverterBase Converter { get; }
+
+ internal IccProfileHeader Header => this.profile.Header;
+
+ internal IccRenderingIntent Intent => this.Header.RenderingIntent;
+
+ internal IccColorSpaceType PcsType => this.Header.ProfileConnectionSpace;
+
+ internal IccVersion Version => this.Header.Version;
+
+ internal bool HasV2PerceptualHandling => this.Intent == IccRenderingIntent.Perceptual && this.Version.Major == 2;
+
+ internal bool HasNoPerceptualHandling => this.Intent == IccRenderingIntent.Perceptual && this.Converter.IsTrc;
+
+ internal bool Is16BitLutEntry => this.Converter.Is16BitLutEntry;
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsRgbCieLab.cs b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsRgbCieLab.cs
index 415dd94c3f..badbcc6831 100644
--- a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsRgbCieLab.cs
+++ b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsRgbCieLab.cs
@@ -12,6 +12,11 @@ public static TTo Convert(this ColorProfileConverter converter, in T
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ return converter.ConvertUsingIccProfile(source);
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS
@@ -33,6 +38,12 @@ public static void Convert(this ColorProfileConverter converter, Rea
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ converter.ConvertUsingIccProfile(source, destination);
+ return;
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS.
diff --git a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsRgbCieXyz.cs b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsRgbCieXyz.cs
index a13f645778..cd7d5e4d65 100644
--- a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsRgbCieXyz.cs
+++ b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsRgbCieXyz.cs
@@ -12,6 +12,11 @@ public static TTo Convert(this ColorProfileConverter converter, in T
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ return converter.ConvertUsingIccProfile(source);
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS
@@ -32,6 +37,12 @@ public static void Convert(this ColorProfileConverter converter, Rea
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ converter.ConvertUsingIccProfile(source, destination);
+ return;
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS.
diff --git a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsRgbRgb.cs b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsRgbRgb.cs
index c1c75dea1b..2a4b64b1ca 100644
--- a/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsRgbRgb.cs
+++ b/src/ImageSharp/ColorProfiles/ColorProfileConverterExtensionsRgbRgb.cs
@@ -12,6 +12,11 @@ public static TTo Convert(this ColorProfileConverter converter, in T
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ return converter.ConvertUsingIccProfile(source);
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS
@@ -33,6 +38,12 @@ public static void Convert(this ColorProfileConverter converter, Rea
where TFrom : struct, IColorProfile
where TTo : struct, IColorProfile
{
+ if (converter.ShouldUseIccProfiles())
+ {
+ converter.ConvertUsingIccProfile(source, destination);
+ return;
+ }
+
ColorConversionOptions options = converter.Options;
// Convert to input PCS.
diff --git a/src/ImageSharp/ColorProfiles/Hsl.cs b/src/ImageSharp/ColorProfiles/Hsl.cs
index 2c98c7df99..7a9365fb75 100644
--- a/src/ImageSharp/ColorProfiles/Hsl.cs
+++ b/src/ImageSharp/ColorProfiles/Hsl.cs
@@ -41,6 +41,16 @@ public Hsl(Vector3 vector)
this.L = vector.Z;
}
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+#pragma warning disable SA1313 // Parameter names should begin with lower-case letter
+ private Hsl(Vector3 vector, bool _)
+#pragma warning restore SA1313 // Parameter names should begin with lower-case letter
+ {
+ this.H = vector.X;
+ this.S = vector.Y;
+ this.L = vector.Z;
+ }
+
///
/// Gets the hue component.
/// A value ranging between 0 and 360.
@@ -83,6 +93,38 @@ public Hsl(Vector3 vector)
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static bool operator !=(Hsl left, Hsl right) => !left.Equals(right);
+ ///
+ public Vector4 ToScaledVector4()
+ => new(this.AsVector3Unsafe() / 360F, 1F);
+
+ ///
+ public static Hsl FromScaledVector4(Vector4 source)
+ => new(source.AsVector3() * 360F, true);
+
+ ///
+ public static void ToScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = source[i].ToScaledVector4();
+ }
+ }
+
+ ///
+ public static void FromScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = FromScaledVector4(source[i]);
+ }
+ }
+
///
public static Hsl FromProfileConnectingSpace(ColorConversionOptions options, in Rgb source)
{
diff --git a/src/ImageSharp/ColorProfiles/Hsv.cs b/src/ImageSharp/ColorProfiles/Hsv.cs
index 7535f2463d..1e013fe1fb 100644
--- a/src/ImageSharp/ColorProfiles/Hsv.cs
+++ b/src/ImageSharp/ColorProfiles/Hsv.cs
@@ -41,6 +41,16 @@ public Hsv(Vector3 vector)
this.V = vector.Z;
}
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+#pragma warning disable SA1313 // Parameter names should begin with lower-case letter
+ private Hsv(Vector3 vector, bool _)
+#pragma warning restore SA1313 // Parameter names should begin with lower-case letter
+ {
+ this.H = vector.X;
+ this.S = vector.Y;
+ this.V = vector.Z;
+ }
+
///
/// Gets the hue component.
/// A value ranging between 0 and 360.
@@ -81,6 +91,38 @@ public Hsv(Vector3 vector)
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static bool operator !=(Hsv left, Hsv right) => !left.Equals(right);
+ ///
+ public Vector4 ToScaledVector4()
+ => new(this.AsVector3Unsafe() / 360F, 1F);
+
+ ///
+ public static Hsv FromScaledVector4(Vector4 source)
+ => new(source.AsVector3() * 360F, true);
+
+ ///
+ public static void ToScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = source[i].ToScaledVector4();
+ }
+ }
+
+ ///
+ public static void FromScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = FromScaledVector4(source[i]);
+ }
+ }
+
///
public static Hsv FromProfileConnectingSpace(ColorConversionOptions options, in Rgb source)
{
diff --git a/src/ImageSharp/ColorProfiles/HunterLab.cs b/src/ImageSharp/ColorProfiles/HunterLab.cs
index 43ad2ac5c0..e978c6de22 100644
--- a/src/ImageSharp/ColorProfiles/HunterLab.cs
+++ b/src/ImageSharp/ColorProfiles/HunterLab.cs
@@ -80,6 +80,49 @@ public HunterLab(Vector3 vector)
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static bool operator !=(HunterLab left, HunterLab right) => !left.Equals(right);
+ ///
+ public Vector4 ToScaledVector4()
+ {
+ Vector3 v3 = default;
+ v3 += this.AsVector3Unsafe();
+ v3 += new Vector3(0, 128F, 128F);
+ v3 /= new Vector3(100F, 255F, 255F);
+ return new Vector4(v3, 1F);
+ }
+
+ ///
+ public static HunterLab FromScaledVector4(Vector4 source)
+ {
+ Vector3 v3 = source.AsVector3();
+ v3 *= new Vector3(100F, 255, 255);
+ v3 -= new Vector3(0, 128F, 128F);
+ return new HunterLab(v3);
+ }
+
+ ///
+ public static void ToScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = source[i].ToScaledVector4();
+ }
+ }
+
+ ///
+ public static void FromScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = FromScaledVector4(source[i]);
+ }
+ }
+
///
public static HunterLab FromProfileConnectingSpace(ColorConversionOptions options, in CieXyz source)
{
@@ -127,7 +170,7 @@ public CieXyz ToProfileConnectingSpace(ColorConversionOptions options)
{
// Conversion algorithm described here:
// http://en.wikipedia.org/wiki/Lab_color_space#Hunter_Lab
- CieXyz whitePoint = options.WhitePoint;
+ CieXyz whitePoint = options.SourceWhitePoint;
float l = this.L, a = this.A, b = this.B;
float xn = whitePoint.X, yn = whitePoint.Y, zn = whitePoint.Z;
diff --git a/src/ImageSharp/ColorProfiles/IColorProfile.cs b/src/ImageSharp/ColorProfiles/IColorProfile.cs
index 6a1b2ee8d0..425e030300 100644
--- a/src/ImageSharp/ColorProfiles/IColorProfile.cs
+++ b/src/ImageSharp/ColorProfiles/IColorProfile.cs
@@ -1,6 +1,8 @@
// Copyright (c) Six Labors.
// Licensed under the Six Labors Split License.
+using System.Numerics;
+
namespace SixLabors.ImageSharp.ColorProfiles;
///
@@ -15,18 +17,60 @@ public interface IColorProfile
public static abstract ChromaticAdaptionWhitePointSource GetChromaticAdaptionWhitePointSource();
}
+///
+/// Defines the contract for all color profiles.
+///
+/// The type of color profile.
+public interface IColorProfile : IColorProfile, IEquatable
+ where TSelf : IColorProfile
+{
+ ///
+ /// Expands the pixel into a generic ("scaled") representation
+ /// with values scaled and clamped between 0 and 1.
+ /// The vector components are typically expanded in least to greatest significance order.
+ ///
+ /// The .
+ public Vector4 ToScaledVector4();
+
+#pragma warning disable CA1000 // Do not declare static members on generic types
+ ///
+ /// Initializes the color instance from a generic a generic ("scaled") representation
+ /// with values scaled and clamped between 0 and 1.
+ ///
+ /// The vector to load the pixel from.
+ /// The .
+ public static abstract TSelf FromScaledVector4(Vector4 source);
+
+ ///
+ /// Converts the span of colors to a generic ("scaled") representation
+ /// with values scaled and clamped between 0 and 1.
+ ///
+ /// The color span to convert from.
+ /// The vector span to write the results to.
+ public static abstract void ToScaledVector4(ReadOnlySpan source, Span destination);
+
+ ///
+ /// Converts the span of colors from a generic ("scaled") representation
+ /// with values scaled and clamped between 0 and 1.
+ ///
+ /// The vector span to convert from.
+ /// The color span to write the results to.
+ public static abstract void FromScaledVector4(ReadOnlySpan source, Span destination);
+#pragma warning restore CA1000 // Do not declare static members on generic types
+}
+
///
/// Defines the contract for all color profiles.
///
/// The type of color profile.
/// The type of color profile connecting space.
-public interface IColorProfile : IColorProfile, IEquatable
+public interface IColorProfile : IColorProfile
where TSelf : IColorProfile
where TProfileSpace : struct, IProfileConnectingSpace
{
#pragma warning disable CA1000 // Do not declare static members on generic types
///
- /// Converts the color from the profile connection space.
+ /// Initializes the color instance from the profile connection space.
///
/// The color profile conversion options.
/// The color profile connecting space.
diff --git a/src/ImageSharp/ColorProfiles/Icc/Calculators/ClutCalculator.cs b/src/ImageSharp/ColorProfiles/Icc/Calculators/ClutCalculator.cs
new file mode 100644
index 0000000000..82d475e578
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/Calculators/ClutCalculator.cs
@@ -0,0 +1,506 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+using System.Numerics;
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Icc.Calculators;
+
+///
+/// Implements interpolation methods for color profile lookup tables.
+/// Adapted from ICC Reference implementation:
+/// https://github.com/InternationalColorConsortium/DemoIccMAX/blob/79ecb74135ad47bac7d42692905a079839b7e105/IccProfLib/IccTagLut.cpp
+///
+internal class ClutCalculator : IVector4Calculator
+{
+ private readonly int inputCount;
+ private readonly int outputCount;
+ private readonly float[] lut;
+ private readonly byte[] gridPointCount;
+ private readonly byte[] maxGridPoint;
+ private readonly int[] indexFactor;
+ private readonly int[] dimSize;
+ private readonly int nodeCount;
+ private readonly float[][] nodes;
+ private readonly float[] g;
+ private readonly uint[] ig;
+ private readonly float[] s;
+ private readonly float[] df;
+ private readonly uint[] nPower;
+ private int n000;
+ private int n001;
+ private int n010;
+ private int n011;
+ private int n100;
+ private int n101;
+ private int n110;
+ private int n111;
+ private int n1000;
+
+ public ClutCalculator(IccClut clut)
+ {
+ Guard.NotNull(clut, nameof(clut));
+ Guard.MustBeGreaterThan(clut.InputChannelCount, 0, nameof(clut.InputChannelCount));
+ Guard.MustBeGreaterThan(clut.OutputChannelCount, 0, nameof(clut.OutputChannelCount));
+
+ this.inputCount = clut.InputChannelCount;
+ this.outputCount = clut.OutputChannelCount;
+ this.g = new float[this.inputCount];
+ this.ig = new uint[this.inputCount];
+ this.s = new float[this.inputCount];
+ this.nPower = new uint[16];
+ this.lut = clut.Values;
+ this.nodeCount = (int)Math.Pow(2, clut.InputChannelCount);
+ this.df = new float[this.nodeCount];
+ this.nodes = new float[this.nodeCount][];
+ this.dimSize = new int[this.inputCount];
+ this.gridPointCount = clut.GridPointCount;
+ this.maxGridPoint = new byte[this.inputCount];
+ for (int i = 0; i < this.inputCount; i++)
+ {
+ this.maxGridPoint[i] = (byte)(this.gridPointCount[i] - 1);
+ }
+
+ this.dimSize[this.inputCount - 1] = this.outputCount;
+ for (int i = this.inputCount - 2; i >= 0; i--)
+ {
+ this.dimSize[i] = this.dimSize[i + 1] * this.gridPointCount[i + 1];
+ }
+
+ this.indexFactor = this.CalculateIndexFactor();
+ }
+
+ public unsafe Vector4 Calculate(Vector4 value)
+ {
+ Vector4 result = default;
+ switch (this.inputCount)
+ {
+ case 1:
+ this.Interpolate1d((float*)&value, (float*)&result);
+ break;
+ case 2:
+ this.Interpolate2d((float*)&value, (float*)&result);
+ break;
+ case 3:
+ this.Interpolate3d((float*)&value, (float*)&result);
+ break;
+ case 4:
+ this.Interpolate4d((float*)&value, (float*)&result);
+ break;
+ default:
+ this.InterpolateNd((float*)&value, (float*)&result);
+ break;
+ }
+
+ return result;
+ }
+
+ private int[] CalculateIndexFactor()
+ {
+ int[] factors = new int[16];
+ switch (this.inputCount)
+ {
+ case 1:
+ factors[0] = this.n000 = 0;
+ factors[1] = this.n001 = this.dimSize[0];
+ break;
+ case 2:
+ factors[0] = this.n000 = 0;
+ factors[1] = this.n001 = this.dimSize[0];
+ factors[2] = this.n010 = this.dimSize[1];
+ factors[3] = this.n011 = this.n001 + this.n010;
+ break;
+ case 3:
+ factors[0] = this.n000 = 0;
+ factors[1] = this.n001 = this.dimSize[0];
+ factors[2] = this.n010 = this.dimSize[1];
+ factors[3] = this.n011 = this.n001 + this.n010;
+ factors[4] = this.n100 = this.dimSize[2];
+ factors[5] = this.n101 = this.n100 + this.n001;
+ factors[6] = this.n110 = this.n100 + this.n010;
+ factors[7] = this.n111 = this.n110 + this.n001;
+ break;
+ case 4:
+ factors[0] = 0;
+ factors[1] = this.n001 = this.dimSize[0];
+ factors[2] = this.n010 = this.dimSize[1];
+ factors[3] = factors[2] + factors[1];
+ factors[4] = this.n100 = this.dimSize[2];
+ factors[5] = factors[4] + factors[1];
+ factors[6] = factors[4] + factors[2];
+ factors[7] = factors[4] + factors[3];
+ factors[8] = this.n1000 = this.dimSize[3];
+ factors[9] = factors[8] + factors[1];
+ factors[10] = factors[8] + factors[2];
+ factors[11] = factors[8] + factors[3];
+ factors[12] = factors[8] + factors[4];
+ factors[13] = factors[8] + factors[5];
+ factors[14] = factors[8] + factors[6];
+ factors[15] = factors[8] + factors[7];
+ break;
+ default:
+ // Initialize ND interpolation variables.
+ factors[0] = 0;
+ int count;
+ for (count = 0; count < this.inputCount; count++)
+ {
+ this.nPower[count] = (uint)(1 << (this.inputCount - 1 - count));
+ }
+
+ uint[] nPower = [0, 1];
+ count = 0;
+ int nFlag = 1;
+ for (uint j = 1; j < this.nodeCount; j++)
+ {
+ if (j == nPower[1])
+ {
+ factors[j] = this.dimSize[count];
+ nPower[0] = (uint)(1 << count);
+ count++;
+ nPower[1] = (uint)(1 << count);
+ nFlag = 1;
+ }
+ else
+ {
+ factors[j] = factors[nPower[0]] + factors[nFlag];
+ nFlag++;
+ }
+ }
+
+ break;
+ }
+
+ return factors;
+ }
+
+ ///
+ /// One dimensional interpolation function.
+ ///
+ /// The input pixel values, which will be interpolated.
+ /// The interpolated output pixels.
+ private unsafe void Interpolate1d(float* srcPixel, float* destPixel)
+ {
+ byte mx = this.maxGridPoint[0];
+
+ float x = UnitClip(srcPixel[0]) * mx;
+
+ uint ix = (uint)x;
+
+ float u = x - ix;
+
+ if (ix == mx)
+ {
+ ix--;
+ u = 1.0f;
+ }
+
+ float nu = (float)(1.0 - u);
+
+ int i;
+ Span p = this.lut.AsSpan((int)(ix * this.n001));
+
+ // Normalize grid units.
+ float dF0 = nu;
+ float dF1 = u;
+
+ int offset = 0;
+ for (i = 0; i < this.outputCount; i++)
+ {
+ destPixel[i] = (float)((p[offset + this.n000] * dF0) + (p[offset + this.n001] * dF1));
+ offset++;
+ }
+ }
+
+ ///
+ /// Two dimensional interpolation function.
+ ///
+ /// The input pixel values, which will be interpolated.
+ /// The interpolated output pixels.
+ private unsafe void Interpolate2d(float* srcPixel, float* destPixel)
+ {
+ byte mx = this.maxGridPoint[0];
+ byte my = this.maxGridPoint[1];
+
+ float x = UnitClip(srcPixel[0]) * mx;
+ float y = UnitClip(srcPixel[1]) * my;
+
+ uint ix = (uint)x;
+ uint iy = (uint)y;
+
+ float u = x - ix;
+ float t = y - iy;
+
+ if (ix == mx)
+ {
+ ix--;
+ u = 1.0f;
+ }
+
+ if (iy == my)
+ {
+ iy--;
+ t = 1.0f;
+ }
+
+ float nt = (float)(1.0 - t);
+ float nu = (float)(1.0 - u);
+
+ int i;
+ Span p = this.lut.AsSpan((int)((ix * this.n001) + (iy * this.n010)));
+
+ // Normalize grid units.
+ float dF0 = nt * nu;
+ float dF1 = nt * u;
+ float dF2 = t * nu;
+ float dF3 = t * u;
+
+ int offset = 0;
+ for (i = 0; i < this.outputCount; i++)
+ {
+ destPixel[i] = (float)((p[offset + this.n000] * dF0) + (p[offset + this.n001] * dF1) + (p[offset + this.n010] * dF2) + (p[offset + this.n011] * dF3));
+ offset++;
+ }
+ }
+
+ ///
+ /// Three dimensional interpolation function.
+ ///
+ /// The input pixel values, which will be interpolated.
+ /// The interpolated output pixels.
+ private unsafe void Interpolate3d(float* srcPixel, float* destPixel)
+ {
+ byte mx = this.maxGridPoint[0];
+ byte my = this.maxGridPoint[1];
+ byte mz = this.maxGridPoint[2];
+
+ float x = UnitClip(srcPixel[0]) * mx;
+ float y = UnitClip(srcPixel[1]) * my;
+ float z = UnitClip(srcPixel[2]) * mz;
+
+ uint ix = (uint)x;
+ uint iy = (uint)y;
+ uint iz = (uint)z;
+
+ float u = x - ix;
+ float t = y - iy;
+ float s = z - iz;
+
+ if (ix == mx)
+ {
+ ix--;
+ u = 1.0f;
+ }
+
+ if (iy == my)
+ {
+ iy--;
+ t = 1.0f;
+ }
+
+ if (iz == mz)
+ {
+ iz--;
+ s = 1.0f;
+ }
+
+ float ns = (float)(1.0 - s);
+ float nt = (float)(1.0 - t);
+ float nu = (float)(1.0 - u);
+
+ Span p = this.lut.AsSpan((int)((ix * this.n001) + (iy * this.n010) + (iz * this.n100)));
+
+ // Normalize grid units
+ float dF0 = ns * nt * nu;
+ float dF1 = ns * nt * u;
+ float dF2 = ns * t * nu;
+ float dF3 = ns * t * u;
+ float dF4 = s * nt * nu;
+ float dF5 = s * nt * u;
+ float dF6 = s * t * nu;
+ float dF7 = s * t * u;
+
+ int offset = 0;
+ for (int i = 0; i < this.outputCount; i++)
+ {
+ destPixel[i] = (float)((p[offset + this.n000] * dF0) +
+ (p[offset + this.n001] * dF1) +
+ (p[offset + this.n010] * dF2) +
+ (p[offset + this.n011] * dF3) +
+ (p[offset + this.n100] * dF4) +
+ (p[offset + this.n101] * dF5) +
+ (p[offset + this.n110] * dF6) +
+ (p[offset + this.n111] * dF7));
+ offset++;
+ }
+ }
+
+ ///
+ /// Four dimensional interpolation function.
+ ///
+ /// The input pixel values, which will be interpolated.
+ /// The interpolated output pixels.
+ private unsafe void Interpolate4d(float* srcPixel, float* destPixel)
+ {
+ byte mw = this.maxGridPoint[0];
+ byte mx = this.maxGridPoint[1];
+ byte my = this.maxGridPoint[2];
+ byte mz = this.maxGridPoint[3];
+
+ float w = UnitClip(srcPixel[0]) * mw;
+ float x = UnitClip(srcPixel[1]) * mx;
+ float y = UnitClip(srcPixel[2]) * my;
+ float z = UnitClip(srcPixel[3]) * mz;
+
+ uint iw = (uint)w;
+ uint ix = (uint)x;
+ uint iy = (uint)y;
+ uint iz = (uint)z;
+
+ float v = w - iw;
+ float u = x - ix;
+ float t = y - iy;
+ float s = z - iz;
+
+ if (iw == mw)
+ {
+ iw--;
+ v = 1.0f;
+ }
+
+ if (ix == mx)
+ {
+ ix--;
+ u = 1.0f;
+ }
+
+ if (iy == my)
+ {
+ iy--;
+ t = 1.0f;
+ }
+
+ if (iz == mz)
+ {
+ iz--;
+ s = 1.0f;
+ }
+
+ float ns = (float)(1.0 - s);
+ float nt = (float)(1.0 - t);
+ float nu = (float)(1.0 - u);
+ float nv = (float)(1.0 - v);
+
+ Span p = this.lut.AsSpan((int)((iw * this.n001) + (ix * this.n010) + (iy * this.n100) + (iz * this.n1000)));
+
+ // Normalize grid units.
+ float[] dF =
+ [
+ ns * nt * nu * nv,
+ ns * nt * nu * v,
+ ns * nt * u * nv,
+ ns * nt * u * v,
+ ns * t * nu * nv,
+ ns * t * nu * v,
+ ns * t * u * nv,
+ ns * t * u * v,
+ s * nt * nu * nv,
+ s * nt * nu * v,
+ s * nt * u * nv,
+ s * nt * u * v,
+ s * t * nu * nv,
+ s * t * nu * v,
+ s * t * u * nv,
+ s * t * u * v,
+ ];
+
+ int offset = 0;
+ for (int i = 0; i < this.outputCount; i++)
+ {
+ float pv = 0.0f;
+ for (int j = 0; j < 16; j++)
+ {
+ pv += p[offset + this.indexFactor[j]] * dF[j];
+ }
+
+ destPixel[i] = pv;
+ offset++;
+ }
+ }
+
+ ///
+ /// Generic N-dimensional interpolation function.
+ ///
+ /// The input pixel values, which will be interpolated.
+ /// The interpolated output pixels.
+ private unsafe void InterpolateNd(float* srcPixel, float* destPixel)
+ {
+ int index = 0;
+ for (int i = 0; i < this.inputCount; i++)
+ {
+ this.g[i] = UnitClip(srcPixel[i]) * this.maxGridPoint[i];
+ this.ig[i] = (uint)this.g[i];
+ this.s[this.inputCount - 1 - i] = this.g[i] - this.ig[i];
+ if (this.ig[i] == this.maxGridPoint[i])
+ {
+ this.ig[i]--;
+ this.s[this.inputCount - 1 - i] = 1.0f;
+ }
+
+ index += (int)this.ig[i] * this.dimSize[i];
+ }
+
+ Span p = this.lut.AsSpan(index);
+ float[] temp = new float[2];
+ bool nFlag = false;
+
+ for (int i = 0; i < this.nodeCount; i++)
+ {
+ this.df[i] = 1.0f;
+ }
+
+ for (int i = 0; i < this.inputCount; i++)
+ {
+ temp[0] = 1.0f - this.s[i];
+ temp[1] = this.s[i];
+ index = (int)this.nPower[i];
+ for (int j = 0; j < this.nodeCount; j++)
+ {
+ this.df[j] *= temp[nFlag ? 1 : 0];
+ if ((j + 1) % index == 0)
+ {
+ nFlag = !nFlag;
+ }
+ }
+
+ nFlag = false;
+ }
+
+ int offset = 0;
+ for (int i = 0; i < this.outputCount; i++)
+ {
+ float pv = 0;
+ for (int j = 0; j < this.nodeCount; j++)
+ {
+ pv += p[offset + this.indexFactor[j]] * this.df[j];
+ }
+
+ destPixel[i] = pv;
+ offset++;
+ }
+ }
+
+ private static float UnitClip(float v)
+ {
+ if (v < 0)
+ {
+ return 0;
+ }
+
+ if (v > 1.0)
+ {
+ return 1.0f;
+ }
+
+ return v;
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/Calculators/ColorTrcCalculator.cs b/src/ImageSharp/ColorProfiles/Icc/Calculators/ColorTrcCalculator.cs
new file mode 100644
index 0000000000..3604642c95
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/Calculators/ColorTrcCalculator.cs
@@ -0,0 +1,65 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+using System.Numerics;
+using System.Runtime.CompilerServices;
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Icc.Calculators;
+
+internal class ColorTrcCalculator : IVector4Calculator
+{
+ private readonly TrcCalculator curveCalculator;
+ private readonly Matrix4x4 matrix;
+ private readonly bool toPcs;
+
+ public ColorTrcCalculator(
+ IccXyzTagDataEntry redMatrixColumn,
+ IccXyzTagDataEntry greenMatrixColumn,
+ IccXyzTagDataEntry blueMatrixColumn,
+ IccTagDataEntry redTrc,
+ IccTagDataEntry greenTrc,
+ IccTagDataEntry blueTrc,
+ bool toPcs)
+ {
+ this.toPcs = toPcs;
+ this.curveCalculator = new TrcCalculator([redTrc, greenTrc, blueTrc], !toPcs);
+
+ Vector3 mr = redMatrixColumn.Data[0];
+ Vector3 mg = greenMatrixColumn.Data[0];
+ Vector3 mb = blueMatrixColumn.Data[0];
+ this.matrix = new Matrix4x4(mr.X, mr.Y, mr.Z, 0, mg.X, mg.Y, mg.Z, 0, mb.X, mb.Y, mb.Z, 0, 0, 0, 0, 1);
+
+ if (!toPcs)
+ {
+ Matrix4x4.Invert(this.matrix, out this.matrix);
+ }
+ }
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ public Vector4 Calculate(Vector4 value)
+ {
+ if (this.toPcs)
+ {
+ // input is always linear RGB
+ value = this.curveCalculator.Calculate(value);
+ CieXyz xyz = new(Vector4.Transform(value, this.matrix).AsVector3());
+
+ // when data to PCS, output from calculator is descaled XYZ
+ // but downstream process requires scaled XYZ
+ // (see DemoMaxICC IccCmm.cpp : CIccXformMatrixTRC::Apply)
+ return xyz.ToScaledVector4();
+ }
+ else
+ {
+ // input is always XYZ
+ Vector4 xyz = Vector4.Transform(value, this.matrix);
+
+ // when data to PCS, upstream process provides scaled XYZ
+ // but input to calculator is descaled XYZ
+ // (see DemoMaxICC IccCmm.cpp : CIccXformMatrixTRC::Apply)
+ xyz = new(CieXyz.FromScaledVector4(xyz).AsVector3Unsafe(), 1);
+ return this.curveCalculator.Calculate(xyz);
+ }
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/Calculators/CurveCalculator.CalculationType.cs b/src/ImageSharp/ColorProfiles/Icc/Calculators/CurveCalculator.CalculationType.cs
new file mode 100644
index 0000000000..d035bd1793
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/Calculators/CurveCalculator.CalculationType.cs
@@ -0,0 +1,14 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+namespace SixLabors.ImageSharp.ColorProfiles.Conversion.Icc;
+
+internal partial class CurveCalculator
+{
+ private enum CalculationType
+ {
+ Identity,
+ Gamma,
+ Lut,
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/Calculators/CurveCalculator.cs b/src/ImageSharp/ColorProfiles/Icc/Calculators/CurveCalculator.cs
new file mode 100644
index 0000000000..c39eaf958f
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/Calculators/CurveCalculator.cs
@@ -0,0 +1,47 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+#nullable disable
+
+using SixLabors.ImageSharp.ColorProfiles.Icc.Calculators;
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Conversion.Icc;
+
+internal partial class CurveCalculator : ISingleCalculator
+{
+ private readonly LutCalculator lutCalculator;
+ private readonly float gamma;
+ private readonly CalculationType type;
+
+ public CurveCalculator(IccCurveTagDataEntry entry, bool inverted)
+ {
+ if (entry.IsIdentityResponse)
+ {
+ this.type = CalculationType.Identity;
+ }
+ else if (entry.IsGamma)
+ {
+ this.gamma = entry.Gamma;
+ if (inverted)
+ {
+ this.gamma = 1f / this.gamma;
+ }
+
+ this.type = CalculationType.Gamma;
+ }
+ else
+ {
+ this.lutCalculator = new LutCalculator(entry.CurveData, inverted);
+ this.type = CalculationType.Lut;
+ }
+ }
+
+ public float Calculate(float value)
+ => this.type switch
+ {
+ CalculationType.Identity => value,
+ CalculationType.Gamma => MathF.Pow(value, this.gamma), // TODO: This could be optimized using a LUT. See SrgbCompanding
+ CalculationType.Lut => this.lutCalculator.Calculate(value),
+ _ => throw new InvalidOperationException("Invalid calculation type"),
+ };
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/Calculators/GrayTrcCalculator.cs b/src/ImageSharp/ColorProfiles/Icc/Calculators/GrayTrcCalculator.cs
new file mode 100644
index 0000000000..8d823c1e95
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/Calculators/GrayTrcCalculator.cs
@@ -0,0 +1,19 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+using System.Numerics;
+using System.Runtime.CompilerServices;
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Icc.Calculators;
+
+internal class GrayTrcCalculator : IVector4Calculator
+{
+ private readonly TrcCalculator calculator;
+
+ public GrayTrcCalculator(IccTagDataEntry grayTrc, bool toPcs)
+ => this.calculator = new TrcCalculator(new IccTagDataEntry[] { grayTrc }, !toPcs);
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ public Vector4 Calculate(Vector4 value) => this.calculator.Calculate(value);
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/Calculators/ISingleCalculator.cs b/src/ImageSharp/ColorProfiles/Icc/Calculators/ISingleCalculator.cs
new file mode 100644
index 0000000000..ce9b7d2f9b
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/Calculators/ISingleCalculator.cs
@@ -0,0 +1,17 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+namespace SixLabors.ImageSharp.ColorProfiles.Icc.Calculators;
+
+///
+/// Represents an ICC calculator with a single floating point value and result
+///
+internal interface ISingleCalculator
+{
+ ///
+ /// Calculates a result from the given value
+ ///
+ /// The input value
+ /// The calculated result
+ float Calculate(float value);
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/Calculators/IVector4Calculator.cs b/src/ImageSharp/ColorProfiles/Icc/Calculators/IVector4Calculator.cs
new file mode 100644
index 0000000000..9beea79503
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/Calculators/IVector4Calculator.cs
@@ -0,0 +1,19 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+using System.Numerics;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Icc.Calculators;
+
+///
+/// Represents an ICC calculator with values and results
+///
+internal interface IVector4Calculator
+{
+ ///
+ /// Calculates a result from the given values
+ ///
+ /// The input values
+ /// The calculated result
+ Vector4 Calculate(Vector4 value);
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/Calculators/LutABCalculator.CalculationType.cs b/src/ImageSharp/ColorProfiles/Icc/Calculators/LutABCalculator.CalculationType.cs
new file mode 100644
index 0000000000..253239cb79
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/Calculators/LutABCalculator.CalculationType.cs
@@ -0,0 +1,18 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+namespace SixLabors.ImageSharp.ColorProfiles.Conversion.Icc;
+
+internal partial class LutABCalculator
+{
+ private enum CalculationType
+ {
+ AtoB = 1 << 3,
+ BtoA = 1 << 4,
+
+ SingleCurve = 1,
+ CurveMatrix = 2,
+ CurveClut = 3,
+ Full = 4,
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/Calculators/LutABCalculator.cs b/src/ImageSharp/ColorProfiles/Icc/Calculators/LutABCalculator.cs
new file mode 100644
index 0000000000..172d806394
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/Calculators/LutABCalculator.cs
@@ -0,0 +1,135 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+#nullable disable
+
+using System.Numerics;
+using SixLabors.ImageSharp.ColorProfiles.Icc.Calculators;
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Conversion.Icc;
+
+internal partial class LutABCalculator : IVector4Calculator
+{
+ private CalculationType type;
+ private TrcCalculator curveACalculator;
+ private TrcCalculator curveBCalculator;
+ private TrcCalculator curveMCalculator;
+ private MatrixCalculator matrixCalculator;
+ private ClutCalculator clutCalculator;
+
+ public LutABCalculator(IccLutAToBTagDataEntry entry)
+ {
+ Guard.NotNull(entry, nameof(entry));
+ this.Init(entry.CurveA, entry.CurveB, entry.CurveM, entry.Matrix3x1, entry.Matrix3x3, entry.ClutValues);
+ this.type |= CalculationType.AtoB;
+ }
+
+ public LutABCalculator(IccLutBToATagDataEntry entry)
+ {
+ Guard.NotNull(entry, nameof(entry));
+ this.Init(entry.CurveA, entry.CurveB, entry.CurveM, entry.Matrix3x1, entry.Matrix3x3, entry.ClutValues);
+ this.type |= CalculationType.BtoA;
+ }
+
+ public Vector4 Calculate(Vector4 value)
+ {
+ switch (this.type)
+ {
+ case CalculationType.Full | CalculationType.AtoB:
+ value = this.curveACalculator.Calculate(value);
+ value = this.clutCalculator.Calculate(value);
+ value = this.curveMCalculator.Calculate(value);
+ value = this.matrixCalculator.Calculate(value);
+ return this.curveBCalculator.Calculate(value);
+
+ case CalculationType.Full | CalculationType.BtoA:
+ value = this.curveBCalculator.Calculate(value);
+ value = this.matrixCalculator.Calculate(value);
+ value = this.curveMCalculator.Calculate(value);
+ value = this.clutCalculator.Calculate(value);
+ return this.curveACalculator.Calculate(value);
+
+ case CalculationType.CurveClut | CalculationType.AtoB:
+ value = this.curveACalculator.Calculate(value);
+ value = this.clutCalculator.Calculate(value);
+ return this.curveBCalculator.Calculate(value);
+
+ case CalculationType.CurveClut | CalculationType.BtoA:
+ value = this.curveBCalculator.Calculate(value);
+ value = this.clutCalculator.Calculate(value);
+ return this.curveACalculator.Calculate(value);
+
+ case CalculationType.CurveMatrix | CalculationType.AtoB:
+ value = this.curveMCalculator.Calculate(value);
+ value = this.matrixCalculator.Calculate(value);
+ return this.curveBCalculator.Calculate(value);
+
+ case CalculationType.CurveMatrix | CalculationType.BtoA:
+ value = this.curveBCalculator.Calculate(value);
+ value = this.matrixCalculator.Calculate(value);
+ return this.curveMCalculator.Calculate(value);
+
+ case CalculationType.SingleCurve | CalculationType.AtoB:
+ case CalculationType.SingleCurve | CalculationType.BtoA:
+ return this.curveBCalculator.Calculate(value);
+
+ default:
+ throw new InvalidOperationException("Invalid calculation type");
+ }
+ }
+
+ private void Init(IccTagDataEntry[] curveA, IccTagDataEntry[] curveB, IccTagDataEntry[] curveM, Vector3? matrix3x1, Matrix4x4? matrix3x3, IccClut clut)
+ {
+ bool hasACurve = curveA != null;
+ bool hasBCurve = curveB != null;
+ bool hasMCurve = curveM != null;
+ bool hasMatrix = matrix3x1 != null && matrix3x3 != null;
+ bool hasClut = clut != null;
+
+ if (hasBCurve && hasMatrix && hasMCurve && hasClut && hasACurve)
+ {
+ this.type = CalculationType.Full;
+ }
+ else if (hasBCurve && hasClut && hasACurve)
+ {
+ this.type = CalculationType.CurveClut;
+ }
+ else if (hasBCurve && hasMatrix && hasMCurve)
+ {
+ this.type = CalculationType.CurveMatrix;
+ }
+ else if (hasBCurve)
+ {
+ this.type = CalculationType.SingleCurve;
+ }
+ else
+ {
+ throw new InvalidIccProfileException("AToB or BToA tag has an invalid configuration");
+ }
+
+ if (hasACurve)
+ {
+ this.curveACalculator = new TrcCalculator(curveA, false);
+ }
+
+ if (hasBCurve)
+ {
+ this.curveBCalculator = new TrcCalculator(curveB, false);
+ }
+
+ if (hasMCurve)
+ {
+ this.curveMCalculator = new TrcCalculator(curveM, false);
+ }
+
+ if (hasMatrix)
+ {
+ this.matrixCalculator = new MatrixCalculator(matrix3x3.Value, matrix3x1.Value);
+ }
+
+ if (hasClut)
+ {
+ this.clutCalculator = new ClutCalculator(clut);
+ }
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/Calculators/LutCalculator.cs b/src/ImageSharp/ColorProfiles/Icc/Calculators/LutCalculator.cs
new file mode 100644
index 0000000000..83704ae214
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/Calculators/LutCalculator.cs
@@ -0,0 +1,77 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+using System.Runtime.CompilerServices;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Icc.Calculators;
+
+internal class LutCalculator : ISingleCalculator
+{
+ private readonly float[] lut;
+ private readonly bool inverse;
+
+ public LutCalculator(float[] lut, bool inverse)
+ {
+ Guard.NotNull(lut, nameof(lut));
+
+ this.lut = lut;
+ this.inverse = inverse;
+ }
+
+ public float Calculate(float value)
+ {
+ if (this.inverse)
+ {
+ return this.LookupInverse(value);
+ }
+
+ return this.Lookup(value);
+ }
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ private float Lookup(float value)
+ {
+ value = Math.Max(value, 0);
+
+ float factor = value * (this.lut.Length - 1);
+ int index = (int)factor;
+ float low = this.lut[index];
+
+ float high = 1F;
+ if (index < this.lut.Length - 1)
+ {
+ high = this.lut[index + 1];
+ }
+
+ return low + ((high - low) * (factor - index));
+ }
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ private float LookupInverse(float value)
+ {
+ int index = Array.BinarySearch(this.lut, value);
+ if (index >= 0)
+ {
+ return index / (float)(this.lut.Length - 1);
+ }
+
+ index = ~index;
+ if (index == 0)
+ {
+ return 0;
+ }
+ else if (index == this.lut.Length)
+ {
+ return 1;
+ }
+
+ float high = this.lut[index];
+ float low = this.lut[index - 1];
+
+ float valuePercent = (value - low) / (high - low);
+ float lutRange = 1 / (float)(this.lut.Length - 1);
+ float lutLow = (index - 1) / (float)(this.lut.Length - 1);
+
+ return lutLow + (valuePercent * lutRange);
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/Calculators/LutEntryCalculator.cs b/src/ImageSharp/ColorProfiles/Icc/Calculators/LutEntryCalculator.cs
new file mode 100644
index 0000000000..c97578ee3f
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/Calculators/LutEntryCalculator.cs
@@ -0,0 +1,80 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+#nullable disable
+
+using System.Numerics;
+using System.Runtime.CompilerServices;
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Icc.Calculators;
+
+internal class LutEntryCalculator : IVector4Calculator
+{
+ private LutCalculator[] inputCurve;
+ private LutCalculator[] outputCurve;
+ private ClutCalculator clutCalculator;
+ private Matrix4x4 matrix;
+ private bool doTransform;
+
+ public LutEntryCalculator(IccLut8TagDataEntry lut)
+ {
+ Guard.NotNull(lut, nameof(lut));
+ this.Init(lut.InputValues, lut.OutputValues, lut.ClutValues, lut.Matrix);
+ this.Is16Bit = false;
+ }
+
+ public LutEntryCalculator(IccLut16TagDataEntry lut)
+ {
+ Guard.NotNull(lut, nameof(lut));
+ this.Init(lut.InputValues, lut.OutputValues, lut.ClutValues, lut.Matrix);
+ this.Is16Bit = true;
+ }
+
+ internal bool Is16Bit { get; }
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ public Vector4 Calculate(Vector4 value)
+ {
+ if (this.doTransform)
+ {
+ value = Vector4.Transform(value, this.matrix);
+ }
+
+ value = CalculateLut(this.inputCurve, value);
+ value = this.clutCalculator.Calculate(value);
+ return CalculateLut(this.outputCurve, value);
+ }
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ private static Vector4 CalculateLut(LutCalculator[] lut, Vector4 value)
+ {
+ ref float f = ref Unsafe.As(ref value);
+ for (int i = 0; i < lut.Length; i++)
+ {
+ Unsafe.Add(ref f, i) = lut[i].Calculate(Unsafe.Add(ref f, i));
+ }
+
+ return value;
+ }
+
+ private void Init(IccLut[] inputCurve, IccLut[] outputCurve, IccClut clut, Matrix4x4 matrix)
+ {
+ this.inputCurve = InitLut(inputCurve);
+ this.outputCurve = InitLut(outputCurve);
+ this.clutCalculator = new ClutCalculator(clut);
+ this.matrix = matrix;
+
+ this.doTransform = !matrix.IsIdentity && inputCurve.Length == 3;
+ }
+
+ private static LutCalculator[] InitLut(IccLut[] curves)
+ {
+ LutCalculator[] calculators = new LutCalculator[curves.Length];
+ for (int i = 0; i < curves.Length; i++)
+ {
+ calculators[i] = new LutCalculator(curves[i].Values, false);
+ }
+
+ return calculators;
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/Calculators/MatrixCalculator.cs b/src/ImageSharp/ColorProfiles/Icc/Calculators/MatrixCalculator.cs
new file mode 100644
index 0000000000..6be1fdbf95
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/Calculators/MatrixCalculator.cs
@@ -0,0 +1,26 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+using System.Numerics;
+using System.Runtime.CompilerServices;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Icc.Calculators;
+
+internal class MatrixCalculator : IVector4Calculator
+{
+ private Matrix4x4 matrix2D;
+ private Vector4 matrix1D;
+
+ public MatrixCalculator(Matrix4x4 matrix3x3, Vector3 matrix3x1)
+ {
+ this.matrix2D = matrix3x3;
+ this.matrix1D = new Vector4(matrix3x1, 0);
+ }
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ public Vector4 Calculate(Vector4 value)
+ {
+ Vector4 transformed = Vector4.Transform(value, this.matrix2D);
+ return Vector4.Add(this.matrix1D, transformed);
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/Calculators/ParametricCurveCalculator.cs b/src/ImageSharp/ColorProfiles/Icc/Calculators/ParametricCurveCalculator.cs
new file mode 100644
index 0000000000..2a3945e270
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/Calculators/ParametricCurveCalculator.cs
@@ -0,0 +1,130 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+using System.Runtime.CompilerServices;
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Icc.Calculators;
+
+internal class ParametricCurveCalculator : ISingleCalculator
+{
+ private readonly IccParametricCurve curve;
+ private readonly IccParametricCurveType type;
+ private const IccParametricCurveType InvertedFlag = (IccParametricCurveType)(1 << 3);
+
+ public ParametricCurveCalculator(IccParametricCurveTagDataEntry entry, bool inverted)
+ {
+ Guard.NotNull(entry, nameof(entry));
+ this.curve = entry.Curve;
+ this.type = entry.Curve.Type;
+
+ if (inverted)
+ {
+ this.type |= InvertedFlag;
+ }
+ }
+
+ public float Calculate(float value)
+ => this.type switch
+ {
+ IccParametricCurveType.Type1 => this.CalculateGamma(value),
+ IccParametricCurveType.Cie122_1996 => this.CalculateCie122(value),
+ IccParametricCurveType.Iec61966_3 => this.CalculateIec61966(value),
+ IccParametricCurveType.SRgb => this.CalculateSRgb(value),
+ IccParametricCurveType.Type5 => this.CalculateType5(value),
+ IccParametricCurveType.Type1 | InvertedFlag => this.CalculateInvertedGamma(value),
+ IccParametricCurveType.Cie122_1996 | InvertedFlag => this.CalculateInvertedCie122(value),
+ IccParametricCurveType.Iec61966_3 | InvertedFlag => this.CalculateInvertedIec61966(value),
+ IccParametricCurveType.SRgb | InvertedFlag => this.CalculateInvertedSRgb(value),
+ IccParametricCurveType.Type5 | InvertedFlag => this.CalculateInvertedType5(value),
+ _ => throw new InvalidIccProfileException("ParametricCurve"),
+ };
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ private float CalculateGamma(float value) => MathF.Pow(value, this.curve.G);
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ private float CalculateCie122(float value)
+ {
+ if (value >= -this.curve.B / this.curve.A)
+ {
+ return MathF.Pow((this.curve.A * value) + this.curve.B, this.curve.G);
+ }
+
+ return 0;
+ }
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ private float CalculateIec61966(float value)
+ {
+ if (value >= -this.curve.B / this.curve.A)
+ {
+ return MathF.Pow((this.curve.A * value) + this.curve.B, this.curve.G) + this.curve.C;
+ }
+
+ return this.curve.C;
+ }
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ private float CalculateSRgb(float value)
+ {
+ if (value >= this.curve.D)
+ {
+ return MathF.Pow((this.curve.A * value) + this.curve.B, this.curve.G);
+ }
+
+ return this.curve.C * value;
+ }
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ private float CalculateType5(float value)
+ {
+ if (value >= this.curve.D)
+ {
+ return MathF.Pow((this.curve.A * value) + this.curve.B, this.curve.G) + this.curve.E;
+ }
+
+ return (this.curve.C * value) + this.curve.F;
+ }
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ private float CalculateInvertedGamma(float value)
+ => MathF.Pow(value, 1 / this.curve.G);
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ private float CalculateInvertedCie122(float value)
+ => (MathF.Pow(value, 1 / this.curve.G) - this.curve.B) / this.curve.A;
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ private float CalculateInvertedIec61966(float value)
+ {
+ if (value >= this.curve.C)
+ {
+ return (MathF.Pow(value - this.curve.C, 1 / this.curve.G) - this.curve.B) / this.curve.A;
+ }
+
+ return -this.curve.B / this.curve.A;
+ }
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ private float CalculateInvertedSRgb(float value)
+ {
+ if (value >= MathF.Pow((this.curve.A * this.curve.D) + this.curve.B, this.curve.G))
+ {
+ return (MathF.Pow(value, 1 / this.curve.G) - this.curve.B) / this.curve.A;
+ }
+
+ return value / this.curve.C;
+ }
+
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ private float CalculateInvertedType5(float value)
+ {
+ if (value >= (this.curve.C * this.curve.D) + this.curve.F)
+ {
+ return (MathF.Pow(value - this.curve.E, 1 / this.curve.G) - this.curve.B) / this.curve.A;
+ }
+
+ return (value - this.curve.F) / this.curve.C;
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/Calculators/TrcCalculator.cs b/src/ImageSharp/ColorProfiles/Icc/Calculators/TrcCalculator.cs
new file mode 100644
index 0000000000..d2fc5d9b55
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/Calculators/TrcCalculator.cs
@@ -0,0 +1,41 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+using System.Numerics;
+using System.Runtime.CompilerServices;
+using SixLabors.ImageSharp.ColorProfiles.Conversion.Icc;
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Icc.Calculators;
+
+internal class TrcCalculator : IVector4Calculator
+{
+ private readonly ISingleCalculator[] calculators;
+
+ public TrcCalculator(IccTagDataEntry[] entries, bool inverted)
+ {
+ Guard.NotNull(entries, nameof(entries));
+
+ this.calculators = new ISingleCalculator[entries.Length];
+ for (int i = 0; i < entries.Length; i++)
+ {
+ this.calculators[i] = entries[i] switch
+ {
+ IccCurveTagDataEntry curve => new CurveCalculator(curve, inverted),
+ IccParametricCurveTagDataEntry parametricCurve => new ParametricCurveCalculator(parametricCurve, inverted),
+ _ => throw new InvalidIccProfileException("Invalid Entry."),
+ };
+ }
+ }
+
+ public unsafe Vector4 Calculate(Vector4 value)
+ {
+ ref float f = ref Unsafe.As(ref value);
+ for (int i = 0; i < this.calculators.Length; i++)
+ {
+ Unsafe.Add(ref f, i) = this.calculators[i].Calculate(Unsafe.Add(ref f, i));
+ }
+
+ return value;
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/IccConverterBase.Checks.cs b/src/ImageSharp/ColorProfiles/Icc/IccConverterBase.Checks.cs
new file mode 100644
index 0000000000..94f906709a
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/IccConverterBase.Checks.cs
@@ -0,0 +1,155 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+#nullable disable
+
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Conversion.Icc;
+
+///
+/// Color converter for ICC profiles
+///
+internal abstract partial class IccConverterBase
+{
+ private static ConversionMethod GetConversionMethod(IccProfile profile, IccRenderingIntent renderingIntent) => profile.Header.Class switch
+ {
+ IccProfileClass.InputDevice or
+ IccProfileClass.DisplayDevice or
+ IccProfileClass.OutputDevice or
+ IccProfileClass.ColorSpace => CheckMethod1(profile, renderingIntent),
+ IccProfileClass.DeviceLink or IccProfileClass.Abstract => CheckMethod2(profile),
+ _ => ConversionMethod.Invalid,
+ };
+
+ private static ConversionMethod CheckMethod1(IccProfile profile, IccRenderingIntent renderingIntent)
+ {
+ ConversionMethod method = CheckMethodD(profile, renderingIntent);
+ if (method != ConversionMethod.Invalid)
+ {
+ return method;
+ }
+
+ method = CheckMethodA(profile, renderingIntent);
+ if (method != ConversionMethod.Invalid)
+ {
+ return method;
+ }
+
+ method = CheckMethodA0(profile);
+ if (method != ConversionMethod.Invalid)
+ {
+ return method;
+ }
+
+ method = CheckMethodTrc(profile);
+ if (method != ConversionMethod.Invalid)
+ {
+ return method;
+ }
+
+ return ConversionMethod.Invalid;
+ }
+
+ private static ConversionMethod CheckMethodD(IccProfile profile, IccRenderingIntent renderingIntent)
+ {
+ if ((HasTag(profile, IccProfileTag.DToB0) || HasTag(profile, IccProfileTag.BToD0))
+ && renderingIntent == IccRenderingIntent.Perceptual)
+ {
+ return ConversionMethod.D0;
+ }
+
+ if ((HasTag(profile, IccProfileTag.DToB1) || HasTag(profile, IccProfileTag.BToD1))
+ && renderingIntent == IccRenderingIntent.MediaRelativeColorimetric)
+ {
+ return ConversionMethod.D1;
+ }
+
+ if ((HasTag(profile, IccProfileTag.DToB2) || HasTag(profile, IccProfileTag.BToD2))
+ && renderingIntent == IccRenderingIntent.Saturation)
+ {
+ return ConversionMethod.D2;
+ }
+
+ if ((HasTag(profile, IccProfileTag.DToB3) || HasTag(profile, IccProfileTag.BToD3))
+ && renderingIntent == IccRenderingIntent.AbsoluteColorimetric)
+ {
+ return ConversionMethod.D3;
+ }
+
+ return ConversionMethod.Invalid;
+ }
+
+ private static ConversionMethod CheckMethodA(IccProfile profile, IccRenderingIntent renderingIntent)
+ {
+ if ((HasTag(profile, IccProfileTag.AToB0) || HasTag(profile, IccProfileTag.BToA0))
+ && renderingIntent == IccRenderingIntent.Perceptual)
+ {
+ return ConversionMethod.A0;
+ }
+
+ if ((HasTag(profile, IccProfileTag.AToB1) || HasTag(profile, IccProfileTag.BToA1))
+ && renderingIntent == IccRenderingIntent.MediaRelativeColorimetric)
+ {
+ return ConversionMethod.A1;
+ }
+
+ if ((HasTag(profile, IccProfileTag.AToB2) || HasTag(profile, IccProfileTag.BToA2))
+ && renderingIntent == IccRenderingIntent.Saturation)
+ {
+ return ConversionMethod.A2;
+ }
+
+ return ConversionMethod.Invalid;
+ }
+
+ private static ConversionMethod CheckMethodA0(IccProfile profile)
+ {
+ bool valid = HasTag(profile, IccProfileTag.AToB0) || HasTag(profile, IccProfileTag.BToA0);
+ return valid ? ConversionMethod.A0 : ConversionMethod.Invalid;
+ }
+
+ private static ConversionMethod CheckMethodTrc(IccProfile profile)
+ {
+ if (HasTag(profile, IccProfileTag.RedMatrixColumn)
+ && HasTag(profile, IccProfileTag.GreenMatrixColumn)
+ && HasTag(profile, IccProfileTag.BlueMatrixColumn)
+ && HasTag(profile, IccProfileTag.RedTrc)
+ && HasTag(profile, IccProfileTag.GreenTrc)
+ && HasTag(profile, IccProfileTag.BlueTrc))
+ {
+ return ConversionMethod.ColorTrc;
+ }
+
+ if (HasTag(profile, IccProfileTag.GrayTrc))
+ {
+ return ConversionMethod.GrayTrc;
+ }
+
+ return ConversionMethod.Invalid;
+ }
+
+ private static ConversionMethod CheckMethod2(IccProfile profile)
+ {
+ if (HasTag(profile, IccProfileTag.DToB0) || HasTag(profile, IccProfileTag.BToD0))
+ {
+ return ConversionMethod.D0;
+ }
+
+ if (HasTag(profile, IccProfileTag.AToB0) || HasTag(profile, IccProfileTag.AToB0))
+ {
+ return ConversionMethod.A0;
+ }
+
+ return ConversionMethod.Invalid;
+ }
+
+ private static bool HasTag(IccProfile profile, IccProfileTag tag)
+ => profile.Entries.Any(t => t.TagSignature == tag);
+
+ private static IccTagDataEntry GetTag(IccProfile profile, IccProfileTag tag)
+ => Array.Find(profile.Entries, t => t.TagSignature == tag);
+
+ private static T GetTag(IccProfile profile, IccProfileTag tag)
+ where T : IccTagDataEntry
+ => profile.Entries.OfType().FirstOrDefault(t => t.TagSignature == tag);
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/IccConverterBase.ConversionMethod.cs b/src/ImageSharp/ColorProfiles/Icc/IccConverterBase.ConversionMethod.cs
new file mode 100644
index 0000000000..43593f0ae9
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/IccConverterBase.ConversionMethod.cs
@@ -0,0 +1,66 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+namespace SixLabors.ImageSharp.ColorProfiles.Conversion.Icc;
+
+///
+/// Color converter for ICC profiles
+///
+internal abstract partial class IccConverterBase
+{
+ ///
+ /// Conversion methods with ICC profiles
+ ///
+ private enum ConversionMethod
+ {
+ ///
+ /// Conversion using anything but Multi Process Elements with perceptual rendering intent
+ ///
+ A0,
+
+ ///
+ /// Conversion using anything but Multi Process Elements with relative colorimetric rendering intent
+ ///
+ A1,
+
+ ///
+ /// Conversion using anything but Multi Process Elements with saturation rendering intent
+ ///
+ A2,
+
+ ///
+ /// Conversion using Multi Process Elements with perceptual rendering intent
+ ///
+ D0,
+
+ ///
+ /// Conversion using Multi Process Elements with relative colorimetric rendering intent
+ ///
+ D1,
+
+ ///
+ /// Conversion using Multi Process Elements with saturation rendering intent
+ ///
+ D2,
+
+ ///
+ /// Conversion using Multi Process Elements with absolute colorimetric rendering intent
+ ///
+ D3,
+
+ ///
+ /// Conversion of more than one channel using tone reproduction curves
+ ///
+ ColorTrc,
+
+ ///
+ /// Conversion of exactly one channel using a tone reproduction curve
+ ///
+ GrayTrc,
+
+ ///
+ /// No valid conversion method available or found
+ ///
+ Invalid,
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/IccConverterbase.Conversions.cs b/src/ImageSharp/ColorProfiles/Icc/IccConverterbase.Conversions.cs
new file mode 100644
index 0000000000..20df08e378
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/IccConverterbase.Conversions.cs
@@ -0,0 +1,109 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+using SixLabors.ImageSharp.ColorProfiles.Icc.Calculators;
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Conversion.Icc;
+
+///
+/// Color converter for ICC profiles
+///
+internal abstract partial class IccConverterBase
+{
+ private IVector4Calculator calculator;
+
+ internal bool Is16BitLutEntry => this.calculator is LutEntryCalculator { Is16Bit: true };
+
+ internal bool IsTrc => this.calculator is ColorTrcCalculator or GrayTrcCalculator;
+
+ ///
+ /// Checks the profile for available conversion methods and gathers all the information's necessary for it.
+ ///
+ /// The profile to use for the conversion.
+ /// True if the conversion is to the Profile Connection Space.
+ /// The wanted rendering intent. Can be ignored if not available.
+ /// Invalid conversion method.
+ protected void Init(IccProfile profile, bool toPcs, IccRenderingIntent renderingIntent)
+ => this.calculator = GetConversionMethod(profile, renderingIntent) switch
+ {
+ ConversionMethod.D0 => toPcs ?
+ InitD(profile, IccProfileTag.DToB0) :
+ InitD(profile, IccProfileTag.BToD0),
+ ConversionMethod.D1 => toPcs ?
+ InitD(profile, IccProfileTag.DToB1) :
+ InitD(profile, IccProfileTag.BToD1),
+ ConversionMethod.D2 => toPcs ?
+ InitD(profile, IccProfileTag.DToB2) :
+ InitD(profile, IccProfileTag.BToD2),
+ ConversionMethod.D3 => toPcs ?
+ InitD(profile, IccProfileTag.DToB3) :
+ InitD(profile, IccProfileTag.BToD3),
+ ConversionMethod.A0 => toPcs ?
+ InitA(profile, IccProfileTag.AToB0) :
+ InitA(profile, IccProfileTag.BToA0),
+ ConversionMethod.A1 => toPcs ?
+ InitA(profile, IccProfileTag.AToB1) :
+ InitA(profile, IccProfileTag.BToA1),
+ ConversionMethod.A2 => toPcs ?
+ InitA(profile, IccProfileTag.AToB2) :
+ InitA(profile, IccProfileTag.BToA2),
+ ConversionMethod.ColorTrc => InitColorTrc(profile, toPcs),
+ ConversionMethod.GrayTrc => InitGrayTrc(profile, toPcs),
+ _ => throw new InvalidIccProfileException("Invalid conversion method."),
+ };
+
+ private static IVector4Calculator InitA(IccProfile profile, IccProfileTag tag)
+ => GetTag(profile, tag) switch
+ {
+ IccLut8TagDataEntry lut8 => new LutEntryCalculator(lut8),
+ IccLut16TagDataEntry lut16 => new LutEntryCalculator(lut16),
+ IccLutAToBTagDataEntry lutAtoB => new LutABCalculator(lutAtoB),
+ IccLutBToATagDataEntry lutBtoA => new LutABCalculator(lutBtoA),
+ _ => throw new InvalidIccProfileException("Invalid entry."),
+ };
+
+ private static IVector4Calculator InitD(IccProfile profile, IccProfileTag tag)
+ {
+ IccMultiProcessElementsTagDataEntry entry = GetTag(profile, tag)
+ ?? throw new InvalidIccProfileException("Entry is null.");
+
+ throw new NotImplementedException("Multi process elements are not supported");
+ }
+
+ private static ColorTrcCalculator InitColorTrc(IccProfile profile, bool toPcs)
+ {
+ IccXyzTagDataEntry redMatrixColumn = GetTag(profile, IccProfileTag.RedMatrixColumn);
+ IccXyzTagDataEntry greenMatrixColumn = GetTag(profile, IccProfileTag.GreenMatrixColumn);
+ IccXyzTagDataEntry blueMatrixColumn = GetTag(profile, IccProfileTag.BlueMatrixColumn);
+
+ IccTagDataEntry redTrc = GetTag(profile, IccProfileTag.RedTrc);
+ IccTagDataEntry greenTrc = GetTag(profile, IccProfileTag.GreenTrc);
+ IccTagDataEntry blueTrc = GetTag(profile, IccProfileTag.BlueTrc);
+
+ if (redMatrixColumn == null ||
+ greenMatrixColumn == null ||
+ blueMatrixColumn == null ||
+ redTrc == null ||
+ greenTrc == null ||
+ blueTrc == null)
+ {
+ throw new InvalidIccProfileException("Missing matrix column or channel.");
+ }
+
+ return new ColorTrcCalculator(
+ redMatrixColumn,
+ greenMatrixColumn,
+ blueMatrixColumn,
+ redTrc,
+ greenTrc,
+ blueTrc,
+ toPcs);
+ }
+
+ private static GrayTrcCalculator InitGrayTrc(IccProfile profile, bool toPcs)
+ {
+ IccTagDataEntry entry = GetTag(profile, IccProfileTag.GrayTrc);
+ return new GrayTrcCalculator(entry, toPcs);
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/IccConverterbase.cs b/src/ImageSharp/ColorProfiles/Icc/IccConverterbase.cs
new file mode 100644
index 0000000000..d9976dc2ac
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/IccConverterbase.cs
@@ -0,0 +1,49 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+#nullable disable
+
+using System.Numerics;
+using System.Runtime.CompilerServices;
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Conversion.Icc;
+
+///
+/// Color converter for ICC profiles
+///
+internal abstract partial class IccConverterBase
+{
+ ///
+ /// Initializes a new instance of the class.
+ ///
+ /// The ICC profile to use for the conversions
+ /// True if the conversion is to the profile connection space (PCS); False if the conversion is to the data space
+ protected IccConverterBase(IccProfile profile, bool toPcs)
+ {
+ Guard.NotNull(profile, nameof(profile));
+ this.Init(profile, toPcs, profile.Header.RenderingIntent);
+ }
+
+ ///
+ /// Converts colors with the initially provided ICC profile
+ ///
+ /// The value to convert
+ /// The converted value
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ public Vector4 Calculate(Vector4 value) => this.calculator.Calculate(value);
+
+ ///
+ /// Converts colors with the initially provided ICC profile
+ ///
+ /// The source colors
+ /// The destination colors
+ public void Calculate(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = this.Calculate(source[i]);
+ }
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/IccDataToDataConverter.cs b/src/ImageSharp/ColorProfiles/Icc/IccDataToDataConverter.cs
new file mode 100644
index 0000000000..cb4d89bb53
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/IccDataToDataConverter.cs
@@ -0,0 +1,22 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+using SixLabors.ImageSharp.ColorProfiles.Conversion.Icc;
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Icc;
+
+///
+/// Color converter for ICC profiles
+///
+internal class IccDataToDataConverter : IccConverterBase
+{
+ ///
+ /// Initializes a new instance of the class.
+ ///
+ /// The ICC profile to use for the conversions
+ public IccDataToDataConverter(IccProfile profile)
+ : base(profile, true) // toPCS is true because in this case the PCS space is also a data space
+ {
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/IccDataToPcsConverter.cs b/src/ImageSharp/ColorProfiles/Icc/IccDataToPcsConverter.cs
new file mode 100644
index 0000000000..6e95d3cb32
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/IccDataToPcsConverter.cs
@@ -0,0 +1,22 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+using SixLabors.ImageSharp.ColorProfiles.Conversion.Icc;
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Icc;
+
+///
+/// Color converter for ICC profiles
+///
+internal class IccDataToPcsConverter : IccConverterBase
+{
+ ///
+ /// Initializes a new instance of the class.
+ ///
+ /// The ICC profile to use for the conversions
+ public IccDataToPcsConverter(IccProfile profile)
+ : base(profile, true)
+ {
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/IccPcsToDataConverter.cs b/src/ImageSharp/ColorProfiles/Icc/IccPcsToDataConverter.cs
new file mode 100644
index 0000000000..d29517fca2
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/IccPcsToDataConverter.cs
@@ -0,0 +1,22 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+using SixLabors.ImageSharp.ColorProfiles.Conversion.Icc;
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Icc;
+
+///
+/// Color converter for ICC profiles
+///
+internal class IccPcsToDataConverter : IccConverterBase
+{
+ ///
+ /// Initializes a new instance of the class.
+ ///
+ /// The ICC profile to use for the conversions
+ public IccPcsToDataConverter(IccProfile profile)
+ : base(profile, false)
+ {
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/IccPcsToPcsConverter.cs b/src/ImageSharp/ColorProfiles/Icc/IccPcsToPcsConverter.cs
new file mode 100644
index 0000000000..30b44ca75c
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/IccPcsToPcsConverter.cs
@@ -0,0 +1,22 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+using SixLabors.ImageSharp.ColorProfiles.Conversion.Icc;
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Icc;
+
+///
+/// Color converter for ICC profiles
+///
+internal class IccPcsToPcsConverter : IccConverterBase
+{
+ ///
+ /// Initializes a new instance of the class.
+ ///
+ /// The ICC profile to use for the conversions
+ public IccPcsToPcsConverter(IccProfile profile)
+ : base(profile, true)
+ {
+ }
+}
diff --git a/src/ImageSharp/ColorProfiles/Icc/SrgbV4Profile.Generated.cs b/src/ImageSharp/ColorProfiles/Icc/SrgbV4Profile.Generated.cs
new file mode 100644
index 0000000000..a4d673488e
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/Icc/SrgbV4Profile.Generated.cs
@@ -0,0 +1,45 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+//
+
+using SixLabors.ImageSharp.Metadata.Profiles.Icc;
+
+namespace SixLabors.ImageSharp.ColorProfiles.Conversion.Icc;
+
+internal static class SrgbV4Profile
+{
+ // Generated using the sRGB-v4.icc profile found at https://github.com/saucecontrol/Compact-ICC-Profiles
+ private static ReadOnlySpan Data => new byte[]
+ {
+ 0, 0, 1, 224, 108, 99, 109, 115, 4, 32, 0, 0, 109, 110, 116, 114, 82, 71, 66, 32, 88, 89, 90, 32, 7, 226, 0, 3, 0,
+ 20, 0, 9, 0, 14, 0, 29, 97, 99, 115, 112, 77, 83, 70, 84, 0, 0, 0, 0, 115, 97, 119, 115, 99, 116, 114, 108, 0, 0,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 246, 214, 0, 1, 0, 0, 0, 0, 211, 45, 104, 97, 110, 100, 163, 178, 171,
+ 223, 92, 167, 3, 18, 168, 85, 164, 236, 53, 122, 209, 243, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 10, 100, 101, 115, 99, 0, 0, 0, 252, 0, 0, 0, 36, 99,
+ 112, 114, 116, 0, 0, 1, 32, 0, 0, 0, 34, 119, 116, 112, 116, 0, 0, 1, 68, 0, 0, 0, 20, 99, 104, 97, 100, 0, 0,
+ 1, 88, 0, 0, 0, 44, 114, 88, 89, 90, 0, 0, 1, 132, 0, 0, 0, 20, 103, 88, 89, 90, 0, 0, 1, 152, 0, 0, 0,
+ 20, 98, 88, 89, 90, 0, 0, 1, 172, 0, 0, 0, 20, 114, 84, 82, 67, 0, 0, 1, 192, 0, 0, 0, 32, 103, 84, 82, 67,
+ 0, 0, 1, 192, 0, 0, 0, 32, 98, 84, 82, 67, 0, 0, 1, 192, 0, 0, 0, 32, 109, 108, 117, 99, 0, 0, 0, 0, 0,
+ 0, 0, 1, 0, 0, 0, 12, 101, 110, 85, 83, 0, 0, 0, 8, 0, 0, 0, 28, 0, 115, 0, 82, 0, 71, 0, 66, 109, 108,
+ 117, 99, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 12, 101, 110, 85, 83, 0, 0, 0, 6, 0, 0, 0, 28, 0, 67, 0,
+ 67, 0, 48, 0, 33, 88, 89, 90, 32, 0, 0, 0, 0, 0, 0, 246, 214, 0, 1, 0, 0, 0, 0, 211, 45, 115, 102, 51, 50,
+ 0, 0, 0, 0, 0, 1, 12, 63, 0, 0, 5, 221, 255, 255, 243, 38, 0, 0, 7, 144, 0, 0, 253, 146, 255, 255, 251, 161, 255,
+ 255, 253, 162, 0, 0, 3, 220, 0, 0, 192, 113, 88, 89, 90, 32, 0, 0, 0, 0, 0, 0, 111, 160, 0, 0, 56, 242, 0, 0,
+ 3, 143, 88, 89, 90, 32, 0, 0, 0, 0, 0, 0, 98, 150, 0, 0, 183, 137, 0, 0, 24, 218, 88, 89, 90, 32, 0, 0, 0,
+ 0, 0, 0, 36, 160, 0, 0, 15, 133, 0, 0, 182, 196, 112, 97, 114, 97, 0, 0, 0, 0, 0, 3, 0, 0, 0, 2, 102, 105,
+ 0, 0, 242, 167, 0, 0, 13, 89, 0, 0, 19, 208, 0, 0, 10, 91,
+ };
+
+ private static readonly Lazy LazyIccProfile = new(() => GetIccProfile());
+
+ public static IccProfile GetProfile() => LazyIccProfile.Value;
+
+ private static IccProfile GetIccProfile()
+ {
+ byte[] buffer = new byte[Data.Length];
+ Data.CopyTo(buffer);
+ return new IccProfile(buffer);
+ }
+}
+
diff --git a/src/ImageSharp/ColorProfiles/KnownYCbCrMatrices.cs b/src/ImageSharp/ColorProfiles/KnownYCbCrMatrices.cs
new file mode 100644
index 0000000000..e2b7bf1026
--- /dev/null
+++ b/src/ImageSharp/ColorProfiles/KnownYCbCrMatrices.cs
@@ -0,0 +1,62 @@
+// Copyright (c) Six Labors.
+// Licensed under the Six Labors Split License.
+
+using System.Numerics;
+
+namespace SixLabors.ImageSharp.ColorProfiles;
+
+///
+/// Provides standard YCbCr matrices for RGB to YCbCr conversion.
+///
+public static class KnownYCbCrMatrices
+{
+#pragma warning disable SA1137 // Elements should have the same indentation
+#pragma warning disable SA1117 // Parameters should be on same line or separate lines
+ ///
+ /// ITU-R BT.601 (SD video standard).
+ ///
+ public static readonly YCbCrMatrix BT601 = new(
+ new Matrix4x4(
+ 0.299000F, 0.587000F, 0.114000F, 0F,
+ -0.168736F, -0.331264F, 0.500000F, 0F,
+ 0.500000F, -0.418688F, -0.081312F, 0F,
+ 0F, 0F, 0F, 1F),
+ new Matrix4x4(
+ 1.000000F, 0.000000F, 1.402000F, 0F,
+ 1.000000F, -0.344136F, -0.714136F, 0F,
+ 1.000000F, 1.772000F, 0.000000F, 0F,
+ 0F, 0F, 0F, 1F),
+ new Vector3(0F, 0.5F, 0.5F));
+
+ ///
+ /// ITU-R BT.709 (HD video, sRGB standard).
+ ///
+ public static readonly YCbCrMatrix BT709 = new(
+ new Matrix4x4(
+ 0.212600F, 0.715200F, 0.072200F, 0F,
+ -0.114572F, -0.385428F, 0.500000F, 0F,
+ 0.500000F, -0.454153F, -0.045847F, 0F,
+ 0F, 0F, 0F, 1F),
+ new Matrix4x4(
+ 1.000000F, 0.000000F, 1.574800F, 0F,
+ 1.000000F, -0.187324F, -0.468124F, 0F,
+ 1.000000F, 1.855600F, 0.000000F, 0F,
+ 0F, 0F, 0F, 1F),
+ new Vector3(0F, 0.5F, 0.5F));
+
+ ///
+ /// ITU-R BT.2020 (UHD/4K video standard).
+ ///
+ public static readonly YCbCrMatrix BT2020 = new(
+ new Matrix4x4(
+ 0.262700F, 0.678000F, 0.059300F, 0F,
+ -0.139630F, -0.360370F, 0.500000F, 0F,
+ 0.500000F, -0.459786F, -0.040214F, 0F,
+ 0F, 0F, 0F, 1F),
+ new Matrix4x4(
+ 1.000000F, 0.000000F, 1.474600F, 0F,
+ 1.000000F, -0.164553F, -0.571353F, 0F,
+ 1.000000F, 1.881400F, 0.000000F, 0F,
+ 0F, 0F, 0F, 1F),
+ new Vector3(0F, 0.5F, 0.5F));
+}
diff --git a/src/ImageSharp/ColorProfiles/Lms.cs b/src/ImageSharp/ColorProfiles/Lms.cs
index 5a6791b2d7..3aa3d72557 100644
--- a/src/ImageSharp/ColorProfiles/Lms.cs
+++ b/src/ImageSharp/ColorProfiles/Lms.cs
@@ -82,20 +82,53 @@ public Lms(Vector3 vector)
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static bool operator !=(Lms left, Lms right) => !left.Equals(right);
- ///
- /// Returns a new representing this instance.
- ///
- /// The .
- [MethodImpl(MethodImplOptions.AggressiveInlining)]
- public Vector3 ToVector3() => new(this.L, this.M, this.S);
+ ///
+ public Vector4 ToScaledVector4()
+ {
+ Vector3 v3 = default;
+ v3 += this.AsVector3Unsafe();
+ v3 += new Vector3(1F);
+ v3 /= 2F;
+ return new Vector4(v3, 1F);
+ }
///
- public static Lms FromProfileConnectingSpace(ColorConversionOptions options, in CieXyz source)
+ public static Lms FromScaledVector4(Vector4 source)
+ {
+ Vector3 v3 = source.AsVector3();
+ v3 *= 2F;
+ v3 -= new Vector3(1F);
+ return new Lms(v3);
+ }
+
+ ///
+ public static void ToScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = source[i].ToScaledVector4();
+ }
+ }
+
+ ///
+ public static void FromScaledVector4(ReadOnlySpan source, Span destination)
{
- Vector3 vector = Vector3.Transform(source.ToVector3(), options.AdaptationMatrix);
- return new Lms(vector);
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = FromScaledVector4(source[i]);
+ }
}
+ ///
+ public static Lms FromProfileConnectingSpace(ColorConversionOptions options, in CieXyz source)
+ => new(Vector3.Transform(source.AsVector3Unsafe(), options.AdaptationMatrix));
+
///
public static void FromProfileConnectionSpace(ColorConversionOptions options, ReadOnlySpan source, Span destination)
{
@@ -110,10 +143,7 @@ public static void FromProfileConnectionSpace(ColorConversionOptions options, Re
///
public CieXyz ToProfileConnectingSpace(ColorConversionOptions options)
- {
- Vector3 vector = Vector3.Transform(this.ToVector3(), options.InverseAdaptationMatrix);
- return new CieXyz(vector);
- }
+ => new(Vector3.Transform(this.AsVector3Unsafe(), options.InverseAdaptationMatrix));
///
public static void ToProfileConnectionSpace(ColorConversionOptions options, ReadOnlySpan source, Span destination)
diff --git a/src/ImageSharp/ColorProfiles/Rgb.cs b/src/ImageSharp/ColorProfiles/Rgb.cs
index 6698e12cb8..9b95278bc9 100644
--- a/src/ImageSharp/ColorProfiles/Rgb.cs
+++ b/src/ImageSharp/ColorProfiles/Rgb.cs
@@ -81,11 +81,54 @@ public Rgb(Vector3 source)
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static bool operator !=(Rgb left, Rgb right) => !left.Equals(right);
+ ///
+ /// Initializes the color instance from a generic scaled .
+ ///
+ /// The vector to load the color from.
+ /// The .
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ public static Rgb FromScaledVector4(Vector4 source)
+ => new(source.AsVector3());
+
+ ///
+ /// Expands the color into a generic ("scaled") representation
+ /// with values scaled and usually clamped between 0 and 1.
+ /// The vector components are typically expanded in least to greatest significance order.
+ ///
+ /// The .
+ [MethodImpl(MethodImplOptions.AggressiveInlining)]
+ public Vector4 ToScaledVector4()
+ => new(this.ToScaledVector3(), 1F);
+
+ ///
+ public static void ToScaledVector4(ReadOnlySpan source, Span destination)
+ {
+ Guard.DestinationShouldNotBeTooShort(source, destination, nameof(destination));
+
+ // TODO: Optimize via SIMD
+ for (int i = 0; i < source.Length; i++)
+ {
+ destination[i] = source[i].ToScaledVector4();
+ }
+ }
+
+ ///