arrays - Rotation of 3D data in XY plane (Matlab) -

-

i have 3d array contains probability values function of spatial position (x , y coordinates). when plot array using imagesc or surf, elliptic profile maximum. rotate data around point of maximum probability in xy plane such y axis becomes parallel major axis of ellipse (which computed using regionprops). don't want rotate image or surface using imrotate or geometric object rotation, since need array containing rotated data. how achieve if know angle of rotation around z axis (obtained regionprops)?

edit: here code used rotate matrix around centroïd. however, when plot result using surf, weirdest result...

% translation , anti-translation matrix translate = [1 0 0 -centroide(1); 0 1 0 -centroide(2); 0 0 1 0; 0 0 0 1]; anti_translate = [1 0 0 centroide(1); 0 1 0 centroide(2); 0 0 1 0; 0 0 0 1];  % rotation matrix rotate = [cosd(-orientation) -sind(-orientation) 0 0; sind(-orientation)   cosd(-orientation) 0 0; 0 0 1 0; 0 0 0 1];  % temporary cell array data transformation data_vector = cell(size(rcb,1),size(rcb,2),size(rcb,3)); % cell array transformed data storage transformed_data = cell(size(rcb,3),1); = 1:(size(rcb,3))     transformed_data{i} = zeros(3, size(rcb,1)*size(rcb,2)); end = 1:size(rcb,3) % temporal dimension     j = 1:size(rcb,2) % x dimension          k = 1:size(rcb,1) % y dimension             % [x y z] vector             data_vector{k,j,i} = [j,k,rcb(k,j,i),1]';             % translation             data_vector{k,j,i} = translate*data_vector{k,j,i};             % rotation             data_vector{k,j,i} = rotate*data_vector{k,j,i};             % anti-translation             data_vector{k,j,i} = anti_translate*data_vector{k,j,i};             data_vector{k,j,i} = data_vector{k,j,i}(1:3);             transformed_data{i}(:,((j-1)*size(rcb,1)+k)) = data_vector{k,j,i};         end     end end  %%%% test %%%% : image @ 100 ms figure, imagesc(transformed_data{100}(2,:),transformed_data{100}(1,:),transformed_data{100}(3,:));

i image:enter image description here

what doing wrong?

have tried rotz rotation matrix rotations around z-axis , check documentation here.

othewise, can multiply matrix following matrix r rotation angle γ :

enter image description here

be careful, function cos , sin in matlab radian based. can use degree cosd.


Comments

Post a Comment

Popular posts from this blog

PHP DOM loadHTML() method unusual warning -

-
i had following code works fine on localhost: $document = new domdocument(); $document->loadhtml($p_result); $form = $document->getelementsbytagname('form')->item(0); // code continues using $form variable after same code updated on outside server, loadhtml() failed , issued warning. warning: domdocument::loadhtml() expects parameter 1 valid path, string given in path/name/to/script.php returned null instead of object code pretty gets fatal error. note contents of $p_result absolutely same on outside server , on localhost. but why displays kind of warning , why not work? doesn't loadhtml() expects argument 1 string in first place? why method expects parameter 1 valid path ? just make clear i'm not calling loadhtmlfile() , i'm calling loadhtml() . thanks. you're affected one of php bugs . issue present in php 5.6.8 , 5.6.9. have affected php version on server, , bug-free version on localhost. the bug forbids null
Read more

python - How to create jsonb index using GIN on SQLAlchemy? -

-
here's current code creating index jsonb. index("mytable_data_idx_id_key", mytable.data['id'].astext, postgresql_using='gin') but got error. sqlalchemy.exc.programmingerror: (psycopg2.programmingerror) data type text has no default operator class access method "gin" hint: must specify operator class index or define default operator class data type. [sql: "create index event_data_idx_id_key on event using gin ((data ->> 'id'))"] is there way create index on sqlalchemy? the postgresql specific sqlalchemy docs @ http://docs.sqlalchemy.org/en/latest/dialects/postgresql.html#operator-classes mention postgresql_ops dictionary provide "operator class" used postgresql, , provide example illustrating use: index('my_index', my_table.c.id, my_table.c.data, postgresql_ops={ 'data': 'text_pattern_ops',
Read more

c# - TransactionScope not rolling back although no complete() is called -

-
i'm using transactionscope rollback transaction fail bool errorreported = false; action<importerrorlog> newerrorcallback = e => { errorreported = true; errorcallback(e); }; using (var transaction = new transactionscope()) { foreach (importtaskdefinition task in taskdefinition) { loader.load(streamfile, newerrorcallback, task.destinationtable, processingtaskid); } if (!errorreported) transaction.complete(); } i'm sure there no transactionscope started ahead or after code. i'm using entity framework insert in db. regardless state of errorreported transaction never rolled in case of error. what missing ? transactionscope sets transaction.current . that's does. wants transacted must @ property. i believe ef each time connection opened reason. probably, connection open when scope installed. open connection inside of scope or enlist manually. ef has nasty "design decision": default opens
Read more