Skip to content

Commit

Permalink
Pseudo True Color for GOES-16
Browse files Browse the repository at this point in the history
  • Loading branch information
@hvanruys
hvanruys committed Jan 17, 2018
1 parent 81ab35e commit fc7051c
Show file tree
Hide file tree
Showing 8 changed files with 81 additions and 68 deletions.
9 changes: 5 additions & 4 deletions core/formgeostationary.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -756,7 +756,7 @@ FormGeostationary::~FormGeostationary()
delete ui;
}

void FormGeostationary::CreateGeoImage(QString type, QVector<QString> spectrumvector, QVector<bool> inversevector, int histogrammethod)
void FormGeostationary::CreateGeoImage(QString type, QVector<QString> spectrumvector, QVector<bool> inversevector, int histogrammethod, bool pseudocolor)
{
// segs->seglgeo[0]->areatype == 1 ==> full
// segs->seglgeo[0]->areatype == 0 ==> europe
Expand Down Expand Up @@ -835,7 +835,7 @@ void FormGeostationary::CreateGeoImage(QString type, QVector<QString> spectrumve
if(opts.geosatellites.at(sl->getGeoSatelliteIndex()).protocol == "HDF" )
CreateGeoImageHDF(sl, type, tex, spectrumvector, inversevector);
else if(opts.geosatellites.at(sl->getGeoSatelliteIndex()).protocol == "netCDF")
CreateGeoImagenetCDF(sl, type, tex, spectrumvector, inversevector, histogrammethod);
CreateGeoImagenetCDF(sl, type, tex, spectrumvector, inversevector, histogrammethod, pseudocolor);
else
CreateGeoImageXRIT(sl, type, tex, spectrumvector, inversevector, histogrammethod);

Expand Down Expand Up @@ -1175,7 +1175,7 @@ void FormGeostationary::CreateGeoImageHDF(SegmentListGeostationary *sl, QString

}

void FormGeostationary::CreateGeoImagenetCDF(SegmentListGeostationary *sl, QString type, QString tex, QVector<QString> spectrumvector, QVector<bool> inversevector, int histogrammethod)
void FormGeostationary::CreateGeoImagenetCDF(SegmentListGeostationary *sl, QString type, QString tex, QVector<QString> spectrumvector, QVector<bool> inversevector, int histogrammethod, bool pseudocolor)
{

QString filetiming;
Expand Down Expand Up @@ -1219,7 +1219,8 @@ void FormGeostationary::CreateGeoImagenetCDF(SegmentListGeostationary *sl, QStri
return;
}
else
sl->ComposeImagenetCDFInThread(llVIS_IR, spectrumvector, inversevector, histogrammethod);
sl->setThreadParameters(llVIS_IR, spectrumvector, inversevector, histogrammethod, pseudocolor);
sl->ComposeImagenetCDFInThread(llVIS_IR, spectrumvector, inversevector, histogrammethod, pseudocolor);
}
}

Expand Down
4 changes: 2 additions & 2 deletions core/formgeostationary.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -38,7 +38,7 @@ class FormGeostationary : public QWidget
void PopulateTreeGeo(int geoindex);
void CreateGeoImageXRIT(SegmentListGeostationary *sl, QString type, QString tex, QVector<QString> spectrumvector, QVector<bool> inversevector, int histogrammethod);
void CreateGeoImageHDF(SegmentListGeostationary *sl, QString type, QString tex, QVector<QString> spectrumvector, QVector<bool> inversevector);
void CreateGeoImagenetCDF(SegmentListGeostationary *sl, QString type, QString tex, QVector<QString> spectrumvector, QVector<bool> inversevector, int histogrammethod);
void CreateGeoImagenetCDF(SegmentListGeostationary *sl, QString type, QString tex, QVector<QString> spectrumvector, QVector<bool> inversevector, int histogrammethod, bool pseudocolor);

Ui::FormGeostationary *ui;
AVHRRSatellite *segs;
Expand All @@ -50,7 +50,7 @@ class FormGeostationary : public QWidget

public slots:
void PopulateTree();
void CreateGeoImage(QString type, QVector<QString> spectrumvector, QVector<bool> inversevector, int histogrammethod);
void CreateGeoImage(QString type, QVector<QString> spectrumvector, QVector<bool> inversevector, int histogrammethod, bool pseudocolor);


private slots:
Expand Down
4 changes: 2 additions & 2 deletions core/formtoolbox.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -2058,7 +2058,7 @@ void FormToolbox::onButtonChannel( QString channel, bool bInverse)
//formimage->adjustPicSize(true);
emit switchstackedwidget(3);

emit getgeosatchannel("VIS_IR", spectrumvector, inversevector, ui->cmbHistogramGeo->currentIndex());
emit getgeosatchannel("VIS_IR", spectrumvector, inversevector, ui->cmbHistogramGeo->currentIndex(), ui->rdbPseudoColor->isChecked());
}

