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

Detach ATT from UR

Browse Source
This commit is contained in:
kradchen
2023-06-15 15:43:18 +08:00
parent 8e3db6877c
commit 6ab314d426
2 changed files with 33 additions and 29 deletions
Download Patch File Download Diff File Expand all files Collapse all files

54
src/reflectionReconstruction/reconstructionSAFT/reconstructionSAFT.cpp
View File

@@ -5,7 +5,7 @@
#include "Function3D.h"
#include "Matrix.h"
#include "SAFTStructs.h"
#include "SAFT_ATT.h"
// #include "SAFT_ATT.h"
#include "SAFT_TOFI.h"
#include "common/dataBlockCreation/removeDataFromArrays.h"
#include "config/config.h"
@@ -95,8 +95,8 @@ namespace Recon {
Matrix_t receiverGeom{nullptr,
(size_t)receiverPosGeom.getDims(),
{(size_t)receiverPosGeom.getDimSize(1),
(size_t)receiverPosGeom.getDimSize(0),
{(size_t)receiverPosGeom.getDimSize(0),
(size_t)receiverPosGeom.getDimSize(1),
(size_t)receiverPosGeom.getDimSize(2)},
receiverPosGeom.getDataSize()};
float *fdata = new float[receiverPosGeom.getDataSize()]{0};
@@ -107,7 +107,7 @@ namespace Recon {
Matrix_t senderGeom{nullptr,
(size_t)senderPosGeom.getDims(),
{(size_t)senderPosGeom.getDimSize(1),
{(size_t)senderPosGeom.getDimSize(0),
(size_t)senderPosGeom.getDimSize(1),
(size_t)senderPosGeom.getDimSize(2)},
senderPosGeom.getDataSize()};
@@ -122,7 +122,7 @@ namespace Recon {
int saftVariant_v[6]{(int)reflectParams::saftVariant[0],(int)reflectParams::saftVariant[1],(int)reflectParams::saftVariant[2],
(int)reflectParams::saftVariant[3],(int)reflectParams::saftVariant[4],(int)reflectParams::saftVariant[5]};
Matrix_t saftVariant{saftVariant_v,2,{1,6,1},3};
Matrix_t saftVariant{saftVariant_v,2,{1,6,1},6};
params.push_back(saftVariant);
Matrix_t SpeedMap3D{nullptr,
@@ -216,20 +216,22 @@ namespace Recon {
Matrix_t otherP{other,2,{1,2,1},2};
params.push_back(otherP);
if (reflectParams::useAscanIndex == 0) {
auto result = SAFT_ATT(params);
delete [] ASdata;
delete [] (float*)imageStartpoint.Data;
delete [] temp1;
delete [] temp2;
delete [] fdata;
delete [] fdata1;
delete [] fdata2;
delete [] fdata3;
delete [] fdata4;
delete [] fdata5;
delete [] gdata;
return Aurora::Matrix::fromRawData((double *)result.Data, result.Dims[0],
result.Dims[1], result.Dims[2]);
// auto result = SAFT_ATT(params);
// delete [] ASdata;
// delete [] (float*)imageStartpoint.Data;
// delete [] temp1;
// delete [] temp2;
// delete [] fdata;
// delete [] fdata1;
// delete [] fdata2;
// delete [] fdata3;
// delete [] fdata4;
// delete [] fdata5;
// delete [] gdata;
// return Aurora::Matrix::fromRawData((double *)result.Data, result.Dims[0],
// result.Dims[1], result.Dims[2]);
return Aurora::Matrix();
} else {
auto result = SAFT_TOFI(params);
delete [] ASdata;
@@ -263,10 +265,10 @@ namespace Recon {
auto senderPos_s = Aurora::transpose(removeDataFromArrays(senderPos, mpIdxs));
auto receiverPos_s = Aurora::transpose(removeDataFromArrays(receiverPos, mpIdxs));
Aurora::Matrix idsSender ;
auto senderPosGeom = Aurora::uniqueByRows(receiverPos_s, idsSender);
auto senderPosGeom = Aurora::transpose(Aurora::uniqueByRows(senderPos_s, idsSender));
idsSender = Aurora::transpose(idsSender);
Aurora::Matrix idsReceiver ;
auto receiverPosGeom = Aurora::uniqueByRows(receiverPos_s, idsReceiver);
auto receiverPosGeom = Aurora::transpose(Aurora::uniqueByRows(receiverPos_s, idsReceiver));
idsReceiver = Aurora::transpose(idsReceiver);
auto numUsedData = countMp;
@@ -280,14 +282,16 @@ namespace Recon {
size_t length = blockIdxs[i + 1] - blockIdxs[i];
size_t begin = blockIdxs[i];
size_t end = blockIdxs[i + 1];
auto _AScans = Aurora::zeros(AScans.getDimSize(0), length);
auto receiverIndex = Aurora::zeros(1, length);
auto senderIndex = Aurora::zeros(1, length);
auto mpIdxs_d = Aurora::zeros(1,mpIdxs.getDataSize());
for (size_t j = begin, k = 0; j < end; j++, k++) {
_AScans(Aurora::$, k) = AScans(Aurora::$, mpIdxs[j]);
receiverIndex[k] = idsReceiver[j];
senderIndex[k] = idsSender[j];
mpIdxs_d[j] =mpIdxs[j];
receiverIndex[k] = idsReceiver[j]+1;
senderIndex[k] = idsSender[j]+1;
}
auto _AScans = removeDataFromArrays(AScans, mpIdxs_d);
Env = callSAFT(_AScans, receiverIndex, senderIndex, receiverPosGeom,
senderPosGeom, SAFT_mode, TimeInterval, transRecos, Env);
}