David Haberthür
369 lines
14 KiB
Matlab
369 lines
14 KiB
Matlab
function displayDumpImage(sourceImage, varargin )
|
|
|
|
newFigure = 0;
|
|
if nargin > 1
|
|
newFigure = varargin{1};
|
|
end
|
|
display = 0;
|
|
if nargin > 2
|
|
display = varargin{2} ;
|
|
end
|
|
|
|
% read the image data
|
|
if isnumeric(sourceImage) | islogical(sourceImage)
|
|
img = sourceImage;
|
|
[height, width] = size(sourceImage);
|
|
sourceImage = 'unknown';
|
|
else
|
|
[img, width, height] = readDumpImage(sourceImage,inf,'float32');
|
|
end
|
|
if nargin > 4
|
|
if varargin{4} < 0
|
|
nimg(:,:) = img(:,width:-1:1);
|
|
img=nimg;
|
|
clear nimg;
|
|
end
|
|
end
|
|
if newFigure
|
|
% display the image in a new figure
|
|
fighnd = figure;
|
|
else
|
|
fighnd = gcf;
|
|
end
|
|
colormap('gray');
|
|
datamin = min(min(img));
|
|
datamax = max(max(img));
|
|
autoscale = 1;
|
|
if ( nargin > 3 ) ...
|
|
&& ( length(varargin{3}) == 2 )
|
|
valuerange = varargin{3};
|
|
minval = valuerange(1);
|
|
maxval = valuerange(2);
|
|
autoscale = 0
|
|
else
|
|
maxval = datamax;
|
|
minval = datamin;
|
|
end
|
|
factor=1/( datamax-datamin );
|
|
shift=datamin;
|
|
if ~display
|
|
if autoscale == 1
|
|
|
|
edthnd = [ findobj(fighnd,'type','uicontrol','style','edit','tag','MIN_EDT') ...
|
|
findobj(fighnd,'type','uicontrol','style','edit','tag','MAX_EDT') ];
|
|
if length( edthnd == 2 )
|
|
value = sscanf(get(edthnd(1),'string'),'%f');;
|
|
minval = value;
|
|
value = sscanf(get(edthnd(2),'string'),'%f');;
|
|
maxval = value;;
|
|
end
|
|
end
|
|
imghnd=imagesc(img,[minval maxval]);
|
|
set(imghnd,'tag','IMAGE_DSP');
|
|
axis equal;
|
|
titlehnd = title(sourceImage);
|
|
set(titlehnd,'tag','TITLE_TXT','HandleVisibility','on','STRING',sourceImage);
|
|
cbarhnd=colorbar('east');
|
|
set(cbarhnd,'tag','COLOR_DSP');
|
|
imgaxhnd=get(imghnd,'parent');
|
|
fighnd=get(imgaxhnd,'parent');
|
|
sldpos=get(cbarhnd,'position');
|
|
sldpos(1) = sldpos(1) + 2 * sldpos(3);
|
|
sldpos(2) = sldpos(2) - 0.015;
|
|
sldpos(3) = 0.025;
|
|
sldpos(4) = sldpos(4) + 0.03;
|
|
sldhnd = [ findobj(fighnd,'type','uicontrol','style','slider','tag','MIN_SLD') ...
|
|
findobj(fighnd,'type','uicontrol','style','slider','tag','MAX_SLD') ...
|
|
findobj(fighnd,'type','uicontrol','style','slider','tag','IMG_SLD') ];
|
|
if length(sldhnd) < 2
|
|
if length(sldhnd) < 1
|
|
sldhnd(1) = uicontrol(fighnd);
|
|
end
|
|
sldhnd(2) = uicontrol(fighnd);
|
|
end
|
|
set(sldhnd(1),'style','slider','units','normalized','position',sldpos);
|
|
set(sldhnd(1),'min',0,'max',1,'sliderstep',[0.0005 0.05],'tag','MIN_SLD');
|
|
set(sldhnd(1),'value',( minval -shift ) * factor,'callback',@callback,'max',1,'min',0);
|
|
sldpos(1) = sldpos(1) + 0.03;
|
|
set(sldhnd(2),'style','slider','units','normalized','position',sldpos);
|
|
set(sldhnd(2),'min',0,'max',1,'sliderstep',[0.0005 0.05],'tag','MAX_SLD');
|
|
set(sldhnd(2),'value',( maxval - shift ) * factor ,'callback',@callback,'min',0,'max',1);
|
|
edtpos = sldpos;
|
|
edtpos(1) = edtpos(1) - 0.075;
|
|
|
|
edtpos(2) = edtpos(2) - 0.075;
|
|
edtpos(3) = 0.1;
|
|
edtpos(4) = 0.05;
|
|
edthnd = [ findobj(fighnd,'type','uicontrol','style','edit','tag','MIN_EDT') ...
|
|
findobj(fighnd,'type','uicontrol','style','edit','tag','MAX_EDT') ];
|
|
%set(gca,'color',[0,0,0],'clim',[minval maxval]);
|
|
set(gca,'color',[0,0,0]);
|
|
if length(edthnd) < 2
|
|
if length(edthnd) < 1
|
|
edthnd(1) = uicontrol(fighnd);
|
|
end
|
|
edthnd(2) = uicontrol(fighnd);
|
|
end
|
|
set(edthnd(1),'style','edit','units','normalized','position',edtpos,'min',min(datamin,minval),'max',max(datamax,maxval));
|
|
set(edthnd(1),'tag','MIN_EDT','string',sprintf('%.4f',minval),'callback',@callback);
|
|
edtpos(2) = sldpos(2) + sldpos(4) + 0.035;
|
|
set(edthnd(2),'style','edit','units','normalized','position',edtpos,'min',datamin,'max',datamax);
|
|
set(edthnd(2),'tag','MAX_EDT','string',sprintf('%.4f',maxval),'callback',@callback);
|
|
if length(sldhnd) < 3
|
|
sldhnd(3) = uicontrol(fighnd);
|
|
end
|
|
pos = get(imgaxhnd,'position');
|
|
set(sldhnd(3),'style','slider','units','normalized','tag','IMG_SLD');
|
|
position = get(sldhnd(3),'position');
|
|
position(1) = pos(1) -0.1;
|
|
position(2) = pos(2);;
|
|
position(3) = 0.03;
|
|
position(4) = pos(4);
|
|
[m dum1 dum2 dum3 sample] = regexp({get(titlehnd,'string') pwd},'(\/|^)([^\s\/]+)\/(rec(_DMP)?|sin|prj2cpr|$;)');
|
|
recdir = '';
|
|
if isempty(m{1})
|
|
sample = sample{2}{1};
|
|
else
|
|
sample = sample{1}{1};
|
|
end
|
|
if isempty(sample{3})
|
|
recdir = ['/' dir([pwd '/rec*']) ]
|
|
else
|
|
recdir = sample{3}
|
|
end;
|
|
count=length(dir([pwd recdir '/' sample{2} '*.DMP']));
|
|
set(sldhnd(3),'position',position,'min',0,'max',1,'sliderstep',[ 1 / ( count - 1 ), 10 / ( count - 1 ) ]);
|
|
[m dum1 dum2 dum3 digits] = regexp(get(titlehnd,'string'),[sample{2} '(\d+).*\.DMP']);
|
|
value = sscanf(digits{1}{1},'%f');
|
|
value = 1 - ( value - 1 ) / ( count - 1 );
|
|
set(sldhnd(3),'value',value,'callback',@callback);
|
|
set(fighnd,'userdata',{length(digits{1}{1}) count pwd} );
|
|
else
|
|
range = maxval - minval;
|
|
title(sourceImage);
|
|
cutoffon = range * 0.1761;
|
|
noiselevel = minval + cutoffon;
|
|
level = ones(2 ,2 ) .* range .* 0.5;
|
|
toplevel = maxval - cutoffon;
|
|
steps = range / 500;
|
|
if display < 2
|
|
mesh(img);
|
|
set(gca,'zlim',[minval maxval]);
|
|
return;
|
|
end
|
|
[ rows cols ] = size(img);
|
|
levelfaces=[];
|
|
subplot(2,2,4);
|
|
set(gca,'units','normalized','visible','off')
|
|
text(0,0.75,0.85,'File:','horizontalalignment','left','verticalalignment','baseline','fontsize',12);
|
|
text(0,0.675,0.85,sourceImage,'horizontalalignment','left','verticalalignment','baseline','fontsize',12);
|
|
text(0,0.5,0.85,'Press j or k to move level up and down','horizontalalignment','left','verticalalignment','baseline','fontsize',12);
|
|
text(0,0.4,0.85,'Press n or t to mark level as nois or signal','horizontalalignment','left','verticalalignment','baseline','fontsize',12);
|
|
text(0,0.3,0.85,'Upper case makes 20 steps, lower case 1','horizontalalignment','left','verticalalignment','baseline','fontsize',12);
|
|
text(0,0.15,0.65,'Press q to quit','horizontalalignment','left','verticalalignment','baseline','fontsize',12);
|
|
set(gca,'position',[ 0.7 0.05 0.25 0.25 ],'view',[ 0 0 ],'cameraposition',[ 10000000 0 0 ],'cameraupvector',[0 1 0]);
|
|
leveltexts = text(0,0.05,0.85,[ 'Signal: ' num2str(toplevel,'%7.4f') ...
|
|
' Noise: ' num2str(noiselevel,'%7.4f') ...
|
|
' SNR: ' num2str(toplevel/noiselevel,'%5.2f') ], ...
|
|
'horizontalalignment','left','verticalalignment','baseline','fontsize',11);
|
|
subplot(2,2,2);
|
|
mesh(img);
|
|
set(gca,'position',[ 0.7 0.3 0.25 0.65 ],'view',[ 0 0 ],'cameraposition',[ 10000000 0 0 ],'cameraupvector',[0 1 0]);
|
|
levelfaces(1,2) = surface( [ 1 cols ; 1 cols ] , [ 1 1 ; rows rows ],level );
|
|
set(gca,'zlim',[ floor(minval * 500) ceil(maxval *500 ) ] ./ 500, 'ylim',[ 0 rows ],'drawmode','fast','color','none');
|
|
set(gca,'box','on');
|
|
rotate3d(gca,'off')
|
|
subplot(2,2,3)
|
|
mesh(img);
|
|
set(gca,'position',[ 0.05 0.05 0.65 0.25 ],'view',[ 0 0 ],'cameraposition',[ 0 -1000000 0 ],'cameraupvector',[0 0 1]);
|
|
set(gca,'zlim',[ floor(minval * 500) ceil(maxval *500 ) ] ./ 500, 'xlim',[ 0 cols ],'drawmode','fast','color','none');
|
|
set(gca,'box','on');
|
|
levelfaces(1,3) = surface( [ 1 cols ; 1 cols ] , [ 1 1 ; rows rows ],level );
|
|
rotate3d(gca,'off')
|
|
subplot(2,2,1)
|
|
mesh(img);
|
|
set(gcf,'units','normalized');
|
|
set(gca,'position',[ 0.05 0.30 0.65 0.65 ],'view',[0 0],'cameraposition',[ 0 0 10 ],'cameraupvector',[0 1 0])
|
|
set(gca,'xlim',[0 cols],'ylim',[0 rows],'visible','off','drawmode','fast');
|
|
set(gca,'color',[0,0,0],'box','on');
|
|
levelfaces(1) = surface( [ 1 cols ; 1 cols ] , [ 1 1 ; rows rows ],level );
|
|
set(levelfaces(1,1),'facecolor',[0.05 0.05 0.05]);
|
|
rotate3d(gca,'off');
|
|
disp 'Press j or k to move level plane up and down';
|
|
disp 'Press n or t to define actual level as noise- or toplevel';
|
|
disp 'Press q to quit';
|
|
disp 'Uppercdase has speed 10 instead of 1';
|
|
while 1
|
|
if ~waitforbuttonpress
|
|
continue;
|
|
end
|
|
key = get(gcf,'currentcharacter');
|
|
switch key
|
|
case {'q' 'Q'}
|
|
break;
|
|
case 'j'
|
|
level = increaseLevel(levelfaces,steps,minval,maxval);
|
|
case 'J'
|
|
level = increaseLevel(levelfaces,steps * 20,minval,maxval);
|
|
case 'k'
|
|
level = increaseLevel(levelfaces,-steps,minval,maxval);
|
|
case 'K'
|
|
level = increaseLevel(levelfaces,-steps * 20,minval,maxval);
|
|
case {'n' 'N' }
|
|
noiselevel = level(1,1);
|
|
set(leveltexts,'string',[ 'Signal: ' num2str(toplevel,'%7.4f') ...
|
|
' Noise: ' num2str(noiselevel,'%7.4f') ...
|
|
' SNR: ' num2str(toplevel/noiselevel,'%5.2f') ]);
|
|
case { 't' 'T' }
|
|
toplevel = level(1,1);
|
|
set(leveltexts,'string',[ 'Signal: ' num2str(toplevel,'%7.4f') ...
|
|
' Noise: ' num2str(noiselevel,'%7.4f') ...
|
|
' SNR: ' num2str(toplevel/noiselevel,'%5.2f') ]);
|
|
otherwise
|
|
end
|
|
end
|
|
end
|
|
clear gcf
|
|
|
|
|
|
function [level] = increaseLevel(hnd,by,minval,maxval)
|
|
levelvalues = get(hnd(1,1),'zdata');
|
|
newlevelvalues = levelvalues + by;
|
|
if ~isempty( find( newlevelvalues < minval ) )...
|
|
|| ~isempty( find( newlevelvalues > maxval ) )
|
|
level = levelvalues;
|
|
return;
|
|
end
|
|
for i=1:length(hnd)
|
|
set(hnd(i),'zdata',newlevelvalues);
|
|
end
|
|
level = newlevelvalues;
|
|
return
|
|
% [r c] = size(img);
|
|
% img(r,:)
|
|
% min(img(r,:))
|
|
% max(img(r,:))
|
|
% a = ind2sub(c,find([img(r,:)] < 0.0 ));
|
|
% img(r,a) = img(r,a) - min(img(r,:));
|
|
% figure
|
|
% plot(img(r,:))
|
|
|
|
%img(r,a) = img(r,a) - min(img(r,:))
|
|
|
|
|
|
|
|
function callback(hnd,data)
|
|
fig = get(hnd,'parent');
|
|
img = findobj(fig,'tag','IMAGE_DSP');
|
|
scb = findobj(fig,'tag','COLOR_DSP');
|
|
axs = get(img,'parent');
|
|
sbd = findobj(scb,'type','image');
|
|
switch get(hnd,'tag');
|
|
case 'MIN_SLD'
|
|
mxsld = findobj(fig,'tag','MAX_SLD');
|
|
mxedt = findobj(fig,'tag','MAX_EDT');
|
|
mnedt = findobj(fig,'tag','MIN_EDT');
|
|
value = get(hnd,'value');
|
|
datmin = get(mxedt,'min');
|
|
datmax = get(mxedt,'max');
|
|
factor = 1/(datmax - datmin);
|
|
shift = datmin;
|
|
value = value / factor + shift;
|
|
limit = get(mxsld,'value') / factor + shift;
|
|
if ( value > limit )
|
|
value = limit - 0.0001;
|
|
set(hnd,'value', ( value - shift ) * factor);
|
|
end
|
|
set(mnedt,'string',sprintf('%.4f',value));
|
|
limits = get(axs,'clim');
|
|
limits(1) = value;
|
|
set(axs,'clim',limits);
|
|
set(sbd,'ydata',limits);
|
|
set(scb,'ylim',limits);
|
|
case 'MAX_SLD'
|
|
mnsld = findobj(fig,'tag','MIN_SLD') ;
|
|
mxedt = findobj(fig,'tag','MAX_EDT') ;
|
|
mnedt = findobj(fig,'tag','MIN_EDT') ;
|
|
value = get(hnd,'value');
|
|
datmin = get(mxedt,'min');
|
|
datmax = get(mxedt,'max');
|
|
factor = 1/(datmax - datmin);
|
|
shift = datmin;
|
|
value = value / factor + shift;
|
|
limit = get(mnsld,'value') / factor + shift;
|
|
if ( value < limit )
|
|
value = limit + 0.0001;
|
|
set(hnd,'value', ( value - shift ) * factor);
|
|
end
|
|
set(mxedt,'string',sprintf('%.4f',value));
|
|
limits = get(axs,'clim')
|
|
limits(2) = value;
|
|
set(axs,'clim',limits);
|
|
set(sbd,'ydata',limits);
|
|
set(scb,'ylim',limits);
|
|
case 'MIN_EDT'
|
|
mnsld = findobj(fig,'tag','MIN_SLD') ;
|
|
mxsld = findobj(fig,'tag','MAX_SLD');
|
|
mxedt = findobj(fig,'tag','MAX_EDT') ;
|
|
mnedt = findobj(fig,'tag','MIN_EDT') ;
|
|
value = sscanf(get(hnd,'string'),'%f');
|
|
datmin = get(mxedt,'min');
|
|
datmax = get(mxedt,'max');
|
|
factor = 1/(datmax - datmin);
|
|
shift = datmin;
|
|
limit = get(mxsld,'value') / factor + shift;
|
|
if ( value > limit )
|
|
value = limit - 0.0001;
|
|
set(hnd,'string', sprintf('%.4f',value ));
|
|
elseif value < datmin
|
|
value = datmin ;
|
|
set(hnd,'string', sprintf('%.4f',value ));
|
|
end
|
|
limits = get(axs,'clim');
|
|
limits(1) = value;
|
|
set(axs,'clim',limits);
|
|
value = (value - shift ) * factor;
|
|
set(mnsld,'value',value);
|
|
set(sbd,'ydata',limits);
|
|
set(scb,'ylim',limits);
|
|
case 'MAX_EDT'
|
|
mnsld = findobj(fig,'tag','MIN_SLD') ;
|
|
mxsld = findobj(fig,'tag','MAX_SLD');
|
|
mxedt = findobj(fig,'tag','MAX_EDT') ;
|
|
mnedt = findobj(fig,'tag','MIN_EDT') ;
|
|
value = sscanf(get(hnd,'string'),'%f');
|
|
datmin = get(mxedt,'min');
|
|
datmax = get(mxedt,'max');
|
|
factor = 1/(datmax - datmin);
|
|
shift = datmin;
|
|
limit = get(mnsld,'value') / factor + shift;
|
|
if ( value < limit )
|
|
value = limit + 0.0001;
|
|
set(hnd,'string', sprintf('%.4f',value));
|
|
elseif value > datmax;
|
|
value = datmax;
|
|
set(hnd,'string', sprintf('%.4f',value));
|
|
end
|
|
limits = get(axs,'clim');
|
|
limits(2) = value;
|
|
set(axs,'clim',limits);
|
|
value = (value - shift ) * factor;
|
|
set(mxsld,'value',value);
|
|
set(sbd,'ydata',limits);
|
|
set(scb,'ylim',limits);
|
|
case 'IMG_SLD'
|
|
digitcount = get(fig,'userdata');
|
|
tithnd = findobj(fig,'tag','TITLE_TXT');
|
|
text = get(tithnd,'string');
|
|
value = round( ( 1 - get(hnd,'value') ) * ( digitcount{2} - 1 ) ) + 1;
|
|
pat = sprintf('\\d{%i,%i}(\\D+\\.DMP)$',digitcount{1},digitcount{1});
|
|
replace = sprintf([ sprintf('%%0.%ii',digitcount{1}) '$1' ],value);
|
|
thisdir = pwd;
|
|
cd(digitcount{3});
|
|
actfig = gcf;
|
|
figure(fig);
|
|
displayDumpImage([regexprep(text,pat,replace) ]);
|
|
figure(actfig);
|
|
cd(thisdir);
|
|
end
|