Changeset 782 for trunk/src/uvmat.m


Ignore:
Timestamp:
Jun 3, 2014, 8:10:03 PM (10 years ago)
Author:
sommeria
Message:

AX and AY changed to Coord_x and Coord_y

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/uvmat.m

    r781 r782  
    8383% velocity vectors. The scalar values in the first case is represented by
    8484% UvData.Field.A, and by UvData.Field.C in the second case. The corresponding set of X
    85 % and Y axes are represented by UvData.Field.AX and UvData.Field.AY, and .X and
     85% and Y axes are represented by UvData.Field.Coord_x and UvData.Field.Coord_y, and .X and
    8686% .Y for C (the same as velocity vectors). If A is a nxxny matrix (scalar
    87 % on a regtular grid), then .AX andf.AY contains only two elements, represneting the
     87% on a regtular grid), then .Coord_x andf.Coord_y contains only two elements, represneting the
    8888% axes of the four image corners. The scalar name is represented by
    8989% the strings .AName and/or .CName.
     
    11451145        set(handles_Fields,'String',{'image'})
    11461146        %set(handles.Coord_x,'Value',1);
    1147         set(handles.Coord_x,'String','AX');
    1148     set(handles.Coord_y,'String','AY');
     1147        set(handles.Coord_x,'String','Coord_x');
     1148    set(handles.Coord_y,'String','Coord_y');
    11491149end
    11501150set(handles.uvmat,'UserData',UvData)
     
    14971497            set(hmask,'CData',imflag)   
    14981498            set(hmask,'AlphaData',flagmask*0.6)
    1499             set(hmask,'XData',MaskField.AX);
    1500             set(hmask,'YData',MaskField.AY);
     1499            set(hmask,'XData',MaskField.Coord_x);
     1500            set(hmask,'YData',MaskField.Coord_y);
    15011501%             uistack(hmask,'top')
    15021502        else
    15031503            axes(handles.PlotAxes)
    15041504            hold on   
    1505             Mask.maskhandle=image(MaskField.AX,MaskField.AY,imflag,'Tag','mask','HitTest','off','AlphaData',0.6*ones(size(flagmask)));
     1505            Mask.maskhandle=image(MaskField.Coord_x,MaskField.Coord_y,imflag,'Tag','mask','HitTest','off','AlphaData',0.6*ones(size(flagmask)));
    15061506            set(handles.CheckMask,'UserData',Mask)
    15071507        end
     
    26912691        if isfield(ObjectData,'NbDim') && isequal(ObjectData.NbDim,2) && isfield(ObjectData,'Mask') && isfield(ObjectData,'A')
    26922692            flag_mask=double(ObjectData.Mask>200);%=0 for masked regions
    2693             AX=ObjectData.AX;%x coordiantes for the scalar field
    2694             AY=ObjectData.AY;%y coordinates for the scalar field
     2693            Coord_x=ObjectData.Coord_x;%x coordiantes for the scalar field
     2694            Coord_y=ObjectData.Coord_y;%y coordinates for the scalar field
    26952695            MaskX=ObjectData.MaskX;%x coordiantes for the mask
    26962696            MaskY=ObjectData.MaskY;%y coordiantes for the mask
    2697             if ~isequal(MaskX,AX)||~isequal(MaskY,AY)
     2697            if ~isequal(MaskX,Coord_x)||~isequal(MaskY,Coord_y)
    26982698                nxy=size(flag_mask);
    26992699                sizpx=(ObjectData.MaskX(end)-ObjectData.MaskX(1))/(nxy(2)-1);%size of a mask pixel
     
    27032703                %project on the positions of the scalar
    27042704                npxy=size(ObjectData.A);
    2705                 dxy(1)=(ObjectData.AY(end)-ObjectData.AY(1))/(npxy(1)-1);%grid mesh in y
    2706                 dxy(2)=(ObjectData.AX(end)-ObjectData.AX(1))/(npxy(2)-1);%grid mesh in x
    2707                 xi=ObjectData.AX(1):dxy(2):ObjectData.AX(end);
    2708                 yi=ObjectData.AY(1):dxy(1):ObjectData.AY(end);
     2705                dxy(1)=(ObjectData.Coord_y(end)-ObjectData.Coord_y(1))/(npxy(1)-1);%grid mesh in y
     2706                dxy(2)=(ObjectData.Coord_x(end)-ObjectData.Coord_x(1))/(npxy(2)-1);%grid mesh in x
     2707                xi=ObjectData.Coord_x(1):dxy(2):ObjectData.Coord_x(end);
     2708                yi=ObjectData.Coord_y(1):dxy(1):ObjectData.Coord_y(end);
    27092709                [XI,YI]=meshgrid(xi,yi);% creates the matrix of regular coordinates
    27102710                flag_mask = interp2(x_mask,y_mask,flag_mask,XI,YI);
     
    36413641        data.DX=UvData.Field.CoordMesh;
    36423642        data.DY=UvData.Field.CoordMesh;
    3643     elseif isfield(UvData.Field,'AX')&isfield(UvData.Field,'AY')& isfield(UvData.Field,'A')%only image
     3643    elseif isfield(UvData.Field,'Coord_x')&isfield(UvData.Field,'Coord_y')& isfield(UvData.Field,'A')%only image
    36443644        np=size(UvData.Field.A);
    3645         meshx=(UvData.Field.AX(end)-UvData.Field.AX(1))/np(2);
    3646         meshy=abs(UvData.Field.AY(end)-UvData.Field.AY(1))/np(1);
     3645        meshx=(UvData.Field.Coord_x(end)-UvData.Field.Coord_x(1))/np(2);
     3646        meshy=abs(UvData.Field.Coord_y(end)-UvData.Field.Coord_y(1))/np(1);
    36473647        data.RangeY=max(meshx,meshy);
    36483648        data.RangeX=max(meshx,meshy);
     
    53815381    figure
    53825382    plot(theta_ref,r_ref)
    5383     azimuth_ima=linspace(DataOut.AY(1),DataOut.AY(2),size(DataOut.A,1));%profile of x index on the transformed image
     5383    azimuth_ima=linspace(DataOut.Coord_y(1),DataOut.Coord_y(2),size(DataOut.A,1));%profile of x index on the transformed image
    53845384    dist_source = interp1(theta_ref,r_ref,azimuth_ima);
    5385     dist_source_pixel=round(size(DataOut.A,2)*(dist_source-DataOut.AX(1))/(DataOut.AX(2)-DataOut.AX(1)));
     5385    dist_source_pixel=round(size(DataOut.A,2)*(dist_source-DataOut.Coord_x(1))/(DataOut.Coord_x(2)-DataOut.Coord_x(1)));
    53865386    line_nan= isnan(dist_source_pixel);
    53875387    dist_source_pixel(line_nan)=1;
Note: See TracChangeset for help on using the changeset viewer.