bryan/hartman-server
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

Format NIcircle.m

Browse Source
This commit is contained in:
Bryan Roessler
2024-07-25 14:49:35 -04:00
parent 9a5c9ffb6a
commit d1b7c87d0b
1 changed files with 26 additions and 29 deletions
Download Patch File Download Diff File Expand all files Collapse all files

55
workflow/templates/easy/NIcircle.m
View File

@@ -3,9 +3,9 @@
%Imaging ToolBox method %Imaging ToolBox method
r=14; r=14;
A=zeros(70,70); %(fIntsc(refPtR:(refPtRExt),refPtC:(refPtCExt))) A=zeros(70,70); %(fIntsc(refPtR:(refPtRExt),refPtC:(refPtCExt)))
m= {40,40}; m={40,40};
A(m{:})=1; A(m{:})=1;
B= imdilate(A,strel('disk',r,0) ); B=imdilate(A,strel('disk',r,0) );
imshow(B) imshow(B)
area=pi*r^2 area=pi*r^2
@@ -18,45 +18,42 @@ clear all
%A=zeros(70,70); %A=zeros(70,70);
%A=zeros(r,r); %A=zeros(r,r);
%P=[40,40]; %P=[40,40];
%center= [refPtR+ round(.5*width), refPtC+ round(.5*width)]; %center=[refPtR+ round(.5*width), refPtC+ round(.5*width)];
%A=zeros(70,70); %A=zeros(70,70);
%--------------------------------------------------------------------- %---------------------------------------------------------------------
%radius=14; %radius=14;
diaExt=2*(radius+1); diaExt=2*(radius+1);
circBoxA=zeros(diaExt,diaExt); circBoxA=zeros(diaExt,diaExt);
center= [radius+2, radius+2]; center=[radius+2, radius+2];
[m, n ]=size(circBoxA); [m,n]=size(circBoxA);
X = bsxfun(@plus,(1:m)', zeros(1,n)); X=bsxfun(@plus,(1:m)', zeros(1,n));
Y = bsxfun(@plus,(1:n), zeros(m,1)); Y=bsxfun(@plus,(1:n), zeros(m,1));
cirMask= sqrt(sum(bsxfun(@minus,cat(3,X,Y),reshape(center,1,1,[])) .^2,3))<=radius; cirMask=sqrt(sum(bsxfun(@minus,cat(3,X,Y),reshape(center,1,1,[])) .^2,3))<=radius;
area=pi*radius^2; area=pi*radius^2;
cirPixA=nnz(cirMask);
cirPixA= nnz(cirMask); optCirMask=double(cirMask);
optCirMask(optCirMask==0)=0.8;
optCirMask= double(cirMask); % Foto Circle Fram(e)
optCirMask(optCirMask==0)= 0.8;
%+++++++++++Foto Circle Fram(e)+++++++++++++++++++++
expansion=2; expansion=2;
radExpan= radius+expansion; radExpan=radius+expansion;
FdiaExt=2*(radExpan); FdiaExt=2*(radExpan);
circBoxA=zeros(FdiaExt,FdiaExt); circBoxA=zeros(FdiaExt,FdiaExt);
center= [radExpan+1, radExpan+1]; center=[radExpan+1, radExpan+1];
[m, n ]=size(circBoxA); [m, n ]=size(circBoxA);
X = bsxfun(@plus,(1:m)', zeros(1,n)); X=bsxfun(@plus,(1:m)', zeros(1,n));
Y = bsxfun(@plus,(1:n), zeros(m,1)); Y=bsxfun(@plus,(1:n), zeros(m,1));
FcirMask= sqrt(sum(bsxfun(@minus,cat(3,X,Y),reshape(center,1,1,[])) .^2,3))<=radExpan; FcirMask=sqrt(sum(bsxfun(@minus,cat(3,X,Y),reshape(center,1,1,[])) .^2,3))<=radExpan;
%FcirPixA= nnz(cirMask); %FcirPixA=nnz(cirMask);
FoptCirMask= double(FcirMask); FoptCirMask=double(FcirMask);
FoptCirMask(FoptCirMask==1)= 2; FoptCirMask(FoptCirMask==1)=2;
%FoptCirMask(FoptCirMask==0)= 1; %FoptCirMask(FoptCirMask==0)=1;
%********Combine Masks to create circular boundry************ %********Combine Masks to create circular boundry************
padOptCirMask= padarray(optCirMask,[expansion-1 expansion-1],0.8); padOptCirMask=padarray(optCirMask,[expansion-1 expansion-1],0.8);
FoptCirMask= FoptCirMask .* padOptCirMask; FoptCirMask=FoptCirMask .* padOptCirMask;
FoptCirMask(FoptCirMask==1.6)= 0.8; FoptCirMask(FoptCirMask==1.6)=0.8;
FoptCirMask(FoptCirMask==0)= 1; FoptCirMask(FoptCirMask==0)=1;
FoptCirMask(FoptCirMask==2)= 1; FoptCirMask(FoptCirMask==2)=1;
%--------------------------------------------------- %---------------------------------------------------
%imagesc(cirMask) %imagesc(cirMask)