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]); } } }