1 | % 'ima_edge_detection': find edges
|
---|
2 |
|
---|
3 | %------------------------------------------------------------------------
|
---|
4 | %%%% Use the general syntax for transform fields with a single input and parameters %%%%
|
---|
5 | % OUTPUT:
|
---|
6 | % DataOut: output field structure
|
---|
7 | %
|
---|
8 | %INPUT:
|
---|
9 | % DataIn: input field structure
|
---|
10 | % Param: matlab structure whose field Param.TransformInput contains the filter parameters
|
---|
11 | %-----------------------------------
|
---|
12 |
|
---|
13 | %=======================================================================
|
---|
14 | % Copyright 2008-2021, LEGI UMR 5519 / CNRS UGA G-INP, Grenoble, France
|
---|
15 | % http://www.legi.grenoble-inp.fr
|
---|
16 | % Joel.Sommeria - Joel.Sommeria (A) legi.cnrs.fr
|
---|
17 | %
|
---|
18 | % This file is part of the toolbox UVMAT.
|
---|
19 | %
|
---|
20 | % UVMAT is free software; you can redistribute it and/or modify
|
---|
21 | % it under the terms of the GNU General Public License as published
|
---|
22 | % by the Free Software Foundation; either version 2 of the license,
|
---|
23 | % or (at your option) any later version.
|
---|
24 | %
|
---|
25 | % UVMAT is distributed in the hope that it will be useful,
|
---|
26 | % but WITHOUT ANY WARRANTY; without even the implied warranty of
|
---|
27 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
---|
28 | % GNU General Public License (see LICENSE.txt) for more details.
|
---|
29 | %=======================================================================
|
---|
30 |
|
---|
31 | function DataOut=ima_patch_detection(DataIn,Param)
|
---|
32 |
|
---|
33 | %% request input parameters
|
---|
34 | if isfield(DataIn,'Action') && isfield(DataIn.Action,'RUN') && isequal(DataIn.Action.RUN,0)
|
---|
35 | prompt = {'npx';'npy';'threshold'};
|
---|
36 | dlg_title = 'get the filter size in x and y';
|
---|
37 | num_lines= 3;
|
---|
38 | def = { '50';'50';'0.3'};
|
---|
39 | if isfield(Param,'TransformInput')&&isfield(Param.TransformInput,'FilterBoxSize_x')&&...
|
---|
40 | isfield(Param.TransformInput,'FilterBoxSize_y')&&isfield(Param.TransformInput,'LumThreshold')
|
---|
41 | def={num2str(Param.TransformInput.FilterBoxSize_x);num2str(Param.TransformInput.FilterBoxSize_y);num2str(Param.TransformInput.LumThreshold)};
|
---|
42 | end
|
---|
43 |
|
---|
44 |
|
---|
45 |
|
---|
46 | %% create the GUI set_threshold
|
---|
47 | set(0,'Units','points')
|
---|
48 | ScreenSize=get(0,'ScreenSize');% get the size of the screen, to put the fig on the upper right
|
---|
49 | Width=350;% fig width in points (1/72 inch)
|
---|
50 | Height=min(0.8*ScreenSize(4),300);
|
---|
51 | Left=ScreenSize(3)- Width-40; %right edge close to the right, with margin=40
|
---|
52 | Bottom=ScreenSize(4)-Height-40; %put fig at top right
|
---|
53 | hfig=findobj(allchild(0),'Tag','set_threshold');
|
---|
54 | if isempty(hfig)
|
---|
55 | hfig=figure('name','set_threshold','tag','set_threshold','MenuBar','none','NumberTitle','off','Units',...
|
---|
56 | 'pixels','Position',[Left,Bottom,Width,Height],'DeleteFcn',@closefcn)%);
|
---|
57 | else
|
---|
58 | figure(hfig)
|
---|
59 | end
|
---|
60 | BackgroundColor=get(hfig,'Color');
|
---|
61 | hh=0.14; % box height (relative)
|
---|
62 | ii=0.01; % gap between uicontrols
|
---|
63 |
|
---|
64 | ww=(1-5*ii)/4; % box width (relative)
|
---|
65 | % first raw of the GUI
|
---|
66 | %Amask=imread('/fsnet/project/coriolis/2015/15MINI_MEDDY/0_REF_FILES/mask_patch.png');
|
---|
67 | uicontrol('Style','text','Units','normalized', 'Position', [ii 0.95-2*ii-0.75*hh ww hh/2],'BackgroundColor',BackgroundColor,...
|
---|
68 | 'String','Threshold','FontUnits','points','FontSize',12,'FontWeight','bold','ForegroundColor','blue','HorizontalAlignment','right');%title
|
---|
69 | % uicontrol('Style','edit','Units','normalized', 'Position', [2*ii+ww 0.95-2*ii-hh ww hh],'tag','num_Z_1','BackgroundColor',[1 1 1],...
|
---|
70 | % 'String',num2str(SliceCoord(1,3)),'FontUnits','points','FontSize',12,'FontWeight','bold','TooltipString','''num_Z_1'': z position of first slice');%edit box
|
---|
71 | hThreshold=uicontrol('Style','edit','Units','normalized', 'Position', [3*ii+2*ww 0.95-2*ii-hh ww hh],'tag','Threshold','BackgroundColor',[1 1 1],...
|
---|
72 | 'FontUnits','points','FontSize',12,'FontWeight','bold');%edit box
|
---|
73 | uicontrol('Style','slider','Units','normalized', 'Position', [0.1 0.1 0.8 0.1],'tag','SetThreshold','BackgroundColor',[1 1 1],...
|
---|
74 | 'FontUnits','points','FontSize',12,'FontWeight','bold','Callback',@(hObject,eventdata)SetThreshold_Callback(hObject,eventdata));%edit box
|
---|
75 | % uiwait(hfig);
|
---|
76 | %
|
---|
77 | %
|
---|
78 | % DataOut.TransformInput.FilterBoxSize_x=str2num(answer{1}); %size of the filtering window
|
---|
79 | % DataOut.TransformInput.FilterBoxSize_y=str2num(answer{2}); %size of the filtering window
|
---|
80 | DataOut.TransformInput.LumThreshold=str2num(get(hThreshold,'String')) %size of the filtering window
|
---|
81 |
|
---|
82 |
|
---|
83 |
|
---|
84 | return
|
---|
85 | end
|
---|
86 |
|
---|
87 | DataOut=DataIn; %default
|
---|
88 | %Amask=imread('/fsnet/project/coriolis/2015/15MINI_MEDDY/0_REF_FILES/mask_patch.png');
|
---|
89 | hfig=findobj(allchild(0),'Tag','set_threshold');
|
---|
90 | hThreshold=findobj(hfig,'Tag','Threshold');
|
---|
91 | Threshold=str2num(get(hThreshold,'String'));
|
---|
92 | plot_mask(DataIn.A,Threshold)
|
---|
93 | % Athreshold=Param.TransformInput.LumThreshold;
|
---|
94 | % %
|
---|
95 | % % DataOut.A=zeros(size(DataIn.A,1),size(DataIn.A,2),3);
|
---|
96 | % DataOut.A=(DataIn.A>Athreshold);%transform to the initial image format
|
---|
97 | % % DataOut.A(:,:,1)=DataIn.A;%transform to the initial image format, red
|
---|
98 | % STATS = regionprops(DataOut.A, 'FilledArea','MinorAxisLength','MajorAxisLength','PixelIdxList');
|
---|
99 | % Area=zeros(size(STATS));
|
---|
100 | % for iobj=1:numel(STATS)
|
---|
101 | % Area(iobj)=STATS(iobj).FilledArea;
|
---|
102 | % end
|
---|
103 | % [Area, main_obj]=max(Area);
|
---|
104 | % MajorAxisLength=STATS(main_obj).MajorAxisLength
|
---|
105 | % MinorAxisLength=STATS(main_obj).MinorAxisLength
|
---|
106 | % for iobj=1:numel(STATS)
|
---|
107 | % if iobj~=main_obj
|
---|
108 | % DataOut.A(STATS(iobj).PixelIdxList)=0;
|
---|
109 | % end
|
---|
110 | % end
|
---|
111 | %
|
---|
112 | % DataOut.A=Amax*DataOut.A;
|
---|
113 |
|
---|
114 |
|
---|
115 | %------------------------------------------------------------------------
|
---|
116 | % function called by selecting CheckRefraction in the GUI set_slices
|
---|
117 | function SetThreshold_Callback(hObject,eventdata)
|
---|
118 | %------------------------------------------------------------------------
|
---|
119 | hfig=get(hObject,'parent');
|
---|
120 | hThreshold=findobj(hfig,'Tag','Threshold');
|
---|
121 | Threshold=get(hObject,'Value');
|
---|
122 |
|
---|
123 | huvmat=findobj(allchild(0),'Tag','uvmat');
|
---|
124 | hhuvmat=guidata(huvmat);
|
---|
125 | UvData=get(huvmat,'UserData');
|
---|
126 | Amax=max(max(UvData.Field.A));
|
---|
127 | Amin=min(min(UvData.Field.A));
|
---|
128 | Threshold=Amin+(Amax-Amin)*Threshold;
|
---|
129 | set(hThreshold,'String',num2str(Threshold))
|
---|
130 |
|
---|
131 | plot_mask(UvData.Field.A,Threshold)
|
---|
132 |
|
---|
133 | UvData.XmlData{1}.TransformInput.LumThreshold=Threshold;
|
---|
134 | set(huvmat,'UserData',UvData)
|
---|
135 |
|
---|
136 |
|
---|
137 |
|
---|
138 |
|
---|
139 |
|
---|
140 |
|
---|
141 | % h_refraction(1)=findobj(hset_slice,'String','surface');
|
---|
142 | % h_refraction(2)=findobj(hset_slice,'Tag','num_H');
|
---|
143 | % h_refraction(3)=findobj(hset_slice,'String','index');
|
---|
144 | % h_refraction(4)=findobj(hset_slice,'Tag','num_RefractionIndex');
|
---|
145 | % if isequal(get(hObject,'Value'),1)
|
---|
146 | % set(h_refraction,'Visible','on')
|
---|
147 | % else
|
---|
148 | % set(h_refraction,'Visible','off')
|
---|
149 | % end
|
---|
150 |
|
---|
151 | function plot_mask(A,Threshold)
|
---|
152 |
|
---|
153 | A=(A>Threshold) ;%& Amask>100;%transform to binary image format
|
---|
154 | try
|
---|
155 | STATS = regionprops(A, 'FilledArea','MinorAxisLength','MajorAxisLength','PixelIdxList');
|
---|
156 | Area=zeros(size(STATS));
|
---|
157 | for iobj=1:numel(STATS)
|
---|
158 | Area(iobj)=STATS(iobj).FilledArea;
|
---|
159 | end
|
---|
160 | [Area, main_obj]=max(Area);
|
---|
161 | MajorAxisLength=STATS(main_obj).MajorAxisLength;
|
---|
162 | MinorAxisLength=STATS(main_obj).MinorAxisLength;
|
---|
163 | for iobj=1:numel(STATS)
|
---|
164 | if iobj~=main_obj
|
---|
165 | A(STATS(iobj).PixelIdxList)=0;
|
---|
166 | end
|
---|
167 | end
|
---|
168 | catch ME
|
---|
169 | disp('image toolbox not available, skipped')
|
---|
170 | end
|
---|
171 | [npy,npx]=size(A);
|
---|
172 | A=cat(3,A,zeros(npy,npx,2));% make a color image
|
---|
173 |
|
---|
174 | huvmat=findobj(allchild(0),'Tag','uvmat');
|
---|
175 | hhuvmat=guidata(huvmat);
|
---|
176 | hmask=findobj(hhuvmat.PlotAxes,'Tag','MaskFig')
|
---|
177 | if isempty(hmask)
|
---|
178 | axes(hhuvmat.PlotAxes)
|
---|
179 | hold on
|
---|
180 | imagesc([0.5 npx],[npy-0.5 0.5],A,'AlphaData',0.2,'Tag','MaskFig')
|
---|
181 | else
|
---|
182 | set(hmask,'CData',A)
|
---|
183 | set(hmask,'AlphaData',0.2)
|
---|
184 | end
|
---|
185 |
|
---|
186 |
|
---|
187 | function closefcn(hObject,eventdata)
|
---|
188 |
|
---|
189 | huvmat=findobj(allchild(0),'Tag','uvmat');
|
---|
190 | hhuvmat=guidata(huvmat);
|
---|
191 | hmask=findobj(hhuvmat.PlotAxes,'Tag','MaskFig')
|
---|
192 | if ~isempty(hmask)
|
---|
193 | delete(hmask)
|
---|
194 | end |
---|