Changeset 148 for aedes_readvnmr.m


Ignore:
Timestamp:
Nov 3, 2010, 4:04:07 PM (9 years ago)
Author:
tjniskan
Message:
  • Added aedes_roifill.m for doing binary flood fill operation. This is

much slower than using the imfill-function, but does not dependend on
Image Processing toolbox. The goal is to make Aedes run with just the
base Matlab install. DICOM support will probably be the last obstacle
in achieving this, because it's going to be a great pain in the a to
code...

  • Fixed overwriting ROIs in aedes_roi_copy_gui.m and added additional

information in the ROI copy button tooltips.

  • Some further iterations to the new readvnmr-function.

M aedes_roi_copy_gui.m
M aedes_getmatlabversion.m
M aedes_readvnmr.m
M aedes.m
M aedes_helpabout.m
M aedes_readfid.m
M aedes_revision.m
A aedes_roifill.m

File:
1 edited

Legend:

Unmodified
Added
Removed
  • aedes_readvnmr.m

    r139 r148  
    11201120% Check for arrayed acquisition
    11211121if isfield(procpar,'array') && ~isempty(procpar.array{1})
    1122   if length(procpar.array)==1 && ~iscell(procpar.array{1}) && ...
    1123       strcmp(procpar.array{1},'pad') && all(procpar.pad==0)
    1124     % Skip the array parsing if the array is a dummy "pad"-array...
    1125     isAcqArrayed = false;
    1126     ArrayLength = 1;
    1127   else
     1122  %if length(procpar.array)==1 && ~iscell(procpar.array{1}) && ...
     1123  %    strcmp(procpar.array{1},'pad') && all(procpar.pad==0)
     1124  %  % Skip the array parsing if the array is a dummy "pad"-array...
     1125  %  isAcqArrayed = false;
     1126  %  ArrayLength = 1;
     1127  %else
    11281128    isAcqArrayed = true;
    11291129    ArrayLength = [];
     
    11381138    end
    11391139    ArrayLength = prod(ArrayLength);
    1140   end
     1140  %end
    11411141else
    11421142  isAcqArrayed = false;
     
    11491149  % Reconstruct 1D data ...
    11501150elseif AcqType==2
     1151 
    11511152  % Reconstruct 2D data
    1152   if strcmpi(seqcon(2:3),'cs')
    1153     kspace = reshape(kspace,[np/2,etl,ArrayLength*ns,nRcvrs,nv/etl]);
    1154     kspace = permute(kspace,[1 2 5 3 4]);
     1153  if strcmpi(seqcon,'ncsnn')
     1154    kspace = reshape(kspace,[np/2,etl,ns,ArrayLength,nv/etl,nRcvrs]);
     1155    kspace = permute(kspace,[1 2 5 3 4 6]);
    11551156    kspace = reshape(kspace,[np/2,nv,ns*ArrayLength,1,nRcvrs]);
    1156   elseif strcmpi(seqcon(2:3),'sc')
    1157     %kspace = reshape(kspace,[np/2,ArrayLength*ns,nv,nRcvrs]);
    1158     kspace = reshape(kspace,[np/2,ns,nRcvrs,ArrayLength,nv]);
    1159     kspace = permute(kspace,[1 5 2 4 3]);
    1160     kspace = reshape(kspace,[np/2,nv,ns*ArrayLength,1,nRcvrs]);
    1161   elseif strcmpi(seqcon(2:3),'cc')
    1162     kspace = reshape(kspace,[np/2,nv,nRcvrs,ns*ArrayLength]);
    1163     kspace = permute(kspace,[1 2 4 3]);
     1157  elseif strcmpi(seqcon,'nscnn')
     1158    if isfield(procpar,'flash_converted')
     1159      kspace = reshape(kspace,[np/2,nRcvrs,ArrayLength,nv,ns]);
     1160      kspace = permute(kspace,[1 4 5 3 2]);
     1161      kspace = reshape(kspace,[np/2,nv,ns*ArrayLength,1,nRcvrs]);
     1162    else
     1163      kspace = reshape(kspace,[np/2,nRcvrs,nv,ArrayLength,ns]);
     1164      kspace = permute(kspace,[1 3 5 4 2]);
     1165      kspace = reshape(kspace,[np/2,nv,ns*ArrayLength,1,nRcvrs]);
     1166    end
     1167  elseif strcmpi(seqcon,'nccnn')
     1168    kspace = reshape(kspace,[np/2,etl,nRcvrs,ns,nv/etl,ArrayLength]);
     1169    kspace = permute(kspace,[1 2 5 4 6 3]);
    11641170    kspace = reshape(kspace,[np/2,nv,ns*ArrayLength,1,nRcvrs]);
    11651171  end
    11661172elseif AcqType==3
    11671173  % Reconstruct 3D data
    1168   if strcmpi(seqcon(2:4),'scs')
     1174  if strcmpi(seqcon,'nscsn')
    11691175    kspace = reshape(kspace,[np/2,etl,nv/etl,nRcvrs,ArrayLength*nv2]);
    11701176    kspace = permute(kspace,[1 2 3 5 4]);
     
    12311237    lpe2 = procpar.lpe2;
    12321238    lro = procpar.lro;
    1233     if lpe2==0
     1239    if AcqType==2
    12341240      % 2D data
    12351241      padSize = [procpar.np/2 ...
    12361242        procpar.np/2*(lpe/lro) ...
    12371243        size(kspace,3)];
    1238     else
     1244    elseif AcqType==3 && lpe2~=0
    12391245      % 3D data
    12401246      padSize = [procpar.np/2 ...
Note: See TracChangeset for help on using the changeset viewer.

Powered by Trac 1.0.9.Copyright © Juha-Pekka Niskanen 2008