Jump to content
  • Advertisement
Sign in to follow this  
dagmare

OpenGL texture help

This topic is 4245 days old which is more than the 365 day threshold we allow for new replies. Please post a new topic.

If you intended to correct an error in the post then please contact us.

Recommended Posts

I wrote an OpenGL program with swiching scenes. I have a problem with one of the scenes- it shows the objects, but the texture put on quad doesnt work. (I have used this textures in previous scenes and everything was ok). I am curious if anybody could check up on my code? this is function switching between the scenes: void __fastcall TFormMain::DrawObjects() { UpDate(); double hour, min, sec; hour=floor(T/60/60); min=floor((T-hour*60*60)/60); sec=floor(T-hour*60*60-min*60); String time, timeh, timem, times; String dates, yrs,mos,ds,hs,mins,ss, timest; timeh=FloatToStr(hour); timem=FloatToStr(min); times=FloatToStr(sec); yrs=IntToStr(int(yr)); mos=IntToStr(int(mo)); ds=IntToStr(int(d)); hs=IntToStr(int(h)); mins=IntToStr(int(m)); ss=IntToStr(int(s)); dates=yrs+'.'+mos+'.'+ds; timest=hs+':'+mins+':'+ss; time=timeh+'h'+timem+'m'+times+'s'; FormMain->Label2->Caption=time; FormMain->Label3->Caption=dates; FormMain->Label4->Caption=timest; switch(object) { case 0: SceneOne(); break; case 1: SceneTwo(); break; case 2: SceneThree(); break; case 3: SceneFour(); break; case 4: SceneFive(); break; case 5: SceneSix(); break; } } this is scene five - the problematic one: void __fastcall TFormMain::SceneFive() { glLoadIdentity(); GroundTrackMap(); if(Child->CheckBox5->Checked) { glLoadIdentity(); DrawGroundTrack(); } glLoadIdentity(); AnimTrack(); } and this the function with texture which is not working: void __fastcall TFormMain::GroundTrackMap() { glLoadIdentity(); glColor3f(1.0,1.0,1.0); glBindTexture(GL_TEXTURE_2D, texture1); glBegin(GL_QUADS); glTexCoord2f(1.0f, 0.0f); glVertex3f(8400,-6400,-15000); glTexCoord2f(0.0f, 0.0f); glVertex3f(-8400,-6400,-15000); glTexCoord2f(0.0f, 1.0f); glVertex3f(-8400,6400,-15000); glTexCoord2f(1.0f, 1.0f); glVertex3f(8400,6400,-15000); glEnd(); } the rest of the functions in scene five is working ok. Could anybody help me with that? actually if there is anybody who can check all my code, than i can pay for that- i am desperate to finish this thing quick :)

Share this post


