@@ -325,12 +325,13 @@ function selectpoint_Callback(hObject, eventdata, handles)
325325segOutline(: ,end )=1 ;
326326segOutline = ~segOutline ;
327327
328+ global numofpixels ;
328329global newsegOutline ;
329330numofpixels= str2num(cell2mat(numberOfPoints(3 )));
330- newsegOutline = Expand(segOutline ,numofpixels );
331+ % newsegOutline = Expand(segOutline,numofpixels);
331332axes(handles .axes );
332333hold on
333- outcome1= evalue + newsegOutline * max(max(evalue ));
334+ outcome1= evalue + segOutline * max(max(evalue ));
334335imshow(outcome1 ,[]); % outline in original picture
335336for i = 1 : length(label )
336337 if myseed(1 ,i )==0
@@ -409,8 +410,7 @@ function selectpoint_Callback(hObject, eventdata, handles)
409410% nonzerovalue = nonzerovalue';
410411% disp(['mean = ',num2str(mean(nonzerovalue)),' v =',num2str(j)]);
411412% %figure;
412- % %plot(nonzerovalue);
413- %
413+ % %plot(nonzerovalue);
414414% if j == 1
415415% numpixel = length(nonzerovalue);
416416% finalmatrix = nonzerovalue;
@@ -423,15 +423,21 @@ function selectpoint_Callback(hObject, eventdata, handles)
423423%
424424% %disp(num(find(temp)));
425425% end
426- %
427- % %disp(finalmatrix);
428426% figure;
429427% imshow(finalmatrix);
430428% for i = 1:c(2,1)
431429% gray(i) = evalue(round(positionX), round(positionY));
432430% end
433431% plot(gray);title(vector);
434432
433+ % D = bwdist(segOutline,'euclidean');
434+ % [px,py] = gradient(D);
435+ % figure;imagesc(D);
436+ % %axis([1,108,90,1]);
437+ % figure;contour(D);
438+
439+
440+
435441
436442% --- Executes on button press in average.
437443function average_Callback(hObject , eventdata , handles )
@@ -635,12 +641,29 @@ function rollback_Callback(hObject, eventdata, handles)
635641% hObject handle to rollback (see GCBO)
636642% eventdata reserved - to be defined in a future version of MATLAB
637643% handles structure with handles and user data (see GUIDATA)
638- global label ;
639- global seed ;
640- global myseed ; % used for plot
641- label = label(2 : end - 1 );
642- seed = seed(2 : end - 1 );
643- myseed = myseed(: ,1 : end - 2 );
644+ % global label;
645+ % global seed;
646+ global myseed ; % used for plot
647+ % label = label(2:end-1);
648+ % seed = seed(2:end-1);
649+ % myseed = myseed(:,1:end-2);
650+ global segOutline ;
651+ global numofpixels ;
652+ global newsegOutline ;
653+ global evalue ;
654+ newsegOutline = Expand(segOutline ,numofpixels );
655+ axes(handles .axes );
656+ hold on
657+ outcome1= evalue + newsegOutline * max(max(evalue ));
658+ imshow(outcome1 ,[]); % outline in original picture
659+ for i = 1 : length(label )
660+ if myseed(1 ,i )==0
661+ plot(myseed(2 ,i ),myseed(3 ,i ),' r.' ,' MarkerSize' ,5 );
662+ else
663+ plot(myseed(2 ,i ),myseed(3 ,i ),' g.' ,' MarkerSize' ,5 );
664+ end
665+ end
666+ hold off
644667
645668
646669% expand the segOutline to 3 pixels
0 commit comments