void Printer::printPos(std::string posFileName) {
std::ofstream posFile(posFileName.c_str(), std::ios::trunc);
if (posFile.is_open()) {
posFile << "NetSize,Strain,YoungMod,pbond,Spr1,Spr2,Spr3" << std::endl;
posFile << netSize << "," << affdel << "," << YOUNGMOD << ","
<< p << ",1,1,1" << std::endl;
for (int i = 0; i < netSize; i++) {
for (int j = 0; j < netSize; j++) {
// Print row, col, position, sprstiff, rlen information to file.
posFile << i << "," << j << "," << pos[(i * netSize + j) * 2] - affposx(i, j)
<< "," << pos[(i * netSize + j) * 2 + 1] - affposy(i);
for (int k = 0; k < 3; k++)
posFile << "," << spr[i][j][k];
posFile << "\n";
}
}
}
posFile.close();
}
void QFESPIMB040CamParamStackConfigWidget2::performStack()
{
if (!(this->use1() || this->use2())) {
QMessageBox::critical(this, tr("B040SPIM: Parameter Stack Acquisition"), tr("Cannot start image acquisition: No camera selected!"));
return;
}
QDateTime startDateTime=QDateTime::currentDateTime();
QList<QFESPIMB040OpticsSetupBase::measuredValues> measured;
bool ok=true;
//////////////////////////////////////////////////////////////////////////////////////
// get array of values
//////////////////////////////////////////////////////////////////////////////////////
QList<double> scanVals=this->stack();
if (scanVals.isEmpty()) {
QMessageBox::critical(this, tr("B040SPIM: Parameter Stack Acquisition"), tr("No values to scan selected!"));
return;
}
int images=scanVals.size()*this->images();
log->log_text(tr("starting parameter stack acquisition:\n"));
log->log_text(tr(" - scan parameter: %1\n").arg(this->stackParameterName()));
log->log_text(tr(" - scan mode: %1\n").arg(this->stackModeName()));
log->log_text(tr(" - scan start: %1\n").arg(this->stackStart()));
log->log_text(tr(" - scan end: %1\n").arg(this->stackEnd()));
log->log_text(tr(" - scan delta: %1\n").arg(this->stackDelta()));
opticsSetup->lockLightpath();
//////////////////////////////////////////////////////////////////////////////////////
// LOCK/INIT CAMERA 1
//////////////////////////////////////////////////////////////////////////////////////
bool useCam1=false;
QFExtension* extension1=NULL;
QFExtensionCamera* ecamera1=NULL;
int camera1=0;
QString acquisitionSettingsFilename1="", previewSettingsFilename1="";
QString acquisitionPrefix1=this->prefix1();
QString TIFFFIlename1=acquisitionPrefix1+".tif";
TIFF* tiff1=NULL;
if (this->use1()) {
if (!(useCam1=opticsSetup->lockCamera(0, &extension1, &ecamera1, &camera1, &previewSettingsFilename1))) {
CAMPARAMSTACK_ERROR(tr("error locking camer 1!\n"));
}
}
if (QFile::exists(this->currentConfigFilename(0))) acquisitionSettingsFilename1=this->currentConfigFilename(0);
//////////////////////////////////////////////////////////////////////////////////////
// LOCK/INIT CAMERA 2
//////////////////////////////////////////////////////////////////////////////////////
bool useCam2=false;
QFExtension* extension2=NULL;
QFExtensionCamera* ecamera2=NULL;
QString acquisitionSettingsFilename2="", previewSettingsFilename2="";
QString acquisitionPrefix2=this->prefix2();
QString TIFFFIlename2=acquisitionPrefix2+".tif";
TIFF* tiff2=NULL;
int camera2=0;
if (this->use2()) {
if(!(useCam2=opticsSetup->lockCamera(1, &extension2, &ecamera2, &camera2, &previewSettingsFilename2))) {
CAMPARAMSTACK_ERROR(tr("error locking camer 2!\n"));
}
}
if (QFile::exists(this->currentConfigFilename(1))) acquisitionSettingsFilename2=this->currentConfigFilename(1);
if (ok && !useCam1 && !useCam2) {
CAMPARAMSTACK_ERROR(tr("Cannot start image acquisition: No camera selected, or both cameras not usable!"));
opticsSetup->unlockLightpath();
ok=false;
}
//////////////////////////////////////////////////////////////////////////////////////
// CHECK WHETHER HARDWRAE SUPPORTS CHANGING THE SELECTED PARAMETER
//////////////////////////////////////////////////////////////////////////////////////
int stackParamI=this->stackParameter();
QFExtensionCamera::CameraSetting stackParam=QFExtensionCamera::CameraSetting(stackParamI);
if (ok) {
if (stackParam>=0 && stackParam<=QFExtensionCamera::CamSetMaxParam) {
if (useCam1 && !ecamera1->isCameraSettingChangable(stackParam)) {
ok=false;
CAMPARAMSTACK_ERROR(tr("camera 1 does not support changing the %1!\n").arg(this->stackParameterName()));
}
if (ok && useCam2 && !ecamera2->isCameraSettingChangable(stackParam)) {
ok=false;
CAMPARAMSTACK_ERROR(tr("camera 2 does not support changing the %1!\n").arg(this->stackParameterName()));
}
} else {
ok=false;
CAMPARAMSTACK_ERROR(tr("invalid stack parameter %1 '%2'\n").arg(stackParam).arg(this->stackParameterName()));
}
}
if (ok) {
if (ok && useCam1) log->log_text(tr(" - storing files with prefix 1: '%1'\n").arg(acquisitionPrefix1));
if (ok && useCam2) log->log_text(tr(" - storing files with prefix 2: '%1'\n").arg(acquisitionPrefix2));
QProgressListDialog progress(tr("Parameter Stack Acquisition"), tr("&Cancel"), 0, 100, this);
progress.addItem(tr("preparation"));
progress.addItem(tr("acquisition"));
progress.addItem(tr("cleanup"));
progress.setWindowModality(Qt::WindowModal);
progress.setValue(0);
progress.start();
//////////////////////////////////////////////////////////////////////////////////////
// SET LIGHTPATH
//////////////////////////////////////////////////////////////////////////////////////
opticsSetup->lockLighpathCombobox();
QString oldLightpath=opticsSetup->getCurrentLightpathFilename();
QString oldLightpathName=opticsSetup->getCurrentLightpath();
QString lightpathName="unknown";
if (this->lightpathActivated()) {
if (!QFile::exists(this->lightpathFilename())) {
CAMPARAMSTACK_ERROR(tr(" - acquisition lighpath configuration '%1' does not exist!\n").arg(this->lightpath()));
opticsSetup->unlockLightpath();
return;
} else {
log->log_text(tr(" - setting acquisition lightpath settings '%1' ...\n").arg(this->lightpath()));
opticsSetup->loadLightpathConfig(this->lightpathFilename(), true, &lightpathName);
log->log_text(tr(" - setting acquisition lightpath settings '%1' ... DONE\n").arg(this->lightpath()));
}
}
//////////////////////////////////////////////////////////////////////////////////////
// switch on light
//////////////////////////////////////////////////////////////////////////////////////
bool formerMainShutterState=opticsSetup->getMainIlluminationShutter();
if (opticsSetup->isMainIlluminationShutterAvailable()){
log->log_text(tr(" - switch main shutter on!\n"));
opticsSetup->setMainIlluminationShutter(true, true);
}
QList<QVariant> realValues1, realValues2;
QTime timAcquisition=QTime::currentTime();
QDateTime timStart;
double duration=0;
progress.nextItem();
if (this->previewMode()) {
log->log_text(tr("acquiring stack in PREVIEW MODE ...\n"));
//////////////////////////////////////////////////////////////////////////////////////
// OPEN OUTPUT TIFF FILES
//////////////////////////////////////////////////////////////////////////////////////
progress.setLabelText(tr("opening output files ..."));
QApplication::processEvents();
if (ok && useCam1) {
QDir().mkpath(QFileInfo(TIFFFIlename1.toLatin1().data()).absolutePath());
tiff1=TIFFOpen(TIFFFIlename1.toLatin1().data(), "w");
if (!tiff1) {
ok=false;
CAMPARAMSTACK_ERROR(tr("error opening TIFF file (camera 1) '%1'!").arg(TIFFFIlename1));
}
}
if (ok && useCam2) {
QDir().mkpath(QFileInfo(TIFFFIlename2.toLatin1().data()).absolutePath());
tiff2=TIFFOpen(TIFFFIlename2.toLatin1().data(), "w");
if (!tiff2) {
ok=false;
CAMPARAMSTACK_ERROR(tr("error opening TIFF file (camera 2) '%1'!").arg(TIFFFIlename2));
}
}
if (progress.wasCanceled()) {
log->log_warning(tr("canceled by user!\n"));
ok=false;
}
//////////////////////////////////////////////////////////////////////////////////////
// CHANGE PARAM, ACQUIRE IMAGE, CHANGE PARAM, ACQUIRE IMAGE, ...
// images are stored in TIFF files using libtiff and they are (possibly) downscaled to 16-bit
//////////////////////////////////////////////////////////////////////////////////////
timAcquisition=QTime::currentTime();
duration=0;
if (ok) {
progress.setLabelText(tr("acquiring images ..."));
bool running=ok;
//double newPos=stageStart;
int stackIdx=0;
int width1=0, height1=0;
uint32_t* buffer1=NULL;
int width2=0, height2=0;
uint32_t* buffer2=NULL;
int imageCnt=0;
measured.append(opticsSetup->getMeasuredValues());
while (running && (stackIdx<scanVals.size())) {
if (ok && useCam1) {
//////////////////////////////////////////////////////////////////////////////////////
// PREPARE CAMERA 1
//////////////////////////////////////////////////////////////////////////////////////
if (ok && useCam1) {
progress.setLabelText(tr("preparing camera 1 ..."));
QApplication::processEvents();
QString tmpName1=QDir::temp().absoluteFilePath("qf3spimb040_cam1tmpsettings.ini");
QFTemporaryFile file1;
if (file1.open()) {
tmpName1=file1.fileName();
}
if (QFile::exists(tmpName1)) QFile::remove(tmpName1);
if (useCam1) QFile::copy(acquisitionSettingsFilename1, tmpName1);
QSettings settings(tmpName1, QSettings::IniFormat);
//QSettings settings(acquisitionSettingsFilename1, QSettings::IniFormat);
ecamera1->changeCameraSetting(settings, stackParam, scanVals[stackIdx]);
ecamera1->useCameraSettings(camera1, settings);
realValues1.append(ecamera1->getCameraCurrentSetting(camera1, stackParam));
//realValues1.append(ecamera1->getExposureTime(camera1));
log->log_text(tr(" - prepared camer 1!\n"));
width1=ecamera1->getCameraImageWidth(camera1);
height1=ecamera1->getCameraImageHeight(camera1);
buffer1=(uint32_t*)qfRealloc(buffer1, width1*height1*sizeof(uint32_t));
if (!buffer1) {
ok=false;
CAMPARAMSTACK_ERROR(tr("could not allocate image buffer for camera 1!\n"));
}
if (QFile::exists(tmpName1)) QFile::remove(tmpName1);
}
}
if (ok && useCam2) {
//////////////////////////////////////////////////////////////////////////////////////
// PREPARE CAMERA 1
//////////////////////////////////////////////////////////////////////////////////////
if (ok && useCam2) {
progress.setLabelText(tr("preparing camera 2 ..."));
QApplication::processEvents();
QString tmpName1=QDir::temp().absoluteFilePath("qf3spimb040_cam2tmpsettings.ini");
QFTemporaryFile file1;
if (file1.open()) {
tmpName1=file1.fileName();
}
if (QFile::exists(tmpName1)) QFile::remove(tmpName1);
if (useCam2) QFile::copy(acquisitionSettingsFilename2, tmpName1);
QSettings settings(tmpName1, QSettings::IniFormat);
ecamera2->changeCameraSetting(settings, stackParam, scanVals[stackIdx]);
ecamera2->useCameraSettings(camera2, settings);
realValues2.append(ecamera2->getCameraCurrentSetting(camera2, stackParam));
//realValues2.append(ecamera2->getExposureTime(camera2));
log->log_text(tr(" - prepared camer 2!\n"));
width2=ecamera2->getCameraImageWidth(camera2);
height2=ecamera2->getCameraImageHeight(camera2);
buffer2=(uint32_t*)qfRealloc(buffer2, width2*height2*sizeof(uint32_t));
if (!buffer2) {
ok=false;
CAMPARAMSTACK_ERROR(tr("could not allocate image buffer for camera 2!\n"));
}
if (QFile::exists(tmpName1)) QFile::remove(tmpName1);
}
}
QApplication::processEvents();
if (progress.wasCanceled()) {
running=false;
log->log_warning(tr(" - acquisition canceled by user!\n"));
} else {
for (int img=0; img<this->images(); img++) {
log->log_text(tr("acquiring images (%1/%2) ...\n").arg(imageCnt+1).arg(images));
progress.setLabelText(tr("acquiring images (%1/%2) ...").arg(imageCnt+1).arg(images));
QApplication::processEvents();
if (progress.wasCanceled()) {
running=false;
log->log_warning(tr(" - acquisition canceled by user!\n"));
break;
}
if (stackIdx==0) {
timAcquisition.start();
timStart=QDateTime::currentDateTime();
}
if (useCam1) {
if (ecamera1->acquireOnCamera(camera1, buffer1)) {
TIFFTWriteUint16from32(tiff1, buffer1, width1, height1, false);
TIFFWriteDirectory(tiff1);
} else {
ok=false;
CAMPARAMSTACK_ERROR(tr("error acquiring image %1/%2 on camera 1!\n").arg(imageCnt+1).arg(images));
}
}
//QApplication::processEvents();
if (useCam2) {
if (ecamera2->acquireOnCamera(camera2, buffer2)) {
TIFFTWriteUint16from32(tiff2, buffer2, width2, height2, false);
TIFFWriteDirectory(tiff2);
} else {
ok=false;
CAMPARAMSTACK_ERROR(tr("error acquiring image %1/%2 on camera 2!\n").arg(imageCnt+1).arg(images));
}
}
imageCnt++;
}
if (this->saveMeasurements()) measured.append(opticsSetup->getMeasuredValues());
//QApplication::processEvents();
}
if (!ok) running=false;
stackIdx++;
progress.setValue((int)round((double)stackIdx/(double)scanVals.size()*100.0));
QApplication::processEvents();
}
duration=timAcquisition.elapsed()/1000.0;
if (buffer1) qfFree(buffer1);
if (buffer2) qfFree(buffer2);
buffer1=buffer2=NULL;
}
} else {
log->log_text(tr("acquiring stack in ACQUISITION MODE ...\n"));
if (progress.wasCanceled()) {
log->log_warning(tr("canceled by user!\n"));
ok=false;
}
QMap<QString, QVariant> acquisitionDescription1;
QList<QFExtensionCamera::CameraAcquititonFileDescription> moreFiles1;
QMap<QString, QVariant> acquisitionDescription2;
QList<QFExtensionCamera::CameraAcquititonFileDescription> moreFiles2;
//////////////////////////////////////////////////////////////////////////////////////
// CHANGE PARAM, ACQUIRE IMAGE, CHANGE PARAM, ACQUIRE IMAGE, ...
// images are stored in TIFF files using libtiff and they are (possibly) downscaled to 16-bit
//////////////////////////////////////////////////////////////////////////////////////
timAcquisition=QTime::currentTime();
duration=0;
if (ok) {
progress.setLabelText(tr("acquiring images ..."));
bool running=ok;
//double newPos=stageStart;
int stackIdx=0;
int width1=0, height1=0;
uint32_t* buffer1=NULL;
int width2=0, height2=0;
uint32_t* buffer2=NULL;
int imageCnt=0;
int frames=this->images();
measured.append(opticsSetup->getMeasuredValues());
while (running && (stackIdx<scanVals.size())) {
QString tmpName1=QDir::temp().absoluteFilePath("qf3spimb040_cam1tmpsettings.ini");
QFTemporaryFile file1;
if (file1.open()) {
tmpName1=file1.fileName();
}
if (QFile::exists(tmpName1)) QFile::remove(tmpName1);
if (useCam1) QFile::copy(acquisitionSettingsFilename1, tmpName1);
QSettings settings1(tmpName1, QSettings::IniFormat);
if (frames>0 && ecamera1->isCameraSettingChangable(QFExtensionCamera::CamSetNumberFrames)) ecamera1->changeCameraSetting(settings1, QFExtensionCamera::CamSetNumberFrames, frames);
ecamera1->changeCameraSetting(settings1, stackParam, scanVals[stackIdx]);
QString tmpName2=QDir::temp().absoluteFilePath("qf3spimb040_cam2tmpsettings.ini");
QFTemporaryFile file2;
if (file2.open()) {
tmpName2=file2.fileName();
}
if (QFile::exists(tmpName2)) QFile::remove(tmpName2);
if (useCam2) QFile::copy(acquisitionSettingsFilename2, tmpName2);
QSettings settings2(tmpName2, QSettings::IniFormat);
if (frames>0 && ecamera2->isCameraSettingChangable(QFExtensionCamera::CamSetNumberFrames)) ecamera2->changeCameraSetting(settings2, QFExtensionCamera::CamSetNumberFrames, frames);
ecamera2->changeCameraSetting(settings2, stackParam, scanVals[stackIdx]);
//////////////////////////////////////////////////////////////////////////////////////
// acquire image series
//////////////////////////////////////////////////////////////////////////////////////
ok = acqTools->acquireSeries(lightpathName, QString("camstack%1").arg(stackIdx,4,10,QLatin1Char('0')), tr("cam parameter stack"), useCam1, extension1, ecamera1, camera1, acquisitionPrefix1+QString("__sidx%1").arg(stackIdx,4,10,QLatin1Char('0')), settings1, acquisitionDescription1, moreFiles1, useCam2, extension2, ecamera2, camera2, acquisitionPrefix2+QString("__sidx%1").arg(stackIdx,4,10,QLatin1Char('0')), settings2, acquisitionDescription2, moreFiles2, &measured, &progress, NULL);
if (!ok) {
if (frames>1) {
CAMPARAMSTACK_ERROR(tr("error acquiring image (%1...%2)/%3 !\n").arg(imageCnt+1).arg(imageCnt+1+frames).arg(images));
} else {
CAMPARAMSTACK_ERROR(tr("error acquiring image %1/%2 !\n").arg(imageCnt+1).arg(images));
}
}
if (QFile::exists(tmpName1)) QFile::remove(tmpName1);
if (QFile::exists(tmpName2)) QFile::remove(tmpName2);
if (!ok) running=false;
stackIdx++;
imageCnt+=frames;
progress.setValue((int)round((double)stackIdx/(double)scanVals.size()*100.0));
QApplication::processEvents();
if (progress.wasCanceled()) {
running=false;
log->log_error(tr("CANCELED BY USER!\n"));
}
}
duration=timAcquisition.elapsed()/1000.0;
if (buffer1) qfFree(buffer1);
if (buffer2) qfFree(buffer2);
buffer1=buffer2=NULL;
}
}
progress.setValue(100);
measured.append(opticsSetup->getMeasuredValues());
progress.nextItem();
//////////////////////////////////////////////////////////////////////////////////////
// collect lightpath data common to all cameras
//////////////////////////////////////////////////////////////////////////////////////
QMap<QString, QVariant> acquisitionDescription;
opticsSetup->saveLightpathConfig(acquisitionDescription, "", "lightpath/");
//////////////////////////////////////////////////////////////////////////////////////
// switch on/off light
//////////////////////////////////////////////////////////////////////////////////////
if (opticsSetup->isMainIlluminationShutterAvailable()){
log->log_text(tr(" - switch main shutter %1!\n").arg((formerMainShutterState)?tr("on"):tr("off")));
opticsSetup->setMainIlluminationShutter(formerMainShutterState, true);
}
//////////////////////////////////////////////////////////////////////////////////////
// RESET LIGHTPATH
//////////////////////////////////////////////////////////////////////////////////////
if (this->lightpathActivated()) {
opticsSetup->loadLightpathConfig(oldLightpath, false);
log->log_text(tr(" - resetting to old lightpath settings (%1) ...\n").arg(oldLightpath));//Name));
}
opticsSetup->unlockLighpathCombobox();
//////////////////////////////////////////////////////////////////////////////////////
// close tiff files and free buffers
//////////////////////////////////////////////////////////////////////////////////////
progress.setLabelText(tr("closing output files ..."));
QApplication::processEvents();
if (tiff1) TIFFClose(tiff1);
if (tiff2) TIFFClose(tiff2);
tiff1=tiff2=NULL;
//////////////////////////////////////////////////////////////////////////////////////
// collect acquisition data common to all cameras
//////////////////////////////////////////////////////////////////////////////////////
QString scanCSV;
if (ok) {
acquisitionDescription["type"]="parameter stack";
acquisitionDescription["stack_parameter"]=this->stackParameterName();
acquisitionDescription["stack_mode"]=this->stackModeName();
acquisitionDescription["stack_start"]=this->stackStart();
acquisitionDescription["stack_delta"]=this->stackDelta();
acquisitionDescription["stack_end"]=this->stackEnd();
acquisitionDescription["sequence_overall_length"]=images;
acquisitionDescription["images_per_step"]=this->images();
acquisitionDescription["stack_length"]=scanVals.size();
acquisitionDescription["start_time"]=timStart;
acquisitionDescription["duration"]=duration;
acquisitionDescription["stack_value_count"]=scanVals.size();
acquisitionDescription["stack_values"]=CDoubleListToQString(scanVals);
if (useCam1) acquisitionDescription["real_stack_values/camera1"]=realValues1;
if (useCam2) acquisitionDescription["real_stack_values/camera2"]=realValues2;
QTextStream pf(&scanCSV);
pf<<"# number, ideal value";
if (useCam1) pf<<", real value camera1";
if (useCam2) pf<<", real value camera2";
pf<<"\n";
for (int i=0; i<qMin(scanVals.size(), qMax(realValues1.size(), realValues2.size())); i++) {
pf<<i<<", "<<CDoubleToQString(scanVals[i]);
if (useCam1 && i<realValues1.size()) pf<<", "<<CDoubleToQString(realValues1[i].toDouble());
if (useCam2 && i<realValues2.size()) pf<<", "<<CDoubleToQString(realValues2[i].toDouble());
pf<<"\n";
}
}
//////////////////////////////////////////////////////////////////////////////////////
// write image stack properties to files, also collects camera specific information
//////////////////////////////////////////////////////////////////////////////////////
if (ok && useCam1) {
QMap<QString, QVariant> acquisitionDescription1=acquisitionDescription;
QList<QFExtensionCamera::CameraAcquititonFileDescription> files;
QFExtensionCamera::CameraAcquititonFileDescription d;
d.name=TIFFFIlename1;
d.description="image stack from camera 1";
d.type="TIFF16";
files.append(d);
QString ParamValuesFilename=acquisitionPrefix1+".param.dat";
QFile posFile(ParamValuesFilename);
if (posFile.open(QIODevice::WriteOnly)) {
posFile.write(scanCSV.toLatin1().data());
posFile.close();
d.name=ParamValuesFilename;
d.description="parameter values";
d.type="CSV";
files.append(d);
} else {
log->log_error(tr(" - could not write parameter values file '%1' for camera 1: %2 ...").arg(ParamValuesFilename).arg(posFile.errorString()));
}
QString MeasurementsFilename=acqTools->saveMeasuredData(acquisitionPrefix1, measured);
if (!MeasurementsFilename.isEmpty() && QFile::exists(MeasurementsFilename)) {
QFExtensionCamera::CameraAcquititonFileDescription d;
d.name=MeasurementsFilename;
d.description="measureable properties of setup";
d.type="CSV";
files.append(d);
}
log->log_text(tr(" - writing acquisition description 1 ..."));
acqTools->savePreviewDescription(0, extension1, ecamera1, camera1, acquisitionPrefix1, acquisitionDescription1, files, startDateTime);
log->log_text(tr(" DONE!\n"));
}
if (ok && useCam2) {
QMap<QString, QVariant> acquisitionDescription2=acquisitionDescription;
QList<QFExtensionCamera::CameraAcquititonFileDescription> files;
QFExtensionCamera::CameraAcquititonFileDescription d;
d.name=TIFFFIlename1;
d.description="image stack from camera 2";
d.type="TIFF16";
files.append(d);
QString ParamValuesFilename=acquisitionPrefix2+".param.dat";
QFile posFile(ParamValuesFilename);
if (posFile.open(QIODevice::WriteOnly)) {
posFile.write(scanCSV.toLatin1().data());
posFile.close();
d.name=ParamValuesFilename;
d.description="parameter values";
d.type="CSV";
files.append(d);
} else {
log->log_error(tr(" - could not write parameter values file '%1' for camera 2: %2 ...").arg(ParamValuesFilename).arg(posFile.errorString()));
}
QString MeasurementsFilename=acqTools->saveMeasuredData(acquisitionPrefix2, measured);
if (!MeasurementsFilename.isEmpty() && QFile::exists(MeasurementsFilename)) {
QFExtensionCamera::CameraAcquititonFileDescription d;
d.name=MeasurementsFilename;
d.description="measureable properties of setup";
d.type="CSV";
files.append(d);
}
log->log_text(tr(" - writing acquisition description 2 ..."));
acqTools->savePreviewDescription(1, extension2, ecamera2, camera2, acquisitionPrefix2, acquisitionDescription2, files, startDateTime);
log->log_text(tr(" DONE!\n"));
}
//////////////////////////////////////////////////////////////////////////////////////
// release cameras
//////////////////////////////////////////////////////////////////////////////////////
if (useCam1) {
opticsSetup->releaseCamera(0);
log->log_text(tr(" - released camera 1!\n"));
}
if (useCam2) {
opticsSetup->releaseCamera(1);
log->log_text(tr(" - released camera 2!\n"));
}
if (ok && useCam1) log->log_text(tr(" - stored files with prefix 1: '%1'\n").arg(acquisitionPrefix1));
if (ok && useCam2) log->log_text(tr(" - stored files with prefix 2: '%1'\n").arg(acquisitionPrefix2));
if (ok) log->log_text(tr("image stack acquisition DONE!\n"));
}
opticsSetup->unlockLightpath();
opticsSetup->ensureLightpath();
}
void StandardStrings::convertLasToRms(QWidget *parent, const QString& lasFileName, const QString& rmsFileName, const QString& posFileName){
//this assumes that MD/Z is the first log in the las file.
QFile lasFile(lasFileName);
QFile rmsFile(rmsFileName);
QFile posFile(posFileName);
if (!lasFile.open(QFile::ReadOnly | QFile::Text)) {
QMessageBox::warning(parent, QString("Unable to open file"), QString("Could not open the file."), QMessageBox::Ok);
return;
}
if (!posFile.open(QFile::ReadOnly | QFile::Text)) {
QMessageBox::warning(parent, QString("Unable to open file"), QString("Could not open the file."), QMessageBox::Ok);
return;
}
if (!rmsFile.open(QFile::WriteOnly | QFile::Text)) {
QMessageBox::warning(parent, QString("Unable to open file"), QString("Could not open the file."), QMessageBox::Ok);
return;
}
QString wellName;
QTextStream lasStream(&lasFile);
QTextStream rmsStream(&rmsFile);
QTextStream posStream(&posFile);
QString line = posStream.readLine();
while(!line.isNull()){
if(!line.isEmpty()){
if(line[0]==QChar('-')){
line = posStream.readLine();
break;
}
}
line = posStream.readLine();
}
QList<double> xList;
QList<double> yList;
QList<double> mdList;
//
while(!line.isNull()){
QStringList list =line.simplified().split(" ");
xList << list[0].toDouble();
yList << list[1].toDouble();
mdList << list[5].toDouble();
line = posStream.readLine();
}
line = lasStream.readLine();
while(!line.isNull()){
if(line.split(" ").at(0)==QString("WELL")){
wellName=line.split("\"").at(1);
while(!wellName[0].isDigit()){
wellName.remove(0,1);
if(wellName.isEmpty()){
QMessageBox::warning(parent, QString("No well name"), QString("Could not find a valid well name."), QMessageBox::Ok);
return;
}
}
wellName=wellName.remove(" ");
break;
}
line = lasStream.readLine();
}
if(wellName.isEmpty()){
QMessageBox::warning(parent, QString("No well name"), QString("Could not find a valid well name."), QMessageBox::Ok);
return;
}
QString wellX;
while(!line.isNull()){
if(line.split(" ").at(0)==QString("XWELL")){
QStringList wellXList=line.split(" ");
for(int i=0;i<wellXList.size();++i){
if(!wellXList[i].isEmpty()){
bool ok;
wellXList[i].toDouble(&ok);
if(ok){
wellX=wellXList[i];
}
}
}
break;
}
line = lasStream.readLine();
}
if(wellX.isEmpty()){
QMessageBox::warning(parent, QString("No well x-coordinate"), QString("Could not find a valid well x-coordinate."), QMessageBox::Ok);
return;
}
QString wellY;
while(!line.isNull()){
if(line.split(" ").at(0)==QString("YWELL")){
QStringList wellYList=line.split(" ");
for(int i=0;i<wellYList.size();++i){
if(!wellYList[i].isEmpty()){
bool ok;
wellYList[i].toDouble(&ok);
if(ok){
wellY=wellYList[i];
}
}
}
break;
}
line = lasStream.readLine();
}
if(wellY.isEmpty()){
QMessageBox::warning(parent, QString("No well x-coordinate"), QString("Could not find a valid well x-coordinate."), QMessageBox::Ok);
return;
}
while(!line.isNull()){
if(!line.isEmpty()){
if(line[0]==QChar('~')){
lasStream.readLine();//read away the comments
lasStream.readLine();
break;
}
}
line = lasStream.readLine();
}
QStringList logNames;
lasStream.readLine();
line = lasStream.readLine();
while(!line.isEmpty()){
logNames << line.split(".")[0].remove(" ");
line = lasStream.readLine();
}
lasStream.readLine();//~ASCII
line = lasStream.readLine();
int minIndex=0;
rmsStream.QTextStream::setRealNumberNotation(QTextStream::FixedNotation);
if(wellX.size()>wellY.size()){
rmsStream.setRealNumberPrecision(wellX.size());
}
else{
rmsStream.setRealNumberPrecision(wellY.size());
}
rmsStream.setRealNumberPrecision(2);
rmsStream << 1.0 << endl;
rmsStream << "undefined" << endl;
rmsStream << wellName << " " << wellX << " " << wellY << endl;
rmsStream << logNames.size() << endl;
for(int i=0;i<logNames.size();++i){
rmsStream << logNames[i] << " " << "UNK lin" << endl;
}
bool first=true;
while(!line.isEmpty()){
if(!first){
rmsStream << endl;
}
else{
first=false;
}
bool printed=false;
QStringList lineList = line.simplified().split(" ");
for(int j=minIndex;j<mdList.size();++j){
if(mdList[j]==lineList.at(0).toDouble()){
rmsStream << xList[j] << " " << yList[j];
for(int i=0;i<lineList.size();++i){
rmsStream << " " << lineList.at(i);
}
printed=true;
break;
} else if(mdList[j]<lineList.at(0).toDouble()){
minIndex=j;
} else{
bool overRule=false;
while(mdList[j]==mdList[minIndex]){
if(j<mdList.size()-1){
++j;
}
else {
overRule=true;
break;
}
}
double fractionDifference=(lineList.at(0).toDouble()-mdList[minIndex])/(mdList[j]-mdList[minIndex]);
if(overRule){
fractionDifference=0;
}
rmsStream << (xList[j]+fractionDifference*(xList[j]-xList[minIndex])) << " ";
rmsStream << (yList[j]+fractionDifference*(yList[j]-yList[minIndex]));
for(int i=0;i<lineList.size();++i){
rmsStream << " " << lineList.at(i);
}
printed=true;
break;
}
}
if(!printed){
rmsStream << (xList[xList.size()-1]) << " " << (yList[yList.size()-1]);
for(int i=0;i<lineList.size();++i){
rmsStream << " " << lineList.at(i);
}
}
line = lasStream.readLine();
}
}