Link to post
Share on other sites
Advertisement
Have you glEnable(GL_TEXTURE_2D)'d the other scenes but not scene five? It might not hurt to toss it into the beginning of the drawing, just to make sure. Also, you seem to have LOTS of glLoadIdentity( )s, which are unnecessary and a performance hit (though not the cause of your problem, just something to note). If that doesn't help, post some more code inside of [source] and [/source] tags (it'll be formatted and pretty, much easier to read).

-jouley

Share this post


Link to post
Share on other sites
yes i did, and even i add this line to scene five, it doesnt help. Is it possible that it is something with my lighting settings? I can send some more code via e-mail if you be so kind to have a look, because it is quite long. Is that ok?

and for now:

//---------------------------------------------------------------------------
#include <vcl.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <windows.h>
#include <fstream>
#include <iostream.h>
#include <condefs.h>
#include <conio.h>
#pragma hdrstop

#include "GlSkel.h"
#include "child1.h"
#include "nSym.h"
#include "incomplete.h"
#include "perturbation.h"

//---------------------------------------------------------------------------
//#pragma hdrstop
#pragma package(smart_init)
#pragma resource "*.dfm"
TFormMain *FormMain;

//FILE *test;

//---------------------------------------------------------------------------
__fastcall TFormMain::TFormMain(TComponent* Owner)
: TForm(Owner)
{ //Child->nowa=FALSE;

//test=fopen("ground.txt","w");
Application->OnIdle = IdleLoop;
//size = 0.50f;
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::IdleLoop(TObject*, bool& done)
{
done = false;
RenderGLScene();
SwapBuffers(hdc);
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::FormCreate(TObject *Sender)
{ StartValue();

hdc = GetDC(Handle);
SetPixelFormatDescriptor();
hrc = wglCreateContext(hdc);
if(hrc == NULL)
ShowMessage(":-)~ hrc == NULL");
if(wglMakeCurrent(hdc, hrc) == false)
ShowMessage("Could not MakeCurrent");
w = ClientWidth;
h = ClientHeight;

glEnable(GL_DEPTH_TEST);
glEnable(GL_CULL_FACE);

glClearColor(0.0f, 0.0f, 0.0f, 0.5f);
SetupLighting();
SetupTextures();

}
//---------------------------------------------------------------------------
void __fastcall TFormMain::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,
16,
0,
0,
PFD_MAIN_PLANE,
0,
0,0,
};
PixelFormat = ChoosePixelFormat(hdc, &pfd);
SetPixelFormat(hdc, PixelFormat, &pfd);
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::FormResize(TObject *Sender)
{
glEnable(GL_TEXTURE_2D); // Enable Texture Mapping
glShadeModel(GL_SMOOTH); // Enable Smooth Shading
glClearColor(0.0f, 0.0f, 0.0f, 0.5f); // Black Background
glClearDepth(1.0f); // Depth Buffer Setup
glEnable(GL_DEPTH_TEST); // Enables Depth Testing
glDepthFunc(GL_LEQUAL); // The Type Of Depth Testing To Do
glHint(GL_PERSPECTIVE_CORRECTION_HINT, GL_NICEST); // Really Nice Perspective Calculations

w = ClientWidth;
h = ClientHeight;

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


glMatrixMode(GL_PROJECTION);
glLoadIdentity();
gluPerspective(45.0f,(GLfloat)w/(GLfloat)h,10000.1f,200000.0f);

glMatrixMode(GL_MODELVIEW);
glLoadIdentity();
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::RenderGLScene()
{
T0=T;T=T+60;
JD=julian(yr, mo, d, h, m, s);
//cieniuj(JD,T);

glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
DrawObjects();
glFlush();
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::DrawObjects()
{ UpDate();

double hour, min, sec;

hour=floor(T/60/60);
min=floor((T-hour*60*60)/60);
sec=floor(T-hour*60*60-min*60);

String time, timeh, timem, times;
String dates, yrs,mos,ds,hs,mins,ss, timest;

timeh=FloatToStr(hour);
timem=FloatToStr(min);
times=FloatToStr(sec);

yrs=IntToStr(int(yr));
mos=IntToStr(int(mo));
ds=IntToStr(int(d));
hs=IntToStr(int(h));
mins=IntToStr(int(m));
ss=IntToStr(int(s));

dates=yrs+'.'+mos+'.'+ds;
timest=hs+':'+mins+':'+ss;

time=timeh+'h'+timem+'m'+times+'s';
FormMain->Label2->Caption=time;
FormMain->Label3->Caption=dates;
FormMain->Label4->Caption=timest;
switch(object)
{
case 0:
SceneOne();
break;
case 1:
SceneTwo();
break;
case 2:
SceneThree();
break;
case 3:
SceneFour();
break;
case 4:
SceneFive();
break;
case 5:
SceneSix();
break;
}
}

//---------------------------------------------------------------------------
void __fastcall TFormMain::SceneOne()
{
glLoadIdentity();
glColor3f(0.9f, 0.5f, 0.5f);
glTranslatef(0.0,0.0,z);
glRotatef(xspeed,1.0f,0.0f,0.0f);
glRotatef(yspeed,0.0f,0.0f,1.0f);
DrawAxis();

glLoadIdentity();
glColor3f(0.9f, 0.5f, 0.5f);
glTranslatef(0.0,0.0,z);
glRotatef(xspeed,1.0f,0.0f,0.0f);
glRotatef(yspeed,0.0f,0.0f,1.0f);
AnimGLSat();

glLoadIdentity();
glTranslatef(0.0,0.0,z);
glRotatef(xspeed,1.0f,0.0f,0.0f);
glRotatef(yspeed,0.0f,0.0f,1.0f);
DrawEarth();
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::SceneTwo()
{

glLoadIdentity();
glColor3f(0.9f, 0.5f, 0.5f);
glTranslatef(0.0,0.0,z);
glRotatef(xspeed,1.0f,0.0f,0.0f);
glRotatef(yspeed,0.0f,0.0f,1.0f);
DrawAxis();

glLoadIdentity();
DrawOrbits();

glLoadIdentity();
glColor3f(0.9f, 0.5f, 0.5f);
glTranslatef(0.0,0.0,z);
glRotatef(xspeed,1.0f,0.0f,0.0f);
glRotatef(yspeed,0.0f,0.0f,1.0f);
AnimGLSat();

glLoadIdentity();
glTranslatef(0.0,0.0,z);
glRotatef(xspeed,1.0f,0.0f,0.0f);
glRotatef(yspeed,0.0f,0.0f,1.0f);
DrawEarth();
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::SceneThree()
{
glLoadIdentity();
glColor3f(0.9f, 0.5f, 0.5f);
glTranslatef(0.0,0.0,z);
glRotatef(xspeed,1.0f,0.0f,0.0f);
glRotatef(yspeed,0.0f,0.0f,1.0f);
DrawAxis();

glLoadIdentity();
AnimGLEcf();

glLoadIdentity();
glTranslatef(0.0,0.0,z);
glRotatef(xspeed,1.0f,0.0f,0.0f);
glRotatef(yspeed,0.0f,0.0f,1.0f);
DrawEarthECF();
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::SceneFour()
{
glLoadIdentity();
glColor3f(0.9f, 0.5f, 0.5f);
glTranslatef(0.0,0.0,z);
glRotatef(xspeed,1.0f,0.0f,0.0f);
glRotatef(yspeed,0.0f,0.0f,1.0f);
DrawAxis();

glLoadIdentity();
AnimGLEcf();

glLoadIdentity();
glTranslatef(0.0,0.0,z);
glRotatef(xspeed,1.0f,0.0f,0.0f);
glRotatef(yspeed,0.0f,0.0f,1.0f);
DrawGLEcf();

glLoadIdentity();
glTranslatef(0.0,0.0,z);
glRotatef(xspeed,1.0f,0.0f,0.0f);
glRotatef(yspeed,0.0f,0.0f,1.0f);
DrawEarthECF();
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::SceneFive()
{
glLoadIdentity();
GroundTrackMap();
if(Child->CheckBox5->Checked)
{
glLoadIdentity();
DrawGroundTrack();
}
glLoadIdentity();
AnimTrack();
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::SceneSix()
{
glLoadIdentity();
FormMain->Label2->Caption=' ';

}
//---------------------------------------------------------------------------
void __fastcall TFormMain::FormPaint(TObject *Sender)
{
glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
glFlush();
DrawObjects();
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::FormDestroy(TObject *Sender)
{
gluDeleteQuadric(Satqua);
gluDeleteQuadric(Earth);
gluDeleteQuadric(Earth1);
delete bitmap;
wglMakeCurrent(NULL, NULL);
wglDeleteContext(hrc);
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::StartValue()
{ for(loop=0;loop<num;loop++)
{ double w[3];
sat[loop].r=rand()%4;
sat[loop].g=rand()%3;
sat[loop].b=rand()%2;
/* while(sat[loop].r<0.3 && sat[loop].g<0.3 && sat[loop].b<0.3)
{
sat[loop].r=rand()%4;
sat[loop].g=rand()%4;
sat[loop].b=rand()%3;
} */

sat[loop].rs=sat[loop].r;
sat[loop].gs=sat[loop].g;
sat[loop].bs=sat[loop].b;
sat[loop].el[5]=loop*10*PI/180.0;//mean anomaly
sat[loop].el[2]=60*PI/180.0;//ink
sat[loop].el[4]=loop*15*PI/180.0;//node
sat[loop].el[0]=20000.0;//semimajor axis
sat[loop].el[1]=0.6; //eccentricity
sat[loop].el[3]=15*PI/180.0; // argument of perigee

secular(sat[loop].el,w);
sat[loop].Ms=(*(w+2));
sat[loop].Nodes=(*(w+1));
sat[loop].oms=(*w);

}
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::DrawAxis()
{
glBegin(GL_LINES);
{
glVertex3f(-100000,0,0);
glVertex3f(100000,0,0);
glVertex3f(0,-100000,0);
glVertex3f(0,100000,0);
glVertex3f(0,0,-100000);
glVertex3f(-0,0,100000);
}
glEnd();
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::DrawEarth()
{

double ziemia;
ziemia=(green+7.29211585530*pow(10,-5)*T)*180/PI;

Earth = gluNewQuadric();
gluQuadricTexture(Earth, GL_TRUE);
gluQuadricDrawStyle(Earth, GLU_FILL);
gluQuadricNormals(Earth, GLU_SMOOTH);
glBindTexture(GL_TEXTURE_2D, texture1);
glColor3f(1.0f,1.0f,1.0f);
glRotatef(-90,0.0f,0.0f,1.0f);
glRotatef(ziemia,0.0f,0.0f,1.0f);
gluSphere(Earth, 6378.137f,320000,320000);
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::DrawEarthECF()
{
Earth1 = gluNewQuadric();
gluQuadricTexture(Earth1, GL_TRUE);
gluQuadricDrawStyle(Earth, GLU_FILL);
gluQuadricNormals(Earth1, GLU_SMOOTH);
glBindTexture(GL_TEXTURE_2D, texture1);
glColor3f(1.0f,1.0f,1.0f);
glRotatef(-90,0.0f,0.0f,1.0f);
gluSphere(Earth1, 6378.137f,320000,320000);

}
//---------------------------------------------------------------------------
void __fastcall TFormMain::DrawSatellite()
{
Satqua = gluNewQuadric();
gluQuadricDrawStyle(Satqua, GLU_FILL);
gluQuadricNormals(Satqua, GLU_SMOOTH);
//glColor3f(sat.r,sat.g,sat.b);
gluSphere(Satqua,200.0f,32,32);

}
//---------------------------------------------------------------------------
void __fastcall TFormMain::AnimGLSat()
{ int i;
double wiekowe[3];

for(i=0;i<num;i++)
{
el[5]=sat.el[5];
el[2]=sat.el[2];
el[4]=sat.el[4];
el[0]=sat.el[0];
el[1]=sat.el[1];
el[3]=sat.el[3];

posvel(T,T0,el,pos,vel); //calculates position and velocity of the satellite

sat.x=pos[0];
sat.y=pos[1];
sat.z=pos[2];

elements(vel,pos,el,T0,T); //calculates orbital elements of the satellite

secular(el,wiekowe);
sat.Ms=(*(wiekowe+2));
sat.Nodes=(*(wiekowe+1));
sat.oms=(*wiekowe);
sat.el[5]=el[5]+(*(wiekowe+2))*(T-T0);
sat.el[2]=el[2];
sat.el[4]=el[4]+(*(wiekowe+1))*(T-T0);
sat.el[0]=el[0];
sat.el[1]=el[1];
sat.el[3]=el[3]+(*wiekowe)*(T-T0);
//fprintf(test,"%.10f %.10f %.10f\n", vel[3], T, sat[0].el[5]);
glLoadIdentity();
glTranslatef(0.0,0.0,z);
glRotatef(xspeed,1.0f,0.0f,0.0f);
glRotatef(yspeed,0.0f,0.0f,1.0f);
glTranslatef(sat.x,sat.y,sat.z);
glColor3f(sat.r,sat.g,sat.b);
DrawSatellite();
}
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::AnimGLEcf()
{ int i;
double wiek[3],rff[4],elf[9],posf[4],velf[4];

for(i=0;i<num;i++)
{
elf[5]=sat.el[5];
elf[2]=sat.el[2];
elf[4]=sat.el[4];
elf[0]=sat.el[0];
elf[1]=sat.el[1];
elf[3]=sat.el[3];

posvel(T,T0,elf,posf,velf); //calculates position and velocity of the satellite
earthfix(posf,rff,T);

elements(velf,posf,elf,T0,T); //calculates orbital elements of the satellite

secular(elf,wiek);
sat.Ms=(*(wiek+2));
sat.Nodes=(*(wiek+1));
sat.oms=(*wiek);
sat.el[5]=elf[5]+(*(wiek+2))*(T-T0);
sat.el[2]=elf[2];
sat.el[4]=elf[4]+(*(wiek+1))*(T-T0);
sat.el[0]=elf[0];
sat.el[1]=elf[1];
sat.el[3]=elf[3]+(*wiek)*(T-T0);

glLoadIdentity();
glTranslatef(0.0,0.0,z);
glRotatef(xspeed,1.0f,0.0f,0.0f);
glRotatef(yspeed,0.0f,0.0f,1.0f);
glTranslatef(rff[0],rff[1],rff[2]);
glColor3f(sat.r,sat.g,sat.b);
DrawSatellite();
}

}
//---------------------------------------------------------------------------
void __fastcall TFormMain::DrawGLEcf()
{ int i,k;
double rf[4],ele[9],polo[4],pre[4];
double Te=0,T0e=0;
double sece[3];
double rep;
int x;

glLoadIdentity();
glTranslatef(0.0,0.0,z);
glRotatef(xspeed,1.0f,0.0f,0.0f);
glRotatef(yspeed,0.0f,0.0f,1.0f);
glColor3f(1.0,0.4,0.4);

if(Child->CheckBox5->Checked)
{
x=1;
rep=Child->repeat;
}
else
{
x=num;
rep=1000;
}

for(i=0;i<x;i++)
{ T0e=T0;Te=T;

ele[5]=sat.el[5];;
ele[2]=sat.el[2];
ele[4]=sat.el[4];
ele[0]=sat.el[0];
ele[1]=sat.el[1];
ele[3]=sat.el[3];
k=0;
glBegin(GL_LINE_STRIP);
while(k<rep)
{
k++;

posvel(Te,T0e,ele,polo,pre); //calculates position and velocity of the satellite
earthfix(polo,rf,Te);

elements(pre,polo,ele,T0e,Te); //calculates orbital elements of the satellite

secular(ele,sece);
ele[5]=ele[5]+(*(sece+2))*(Te-T0e);
ele[4]=ele[4]+(*(sece+1))*(Te-T0e);
ele[3]=ele[3]+(*sece)*(Te-T0e);

glLoadIdentity();
glVertex3f(rf[0],rf[1],rf[2]);
T0e=Te;
Te=Te+60;
} glEnd();
}
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::DrawOrbits()
{
int l;
double ele[9],pre[4],polo[4];
double To=0,T0o=0;
double P,nn;
double GM=398600.64;
double sec[3];

To=T; T0o=T0;

for(l=0;l<num;l++)
{
To=0,T0o=0;

ele[0]=sat[l].el[0];
ele[1]=sat[l].el[1];
ele[2]=sat[l].el[2];
ele[3]=sat[l].el[3];
ele[4]=sat[l].el[4];
ele[5]=sat[l].el[5];

nn=sqrt(GM/pow(ele[0],3));
P=2*PI/nn;

glLoadIdentity();
glColor3f(1.0,0.4,0.4);
glTranslatef(0.0,0.0,z);
glRotatef(xspeed,1.0f,0.0f,0.0f);
glRotatef(yspeed,0.0f,0.0f,1.0f);

glBegin(GL_LINE_STRIP);
{
while(T0o<P+100)
{ T0o=To;
To=To+60;
posvel(To,T0o,ele,polo,pre);
elements(pre,polo,ele,T0o,To);
secular(ele,sec);
ele[5]=ele[5]+(*(sec+2))*(To-T0o);
ele[4]=ele[4]+(*(sec+1))*(To-T0o);
ele[3]=ele[3]+(*sec)*(To-T0o);
glVertex3f(polo[0],polo[1],polo[2]);

}
}glEnd();

}
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::GroundTrackMap()
{
glLoadIdentity();
glColor3f(1.0,1.0,1.0);
glEnable(GL_TEXTURE_2D);
glBindTexture(GL_TEXTURE_2D, texture1);

glBegin(GL_QUADS);
glTexCoord2f(1.0f, 0.0f); glVertex3f(8400,-6400,-15000);
glTexCoord2f(0.0f, 0.0f); glVertex3f(-8400,-6400,-15000);
glTexCoord2f(0.0f, 1.0f); glVertex3f(-8400,6400,-15000);
glTexCoord2f(1.0f, 1.0f); glVertex3f(8400,6400,-15000);
glEnd();
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::AnimTrack()
{
int l,k;
double elt[9],velt[4],post[4];
//double Tt=0,T0t=0;
double ratio1, ratio2;
double sect[3],r[4];
double earthrot=7.29211585530*pow(10,-5);
double teta, latitude, longitude,teta0=0;
teta0=green;
ratio1=2*8400/(2*PI);
ratio2=2*6400/PI;
glColor3f(0.5,0.5,0.5);

glLoadIdentity();
glPointSize(5);
glBegin(GL_POINTS);
for(l=0;l<num;l++)
{ //double T0t=T0, Tt=T;
elt[0]=sat[l].el[0];
elt[1]=sat[l].el[1];
elt[2]=sat[l].el[2];
elt[3]=sat[l].el[3];
elt[4]=sat[l].el[4];
elt[5]=sat[l].el[5];
// glPointSize(0.5);
// glBegin(GL_POINTS);
// for(k=0;k<2000;k++)
// { T0t=Tt;
// Tt=Tt+30;
posvel(T,T0,elt,post,velt);
earthfix(post,r,T);
// sat[l].x=polo[0];
// sat[l].y=polo[1];
// sat[l].z=polo[2];

teta=teta0+earthrot*T;
// groundtrack(post,teta,latitude, longitude);
// fprintf(test,"%f %f\n",latitude,longitude);
latitude=asin(post[2]/post[3]);

L99: if(latitude<-PI/2)
latitude=latitude+PI;
if(latitude<-PI/2) goto L99;
// glColor3f(1,0,0);

longitude=atan2(post[1],post[0]);
longitude=longitude-teta;
L100: if(longitude<0)
{
longitude=longitude+2*PI;
}
if(longitude<0) goto L100;

elements(velt,post,elt,T0,T);
secular(elt,sect);
elt[5]=elt[5]+(*(sect+2))*(T-T0);
elt[4]=elt[4]+(*(sect+1))*(T-T0);
elt[3]=elt[3]+(*sect)*(T-T0);

sat[l].oms=(*sect);
sat[l].Nodes=(*(sect+1));
sat[l].Ms=(*(sect+2));

sat[l].el[0]=elt[0];
sat[l].el[1]=elt[1];
sat[l].el[2]=elt[2];
sat[l].el[3]=elt[3];
sat[l].el[4]=elt[4];
sat[l].el[5]=elt[5];
glLoadIdentity();
glColor3f(sat[l].r,sat[l].g,sat[l].b);
glVertex3f(ratio1*longitude-8400,ratio2*latitude,-15000);
}
glEnd(); //Sleep(40);
// } glEnd();
}

//---------------------------------------------------------------------------
//---------------------------------------------------------------------------
void __fastcall TFormMain::DrawGroundTrack()
{
int l,k;
double elt[9],velt[4],post[4];
double Tt=0,T0t=0;
double ratio1, ratio2;
double sect[3],r[4];
double earthrot=7.29211585530*pow(10,-5);
double teta, latitude, longitude,teta0=0;


ratio1=2*8400/(2*PI);
ratio2=2*6400/PI;
glColor3f(1.0,1.0,1.0);
teta0=green;
glLoadIdentity();
glPointSize(1);
// elt[5]=sat[l].el[5];
//glBegin(GL_POINTS);
for(l=0;l<1;l++)
{ T0t=T0, Tt=T;
elt[0]=sat[l].el[0];
elt[1]=sat[l].el[1];
elt[2]=sat[l].el[2];
elt[3]=sat[l].el[3];
elt[4]=sat[l].el[4];
elt[5]=sat[l].el[5];


glPointSize(1.0);
glBegin(GL_POINTS);
for(k=0;k<2000;k++)
{ T0t=Tt;
Tt=Tt+600;
posvel(Tt,T0t,elt,post,velt);
earthfix(post,r,Tt);
// sat[l].x=polo[0];
// sat[l].y=polo[1];
// sat[l].z=polo[2];

teta=teta0+earthrot*Tt;
// groundtrack(post,teta,latitude, longitude);
// fprintf(test,"%f %f\n",latitude,longitude);
latitude=asin(post[2]/post[3]);

L99: if(latitude<-PI/2)
latitude=latitude+PI;
if(latitude<-PI/2) goto L99;
// glColor3f(1,0,0);

longitude=atan2(post[1],post[0]);
longitude=longitude-teta;
L100: if(longitude<0)
{
longitude=longitude+2*PI;
}
if(longitude<0) goto L100;

elements(velt,post,elt,T0t,Tt);
secular(elt,sect);
elt[5]=elt[5]+(*(sect+2))*(Tt-T0t);
elt[4]=elt[4]+(*(sect+1))*(Tt-T0t);
elt[3]=elt[3]+(*sect)*(Tt-T0t);
/*
sat[l].oms=(*sect);
sat[l].Nodes=(*(sect+1));
sat[l].Ms=(*(sect+2));
*/
/*
sat[l].el[0]=elt[0];
sat[l].el[1]=elt[1];
sat[l].el[2]=elt[2];
sat[l].el[3]=elt[3];
sat[l].el[4]=elt[4];
sat[l].el[5]=elt[5];*/

glLoadIdentity();
glColor3f(1.0,1.0,1.0);
//glColor3f(0,0,0);
glVertex3f(ratio1*longitude-8400,ratio2*latitude,-15000);
}
glEnd();// Sleep(40);
} //glEnd();
}

//----------------------------------------------------------------------------
void __fastcall TFormMain::FormKeyDown(TObject *Sender, WORD &Key,
TShiftState Shift)
{
if(Key == VK_UP)
xspeed-=0.5f;
if(Key == VK_DOWN)
xspeed+=0.5f;
if(Key == VK_LEFT)
yspeed-=0.5f;
if(Key == VK_RIGHT)
yspeed+=0.5f;
if(Key == VK_PRIOR)
z-=200.5f;
if(Key == VK_NEXT)
z+=200.5f;
if (Key==' ')
{
object++;
if(object>4)
object=0;
}
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::SetupLighting()
{
GLfloat MaterialAmbient[] = {0.5, 0.5, 0.5, 1.0};
GLfloat MaterialDiffuse[] = {1.0, 1.0, 1.0, 1.0};
GLfloat MaterialSpecular[] = {1.0, 1.0, 1.0, 1.0};
GLfloat MaterialShininess[] = {50.0};
GLfloat AmbientLightPosition[] = {1.0, 1.0, 1.0, 0.0};
GLfloat LightAmbient[] = {0.5, 0.5, 0.5, 1.0};

glMaterialfv(GL_FRONT, GL_AMBIENT, MaterialAmbient);
glMaterialfv(GL_FRONT, GL_DIFFUSE, MaterialDiffuse);
glMaterialfv(GL_FRONT, GL_SPECULAR, MaterialSpecular);
glMaterialfv(GL_FRONT, GL_SHININESS, MaterialShininess);
glLightfv(GL_LIGHT0,GL_POSITION, AmbientLightPosition);
glLightModelfv(GL_LIGHT_MODEL_AMBIENT, LightAmbient);

glEnable(GL_LIGHTING);
glEnable(GL_LIGHT0);
glEnable(GL_COLOR_MATERIAL);
glColorMaterial(GL_FRONT, GL_AMBIENT_AND_DIFFUSE);
glShadeModel(GL_SMOOTH);
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::SetupTextures()
{
bitmap = new Graphics::TBitmap;
bitmap->LoadFromFile("Data/earth.bmp");
GLubyte bits[256][256][4];
for(int i = 0; i < 256; i++)
{
for(int j = 0; j < 256; j++)
{
bits[j][0]= (GLbyte)GetRValue(bitmap->Canvas->Pixels[j]);
bits[j][1]= (GLbyte)GetGValue(bitmap->Canvas->Pixels[j]);
bits[j][2]= (GLbyte)GetBValue(bitmap->Canvas->Pixels[j]);
bits[j][3]= (GLbyte)255;
}
}
glPixelStorei(GL_UNPACK_ALIGNMENT, 4);
glGenTextures(1, &texture1);
glBindTexture(GL_TEXTURE_2D, texture1);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_CLAMP);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_CLAMP);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
glTexImage2D(GL_TEXTURE_2D, 0, GL_RGBA, 256, 256, 0, GL_RGBA, GL_UNSIGNED_BYTE, bits);
/*
bitmap->LoadFromFile("Data/earth.bmp");
for(int i = 0; i < 64; i++)
{
for(int j = 0; j < 64; j++)
{
bits[j][0]= (GLbyte)GetRValue(bitmap->Canvas->Pixels[j]);
bits[j][1]= (GLbyte)GetGValue(bitmap->Canvas->Pixels[j]);
bits[j][2]= (GLbyte)GetBValue(bitmap->Canvas->Pixels[j]);
bits[j][3]= (GLbyte)255;
}
}
glPixelStorei(GL_UNPACK_ALIGNMENT, 4);
glGenTextures(1, &texture2);
glBindTexture(GL_TEXTURE_2D, texture2);
//glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_CLAMP);
//glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_CLAMP);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_LINEAR);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR);
//glTexParameteri(GL_TEXTURE_2D,GL_TEXTURE_MIN_FILTER,GL_LINEAR_MIPMAP_NEAREST);
glTexImage2D(GL_TEXTURE_2D, 0, 3, 64, 64, 0, GL_RGB, GL_UNSIGNED_BYTE, bits);

bitmap->LoadFromFile("Data/earth.bmp");
for(int i = 0; i < 64; i++)
{
for(int j = 0; j < 64; j++)
{
bits[j][0]= (GLbyte)GetRValue(bitmap->Canvas->Pixels[j]);
bits[j][1]= (GLbyte)GetGValue(bitmap->Canvas->Pixels[j]);
bits[j][2]= (GLbyte)GetBValue(bitmap->Canvas->Pixels[j]);
bits[j][3]= (GLbyte)255;
}
}



glPixelStorei(GL_UNPACK_ALIGNMENT, 4);
glGenTextures(1, &texture3);
glBindTexture(GL_TEXTURE_2D, texture3);
//glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_CLAMP);
// glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_CLAMP);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_LINEAR);
glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR);
// glTexParameteri(GL_TEXTURE_2D,GL_TEXTURE_MIN_FILTER,GL_LINEAR_MIPMAP_NEAREST);
glTexImage2D(GL_TEXTURE_2D, 0, 3, 64, 64, 0, GL_RGB, GL_UNSIGNED_BYTE, bits);

glEnable(GL_TEXTURE_2D);
//glTexEnvf(GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_DECAL); */

}
//--------------------------------------------------------------------------
// ASTRODYNAMICS
//--------------------------------------------------------------------------

void __fastcall TFormMain:: posvel(double T, double T0, double EL[9], double R[4], double PR[4])
{
double GM=398600.64; //km^3*s^-2
double AR,VR,VN,P,EA;
double AKE,A,EC,AI,AM,V,U,SOM,COM,SI,CI,SU,CU; //ED;
AKE=sqrt(GM);
A=EL[0];
EC=EL[1];
AI=EL[2];
// A=20000,0;
EL[8]=AKE/(pow(A,1.5));//mean motion
AM=EL[8]*(T-T0)+EL[5];
// printf("AM %.15f \n",AM);
//ED=AM;
EA=KEPL(AM,EC);

EL[6]=EA;
V=2.0*atan((sqrt((1.0+EC)/(1.0-EC)))*tan(EA/2.0));
if(V<0) V=V+2*PI;
EL[7]=V;
U=V+EL[3];
AR=A*(1.0-EC*EC)/(1.0+EC*cos(V));
SOM=sin(EL[4]);
COM=cos(EL[4]);
SU=sin(U);
CU=cos(U);
SI=sin(AI);
CI=cos(AI);

R[0]=AR*(CU*COM-SU*SOM*CI);
R[1]=AR*(CU*SOM+SU*COM*CI);
R[2]=AR*SU*SI;
R[3]=AR;
P=A*(1.0-EC*EC);
VR=AKE*EC*sin(V)/sqrt(P);
VN=AKE*(1.0+EC*cos(V))/(sqrt(P));

PR[0]=R[0]*VR/AR+(-SU*COM-CU*SOM*CI)*VN;
PR[1]=R[1]*VR/AR+(-SU*SOM+CU*COM*CI)*VN;
PR[2]=R[2]*VR/AR+CU*SI*VN;
PR[3]=sqrt(PR[0]*PR[0]+PR[1]*PR[1]+PR[2]*PR[2]);
// fprintf(test, "%.10f %f.10\n",V*180/M_PI, PR[3]);
}
//--------------------------------------------------------------------------
double KEPL(double AM, double EC)
{
double E1,ED;
E1=AM;

L4: ED=AM+EC*sin(E1);
if(sqrt(pow(ED-E1,2))>0.00000000000001)
{
E1=ED;
goto L4;
}

return ED;
}
//---------------------------------------------------------------------------
double cross(double vec1[3],double vec2[3], double vec4[3])
{
vec4[0]=vec1[1]*vec2[2]-vec1[2]*vec2[1];
vec4[1]=vec1[2]*vec2[0]-vec1[0]*vec2[2];
vec4[2]=vec1[0]*vec2[1]-vec1[1]*vec2[0];
return 0;
}
//---------------------------------------------------------------------------
void __fastcall TFormMain:: elements(double vel[4], double pos[4], double el[9], double T0, double T)
{
double h[3],dlh,dlv,dlr,en, node, omega,EA,M;
double PI=atan(1.0)*4.0;
double GM=398600.64; //km^3*s^-2
double a,ec,ink,b,c,u,p,v;
cross(pos,vel,h);
dlh=lenght(h);
dlv=lenght(vel);
dlr=lenght(pos);
en=pow(dlv,2)/2.0-GM/dlr;
a=-GM/(2.0*en);
el[0]=a;

if(dlh*dlh/(GM*a)>=1.0) ec=el[1];
else ec=sqrt(1.0-dlh*dlh/(GM*a));

// ec=sqrt(1.0-dlh*dlh/(GM*a));
el[1]=ec;
ink=acos(h[2]/dlh);
if (ink<0.0) ink=ink+PI;
el[2]=ink;
if(ink==0.0000) node=el[4];
else
node=atan2(h[0],-h[1]);
if(node<0.0) node=node+2*PI;
el[4]=node;
//if (ink==0.0) ink=0.0001;
//else
b=pos[2]/(sin(ink));
c=pos[0]*cos(node)+pos[1]*sin(node);
u=atan2(b,c);
if(ink==0.0000) u=0.0;
if(u<0.0) u=u+2.0*PI;
p=a*(1.0-pow(ec,2));
b=sqrt(p/GM)*scalar(vel,pos);
c=p-dlr;
v=atan2(b,c);
if(v<0.0) v=v+2*PI;
omega=u-v;
if(ink==0.0000) omega=el[3];
if(omega<0.0) omega=omega+2*PI;
el[3]=omega;
el[8]=sqrt(GM)/(pow(a,1.5));//mean motion
EA=2.0*atan(sqrt((1-ec)/(1+ec))*tan(v/2.0));
M=EA-ec*sin(EA);
//M=el[5]+el[8]*(T-T0);
if(M<0.0) M=M+2*PI;
el[5]=M;
}
//---------------------------------------------------------------------------
double lenght(double vec[3])
{
double lenght;
lenght=sqrt(pow(vec[0],2)+pow(vec[1],2)+pow(vec[2],2));
return lenght;
}
//---------------------------------------------------------------------------
double scalar(double vec1[3],double vec2[3])
{
double scalar;
scalar=vec1[0]*vec2[0]+vec1[1]*vec2[1]+vec1[2]*vec2[2];
return scalar;
}
//---------------------------------------------------------------------------
double *secular(double el[9], double wiek[3])
{
double Re=6378.137,GM=398600.64;//km
double J2=-0.00108262668355;
double nm, ink, p, ec,aa, ksi;
double omegas,Ms,nodes;

ink=el[2];
ec=el[1];
aa=el[0];
p=aa*(1.0-pow(ec,2));
nm=sqrt(GM/(pow(aa,3)));

ksi=3*Re*Re*J2/(4*p*p);
//printf("ksi %.15f\n",ksi);
omegas=ksi*nm*(4-5*sin(ink)*sin(ink));
nodes=-2*ksi*nm*cos(ink);
Ms=ksi*nm*sqrt(1-ec*ec)*(2-3*sin(ink)*sin(ink));
wiek[0]=omegas;
wiek[1]=nodes;
wiek[2]=Ms;

return wiek;
// printf("secular %.15f %.15f %.15f\n",omegas, nodes, Ms);
}
//---------------------------------------------------------------------------
double earthfix(double r[4], double rf[4], double t)
{
//double PI=atan(1.0)*4.0;
double teta0ecf=0,tetaecf=0;
teta0ecf=green;
//double teta=1.00273790926493*2.0*PI;
tetaecf=teta0ecf+7.29211585530*pow(10,-5)*t;
rf[0]=cos(tetaecf)*r[0]+sin(tetaecf)*r[1];
rf[1]=-sin(tetaecf)*r[0]+cos(tetaecf)*r[1];
rf[2]=r[2];
return 0;
}
//----------------------------------------------------------------------------
void __fastcall TFormMain::ECIframe1Click(TObject *Sender)
{
object=0;
}
//---------------------------------------------------------------------------

void __fastcall TFormMain::ECIorbits1Click(TObject *Sender)
{
object=1;
}
//---------------------------------------------------------------------------

void __fastcall TFormMain::ECFFrame1Click(TObject *Sender)
{
object=2;
}
//---------------------------------------------------------------------------

void __fastcall TFormMain::ECFandTrajectory1Click(TObject *Sender)
{
object=3;
}
//---------------------------------------------------------------------------

void __fastcall TFormMain::GroundTrack1Click(TObject *Sender)
{
object=4;
}
//---------------------------------------------------------------------------

void __fastcall TFormMain::CreateNewConstelatio1Click(TObject *Sender)
{

Child->Width=350;
Child->Height=600;
Child->Visible=true;

Child->Edit7->Text=IntToStr(num);
Child->Edit1->Text=FloatToStr(sat[0].el[0]);
Child->Edit2->Text=FloatToStr(sat[0].el[1]);
Child->Edit3->Text=FloatToStr(sat[0].el[2]);
Child->Edit4->Text=FloatToStr(sat[0].el[3]);
Child->Edit5->Text=FloatToStr(sat[0].el[4]);
Child->Edit6->Text=FloatToStr(sat[0].el[5]);


}
//---------------------------------------------------------------------------
void __fastcall TFormMain::UpDate()
{

int i;
// double deltaOMEGA,
double deltaM;
int max;

if(Child->CheckBox5->Checked)
{
Child->CheckBox6->Enabled=FALSE;
}
else if(Child->CheckBox6->Checked)
{
Child->CheckBox5->Enabled=FALSE;
}
else if(!(Child->CheckBox5->Checked) && !(Child->CheckBox6->Checked))
{
Child->CheckBox6->Enabled=TRUE;
Child->CheckBox5->Enabled=TRUE;
}

if(Child->CheckBox1->Checked)
{
Child->CheckBox2->Enabled=FALSE;
Child->CheckBox3->Enabled=FALSE;
Child->CheckBox4->Enabled=FALSE;
Child->Edit12->Enabled=FALSE;
}

else if(Child->CheckBox2->Checked)
{
Child->Edit12->Enabled=TRUE;
Child->CheckBox1->Enabled=FALSE;
Child->CheckBox3->Enabled=FALSE;
Child->CheckBox4->Enabled=FALSE;
}

else if(Child->CheckBox3->Checked)
{
Child->CheckBox2->Enabled=FALSE;
Child->CheckBox1->Enabled=FALSE;
Child->CheckBox4->Enabled=FALSE;
Child->Edit12->Enabled=FALSE;
}

else if(Child->CheckBox4->Checked)
{
Child->CheckBox2->Enabled=FALSE;
Child->CheckBox3->Enabled=FALSE;
Child->CheckBox1->Enabled=FALSE;
Child->Edit12->Enabled=FALSE;
}
else if(
!(Child->CheckBox4->Checked) &&
!(Child->CheckBox3->Checked) &&
!(Child->CheckBox2->Checked) &&
!(Child->CheckBox1->Checked))
{
Child->CheckBox2->Enabled=TRUE;
Child->CheckBox3->Enabled=TRUE;
Child->CheckBox1->Enabled=TRUE;
Child->CheckBox4->Enabled=TRUE;
Child->Edit12->Enabled=TRUE;
}


if(Child->nowa)
{ T=0.0; T0=0.0;
yr=Child->Edit13->Text.ToDouble();
mo=Child->Edit14->Text.ToDouble();
d=Child->Edit15->Text.ToDouble();
h=Child->Edit16->Text.ToDouble();
m=Child->Edit17->Text.ToDouble();
s=Child->Edit18->Text.ToDouble();
JD=julian(yr, mo, d, h, m, s);
green=sidereal(JD,h);

if(Child->CheckBox5->Checked)
{
FormMain->semi();
}

num=Child->Edit7->Text.ToInt();
sat[0].el[0]=Child->Edit1->Text.ToDouble();
sat[0].el[1]=Child->Edit2->Text.ToDouble();
if(sat[0].el[1]<=0) sat[0].el[1]=0,0001;
sat[0].el[2]=Child->Edit3->Text.ToDouble();
sat[0].el[3]=Child->Edit4->Text.ToDouble();
sat[0].el[4]=Child->Edit5->Text.ToDouble();
sat[0].el[5]=Child->Edit6->Text.ToDouble();

for(i=0;i<num;i++)
{
sat.r=rand()%4;
sat.g=rand()%2;
sat.b=rand()%3;
}



for(i=1;i<num;i++)
{
sat.el[5]=i*20*PI/180+sat[0].el[5];//loop*10*PI/180.0;//mean anomaly
sat.el[2]=sat[0].el[2];//60*PI/180.0;//ink
sat.el[4]=sat[0].el[4]+i*25*PI/180.0;//loop*15*PI/180.0;//node
sat.el[0]=sat[0].el[0];//20000.0;//semimajor axis
sat.el[1]=sat[0].el[1];//0.6; //eccentricity
sat.el[3]=sat[0].el[3];//15*PI/180.0; // argument of perigee
}

if(Child->CheckBox1->Checked && Child->CheckBox1->Enabled)
{
Fn=Child->Edit8->Text.ToInt();
Fd=Child->Edit9->Text.ToInt();
Np=Child->Edit10->Text.ToInt();
Nd=Child->Edit11->Text.ToInt();
num=Child->Edit7->Text.ToInt();
sat[0].el[0]=Child->Edit1->Text.ToDouble();
sat[0].el[1]=Child->Edit2->Text.ToDouble();


if(sat[0].el[1]<=0) sat[0].el[1]=0,0001;
sat[0].el[2]=Child->Edit3->Text.ToDouble();
//deltaOMEGA=-2.0*PI*Fn/Fd;
max=Fd*Nd;
//if(max>1000) max=1000;
if(num>max) num=max;
Child->Edit7->Text=IntToStr(num);

for(i=1;i<num;i++)
{ sat.el[0]=Child->Edit1->Text.ToDouble();
sat.el[1]=Child->Edit2->Text.ToDouble();
sat.el[2]=Child->Edit3->Text.ToDouble();
sat.el[4]= sat[i-1].el[4]-2*PI*Fn/Fd;
sat.el[5]= sat[i-1].el[5]+2*PI*Fn/Fd*Np/Nd;
}
Child->nowa=FALSE;
}
else if(Child->CheckBox2->Checked && Child->CheckBox2->Enabled)
{

dOMEGA=Child->Edit12->Text.ToDouble();

Fn=Child->Edit8->Text.ToInt();
Fd=Child->Edit9->Text.ToInt();
Np=Child->Edit10->Text.ToInt();
Nd=Child->Edit11->Text.ToInt();
num=Fd+1;
if(num>1000) num=1000;
Child->Edit7->Text=IntToStr(num);
for(i=1;i<num;i++)
{
sat.el[0]=Child->Edit1->Text.ToDouble();
sat.el[1]=Child->Edit2->Text.ToDouble();
sat.el[2]=Child->Edit3->Text.ToDouble();
sat.el[4]= sat[i-1].el[4]-dOMEGA*Fn/Fd;
//if(sat.el[4]>sat[0].el[4]+dOMEGA || sat.el[4]<-dOMEGA+sat[0].el[4])
// sat.el[4]=sat[0].el[4];
sat.el[5]= sat[i-1].el[5]+dOMEGA*Fn/Fd*Np/Nd;
}
Child->nowa=FALSE;
}
else if(Child->CheckBox3->Checked && Child->CheckBox3->Enabled)
{
NonSym->Visible=TRUE;

}
if(Child->CheckBox4->Checked && Child->CheckBox4->Enabled)
{
incom->Visible=TRUE;
}
Child->nowa=FALSE;
}

if(NonSym->zrobione)
{
double startOM;
double deltaOM[10];
Fn=Child->Edit8->Text.ToInt();
Fd=Child->Edit9->Text.ToInt();
Np=Child->Edit10->Text.ToInt();
Nd=Child->Edit11->Text.ToInt();
startOM=NonSym->Edit1->Text.ToDouble();
deltaOM[0]=NonSym->Edit2->Text.ToDouble();
deltaOM[1]=NonSym->Edit3->Text.ToDouble();
deltaOM[2]=NonSym->Edit4->Text.ToDouble();
deltaOM[3]=NonSym->Edit5->Text.ToDouble();
deltaOM[4]=NonSym->Edit6->Text.ToDouble();
deltaOM[5]=NonSym->Edit7->Text.ToDouble();
deltaOM[6]=NonSym->Edit8->Text.ToDouble();
deltaOM[7]=NonSym->Edit9->Text.ToDouble();
deltaOM[8]=NonSym->Edit10->Text.ToDouble();

if(num>10) num=10;
Child->Edit7->Text=IntToStr(num);

sat[0].el[5]==-1.0*startOM*Nd/Np;
sat[0].el[4]=startOM;

for(loop=1;loop<num;loop++)
{
sat[loop].el[4]=sat[loop-1].el[4]-deltaOM[loop-1];
deltaM=-1.0*deltaOM[loop-1]*Nd/Np;
sat[loop].el[5]= sat[loop-1].el[5]+deltaM;
}
T=0.0; T0=0.0;
NonSym->Visible=FALSE;
}

if(incom->yes)
{
int j,i,k=0,x=0;
int rem=0;
int which=0;
rem=0;
rem=incom->Edit1->Text.ToInt();

for(j=0;j<100;j++)
incom->remove[j]=0;

for(j=0;j<rem;j++)
incom->remove[j]=incom->Memo1->Lines->Strings[j].ToInt();

while(k<rem)
{
for(j=0;j<rem;j++)
{
if(incom->remove[j]>x)
{
x=incom->remove[j];
which=j;
}

}

if(x>-1)
{
sat[x].x=sat[num-1].x;
sat[x].y=sat[num-1].y;
sat[x].z=sat[num-1].z;

for(i=0;i<9;i++)
sat[x].el=sat[num-1].el;
x=0;
}
k++;
incom->remove[which]=-1;
num--;


}
Child->Edit7->Text=IntToStr(num);
incom->Memo1->Clear();
incom->yes=FALSE;
}





Child->nowa=FALSE;
NonSym->zrobione=FALSE;
}
void __fastcall TFormMain::Exit1Click(TObject *Sender)
{
Application->Terminate();
}
//---------------------------------------------------------------------------


void __fastcall TFormMain::SaveConstellationas1Click(TObject *Sender)
{ int one, two, three, flower;
if(Child->CheckBox5->Checked) flower=1; else flower=0;
if(Child->CheckBox1->Checked) one=1; else one=0;
if(Child->CheckBox2->Checked) two=1; else two=0;
if(Child->CheckBox3->Checked) three=1; else three=0;

if(!SaveDialog1->Execute()) return;
String Name=SaveDialog1->FileName;

if(SaveDialog1->FileName!="")
{
ofstream Students(SaveDialog1->FileName.c_str(), ios::out);
Students <<flower<<"\n"<<num<<"\n"<<(FloatToStr(sat[0].el[0])).c_str()<<"\n"
<<(FloatToStr(sat[0].el[1])).c_str()<<"\n"
<<(FloatToStr(sat[0].el[2])).c_str()<<"\n"
<<(FloatToStr(sat[0].el[3])).c_str()<<"\n"
<<(FloatToStr(sat[0].el[4])).c_str()<<"\n"
<<(FloatToStr(sat[0].el[5])).c_str()<<"\n"
<<one<<"\n"<<two<<"\n"<<three<<"\n"
<<Fn<<"\n"<<Fd<<"\n"<<Np<<"\n"<<Nd<<"\n"
<<(FloatToStr(dOMEGA)).c_str();
}
}
//---------------------------------------------------------------------------

void __fastcall TFormMain::OpenConstellation1Click(TObject *Sender)
{
int k=0;
if(!OpenDialog1->Execute()) return;
String Name=OpenDialog1->FileName;

std::ifstream is(Name.c_str()); if(!is.is_open()) return;
std::string s;

while(!is.eof())
{

std::getline(is,s);
String S=s.c_str(); // convert back to Borland AnsiString
if(k==0) //check if it is flower constellation
{
if(S.ToInt()==1) //flower constellation?//yes
{
Child->CheckBox5->Checked=true;
Child->CheckBox6->Checked=false;
}
else //no
{
Child->CheckBox6->Checked=true;
Child->CheckBox5->Checked=false;

}
}

if(k==1) //load number of satellites
{
num=S.ToInt();
Child->Edit7->Text=S;
}
if(Child->CheckBox5->Checked)
{
if(k==2) //load semi-major axis
{
sat[0].el[0]=S.ToDouble();
Child->Edit1->Text=S;
}

if(k==3) //load eccentricity
{
sat[0].el[1]=S.ToDouble();
Child->Edit2->Text=S;
}

if(k==4) //load inclination
{
sat[0].el[2]=S.ToDouble();
Child->Edit3->Text=S;
}

if(k==5)
{
sat[0].el[3]=S.ToDouble();
Child->Edit4->Text=S;
}

if(k==6)
{
sat[0].el[4]=S.ToDouble();
Child->Edit5->Text=S;
}

if(k==7)
{
sat[0].el[5]=S.ToDouble();
Child->Edit6->Text=S;
}

if(k==8)
{
if(S.ToInt()==1)
{
Child->CheckBox1->Checked=true;
Child->CheckBox1->Enabled=true;
Child->CheckBox2->Checked=false;
Child->CheckBox3->Checked=false;
Child->CheckBox4->Checked=false;
}
}

if(k==9)
{
if(S.ToInt()==1)
{
Child->CheckBox2->Checked=true;
Child->CheckBox2->Enabled=true;
Child->CheckBox1->Checked=false;
Child->CheckBox3->Checked=false;
Child->CheckBox4->Checked=false;
}
}

if(k==10)
{
if(S.ToInt()==1)
{
Child->CheckBox3->Checked=true;
Child->CheckBox3->Enabled=true;
Child->CheckBox1->Checked=false;
Child->CheckBox2->Checked=false;
Child->CheckBox4->Checked=false;
}
}

if(k==11)
{
Fn=S.ToInt();
Child->Edit8->Text=S;
}

if(k==12)
{
Fd=S.ToInt();
Child->Edit9->Text=S;
}

if(k==13)
{
Np=S.ToInt();
Child->Edit10->Text=S;
}

if(k==14)
{
Nd=S.ToInt();
Child->Edit11->Text=S;
}

if(k==15)
{
dOMEGA=S.ToDouble();
Child->Edit12->Text=S;
}
}
if(Child->CheckBox6->Checked)
{
//int l=0;
if(k>15 && k<=15+num)
{
sat[k-15].el[0]=S.ToDouble();
Child->Edit1->Text=S;
}

if(k>15+num && k<=15+2*num)
{
sat[k-num-15].el[1]=S.ToDouble();
Child->Edit2->Text=S;
}

if(k>15+2*num && k<=15+3*num)
{
sat[k-2*num-15].el[2]=S.ToDouble();
Child->Edit3->Text=S;
}

if(k>15+3*num && k<=15+4*num)
{
sat[k-3*num-15].el[3]=S.ToDouble();
Child->Edit4->Text=S;
}

if(k>15+4*num && k<=15+5*num)
{
sat[k-4*num-15].el[4]=S.ToDouble();
Child->Edit5->Text=S;
}

if(k>15+5*num && k<=15+6*num)
{
sat[k-5*num-15].el[5]=S.ToDouble();
Child->Edit6->Text=S;
}
Child->nowa=false;
}
k++;

object=0;
}

if(Child->CheckBox5->Checked)
Child->nowa=true;
}
//---------------------------------------------------------------------------
void __fastcall TFormMain::semi(void)
{

double rotz; // earth rotation
double ksi; //
double J2; // J2
double Rz; // eart radius
double GM; // gravitational constant
double Nd;
int Np, i;
double a, a0 ,p, hp, ink, n, ec, roznica;
a0=Child->Edit1->Text.ToDouble();
ec=Child->Edit2->Text.ToDouble();
ink=Child->Edit3->Text.ToDouble();
Np=Child->Edit10->Text.ToInt();
Nd=Child->Edit11->Text.ToInt();
//constants
rotz=7.29211585530*pow(10.0,-5.0);
J2=-1.08262668355*pow(10.0,-3.0);
Rz=6378.137;
GM=398600.0;

hp=a0-Rz;
roznica=1.0;

while(sqrt(pow(roznica,2))>0.0000000000000001)
{
hp=a0*(1-ec);
ec=-1*hp/a0+1;
p=a0*(1.0-ec*ec);
n=sqrt(GM/pow(a0,3.0));
ksi=3.0*Rz*Rz*J2/(4.0*p*p);
a=pow(GM,1.0/3.0)*pow((pow((1.0/rotz*Nd/Np*pow((1.0+2.0*ksi*n/rotz*cos(ink)),-1.0)*(1.0+ksi*(4.0+2.0*sqrt(1.0-ec*ec)-(5.0+3.0*sqrt(1.0-ec*ec))*sin(ink)*sin(ink)))),2.0)),1.0/3.0);
roznica=a-a0;
a0=a;
}

Child->Edit1->Text=FloatToStr(a);
sat[0].el[0]=a;
Child->repeat=2*PI/rotz*Nd/Np*pow((1.0+2.0*ksi*n/rotz*cos(ink)),-1.0)*(1.0+ksi*(4.0+2.0*sqrt(1.0-ec*ec)-(5.0+3.0*sqrt(1.0-ec*ec))*sin(ink)*sin(ink)));

}

void groundtrack(double r[4], double tetaGMST, double phi, double lambda)
{
double rsat,a,b;
double Re=6378.137;
double EE,alpha,F,P,Q,D,ni,G,t;
double ec=0.081819221456;
rsat=sqrt(r[0]*r[0]+r[1]*r[1]);

a=Re;
if(r[2]<0)
b=-6356.75160056;
else if(r[2]>0)
b=6356.75160056;
else b=0;

EE=(b*r[2]-(a*a-b*b))/(a*rsat);
alpha=atan2(r[1]/rsat,r[0]/rsat);
lambda=alpha-tetaGMST;

while(lambda<0.0)
lambda=lambda+2.0*PI;
// if(lambda<0.0) goto G1;

while(lambda>2.0*PI)
lambda=lambda-2.0*PI;
//if(lambda>2.0*PI)
//goto G2;


F=(b*r[2]+(a*a-b*b))/(a*rsat);
P=4.0*(EE*F+1)/3.0;
Q=2.0*(EE*EE-F*F);
D=pow(P,3.0)+Q*Q;

if(D>0)
ni=pow((sqrt(D)-Q),(1.0/3.0))-pow((sqrt(D)+Q),(1.0/3.0));
if(D<0)
ni=2*sqrt(-P)*cos(1.0/3.0*pow(cos(Q/(P*sqrt(-P))),-1.0));
else ni=0;
G=0.5*(sqrt(EE*EE+ni)+EE);
t=sqrt(G*G+(F-ni*G)/(2.0*G-EE))-G;
phi=atan(a*(1.0-t*t)/(2.0*b*t));
phi=atan((1-ec*ec)*tan(phi)); /*
G3: if(phi<-PI/2.0)
phi=phi+PI;
if(phi<-PI/2.0) goto G3;

G4: if(phi>PI/2.0)
phi=phi-PI;
if(phi>PI/2.0) goto G4;
*/

//fprintf(test,"phi=%f lambda=%f\n",phi*180.0/PI,lambda*180.0/PI);
}
//----------------------------------------------------------------------------
double __fastcall TFormMain:: julian(double yr, double mo, double d, double h, double min, double s)
{
double juliand;
double one;

one=7.0*(yr+int((mo+9.0)/12.0))/4.0;
juliand=367.0*yr-int(one)+int(275*mo/9.0)+d+1721013.5+((s/60.0+min)/(60.0)+h)/(24.0);
return juliand;

}

//-----------------------------------------------------------------------------
bool __fastcall TFormMain:: shadow(double rsat[3], double JD)
{
double TUT1=0;

TUT1=(JD-2451545.0)/(36525.0);
double TTDB, Lsun, Msun, Lecl;
double rs, epsilon, rsun[3];

Lsun=280.460+36000.770*TUT1;
while(Lsun<360.0)
Lsun=Lsun+360.0;

while(Lsun>360)
Lsun=Lsun-360;

TTDB=TUT1;
Msun=357.5277233+35999.05034*TTDB;

while(Msun<360.0)
Msun=Msun+360.0;

while(Msun>360)
Msun=Msun-360;

Lecl=Lsun+1.914666471*sin(Msun*M_PI/180.0)+0.019994643*sin(2.0*Msun*M_PI/180.0);
rs=1.000140612-0.016708617*cos(Msun*M_PI/180.0)-0.000139589*cos(2.0*Msun*M_PI/180.0);
epsilon=23.439291-0.0130042*TTDB;

double AU=149597870.691;
rsun[0]=rs*cos(Lecl*M_PI/180)*AU;
rsun[1]=rs*cos(epsilon*M_PI/180)*sin(Lecl*M_PI/180)*AU;
rsun[2]=rs*sin(epsilon*M_PI/180)*sin(Lecl*M_PI/180)*AU;

double deltas, alphas, a1,a2;

deltas=asin(sin(epsilon*M_PI/180.0)*sin(Lecl*M_PI/180.0));

a1=cos(epsilon*M_PI/180.0)*sin(Lecl*M_PI/180.0)/(cos(deltas));
a2=cos(Lecl*M_PI/180.0)/(cos(deltas));

alphas=atan2(a1,a2);

double z;

double Rz=6378.137;
double ksi;
double horiz, vert;
double pen, alphapen, umb;
bool umbra=false, preumbra=false;
alphapen=0.26900424*M_PI/180.0;
double alphaumb=0.26411888*M_PI/180.0;
z=scalar(rsun,rsat);

if(z<0.0)
{
ksi=acos(z/(lenght(rsun)*lenght(rsat)));
horiz=lenght(rsat)*cos(ksi);
vert=lenght(rsat)*sin(ksi);
pen=Rz+tan(alphapen)*horiz;
if(vert<=pen)
{
preumbra=true;
umb=Rz-tan(alphaumb)*horiz;
if(vert<=umb)
{
umbra=true;
}

}
}
if(preumbra || umbra)
return true;
else
return false;
}
//----------------------------------------------------------------------------
double sidereal(double JD, double HG)
{
double JD0,D,D0,GST,TG;
//julian date of previous midnight
if((JD-int(JD))>0.5)
JD0=int(JD)+0.5;
else
JD0=int(JD)-0.5;

//number of days and fraction from 2000 January
D=JD-2451545.0;
D0=JD0-2451545.0;

TG = D/36525;

//greenwich sidereal time in hours:
GST = 6.697374558+0.06570982441908*D0+1.00273790935*HG+0.000026*TG*TG;
while(GST>2*M_PI)
GST=GST-2*M_PI;
while(GST<2*M_PI)
GST=GST+2*M_PI;

return GST;
}


//------------------------------------------------------------------------
bool checkshadow(double rsat[3], double JD, double T)
{ bool cien,JDn;
double dsh,hsh,minsh,ssh;
double dni, godz, minuty, sekundy;

dni=floor(T/(24.0*3600.0));
godz=floor((T-dni*(24.0*3600.0))/3600.0);
minuty=floor(((T-dni*(24.0*3600.0))-godz*3600.0)/60.0);
sekundy=floor((((T-dni*(24.0*3600.0))-godz*3600.0)-minuty*60.0));

dsh=d+dni;
hsh=h+godz;
minsh=m+minuty;
ssh=s+sekundy;

JDn=FormMain->julian(yr,mo,dsh,hsh,minsh,ssh);
cien=FormMain->shadow(rsat,JDn);
return cien;

}

//----------------------------------------------------------------------

void cieniuj(double JD, double T)
{
double polsat[3];
int p;
for(p=0;p<num;p++)
{
polsat[0]=sat

.x;
polsat[1]=sat

.y;
polsat[2]=sat

.z;

if(checkshadow(polsat,JD,T))
{ /*
sat

.r=0.2;
sat

.g=0.2;
sat

.b=0.2;*/
}
else
{
sat

.r=sat

.rs;
sat

.g=sat

.gs;
sat

.b=sat

.bs;
}
}
}
void __fastcall TFormMain::Perturbation1Click(TObject *Sender)
{
pert->Visible=true;
}
//---------------------------------------------------------------------------



Share this post


Link to post
Share on other sites
Sign in to follow this  

  • Advertisement
×

Important Information

By using GameDev.net, you agree to our community Guidelines, Terms of Use, and Privacy Policy.

We are the game development community.

Whether you are an indie, hobbyist, AAA developer, or just trying to learn, GameDev.net is the place for you to learn, share, and connect with the games industry. Learn more About Us or sign up!

Sign me up!