在定量遥感中,通常需要计算地方时和太阳高度角,本文采用C++语言实现。
#include <cmath>
#include <iostream>
#include <fstream>
using namespace std;
void main()
{
int JD,NF,Y,R,s[5],F[5];
float JF,WD;
float N0;
ifstream data1("d:\\result\\data-input.txt");
ofstream out1("d:\\result\\6gaodujiao.txt");
data1>>JD>>JF;
data1>>WD;
data1>>NF>>Y>>R;
out1<<JD<<" "<<JF<<endl;
out1<<WD<<endl;
out1<<NF<<" "<<Y<<" "<<R<<endl;
out1<<"太阳方位角在正南为0,由南向西为0——+180,由南向东为0——-180。午前取负值,午后取正值。"<<endl;
N0=79.6764+0.2422*(NF-1985)-int((NF-1985)/4);
float A,B,G,L;
A=NF/4;
double C;
C=32.8;
B=A-int(A);
if(Y<=2)
C=30.6;
if(B=0 && Y>2 )
C=31.8;
G=int(30.6*Y-C+0.5)+R;
L=(JD+JF/60)/15;
float H[5], N[5], theta[5],Er[5],Ed[5],Et[5];
double Sd[5];
float St[5], shj