Skip to content

GitLab

Commit eff03c5d authored by smaretz's avatar smaretz
Browse files

debugged

parent d65598d5
Pipeline #100823 passed with stage
in 1 minute and 22 seconds
Hide whitespace changes
Inline Side-by-side
functions/tomography/astraWrappers/astraFDK.m
View file @ eff03c5d
...@@ -158,6 +158,9 @@ numChunks = ceil(settings.outputSize.^2 * settings.numSlices / MAX_VOLUME_SIZE); ...@@ -158,6 +158,9 @@ numChunks = ceil(settings.outputSize.^2 * settings.numSlices / MAX_VOLUME_SIZE);
if numChunks > 1 if numChunks > 1
vol = zeros([settings.outputSize*[1,1], settings.numSlices], 'single'); vol = zeros([settings.outputSize*[1,1], settings.numSlices], 'single');
slicesPerChunk = ceil(settings.numSlices/numChunks); slicesPerChunk = ceil(settings.numSlices/numChunks);
% convert input to single precision to avoid copying in each call of astraFDK in the loop below
projs = single(projs);
for chunkIdx = 1:numChunks for chunkIdx = 1:numChunks
firstSlice = 1+(chunkIdx-1)*slicesPerChunk; firstSlice = 1+(chunkIdx-1)*slicesPerChunk;
...@@ -166,7 +169,7 @@ if numChunks > 1 ...@@ -166,7 +169,7 @@ if numChunks > 1
settingsChunk = settings; settingsChunk = settings;
settingsChunk.numSlices = lastSlice - firstSlice + 1; settingsChunk.numSlices = lastSlice - firstSlice + 1;
settingsChunk.offset = settings.offset + (midSlice-1) - (numSlices-1)/2; settingsChunk.offset = settings.offset + (midSlice-1) - (settings.numSlices-1)/2;
vol(:,:,firstSlice:lastSlice) = astraFDK(projs, tomoAngles, z01, z02, dx, settingsChunk); vol(:,:,firstSlice:lastSlice) = astraFDK(projs, tomoAngles, z01, z02, dx, settingsChunk);
end end
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment