img = ColorConvert[ExampleData[{"TestImage", "Lena"}], "Grayscale"];
fft = Fourier[ImageData[img]];
fft = RotateLeft[fft, Floor[Dimensions[fft]/2]];
fftAbsData = Log[Abs[fft]] - Min[Log[Abs[fft]]];Abs[fft];
{fftAbsMin, fftAbsMax} = MinMax[fftAbsData];
myColorTable =
Flatten@{Table[{Blend[{Blue, Green, Yellow, Orange, Red}, x]}, {x, 1/256, 1, 1/256}]};
g = Colorize[Image[(fftAbsData - fftAbsMin)/(fftAbsMax - fftAbsMin)]Colorize[Image[fftAbsData],
ColorFunction -> (Blend[myColorTable, #] &)];
myLegend =
BarLegend[
{myColorTable, {Exp[fftAbsMin + Min[Log[Abs[fft]]]]fftAbsMin,
Exp[fftAbsMax + Min[Log[Abs[fft]]]]fftAbsMax}},
LegendLabel -> "amplitude", LegendMarkerSize -> {40, 300},
LabelStyle -> {FontFamily -> "Calibri", FontSize -> 15}
];
Legended[
Graphics[
Inset[g, Scaled[{.5, .5}], Automatic, Scaled[1]], Frame -> True,
FrameLabel -> {{"cycles/pixel", ""}, {"cycles/pixel", ""}},
PlotRange -> {{-0.5, 0.5}, {-0.5, 0.5}},
AspectRatio -> ImageAspectRatio@g
]
, myLegend
]
With a modified color table:
myColorTable =
Flatten@{Table[{Blend[{Blue, Green, Yellow, Orange}, x]}, {x, 1/50,
1, 1/50}],
Table[{Blend[{Orange, Red, White}, x]}, {x, 1/206, 1, 1/206}]};
one gets:
The same way I getthe FFT for the phases is:
fftPhasesData = Arg@fft/Pi;
{fftPhasesMin, fftPhasesMax} = MinMax[fftPhasesData];
myColorTable =
Flatten@{Table[{Blend[{Blue, Green, Yellow, Orange, Red}, x]}, {x, 1/256, 1, 1/256}]};
g = Colorize[Image[(fftPhasesData - fftPhasesMin)/(fftPhasesMax -
fftPhasesMin)], ColorFunction -> (Blend[myColorTable, #] &)];
myLegend =
BarLegend[
{myColorTable, {fftPhasesMin, fftPhasesMax}},
LegendLabel -> "phase (\[Pi])", LegendMarkerSize -> {40, 300},
LabelStyle -> {FontFamily -> "Calibri", FontSize -> 15}
];
Legended[
Graphics[
Inset[g, Scaled[{.5, .5}], Automatic, Scaled[1]], Frame -> True,
FrameLabel -> {{"cycles/pixel", ""}, {"cycles/pixel", ""}},
PlotRange -> {{-0.5, 0.5}, {-0.5, 0.5}},
AspectRatio -> ImageAspectRatio@g
]
, myLegend
]