void FormToolbox::on_btnGeoColor_clicked()
Expand Down Expand Up @@ -2479,7 +2479,7 @@ void FormToolbox::onButtonColorHRV(QString type)
}

emit switchstackedwidget(3);
emit getgeosatchannel(type, spectrumvector, inversevector, ui->cmbHistogramGeo->currentIndex());
emit getgeosatchannel(type, spectrumvector, inversevector, ui->cmbHistogramGeo->currentIndex(), ui->rdbPseudoColor->isChecked());
}

void FormToolbox::on_btnTextureMet_clicked()
Expand Down
2 changes: 1 addition & 1 deletion core/formtoolbox.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -118,7 +118,7 @@ public slots:
void slotDisplayDNBGraph();

signals:
void getgeosatchannel(QString, QVector<QString>, QVector<bool>, int);
void getgeosatchannel(QString, QVector<QString>, QVector<bool>, int, bool);
void overlaycorrection(int,int);
void switchstackedwidget(int);

Expand Down
9 changes: 8 additions & 1 deletion core/formtoolbox.ui
View file
Original file line number Diff line number Diff line change
Expand Up @@ -32,7 +32,7 @@
<item>
<widget class="QTabWidget" name="tabWidget">
<property name="currentIndex">
<number>1</number>
<number>3</number>
</property>
<property name="tabBarAutoHide">
<bool>false</bool>
Expand Down Expand Up @@ -4942,6 +4942,13 @@
</item>
</layout>
</item>
<item>
<widget class="QRadioButton" name="rdbPseudoColor">
<property name="text">
<string>Calculate pseudo True Color for GOES-16</string>
</property>
</widget>
</item>
<item>
<layout class="QHBoxLayout" name="horizontalLayout_41">
<item>
Expand Down
2 changes: 1 addition & 1 deletion core/mainwindow.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -166,7 +166,7 @@ MainWindow::MainWindow(QWidget *parent) :
connect( formgeostationary, SIGNAL(geostationarysegmentschosen(int, QStringList)), formtoolbox, SLOT(geostationarysegmentsChosen(int, QStringList)));
connect( formgeostationary, SIGNAL(setbuttonlabels(int, bool)), formtoolbox, SLOT(setButtons(int, bool)));

connect( formtoolbox, SIGNAL(getgeosatchannel(QString, QVector<QString>, QVector<bool>, int)), formgeostationary, SLOT(CreateGeoImage(QString, QVector<QString>, QVector<bool>, int)));
connect( formtoolbox, SIGNAL(getgeosatchannel(QString, QVector<QString>, QVector<bool>, int, bool)), formgeostationary, SLOT(CreateGeoImage(QString, QVector<QString>, QVector<bool>, int, bool)));
connect( formtoolbox, SIGNAL(switchstackedwidget(int)), this, SLOT(slotSwitchStackedWindow(int)));

connect( formgeostationary, SIGNAL(enabletoolboxbuttons(bool)), formtoolbox, SLOT(setToolboxButtons(bool)));
Expand Down
95 changes: 48 additions & 47 deletions core/segmentlistgeostationary.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -73,12 +73,12 @@ extern gshhsData *gshhsdata;
// Electro_n1
// Height = 2784 / 6 = 464 Width = 2784

void doComposeGeostationaryXRIT(SegmentListGeostationary *sm, QString segment_path, int channelindex, QVector<QString> spectrumvector, QVector<bool> inversevector)
void SegmentListGeostationary::doComposeGeostationaryXRIT(SegmentListGeostationary *sm, QString segment_path, int channelindex, QVector<QString> spectrumvector, QVector<bool> inversevector)
{
sm->ComposeSegmentImageXRIT(segment_path, channelindex, spectrumvector, inversevector);
}

