---

MDISlice.cpp

//----------------------------------------------------------------------------
//Borland C++Builder
//Copyright (c) 1987, 1998-2002 Borland International Inc. All Rights Reserved.
//----------------------------------------------------------------------------
//---------------------------------------------------------------------------
#include <vcl.h>
#include <math.h>
#pragma hdrstop

#include "MDISlice.h"
#include "MDIEdit.h"
#include "MDIFrame.h"
#include "DContrast.h"
//---------------------------------------------------------------------------
#pragma resource "*.dfm"
TSliceForm *SliceForm;

rgba *sliceview;

//---------------------------------------------------------------------------
__fastcall TSliceForm::TSliceForm(TComponent* Owner)
    : TForm(Owner)
{

}
//---------------------------------------------------------------------------
void __fastcall TSliceForm::SetPixelFormatDescriptor()
{
    PIXELFORMATDESCRIPTOR pfd = {
        sizeof(PIXELFORMATDESCRIPTOR),
        1,
        PFD_DRAW_TO_WINDOW | PFD_SUPPORT_OPENGL | PFD_DOUBLEBUFFER,
        PFD_TYPE_RGBA,
        24,
        0,0,0,0,0,0,
        0,0,
        0,0,0,0,0,
        32,
        0,
        0,
        PFD_MAIN_PLANE,
        0,
        0,0,
    };
    PixelFormat = ChoosePixelFormat(hdc, &pfd);
    SetPixelFormat(hdc, PixelFormat, &pfd);
}
//---------------------------------------------------------------------
void __fastcall TSliceForm::FormClose(TObject *Sender, TCloseAction &Action)
{
Action=caFree;
FrameForm->Slices->Checked = false;
delete [] sliceview;
memusage = memusage - sizeof(rgba)*512*512;
FrameForm->memupdate();
wglDeleteContext(hrc);
ReleaseDC(Handle,hdc);
}
//---------------------------------------------------------------------
void __fastcall TSliceForm::FormCreate(TObject *Sender)
{
    int i,a;

    sliceview = new rgba[512*512];
    memusage = memusage + sizeof(rgba)*512*512;
    FrameForm->memupdate();

    hdc = GetDC(Handle);
    SetPixelFormatDescriptor();
    hrc = wglCreateContext(hdc);
    if(hrc == NULL) ShowMessage("Could not create DC!");
    wglMakeCurrent(hdc, hrc);

    w = ClientWidth;
    h = ClientHeight;

    glEnable(GL_DEPTH_TEST);
    glEnable(GL_CULL_FACE);

    glClearColor(0.2f, 0.2f, 0.2f, 1.0f);
    glOrtho (-50.0,50.0,-50.0,50.0,-50.0,50.0);

    glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
    glFlush();
    SwapBuffers(hdc);
    FormResize(Sender);
    pagesize = sizeX*sizeY;
        for(int a = 0; a<sizeY; a+=1) {
        for(int i = 0; i<sizeX; i+=1) {
                sliceview[a*sizeX+i].r =
                        sliceview[a*sizeX+i].g =
                        sliceview[a*sizeX+i].b =
                                contr->getVal(data[pagesize*TrackBar->Position+a*sizeX+i]);
        }
        }
    TrackBar->Max = sizeZ-1;
}
//---------------------------------------------------------------------
void __fastcall TSliceForm::FormResize(TObject *Sender)
{
        tPanel->Height = ClientHeight-34;
        TrackBar->Height = ClientHeight-36;
        kontrast->Height = ClientHeight-188;
        kontrastpanel->Height = ClientHeight-235;
        CbTransfer->Top = ClientHeight-211;
        PrepareGLWindow();
        PaintGL();
}
//---------------------------------------------------------------------------
void __fastcall TSliceForm::FormPaint(TObject *Sender)
{
        //histo->drawTo(HistoArea);
        PaintGL();
}
//---------------------------------------------------------------------------
void __fastcall TSliceForm::PaintGL()
{
        float i;

        glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
        PrepareGLWindow();
        glColor3f(0.6f,0.6f,0.6f);
        glLineWidth(3.0f);
        glBegin(GL_LINE_STRIP);
        glVertex3f(-50,-50,0);
        glVertex3f(50,-50,0);
        glVertex3f(50,50,0);
        glVertex3f(-50,50,0);
        glVertex3f(-50,-50,0);
        glEnd();


//        glDrawPixels(sizeX, sizeY, GL_RGB, GL_UNSIGNED_BYTE,
//                     color+(sizeX*sizeY*TrackBar->Position));
if (RadioButton1->Checked) glDrawPixels(sizeX, sizeY, GL_RGBA, GL_UNSIGNED_BYTE, sliceview);
if (RadioButton2->Checked) glDrawPixels(sizeX, sizeZ, GL_RGBA, GL_UNSIGNED_BYTE, sliceview);
if (RadioButton3->Checked) glDrawPixels(sizeY, sizeZ, GL_RGBA, GL_UNSIGNED_BYTE, sliceview);
                        // we use unsigned byte (char) to meet speed issues.
                        // 16 bit data (__int16) slows down drawing to much.

        glFlush();
        SwapBuffers(hdc);
}

