diff options
Diffstat (limited to 'modules/graphics/macros/polarplot.sci')
| -rwxr-xr-x | modules/graphics/macros/polarplot.sci | 234 |
1 files changed, 234 insertions, 0 deletions
diff --git a/modules/graphics/macros/polarplot.sci b/modules/graphics/macros/polarplot.sci new file mode 100755 index 000000000..89390c160 --- /dev/null +++ b/modules/graphics/macros/polarplot.sci @@ -0,0 +1,234 @@ +// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab +// Copyright (C) INRIA +// Copyright (C) 2010 - DIGITEO - Manuel Juliachs +// This file must be used under the terms of the CeCILL. +// This source file is licensed as described in the file COPYING, which +// you should have received as part of this distribution. The terms +// are also available at +// http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt + +function polarplot(theta,rho,style,strf,leg,rect) + [lhs,rhs]=argn(0) + if rhs<=0 then + theta=0:.01:2*%pi; + rho=sin(2*theta).*cos(2*theta) + clf(); + polarplot(theta,rho) + return + end + if size(theta,1)==1 then + theta=theta(:), + end + if size(rho,1)==1 then + rho=rho(:), + end + rm=max(abs(rho)) + x=rho.*cos(theta) + y=rho.*sin(theta) + + opts=[] + isstrf=%f; + isframeflag=%f; + isrect=%f; + if exists("style","local")==1 then + opts=[opts,"style=style"] + end + if exists("strf","local")==1 then + opts=[opts,"strf=strf"] + isstrf=%t + end + if exists("leg","local")==1 then + opts=[opts,"leg=leg"] + end + if exists("rect","local")==1 then + opts=[opts,"rect=rect"] + isrect=%t + end + if exists("frameflag","local")==1 then + opts=[opts,"frameflag=frameflag"] + isframeflag=%t + end + + if size(opts,2)<rhs-2 then + error(msprintf(gettext("%s: Wrong value for input argument: ''%s'', ''%s'', ''%s'', ''%s'' or ''%s'' expected.\n"),"polarplot","style","strf","leg","rect","frameflag")); + end + + // Some default values: + Amin=0 // starting angle for the frame + dA=360 // span of the angular frame + nn=4 // number of quadrants to be drawn + + xmin=min(x); + xmax=max(x); + L=(xmax-xmin)*1.07; + ymin=min(y); + ymax=max(y); + H=(ymax-ymin)*1.07; + // Angle at which Radial labels will be displayed + A=round(atan((ymin+ymax)/2,(xmin+xmax)/2)/%pi*180/45)*45; + dx=-0.5, dy=-0.5 // H & V shifts in string-width and string-height units + + // Case without rect= + if ~isrect then + // Determines quadrant(s) to be drawn + Q=[%T %T %T %T]; + e=rm/500; + + if min(x)<-e then + xmin=-rm; + else + xmin=0; Q([2 3])=%F; + end + + if max(x)>e then + xmax= rm; + else + xmax=0; Q([1 4])=%F; + end + + if min(y)<-e then + ymin=-rm; + else + ymin=0; Q([3 4])=%F; + end + + if max(y)>e then + ymax= rm; + else + ymax=0; Q([1 2])=%F; + end + + L=(xmax-xmin)*1.1; if L==0, L=2*rm*1.1; end + H=(ymax-ymin)*1.1; if H==0, H=2*rm*1.1; end + x0=(xmin+xmax)/2; y0=(ymin+ymax)/2; + rect=[x0-L/2 y0-H/2 x0+L/2 y0+H/2] + + // Special case: data aligned on the x or y axis + if Q==[%F %F %F %F], + if (ymax-ymin)<2*e, // on x axis + if xmin<-e then + Q([2 3])=%T + end + if xmax> e then + Q([1 4])=%T + end + else // on y axis + if ymin<-e then + Q([3 4])=%T + end + if ymax> e then + Q([1 2])=%T + end + end + end + + n=find(Q); // id numbers of quadrants to be drawn + nn=length(n) // number of quadrants to be drawn + Amin=(n(1)-1)*90 + + select nn + case 1, + dA=90; + if n==1, A=90, dx=-1.1, dy=-0.5 + elseif n==2, A=90, dx=0.2, dy=-0.5 + elseif n==3, A=270, dx=0.2, dy=-0.5 + else A=270, dx=-1.1, dy=-0.5 + end + case 2 + dA=180; + if n(1)==1 + if n(2)==2, //A=90, dx=0.1, dy=-0.5 + else Amin=-90, A=90, dx=-1.2, dy=-0.5, end + elseif n(1)==2, A=90, dx=0.2, dy=-0.5 + else A=0, dx=-0.5, dy=0.2 + end + else + Amin=0, dA=360 + end + opts=[opts,"rect=rect"] + end // if ~isrect + + if isstrf& isframeflag then + error(msprintf(gettext("%s: ''%s'' cannot be used with ''%s''.\n"),"polarplot","frameflag","strf")); + end + if ~(isstrf) then + axesflag=0 + opts=[opts,"axesflag=axesflag"], + end + if ~(isstrf|isframeflag) then + frameflag=4 + opts=[opts,"frameflag=frameflag"], + end + drawlater() + execstr("plot2d(x,y,"+strcat(opts,",")+")") + + fcolor=color("grey70"); + xset("dashes",1) + + // CIRCULAR FRAME: + // Radial values for the frame: + fmt_in=format(), format("v",9) + // Tunning for smart values: + p=floor(log10(abs(rm))); + m=rm/10^p; + if m<1.3, dm=0.2 + elseif m<=2, dm=0.3 + elseif m<4, dm=0.5 + else dm=1, + end + k=fix(m/dm) + if m-k*dm < dm/5, k=k-1, end + R=[(1:k)*dm*10^p ] + // Tuning for smart 10^ display using LaTeX instead of D+## exponential display + if abs(p)<4, + Rtxt=string(R) + [v,k]=max(length(Rtxt)), tmp=xstringl(0,0,Rtxt(k)) + else + if dm<1, dm=dm*10, p=p-1, end + tmp=string(R/10^p)+"108" + [v,k]=max(length(tmp)), tmp=xstringl(0,0,tmp(k)) + Rtxt="$\scriptstyle "+string(R/10^p)+"\:.10^{"+string(p)+"}$"; + end + w=tmp(3); h=tmp(4); + format(fmt_in(2),fmt_in(1)) // Restoring entrance format + R = [ R rm ] + + // Drawing & labelling the radial frame + kM=size(R,"*"); + for k=1:kM + r=R(k) + xarc(-r,r,2*r,2*r,Amin*64,dA*64) + e = gce();,e.line_style=3 + e.foreground=fcolor; + if k==kM + e.line_style=1; // solid outer arc + else + xstring(r*cosd(A)+w*dx, r*sind(A)+h*dy, Rtxt(k)) + e = gce(); + e.clip_state = "off"; + end + end + + // ANGULAR FRAME: + if nn<3, eA=10, else eA=30; end // adaptative angular sampling + an=linspace(Amin,Amin+dA,round(dA/eA)+1); + // avoiding 360 == 0 + if nn>2, tmp=find(abs(an-360)<eA/10); an(tmp)=[]; end + // Adjusting H-shifts of angular labels + tmp=xstringl(0,0,"360"); + w=tmp(3); h=tmp(4); + rL=rm*1.03; // Radius of angular labels + for k=an // draws and labels angular rays + xsegs([0;rm*cosd(k)],[0;rm*sind(k)]) + e = gce(); e.segs_color=fcolor; e.line_style=3; + xstring((rL+w/2)*cosd(k)-w/2, (rL+h/2)*sind(k)-h/2, string(k)) + e = gce(); + e.clip_state = "off"; + end + + a=gca(); + a.data_bounds=[rect(1:2);rect(3:4)] + a.margins=[0.07 0.07 0.12 0.07] + + drawnow() +endfunction |