void doComposeGeostationaryXRITHimawari(SegmentListGeostationary *sm, QString segment_path, int channelindex, QVector<QString> spectrumvector, QVector<bool> inversevector)
void SegmentListGeostationary::doComposeGeostationaryXRITHimawari(SegmentListGeostationary *sm, QString segment_path, int channelindex, QVector<QString> spectrumvector, QVector<bool> inversevector)
{
sm->ComposeSegmentImageXRITHimawari(segment_path, channelindex, spectrumvector, inversevector);
}
Expand All @@ -90,14 +90,14 @@ void doComposeGeostationaryHDF(SegmentListGeostationary *sm, QString segment_pat
}
*/

void doComposeGeostationaryHDFInThread(SegmentListGeostationary *sm, QStringList filelist, QVector<QString> spectrumvector, QVector<bool> inversevector)
void SegmentListGeostationary::doComposeGeostationaryHDFInThread(SegmentListGeostationary *sm, QStringList filelist, QVector<QString> spectrumvector, QVector<bool> inversevector)
{
sm->ComposeSegmentImageHDFInThread(filelist, spectrumvector, inversevector);
}

void doComposeGeostationarynetCDFInThread(SegmentListGeostationary *sm, QStringList strlist, QVector<QString> spectrumvector, QVector<bool> inversevector, int histogrammethod)
void SegmentListGeostationary::doComposeGeostationarynetCDFInThread(SegmentListGeostationary *sm) //, QStringList strlist, QVector<QString> spectrumvector, QVector<bool> inversevector, int histogrammethod)
{
sm->ComposeSegmentImagenetCDFInThread(strlist, spectrumvector, inversevector, histogrammethod);
sm->ComposeSegmentImagenetCDFInThread(); //strlist, spectrumvector, inversevector, histogrammethod);
}


Expand Down Expand Up @@ -345,13 +345,25 @@ bool SegmentListGeostationary::ComposeImageHDFInThread(QStringList strlist, QVec
return true;
}

bool SegmentListGeostationary::ComposeImagenetCDFInThread(QStringList strlist, QVector<QString> spectrumvector, QVector<bool> inversevector, int histogrammethod)
void SegmentListGeostationary::setThreadParameters(QStringList strlist, QVector<QString> spectrumvector, QVector<bool> inversevector, int histogrammethod, bool pseudocolor)
{
this->segmentfilelist = strlist;
this->spectrumvector = spectrumvector;
this->inversevector = inversevector;
this->histogrammethod = histogrammethod;
this->pseudocolor = pseudocolor;

}

bool SegmentListGeostationary::ComposeImagenetCDFInThread(QStringList strlist, QVector<QString> spectrumvector, QVector<bool> inversevector, int histogrammethod, bool pseudocolor)
{

qDebug() << QString("SegmentListGeostationary::ComposeImagenetCDFInThread spectrumvector = %2 %3 %4").arg(spectrumvector.at(0)).arg(spectrumvector.at(1)).arg(spectrumvector.at(2));


QApplication::setOverrideCursor(( Qt::WaitCursor));
QtConcurrent::run(doComposeGeostationarynetCDFInThread, this, strlist, spectrumvector, inversevector, histogrammethod);
setThreadParameters(strlist, spectrumvector, inversevector, histogrammethod, pseudocolor);
QtConcurrent::run(doComposeGeostationarynetCDFInThread, this); //, strlist, spectrumvector, inversevector, histogrammethod);
return true;
}

Expand Down Expand Up @@ -1517,7 +1529,7 @@ void SegmentListGeostationary::ComposeSegmentImageHDFInThread(QStringList fileli

}