void __fastcall TSliceForm::PrepareGLWindow()
{
        wglMakeCurrent(hdc, hrc);

        GLfloat nRange = 50.0;
        GLfloat xRange, yRange;
        float yrat;

        w = ClientWidth-200;
        h = ClientHeight;

        if(h == 0) h = 1;
        if(w == 0) w = 1;

        glViewport(200, 0, w, h);

        glMatrixMode(GL_PROJECTION);
        glLoadIdentity();
if (RadioButton1->Checked) {
        yrat = (float)sizeX/sizeY;
        if (w <= h*yrat)
        {        glOrtho (-nRange, nRange, -nRange*h*yrat/w, nRange*h*yrat/w, -nRange, nRange);
                glPixelZoom((float)w/(float)sizeX, (float)w/((float)sizeY*yrat));}
         else
        {        glOrtho (-nRange*w/(h*yrat), nRange*w/(h*yrat), -nRange, nRange, -nRange, nRange);
                glPixelZoom((float)h*yrat/(float)sizeX, (float)h/((float)sizeY));}
}
if (RadioButton2->Checked) {
        yrat = (float)sizeX/sizeZ;
        if (w <= h*yrat)
        {        glOrtho (-nRange, nRange, -nRange*h*yrat/w, nRange*h*yrat/w, -nRange, nRange);
                glPixelZoom((float)w/(float)sizeX, (float)w/((float)sizeZ*yrat));}
         else
        {        glOrtho (-nRange*w/(h*yrat), nRange*w/(h*yrat), -nRange, nRange, -nRange, nRange);
                glPixelZoom((float)h*yrat/(float)sizeX, (float)h/(float)sizeZ);}
}
if (RadioButton3->Checked) {
        yrat = (float)sizeY/sizeZ;
        if (w <= h*yrat)
        {        glOrtho (-nRange, nRange, -nRange*h*yrat/w, nRange*h*yrat/w, -nRange, nRange);
                glPixelZoom((float)w/(float)sizeY, (float)w/((float)sizeZ*yrat));}
         else
        {        glOrtho (-nRange*w/(h*yrat), nRange*w/(h*yrat), -nRange, nRange, -nRange, nRange);
                glPixelZoom((float)h*yrat/(float)sizeY, (float)h/(float)sizeZ);}
}

        glRasterPos2d(-50,-50);
        glMatrixMode(GL_MODELVIEW);
        glLoadIdentity();
}
//---------------------------------------------------------------------------
void __fastcall TSliceForm::TrackBarChange(TObject *Sender)
{
if (CbTransfer->Checked)
{

if (RadioButton1->Checked)
{
        //float ratio = (float)0xFF / (float)4096;
        for(int a = 0; a<sizeY; a+=1) {
        for(int i = 0; i<sizeX; i+=1) {
                sliceview[a*sizeX+i] =
                        transe->getColor(data[pagesize*TrackBar->Position+a*sizeX+i]);
                sliceview[a*sizeX+i].r =
                        (int)((float)sliceview[a*sizeX+i].r*(float)sliceview[a*sizeX+i].a/(float)0xff);
                sliceview[a*sizeX+i].g =
                        (int)((float)sliceview[a*sizeX+i].g*(float)sliceview[a*sizeX+i].a/(float)0xff);
                sliceview[a*sizeX+i].b =
                        (int)((float)sliceview[a*sizeX+i].b*(float)sliceview[a*sizeX+i].a/(float)0xff);
        }
        }
}
if (RadioButton2->Checked)
{
        for(int a = 0; a<sizeZ; a+=1) {
        for(int i = 0; i<sizeX; i+=1) {
                //sliceview[a*sizeX+i].r = sliceview[a*sizeX+i].g = sliceview[a*sizeX+i].b = color[pagesize*a+TrackBar->Position*sizeX+i].b;
                sliceview[a*sizeX+i] =
                        transe->getColor(data[pagesize*a+TrackBar->Position*sizeX+i]);
                sliceview[a*sizeX+i].r =
                        (int)((float)sliceview[a*sizeX+i].r*(float)sliceview[a*sizeX+i].a/(float)0xff);
                sliceview[a*sizeX+i].g =
                        (int)((float)sliceview[a*sizeX+i].g*(float)sliceview[a*sizeX+i].a/(float)0xff);
                sliceview[a*sizeX+i].b =
                        (int)((float)sliceview[a*sizeX+i].b*(float)sliceview[a*sizeX+i].a/(float)0xff);
        }
        }
}
if (RadioButton3->Checked)
{
        for(int a = 0; a<sizeZ; a+=1) {
        for(int i = 0; i<sizeY; i+=1) {
                //sliceview[a*sizeY+i].r = sliceview[a*sizeY+i].g = sliceview[a*sizeY+i].b = color[TrackBar->Position+i*sizeX+a*pagesize].b;
                sliceview[a*sizeY+i] = transe->getColor(data[TrackBar->Position+i*sizeX+a*pagesize]);
                sliceview[a*sizeY+i].r =
                        (int)((float)sliceview[a*sizeY+i].r*(float)sliceview[a*sizeY+i].a/(float)0xff);
                sliceview[a*sizeY+i].g =
                        (int)((float)sliceview[a*sizeY+i].g*(float)sliceview[a*sizeY+i].a/(float)0xff);
                sliceview[a*sizeY+i].b =
                        (int)((float)sliceview[a*sizeY+i].b*(float)sliceview[a*sizeY+i].a/(float)0xff);
        }
        }
}

} else {

if (RadioButton1->Checked)
{
        //float ratio = (float)0xFF / (float)4096;
        for(int a = 0; a<sizeY; a+=1) {
        for(int i = 0; i<sizeX; i+=1) {
                //sliceview[a*sizeX+i].r = sliceview[a*sizeX+i].g = sliceview[a*sizeX+i].b = color[pagesize*TrackBar->Position+a*sizeX+i].b;
                sliceview[a*sizeX+i].r =
                        sliceview[a*sizeX+i].g =
                        sliceview[a*sizeX+i].b =
                                contr->getVal(data[pagesize*TrackBar->Position+a*sizeX+i]);
        }
        }
}
if (RadioButton2->Checked)
{
        for(int a = 0; a<sizeZ; a+=1) {
        for(int i = 0; i<sizeX; i+=1) {
                //sliceview[a*sizeX+i].r = sliceview[a*sizeX+i].g = sliceview[a*sizeX+i].b = color[pagesize*a+TrackBar->Position*sizeX+i].b;
                sliceview[a*sizeX+i].r =
                        sliceview[a*sizeX+i].g =
                        sliceview[a*sizeX+i].b =
                                contr->getVal(data[pagesize*a+TrackBar->Position*sizeX+i]);
        }
        }
}
if (RadioButton3->Checked)
{
        for(int a = 0; a<sizeZ; a+=1) {
        for(int i = 0; i<sizeY; i+=1) {
                //sliceview[a*sizeY+i].r = sliceview[a*sizeY+i].g = sliceview[a*sizeY+i].b = color[TrackBar->Position+i*sizeX+a*pagesize].b;
                sliceview[a*sizeY+i].r =
                        sliceview[a*sizeY+i].g =
                        sliceview[a*sizeY+i].b =
                                contr->getVal(data[TrackBar->Position+i*sizeX+a*pagesize]);
        }
        }
}

}
lslice->Caption = IntToStr(TrackBar->Position) + " / " + IntToStr(TrackBar->Max);
PaintGL();
}
//---------------------------------------------------------------------------

