protected override void OnPaint( PaintEventArgs e )
{
base.OnPaint( e );
if ( m_Slice != null )
e.Graphics.DrawImage( m_Slice, 0, 0, Width, Height );
if ( ShowIsoLines )
{
PointF P0 = new PointF(), P1 = new PointF();
double ThetaHalf, PhiHalf, ThetaDiff, PhiDiff;
double ThetaHalfN, PhiHalfN, ThetaDiffN, PhiDiffN;
for ( int IsolineIndex=0; IsolineIndex < 20; IsolineIndex++ )
{
double ThetaI = (1+IsolineIndex) / 20.0 * 0.495 * Math.PI;
// double ThetaI = Math.Sqrt( (1+IsolineIndex) / 12.0 ) * 0.49 * Math.PI;
for ( int i=0; i < 20; i++ )
{
ThetaI = (1.0+IsolineIndex) / 20.0 * 0.495 * Math.PI;
double HalfPhi = 0.5 * i * Math.PI / 20;
DisplayForm.std_coords_to_half_diff_coords( ThetaI, HalfPhi, ThetaI, -HalfPhi, out ThetaHalf, out PhiHalf, out ThetaDiff, out PhiDiff );
//Convert( ThetaI, Phi, out ThetaHalf, out ThetaDiff );
// Change Phi_d
double TestThetaI, TestPhiI, TestThetaO, TestPhiO;
DisplayForm.half_diff_coords_to_std_coords( ThetaHalf, PhiHalf, ThetaDiff, PhiD, out TestThetaI, out TestPhiI, out TestThetaO, out TestPhiO );
DisplayForm.std_coords_to_half_diff_coords( TestThetaI, TestPhiI, TestThetaO, TestPhiO, out ThetaHalf, out PhiHalf, out ThetaDiff, out PhiDiff );
ThetaI = (1.2+IsolineIndex) / 20.0 * 0.495 * Math.PI;
HalfPhi = 0.5 * (i+1) * Math.PI / 20;
DisplayForm.std_coords_to_half_diff_coords( ThetaI, HalfPhi, ThetaI, -HalfPhi, out ThetaHalfN, out PhiHalfN, out ThetaDiffN, out PhiDiffN );
//Convert( ThetaI, Phi, out ThetaHalfN, out ThetaDiffN );
// Change Phi_d
DisplayForm.half_diff_coords_to_std_coords( ThetaHalfN, PhiHalfN, ThetaDiffN, PhiD, out TestThetaI, out TestPhiI, out TestThetaO, out TestPhiO );
DisplayForm.std_coords_to_half_diff_coords( TestThetaI, TestPhiI, TestThetaO, TestPhiO, out ThetaHalfN, out PhiHalfN, out ThetaDiffN, out PhiDiffN );
P0.X = (float) (Width * ThetaHalf * 0.63661977236758134307553505349006); // divided by PI/2
P0.Y = (float) (Height * (1.0f - ThetaDiff * 0.63661977236758134307553505349006)); // divided by PI/2
P1.X = (float) (Width * ThetaHalfN * 0.63661977236758134307553505349006); // divided by PI/2
P1.Y = (float) (Height * (1.0f - ThetaDiffN * 0.63661977236758134307553505349006)); // divided by PI/2
e.Graphics.DrawLine( m_Pen, P0, P1 );
}
//Convert( ThetaI, 0.8 * Math.PI, out ThetaHalf, out ThetaDiff );
{
double ThetaDH = Math.Acos( Math.Sqrt( Math.Cos( ThetaI ) ) );
double PhiI = 2.0 * Math.Atan( 1.0 / Math.Cos( ThetaDH ) );
// DisplayForm.std_coords_to_half_diff_coords( ThetaI, 0, ThetaI, PhiI, out ThetaHalf, out PhiHalf, out ThetaDiff, out PhiDiff );
Convert( ThetaI, PhiI, out ThetaHalf, out ThetaDiff );
}
P0.X = (float) (Width * ThetaHalf * 0.63661977236758134307553505349006) - 16.0f; // divided by PI/2
P0.Y = (float) (Height * (1.0f - ThetaDiff * 0.63661977236758134307553505349006)) - 0*16.0f; // divided by PI/2
e.Graphics.DrawString( "" + (int) (ThetaI * 180 / Math.PI) + "°", Font, Brushes.Black, P0 );
}
}
if ( ShowIsoLines )
{
PointF P0 = new PointF(), P1 = new PointF();
double ThetaHalf, PhiHalf, ThetaDiff, PhiDiff;
double ThetaHalfN, PhiHalfN, ThetaDiffN, PhiDiffN;
for ( int IsolineIndex=0; IsolineIndex < 19; IsolineIndex++ )
{
double Phi = (1.0+IsolineIndex) / 20.0 * Math.PI;
// double Phi = Math.Sqrt( (0.5+IsolineIndex) / 12.0 ) * Math.PI;
for ( int i=0; i < 40; i++ )
{
double Angle = i * 0.50 * Math.PI / 40;
DisplayForm.std_coords_to_half_diff_coords( Angle, 0, Angle, Phi, out ThetaHalf, out PhiHalf, out ThetaDiff, out PhiDiff );
Convert( Angle, Phi, out ThetaHalf, out ThetaDiff );
Angle = (i+1) * 0.50 * Math.PI / 40;
DisplayForm.std_coords_to_half_diff_coords( Angle, 0, Angle, Phi, out ThetaHalfN, out PhiHalfN, out ThetaDiffN, out PhiDiffN );
Convert( Angle, Phi, out ThetaHalfN, out ThetaDiffN );
P0.X = (float) (Width * ThetaHalf * 0.63661977236758134307553505349006); // divided by PI/2
P0.Y = (float) (Height * (1.0f - ThetaDiff * 0.63661977236758134307553505349006)); // divided by PI/2
P1.X = (float) (Width * ThetaHalfN * 0.63661977236758134307553505349006); // divided by PI/2
P1.Y = (float) (Height * (1.0f - ThetaDiffN * 0.63661977236758134307553505349006)); // divided by PI/2
e.Graphics.DrawLine( m_Pen, P0, P1 );
}
DisplayForm.std_coords_to_half_diff_coords( 0.5*Math.PI, 0, 0.5*Math.PI, Phi, out ThetaHalf, out PhiHalf, out ThetaDiff, out PhiDiff );
P0.X = (float) (Width * ThetaHalf * 0.63661977236758134307553505349006) - 30.0f; // divided by PI/2
P0.Y = (float) (Height * (1.0f - ThetaDiff * 0.63661977236758134307553505349006)) + 2.0f; // divided by PI/2
e.Graphics.DrawString( "" + (int) (Phi * 180 / Math.PI) + "°", Font, Brushes.Black, P0 );
}
}
}