void SegmentListGeostationary::ComposeSegmentImagenetCDFInThread(QStringList filelist, QVector<QString> spectrumvector, QVector<bool> inversevector , int histogrammethod)
void SegmentListGeostationary::ComposeSegmentImagenetCDFInThread() //(QStringList filelist, QVector<QString> spectrumvector, QVector<bool> inversevector , int histogrammethod)
{

QString ncfile[3];
Expand Down Expand Up @@ -1546,12 +1558,12 @@ void SegmentListGeostationary::ComposeSegmentImagenetCDFInThread(QStringList fil

QStringList filelistsorted;

for(int i = 0; i < spectrumvector.count(); i++)
for(int i = 0; i < this->spectrumvector.count(); i++)
{
for(int j = 0; j < filelist.count(); j++)
for(int j = 0; j < this->segmentfilelist.count(); j++)
{
if(filelist[j].mid(opts.geosatellites.at(geoindex).indexspectrum, opts.geosatellites.at(geoindex).spectrumlist.at(0).length()) == spectrumvector.at(i))
filelistsorted.append(filelist[j]);
if(this->segmentfilelist[j].mid(opts.geosatellites.at(geoindex).indexspectrum, opts.geosatellites.at(geoindex).spectrumlist.at(0).length()) == this->spectrumvector.at(i))
filelistsorted.append(this->segmentfilelist[j]);
}
}

Expand All @@ -1564,7 +1576,7 @@ void SegmentListGeostationary::ComposeSegmentImagenetCDFInThread(QStringList fil

qDebug() << "Starting netCDF file " + ncfile[j];
retval = nc_open(pncfile[j], NC_NOWRITE, &ncfileid[j]);
if(retval != NC_NOERR) qDebug() << "error opening netCDF file " << filelist.at(j);
if(retval != NC_NOERR) qDebug() << "error opening netCDF file " << this->segmentfilelist.at(j);

retval = nc_inq_varid(ncfileid[j], "max_radiance_value_of_valid_pixels", &varid);
if (retval != NC_NOERR) qDebug() << "error reading max_radiance_value_of_valid_pixels id";
Expand Down Expand Up @@ -1754,20 +1766,23 @@ void SegmentListGeostationary::ComposeSegmentImagenetCDFInThread(QStringList fil
emit this->progressCounter(progcounter);


for(int line = 0; line < ydim; line++)
if(this->pseudocolor)
{
for (int pixelx = 0 ; pixelx < xdim; pixelx++)
for(int line = 0; line < ydim; line++)
{
if(kindofimage == "VIS_IR Color")
for (int pixelx = 0 ; pixelx < xdim; pixelx++)
{
rc = *(imageptrs->ptrRed[0] + line * xdim + pixelx);
gc = *(imageptrs->ptrGreen[0] + line * xdim + pixelx);
bc = *(imageptrs->ptrBlue[0] + line * xdim + pixelx);

if(rc != fillvalue[0] && gc != fillvalue[1] && bc != fillvalue[2] )
if(kindofimage == "VIS_IR Color")
{
true_green = (quint16)qMin(qMax(qRound(0.48358168 * (float)rc + 0.45706946 * (float)bc + 0.06038137 * (float)gc), 0), 1023);
*(imageptrs->ptrGreen[0] + line * xdim + pixelx) = true_green;
rc = *(imageptrs->ptrRed[0] + line * xdim + pixelx);
gc = *(imageptrs->ptrGreen[0] + line * xdim + pixelx);
bc = *(imageptrs->ptrBlue[0] + line * xdim + pixelx);

if(rc != fillvalue[0] && gc != fillvalue[1] && bc != fillvalue[2] )
{
true_green = (quint16)qMin(qMax(qRound(0.48358168 * (float)rc + 0.45706946 * (float)bc + 0.06038137 * (float)gc), 0), 1023);
*(imageptrs->ptrGreen[0] + line * xdim + pixelx) = true_green;
}
}
}
}
Expand All @@ -1790,7 +1805,7 @@ void SegmentListGeostationary::ComposeSegmentImagenetCDFInThread(QStringList fil
dqfvalue[1] = *(imageptrs->ptrDQF[1] + line * xdim + pixelx);
dqfvalue[2] = *(imageptrs->ptrDQF[2] + line * xdim + pixelx);

if(histogrammethod == CMB_HISTO_NONE_95)
if(this->histogrammethod == CMB_HISTO_NONE_95)
{
if(rc != fillvalue[0])
indexoutrc = (quint16)qMin(qMax(qRound(1023.0 * (float)(rc - imageptrs->minRadianceIndex[0] ) / (float)(imageptrs->maxRadianceIndex[0] - imageptrs->minRadianceIndex[0])), 0), 1023);
Expand All @@ -1799,7 +1814,7 @@ void SegmentListGeostationary::ComposeSegmentImagenetCDFInThread(QStringList fil
if(bc != fillvalue[2])
indexoutbc = (quint16)qMin(qMax(qRound(1023.0 * (float)(bc - imageptrs->minRadianceIndex[2] ) / (float)(imageptrs->maxRadianceIndex[2] - imageptrs->minRadianceIndex[2])), 0), 1023);
}
else if(histogrammethod == CMB_HISTO_NONE_100)
else if(this->histogrammethod == CMB_HISTO_NONE_100)
{
if(rc != fillvalue[0])
indexoutrc = rc;
Expand All @@ -1808,27 +1823,13 @@ void SegmentListGeostationary::ComposeSegmentImagenetCDFInThread(QStringList fil
if(bc != fillvalue[2])
indexoutbc = bc;
}
else if(histogrammethod == CMB_HISTO_EQUALIZE)
else if(this->histogrammethod == CMB_HISTO_EQUALIZE)
{
if( rc != fillvalue[0]) indexoutrc = (quint16)qMin(qMax(qRound((float)imageptrs->lut_ch[0][rc]), 0), 1023);
if( gc != fillvalue[1]) indexoutgc = (quint16)qMin(qMax(qRound((float)imageptrs->lut_ch[1][gc]), 0), 1023);
if( bc != fillvalue[2]) indexoutbc = (quint16)qMin(qMax(qRound((float)imageptrs->lut_ch[2][bc]), 0), 1023);
}

if( (rc == fillvalue[0] && dqfvalue[0] == -1) || (gc == fillvalue[1] && dqfvalue[1] == -1) || (bc == fillvalue[2] && dqfvalue[2] == -1))
{
r = 0;
g = 0;
b = 0;
}
else if( (rc == fillvalue[0] && dqfvalue[0] == 2) || (gc == fillvalue[1] && dqfvalue[1] == 2) || (bc == fillvalue[2] && dqfvalue[2] == 2))
{
r = 255;
g = 255;
b = 255;

}

if( (rc == fillvalue[0] && dqfvalue[0] == -1) || (gc == fillvalue[1] && dqfvalue[1] == -1) || (bc == fillvalue[2] && dqfvalue[2] == -1))
{
r = 0;
Expand All @@ -1837,9 +1838,9 @@ void SegmentListGeostationary::ComposeSegmentImagenetCDFInThread(QStringList fil
}
else
{
r = quint16(inversevector[0] ? (1023 - indexoutrc)/4 : indexoutrc/4);
g = quint16(inversevector[1] ? (1023 - indexoutgc)/4 : indexoutgc/4);
b = quint16(inversevector[2] ? (1023 - indexoutbc)/4 : indexoutbc/4);
r = quint16(this->inversevector[0] ? (1023 - indexoutrc)/4 : indexoutrc/4);
g = quint16(this->inversevector[1] ? (1023 - indexoutgc)/4 : indexoutgc/4);
b = quint16(this->inversevector[2] ? (1023 - indexoutbc)/4 : indexoutbc/4);
}

}
Expand All @@ -1848,17 +1849,17 @@ void SegmentListGeostationary::ComposeSegmentImagenetCDFInThread(QStringList fil
rc = *(imageptrs->ptrRed[0] + line * xdim + pixelx);
dqfvalue[0] = *(imageptrs->ptrDQF[0] + line * xdim + pixelx);

if(histogrammethod == CMB_HISTO_NONE_95)
if(this->histogrammethod == CMB_HISTO_NONE_95)
{
if(rc != fillvalue[0])
indexoutrc = (quint16)qMin(qMax(qRound(1023.0 * (float)(rc - imageptrs->minRadianceIndex[0] ) / (float)(imageptrs->maxRadianceIndex[0] - imageptrs->minRadianceIndex[0])), 0), 1023);
}
else if(histogrammethod == CMB_HISTO_NONE_100)
else if(this->histogrammethod == CMB_HISTO_NONE_100)
{
if(rc != fillvalue[0])
indexoutrc = rc;
}
else if(histogrammethod == CMB_HISTO_EQUALIZE)
else if(this->histogrammethod == CMB_HISTO_EQUALIZE)
{
if( rc != fillvalue[0]) indexoutrc = (quint16)qMin(qMax(qRound((float)imageptrs->lut_ch[0][rc]), 0), 1023);
}
Expand All @@ -1878,7 +1879,7 @@ void SegmentListGeostationary::ComposeSegmentImagenetCDFInThread(QStringList fil
}
else
{
r = quint16(inversevector[0] ? (1023 - indexoutrc)/4 : indexoutrc/4);
r = quint16(this->inversevector[0] ? (1023 - indexoutrc)/4 : indexoutrc/4);
g = r;
b = r;
}
Expand Down
Loading

0 comments on commit fc7051c

Please sign in to comment.