void __fastcall TSliceForm::RadioButton1Click(TObject *Sender)
{
TrackBar->Max = sizeZ-1;
TrackBarChange(this);
}
//---------------------------------------------------------------------------

void __fastcall TSliceForm::RadioButton2Click(TObject *Sender)
{
TrackBar->Max = sizeY-1;
TrackBarChange(this);
}
//---------------------------------------------------------------------------

void __fastcall TSliceForm::RadioButton3Click(TObject *Sender)
{
TrackBar->Max = sizeX-1;
TrackBarChange(this);
}
//---------------------------------------------------------------------------

void __fastcall TSliceForm::FormActivate(TObject *Sender)
{
        wglMakeCurrent(hdc, hrc);
        PaintGL();
}
//---------------------------------------------------------------------------
void __fastcall TSliceForm::HistoAreaDblClick(TObject *Sender)
{
        ContrastForm = new TContrastForm(this, histo, contr);
        ContrastForm->Show();
}
//---------------------------------------------------------------------------

void __fastcall TSliceForm::HistoAreaPaint(TObject *Sender)
{
        histo->drawTo(HistoArea);
        int l = (int)(contr->Lower*(float)histo->size/histo->getHPP(HistoArea));
        int u = (int)(contr->Upper*(float)histo->size/histo->getHPP(HistoArea));
        HistoArea->Canvas->Pen->Color = clBlack;
        if (histo->horizontal) {
                HistoArea->Canvas->MoveTo(l,0);
                HistoArea->Canvas->LineTo(l,HistoArea->Height);
                HistoArea->Canvas->MoveTo(u,0);
                HistoArea->Canvas->LineTo(u,HistoArea->Height);
        } else {
                HistoArea->Canvas->MoveTo(0,l);
                HistoArea->Canvas->LineTo(HistoArea->Width,l);
                HistoArea->Canvas->MoveTo(0,u);
                HistoArea->Canvas->LineTo(HistoArea->Width,u);
        }
}
//---------------------------------------------------------------------------

---