Images don't show properly with imshow after certain calculation
4 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Walker
el 10 de Dic. de 2019
Comentada: Walker
el 16 de Dic. de 2019
Hello all,
The question is I need to calculate the matrix of a uint8 image. After the calculation (one order difference of the field) and normalization, I tried to figure out the image with imshow, however, the reults is a white blank or no longer smooth. The figure1.jpg can be any picture as you can find.
I should clarify that the values are uint8 after imread and I used "im2double" to transfer the type to double. The final imshow image is black and the final image will looks like some discrete dots when the matrix times 255.
clear all
deltax = 1;
deltay = 1;
AmpDistribution = imread('figure1.jpg')
AmpDistributionDouble = im2double(AmpDistribution);
AmpDistributionDouble(isnan(AmpDistributionDouble))=0;
[MSize, NSize, VSize] = size(AmpDistributionDouble);
DiffoverA = zeros(MSize,NSize,VSize);
for i = 2: MSize-1
for j = 2:NSize-1
for k = 1:VSize
DiffoverA(i,j,k) = (( AmpDistributionDouble(i+1,j,k) -AmpDistributionDouble(i,j,k))./(deltax) + (AmpDistributionDouble(i,j+1,k) -AmpDistributionDouble(i,j,k))./(deltay^2 ));
end
end
end
DiffoverA(isnan(DiffoverA))=0;
DiffoverA(DiffoverA==inf)=0;
MinValue = min(DiffoverA,[],'all');
MaxValue = max(DiffoverA,[],'all');
NormDelta = 255*(DiffoverA -abs(MinValue))/(MaxValue + abs(MinValue));
figure(1)
image(AmpDistribution);
figure(2)
imshow(AmpDistributionDouble,[]);
figure(3)
imshow(NormDelta,[]);
shading interp
0 comentarios
Respuesta aceptada
Bjorn Gustavsson
el 10 de Dic. de 2019
When you do things like this the errors are always (in my experience) in different normalizations and type-castings. To work these out I always go to the stage of looking at the individual R/G/B-planes of my images and check where they lie in intensity. Typically something like this:
subplot(2,2,1)
imagesc(NormDelta(:,:,1)),colorbar
subplot(2,2,2)
imagesc(NormDelta(:,:,2)),colorbar
subplot(2,2,3)
imagesc(NormDelta(:,:,3)),colorbar
subplot(2,2,4)
manualScale = @(I) (I-min(I(:)))/(max(I(:))-min(I(:)));
imagesc(manualScale(NormDelta))
Then you'll see where your normalised NormDelta are in intensity.
Finally: Your normalization step looks "very peculiar" to me, but I cant say if that is what you want.
HTH
14 comentarios
Image Analyst
el 13 de Dic. de 2019
Instead of
NormDelta = 255*(DiffoverA -abs(MinValue))/(MaxValue + abs(MinValue));
try using mat2gray(), rescale(), or imadjust(). And I would use imshow() with [] rather than image() or imagesc(). If the image still looks uniform, I'd examine the image in the workspace variable editor or look at its histogram. Maybe it really is largely the same value except for a very few pixels that are ruining the normalization.
Más respuestas (0)
Ver también
Productos
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!