Images don't show properly with imshow after certain calculation
2 次查看(过去 30 天)
显示 更早的评论
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 个评论
采纳的回答
Bjorn Gustavsson
2019-12-10
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 个评论
Image Analyst
2019-12-13
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.
更多回答(0 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!