first commnit
This commit is contained in:
246
othersoft/calibration_console/Source_Files/calibration.cpp
Normal file
246
othersoft/calibration_console/Source_Files/calibration.cpp
Normal file
@ -0,0 +1,246 @@
|
||||
#include "Header_Files/calibration.h"
|
||||
|
||||
CalibrationAlgorithm::CalibrationAlgorithm()
|
||||
{
|
||||
|
||||
}
|
||||
|
||||
CalibrationAlgorithm::~CalibrationAlgorithm()
|
||||
{
|
||||
|
||||
}
|
||||
|
||||
void CalibrationAlgorithm::readAndResample_StandardLightFile(QString filePath, int integratingSphereDetectorValue, DeviceAttribute deviceAttribute, DeviceInfo deviceInfo)
|
||||
{
|
||||
|
||||
QFile file(filePath);
|
||||
if (!file.open(QIODevice::ReadOnly | QIODevice::Text))
|
||||
{
|
||||
std::cout << "<EFBFBD>ļ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ڣ<EFBFBD>" << std::endl;
|
||||
|
||||
return;
|
||||
}
|
||||
|
||||
//<2F><>ȡ<EFBFBD><C8A1><EFBFBD><D7BC><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
||||
int lineCount = 0;
|
||||
while (!file.atEnd())
|
||||
{
|
||||
QByteArray tmp = file.readLine();
|
||||
lineCount++;
|
||||
}
|
||||
double * StandardLightWavelength_tmp = new double[lineCount - 1];
|
||||
double * StandardLightData_tmp = new double[lineCount - 1];
|
||||
|
||||
|
||||
file.seek(0);
|
||||
for (size_t i = 0; i < lineCount; i++)
|
||||
{
|
||||
QByteArray line = file.readLine();
|
||||
QString str(line);
|
||||
//cout << str.section('\t', 1).trimmed().toStdString() << endl;
|
||||
|
||||
if (i == 0)
|
||||
{
|
||||
QString first = str.section('\t', 0, 0);
|
||||
m_dStandardLightDataBase = first.toDouble();
|
||||
}
|
||||
else
|
||||
{
|
||||
QString first = str.section('\t', 0, 0);
|
||||
QString second = str.section('\t', 1, 1);
|
||||
StandardLightWavelength_tmp[i - 1] = first.toDouble();
|
||||
StandardLightData_tmp[i - 1] = second.toDouble();
|
||||
|
||||
//if (i== lineCount-1)//<2F>鿴<EFBFBD><E9BFB4><EFBFBD><EFBFBD>һ<EFBFBD><D2BB><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>Ƿ<EFBFBD><C7B7><EFBFBD>ȷ
|
||||
//{
|
||||
// double xx = first.toDouble();
|
||||
// double yy = second.toDouble();
|
||||
// std::cout << "xx<78><78>" << xx <<std::endl;
|
||||
// std::cout << "yy<79><79>" << yy << std::endl;
|
||||
//}
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
|
||||
//<2F><>ȡ<EFBFBD><C8A1><EFBFBD><D7BC><EFBFBD><EFBFBD><EFBFBD>ݵ<EFBFBD><DDB5><EFBFBD>Ч<EFBFBD><D0A7>
|
||||
int startPos, endPos;
|
||||
int buffer = 3;//<2F><>Ҫ<EFBFBD><D2AA><EFBFBD>DZ<EFBFBD><EFBFBD><D7BC>ԭʼ<D4AD>ļ<EFBFBD><C4BC><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
||||
//QString biaozhundengfanwei = QString::number(StandardLightWavelength_tmp[0]) + "--" + QString::number(StandardLightWavelength_tmp[lineCount - 2]);
|
||||
|
||||
if (deviceAttribute.fWaveLengthInNM[0] < StandardLightWavelength_tmp[0])//<2F><><EFBFBD><D7BC><EFBFBD>ļ<EFBFBD><C4BC><EFBFBD>Χδ<CEA7><CEB4><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>С<EFBFBD><D0A1><EFBFBD><EFBFBD>
|
||||
{
|
||||
startPos = 0;
|
||||
}
|
||||
else
|
||||
{
|
||||
for (size_t i = 0; i < lineCount - 1; i++)
|
||||
{
|
||||
if (deviceAttribute.fWaveLengthInNM[0] < StandardLightWavelength_tmp[i])
|
||||
{
|
||||
startPos = i - buffer;
|
||||
|
||||
break;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
if (deviceAttribute.fWaveLengthInNM[deviceAttribute.iPixels - 1] > StandardLightWavelength_tmp[lineCount - 2])//<2F><><EFBFBD><D7BC><EFBFBD>ļ<EFBFBD><C4BC><EFBFBD>Χδ<CEA7><CEB4><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
||||
{
|
||||
endPos = lineCount - 2;
|
||||
}
|
||||
else
|
||||
{
|
||||
for (size_t i = 0; i < lineCount - 1; i++)
|
||||
{
|
||||
if (deviceAttribute.fWaveLengthInNM[deviceAttribute.iPixels - 1] < StandardLightWavelength_tmp[i])
|
||||
{
|
||||
endPos = i + buffer;//??
|
||||
|
||||
break;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
m_dStandardLightWavelength = new double[endPos - startPos];
|
||||
m_dStandardLightData = new double[endPos - startPos];
|
||||
|
||||
for (size_t i = 0; i < endPos - startPos; i++)
|
||||
{
|
||||
m_dStandardLightWavelength[i] = StandardLightWavelength_tmp[i + startPos];
|
||||
m_dStandardLightData[i] = StandardLightData_tmp[i + startPos];
|
||||
}
|
||||
|
||||
|
||||
|
||||
// //<2F><><EFBFBD>ضϱ<D8B6><EFBFBD><D7BC><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>csv
|
||||
// QFileInfo fileInfo(filePath);
|
||||
// QString standardLightFileFolder = fileInfo.path();
|
||||
// QString standardLightFileName = fileInfo.fileName();
|
||||
//
|
||||
// QDateTime curDateTime = QDateTime::currentDateTime();
|
||||
// QString currentTime = curDateTime.toString("yyyy_MM_dd_hh_mm_ss");
|
||||
//
|
||||
// QString tmp = standardLightFileFolder + "/" + currentTime + "_" + QString::fromStdString(deviceInfo.strSN) +"_"+ standardLightFileName + "_truncation.csv";
|
||||
//
|
||||
// std::ofstream outfile1(tmp.toStdString().c_str());
|
||||
//
|
||||
// for (size_t i = 0; i < endPos - startPos; i++)
|
||||
// {
|
||||
// if (i == 0)
|
||||
// {
|
||||
// outfile1 << m_dStandardLightDataBase << std::endl;
|
||||
// }
|
||||
//
|
||||
// outfile1 << m_dStandardLightWavelength[i] << "," << m_dStandardLightData[i] << std::endl;
|
||||
// }
|
||||
// outfile1.close();
|
||||
|
||||
|
||||
|
||||
//<2F>ز<EFBFBD><D8B2><EFBFBD><EFBFBD><EFBFBD><EFBFBD><D7BC><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
||||
Eigen::VectorXd vx = Eigen::VectorXd::Map(m_dStandardLightWavelength, endPos - startPos/*x_vec.size()*/);
|
||||
Eigen::VectorXd vy = Eigen::VectorXd::Map(m_dStandardLightData, endPos - startPos/*y_vec.size()*/);
|
||||
|
||||
using namespace ZZ_MATH::SplineFit;
|
||||
SplineInterpolation m_sfLine(vx, vy);
|
||||
|
||||
m_dStandardLightWavelengthResampled = new double[deviceAttribute.iPixels];
|
||||
m_dStandardLightDataResampled = new double[deviceAttribute.iPixels];
|
||||
|
||||
double dTemp,scaleFactor;
|
||||
if(integratingSphereDetectorValue <= 0)
|
||||
scaleFactor=1;
|
||||
else
|
||||
scaleFactor=integratingSphereDetectorValue/m_dStandardLightDataBase;
|
||||
|
||||
for (size_t i = 0; i < deviceAttribute.iPixels; i++)
|
||||
{
|
||||
if (deviceAttribute.fWaveLengthInNM[i] < StandardLightWavelength_tmp[0])//<2F>˲<EFBFBD><CBB2><EFBFBD> < <20><><EFBFBD><D7BC><EFBFBD>ļ<EFBFBD><C4BC><EFBFBD>С<EFBFBD><D0A1><EFBFBD><EFBFBD>
|
||||
{
|
||||
dTemp = m_sfLine(StandardLightWavelength_tmp[0]);//???????????????
|
||||
}
|
||||
else if (StandardLightWavelength_tmp[lineCount - 2] < deviceAttribute.fWaveLengthInNM[i])//<2F><><EFBFBD><D7BC><EFBFBD>ļ<EFBFBD><C4BC><EFBFBD><EFBFBD><EFBFBD> < <20>˲<EFBFBD><CBB2><EFBFBD>
|
||||
{
|
||||
dTemp = m_sfLine(StandardLightWavelength_tmp[lineCount - 2]);//???????????????
|
||||
}
|
||||
else
|
||||
{
|
||||
dTemp = m_sfLine(deviceAttribute.fWaveLengthInNM[i]);//?????
|
||||
}
|
||||
|
||||
//double dTemp2 = m_sfLine(deviceAttribute.fWaveLengthInNM[i]);//<2F><><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>û<EFBFBD><C3BB>ʹ<EFBFBD><CAB9>
|
||||
|
||||
m_dStandardLightWavelengthResampled[i] = deviceAttribute.fWaveLengthInNM[i];
|
||||
m_dStandardLightDataResampled[i] = dTemp * scaleFactor;
|
||||
}
|
||||
|
||||
|
||||
// //<2F><><EFBFBD>ز<EFBFBD><D8B2><EFBFBD><EFBFBD>ı<EFBFBD><EFBFBD><D7BC><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>csv
|
||||
// QString outputName = standardLightFileFolder + "/" + currentTime + "_" + QString::fromStdString(deviceInfo.strSN) + "_" + standardLightFileName + "_resample.csv";
|
||||
// std::ofstream outfile2(outputName.toStdString().c_str());
|
||||
// for (size_t i = 0; i < deviceAttribute.iPixels; i++)
|
||||
// {
|
||||
// if (i == 0)
|
||||
// {
|
||||
// outfile2 << m_dStandardLightDataBase << std::endl;
|
||||
// }
|
||||
//
|
||||
// outfile2 << deviceAttribute.fWaveLengthInNM[i] << "," << m_dStandardLightDataResampled[i] << std::endl;
|
||||
// }
|
||||
// outfile2.close();
|
||||
}
|
||||
|
||||
void CalibrationAlgorithm::produceCalfile(QString calFilePath, DeviceAttribute deviceAttribute, DataFrame integratingSphereData, DataFrame darkData)
|
||||
{
|
||||
using namespace ZZ_MISCDEF;//ZZ_U32
|
||||
|
||||
int errorCode;
|
||||
size_t writeCounter;
|
||||
|
||||
double* m_gain = new double[deviceAttribute.iPixels];//double*
|
||||
double* m_offset = new double[deviceAttribute.iPixels];//double*
|
||||
|
||||
for (size_t i = 0; i < deviceAttribute.iPixels; i++)
|
||||
{
|
||||
if (integratingSphereData.lData[i] - darkData.lData[i] == 0)//<2F><><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ĸΪ<C4B8><CEAA>
|
||||
{
|
||||
m_gain[i] = 0;
|
||||
}
|
||||
else
|
||||
{
|
||||
m_gain[i] = m_dStandardLightDataResampled[i] / (integratingSphereData.lData[i] - darkData.lData[i]);
|
||||
}
|
||||
|
||||
m_offset[i] = 0;
|
||||
}
|
||||
|
||||
//д<>뵽<EFBFBD><EBB5BD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ļ<EFBFBD>
|
||||
FILE * calFileHandle = fopen(calFilePath.toStdString().c_str(), "wb");
|
||||
|
||||
writeCounter = fwrite(&integratingSphereData.usExposureTimeInMS,sizeof(ZZ_U32), 1, calFileHandle);//<2F>ع<EFBFBD>ʱ<EFBFBD><CAB1>
|
||||
writeCounter = fwrite(&integratingSphereData.fTemperature, sizeof(float), 1, calFileHandle);//<2F>¶<EFBFBD>
|
||||
writeCounter = fwrite(&deviceAttribute.iPixels, sizeof(int), 1, calFileHandle);//<2F><><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
||||
writeCounter = fwrite(&deviceAttribute.fWaveLengthInNM, sizeof(float), deviceAttribute.iPixels, calFileHandle);//<2F><><EFBFBD><EFBFBD>
|
||||
writeCounter = fwrite(m_gain, sizeof(double), deviceAttribute.iPixels, calFileHandle);//gain
|
||||
writeCounter = fwrite(m_offset, sizeof(double), deviceAttribute.iPixels, calFileHandle);//offset
|
||||
|
||||
fclose(calFileHandle);
|
||||
|
||||
|
||||
//д<>뵽CSV<53>ļ<EFBFBD>
|
||||
QString calFile_csv = calFilePath.split(".")[0] + ".csv";
|
||||
std::ofstream outfile(calFile_csv.toStdString().c_str());
|
||||
for (int i = 0; i < deviceAttribute.iPixels; i++)
|
||||
{
|
||||
if (i==0)
|
||||
{
|
||||
outfile << integratingSphereData.usExposureTimeInMS << std::endl;
|
||||
}
|
||||
|
||||
outfile << deviceAttribute.fWaveLengthInNM[i] << "," << m_gain[i] << std::endl;
|
||||
}
|
||||
outfile.close();
|
||||
|
||||
delete[] m_gain;
|
||||
}
|
||||
Reference in New Issue
Block a user