sunwen/UR
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Fix discretizePositions

Browse Source
This commit is contained in:
kradchen
2023-05-31 14:43:09 +08:00
parent 1170b72d04
commit 30ea494d51
3 changed files with 17 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
src/transmissionReconstruction/reconstruction/reconstruction.cpp
View File

@@ -83,10 +83,22 @@ namespace Recon {
result.res = calculateResolution(result.ddims, result.dims); result.res = calculateResolution(result.ddims, result.dims);
int numDim = aVSenderCoordList.getDimSize(0); int numDim = aVSenderCoordList.getDimSize(0);
auto tempDdims = result.ddims.block(0,0,numDim-1); auto tempDdims = result.ddims.block(0,0,numDim-1);
auto tempTransRes = repmat(transpose(result.res),1,2); auto tempTransRes = repmat(transpose(result.res),1,aVSenderCoordList.getDimSize(1));
auto tempDivRes = repmat(transpose(tempDdims / result.res),1,2); auto tempDivRes = repmat(transpose(tempDdims / result.res),1,aVSenderCoordList.getDimSize(1));
aVSenderCoordList = (aVSenderCoordList / tempTransRes) - tempDivRes + 1; aVSenderCoordList = (aVSenderCoordList / tempTransRes) - tempDivRes + 1;
aVReceiverCoordList = (aVReceiverCoordList / tempTransRes) - tempDivRes + 1; aVReceiverCoordList = (aVReceiverCoordList / tempTransRes) - tempDivRes + 1;
{
float* temp = new float[aVReceiverCoordList.getDataSize()]{0};
std::copy(aVReceiverCoordList.getData(),aVReceiverCoordList.getData()+aVReceiverCoordList.getDataSize(),temp);
std::copy(temp,temp+aVReceiverCoordList.getDataSize(),aVReceiverCoordList.getData());
delete [] temp;
}
{
float* temp = new float[aVSenderCoordList.getDataSize()]{0};
std::copy(aVSenderCoordList.getData(),aVSenderCoordList.getData()+aVSenderCoordList.getDataSize(),temp);
std::copy(temp,temp+aVSenderCoordList.getDataSize(),aVSenderCoordList.getData());
delete [] temp;
}
result.receiverCoordList = aVReceiverCoordList; result.receiverCoordList = aVReceiverCoordList;
result.senderCoordList = aVSenderCoordList; result.senderCoordList = aVSenderCoordList;
return result; return result;

3
test/Reconstruction_Test.cpp
View File

@@ -71,12 +71,13 @@ TEST_F(Reconstruction_Test, discretizePositions) {
MatlabReader m("/home/sun/testData/discretizePositions.mat"); MatlabReader m("/home/sun/testData/discretizePositions.mat");
auto senderList = m.read("senderList"); auto senderList = m.read("senderList");
auto receiverList = m.read("receiverList"); auto receiverList = m.read("receiverList");
auto result = Recon::discretizePositions(senderList,receiverList,Recon::transParams::numPixelXY);
auto senderListResult = m.read("senderListResult"); auto senderListResult = m.read("senderListResult");
auto receiverListResult = m.read("receiverListResult"); auto receiverListResult = m.read("receiverListResult");
auto dims = m.read("dims"); auto dims = m.read("dims");
auto ddims = m.read("ddims"); auto ddims = m.read("ddims");
auto res = m.read("res"); auto res = m.read("res");
auto result = Recon::discretizePositions(senderList,receiverList,Recon::transParams::numPixelXY);
EXPECT_DOUBLE_AE(senderListResult.getDataSize(), result.senderCoordList.getDataSize()); EXPECT_DOUBLE_AE(senderListResult.getDataSize(), result.senderCoordList.getDataSize());
for(size_t i=0; i<senderListResult.getDataSize(); ++i) for(size_t i=0; i<senderListResult.getDataSize(); ++i)