Changeset 148 for aedes_readvnmr.m
 Timestamp:
 Nov 3, 2010, 4:04:07 PM (9 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

aedes_readvnmr.m
r139 r148 1120 1120 % Check for arrayed acquisition 1121 1121 if 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 else1122 %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 1128 1128 isAcqArrayed = true; 1129 1129 ArrayLength = []; … … 1138 1138 end 1139 1139 ArrayLength = prod(ArrayLength); 1140 end1140 %end 1141 1141 else 1142 1142 isAcqArrayed = false; … … 1149 1149 % Reconstruct 1D data ... 1150 1150 elseif AcqType==2 1151 1151 1152 % 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]); 1155 1156 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]); 1164 1170 kspace = reshape(kspace,[np/2,nv,ns*ArrayLength,1,nRcvrs]); 1165 1171 end 1166 1172 elseif AcqType==3 1167 1173 % Reconstruct 3D data 1168 if strcmpi(seqcon (2:4),'scs')1174 if strcmpi(seqcon,'nscsn') 1169 1175 kspace = reshape(kspace,[np/2,etl,nv/etl,nRcvrs,ArrayLength*nv2]); 1170 1176 kspace = permute(kspace,[1 2 3 5 4]); … … 1231 1237 lpe2 = procpar.lpe2; 1232 1238 lro = procpar.lro; 1233 if lpe2==01239 if AcqType==2 1234 1240 % 2D data 1235 1241 padSize = [procpar.np/2 ... 1236 1242 procpar.np/2*(lpe/lro) ... 1237 1243 size(kspace,3)]; 1238 else 1244 elseif AcqType==3 && lpe2~=0 1239 1245 % 3D data 1240 1246 padSize = [procpar.np/2 ...
Note: See TracChangeset
for help on using the changeset viewer.