void SegmentListVIIRSM::CalculateProjectionLUT()
{
qDebug() << "start SegmentListVIIRSM::CalculateProjectionLUT()";
int earth_views;
long stats_ch[3][256];
long cnt_active_pixels = 0;
for(int k = 0; k < 3; k++)
{
for (int j = 0; j < 256; j++)
{
stats_ch[k][j] = 0;
}
}
for (int i=0; i < 3; i++)
{
for (int j=0; j < 1024; j++)
{
imageptrs->lut_ch[i][j] = 0;
}
}
for(int k = 0; k < 3; k++)
{
imageptrs->stat_max_ch[k] = 0;
imageptrs->stat_min_ch[k] = 9999999;
this->stat_max_ch[k] = 0;
this->stat_min_ch[k] = 9999999;
}
bool composecolor;
int x, y;
QList<Segment *>::iterator segsel = segsselected.begin();
while ( segsel != segsselected.end() )
{
SegmentVIIRSM *segm = (SegmentVIIRSM *)(*segsel);
segm->recalculateStatsInProjection();
++segsel;
}
segsel = segsselected.begin();
while ( segsel != segsselected.end() )
{
SegmentVIIRSM *segm = (SegmentVIIRSM *)(*segsel);
composecolor = segm->composeColorImage();
for(int i = 0; i < (composecolor ? 3 : 1); i++)
{
if( segm->stat_max_projection[i] > this->stat_max_ch[i])
this->stat_max_ch[i] = segm->stat_max_projection[i];
if( segm->stat_min_projection[i] < this->stat_min_ch[i])
this->stat_min_ch[i] = segm->stat_min_projection[i];
}
cnt_active_pixels += segm->active_pixels[0];
++segsel;
}
for(int i = 0; i < (composecolor ? 3 : 1); i++)
{
imageptrs->stat_max_ch[i] = this->stat_max_ch[i];
imageptrs->stat_min_ch[i] = this->stat_min_ch[i];
}
imageptrs->active_pixels = cnt_active_pixels;
segsel = segsselected.begin();
while ( segsel != segsselected.end() )
{
SegmentVIIRSM *segm = (SegmentVIIRSM *)(*segsel);
composecolor = segm->composeColorImage();
earth_views = segm->earth_views_per_scanline;
for(int k = 0; k < (composecolor ? 3 : 1); k++)
{
for (int line = 0; line < segm->NbrOfLines; line++)
{
for (int pixelx = 0; pixelx < segm->earth_views_per_scanline; pixelx++)
{
x = segm->getProjectionX(line, pixelx);
y = segm->getProjectionY(line, pixelx);
if(x >= 0 && x < imageptrs->ptrimageProjection->width() && y >= 0 && y < imageptrs->ptrimageProjection->height())
{
int pixel = *(segm->ptrbaVIIRS[k].data() + line * segm->earth_views_per_scanline + pixelx);
int pixcalc = 256 * (pixel - imageptrs->stat_min_ch[k]) / (imageptrs->stat_max_ch[k] - imageptrs->stat_min_ch[k]);
pixcalc = ( pixcalc < 0 ? 0 : pixcalc);
pixcalc = ( pixcalc > 255 ? 255 : pixcalc );
stats_ch[k][pixcalc]++;
}
}
}
}
++segsel;
}
float scale = 256.0 / (float)imageptrs->active_pixels;
unsigned long long sum_ch[3];
for (int i=0; i < 3; i++)
{
sum_ch[i] = 0;
}
for( int i = 0; i < 256; i++)
{
for(int k = 0; k < (composecolor ? 3 : 1); k++)
{
sum_ch[k] += stats_ch[k][i];
imageptrs->lut_ch[k][i] = (quint16)(sum_ch[k] * scale);
imageptrs->lut_ch[k][i] = ( imageptrs->lut_ch[k][i] > 255 ? 255 : imageptrs->lut_ch[k][i]);
}
}
}
