Changeset 402 for trunk/src/proj_field.m


Ignore:
Timestamp:
Apr 29, 2012, 11:29:24 AM (12 years ago)
Author:
sommeria
Message:

bugs corrected and improved procedure for object projection

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/proj_field.m

    r399 r402  
    8282function [ProjData,errormsg]=proj_field(FieldData,ObjectData)
    8383errormsg='';%default
    84 if ~exist('FieldName','var')
    85     FieldName='';
    86 end
     84% if ~exist('FieldName','var')
     85%     FieldName='';
     86% end
    8787%% case of no projection (object is used only as graph display)
    8888if isfield(ObjectData,'ProjMode') && (isequal(ObjectData.ProjMode,'none')||isequal(ObjectData.ProjMode,'mask_inside')||isequal(ObjectData.ProjMode,'mask_outside'))
     
    11951195            case 'filter'
    11961196                if ~isempty(VarType.coord_tps)
    1197                     VarType.coord_tps
     1197                     %Coord_tps=FieldData.ListVarName{VarType.coord_tps};
     1198                     %TODO: possibly translate and rotate coordiantes translate  initial coordinates
    11981199                    coord_x_proj=XMin:DX:XMax;
    11991200                    coord_y_proj=YMin:DY:YMax;
     
    12011202                    np_y=numel(coord_y_proj);
    12021203                    [XI,YI]=meshgrid(coord_x_proj,coord_y_proj');
    1203                     XI=reshape(XI,[],1);
    1204                     YI=reshape(YI,[],1);
     1204                    XI=reshape(XI,[],1)+ObjectData.Coord(1,1);
     1205                    YI=reshape(YI,[],1)+ObjectData.Coord(1,2);
    12051206                    ProjData=calc_field(FieldData.FieldList,FieldData,[XI YI]);
    12061207                    for ilist=3:length(ProjData.ListVarName)% reshape data, excluding coordinates (ilist=1-2), TODO: rationalise
Note: See TracChangeset for help on using the changeset viewer.