/*?本程序為日本Himawari-8衛星hsd格式數據經投影轉換為LDF格式數據，并可用swap打開
?*?主要功能：能對圓盤分幅數據聯合要求所有分幅文件在同一目錄下，并將坐標轉換為經緯度坐標
?*?如果要地形矯正將地形文件china.img放在執行時的當前目錄，
?*?默認情況下投影區域為?經度70-150，緯度0-60范圍，格點精度0.01度，默認處理通道為所提供文件的通道值，
?*?默認輸出文件為將hsd文件的.DAT替換為.ldf
???用法：H8proj.exe?[-dstwin?LULON?LULAT?DX?DY]?[-size?nx?ny]?[-bands?xyz...]?[-addsun]?[-addsat]?srcfile?[dstfile]
?*?例子：1.?H8proj.exe?-dstwin?70?60?0.01?-0.01?-size?8000?6000?-bands?123?HS_H08_20150715_0130_B03_FLDK_R05_S0210.DAT?H08_20150715_0130.img
????????2.?decodeHSD.exe?HS_H08_20150715_0130_B03_FLDK_R05_S0210.DAT

?*?內存要求，因為程序在處理投影時，會自動打開投影區域內涉及的所有分幅文件，當選擇區域較大時系統內存要較大
?*?浙江氣候中心??張小偉??qq：51903418
?*
?*
?*?2015.7.16?修正計算亮溫時的錯誤，還未測試，亮溫計算是否正確
?*?2015.7.17?增加ldf格式輸出，使輸出結果可以在SWAP軟件下打開
?*?2015.7.20?使用官方swap?ldf文件頭修正計算亮溫時，公式中的波長單位錯誤
?*?2016.1.15?東華王新剛增加多通道處理，衛星天頂角層?，多字節字符處理等功能
?*?2016.4.29?和衛星中心陳杰討論，增加地形校準功能
?*?2016.5.3??將衛星天頂角的計算從原來的讀取插值法換成直接計算法不同系統不同路徑分割
?*?2016.8.2??修改路徑處理中的bug，增加pack選項來適應hsd數據結構的對齊問題
?*?2016.10.17?增加太陽高度角計算部分，利用了hsd文件中的太陽位置，其值與普通天文計算公式位置有1-2%差異，
?*??還未確切差異的原因，歡迎賜教和討論。
?*/

#include?
#include?
#include?
#include?
#include?
#include?
#include?
#include?
#include?
#include?

#ifdef?_MSC_VER
#include?“stdafx.h“
#endif

using?namespace?std;

#if?defined（WIN32）?||?defined（WIN64）
#define?pathsplit?“\\“
#else
#define?pathsplit?“/“
#endif

const?double?omega_e=1.00273790934;?//solar?day?over?sidereal?day
const?double?pi=3.14159265358979323846;
const?double?r2d=180/pi;
const?double?d2r=pi/180;
const?double?Re=6378.135;
const?double?AU=1.49597870691E8;	/*?Astronomical?unit?-?km?（IAU?76）?*/

const?double?topdx=0.008333333333333333topdy=0.008333333333333333;
const?double?topulat=70.0topulon=70.0;
const?int?topnx=8400topny=8400;
const?long?topnn=topnx*topny;
int?topflag=0;
short?topdata[topnn];

double?getHeight（const?char?*filedouble?londouble?lat）{
	FILE?*fp;
	int?ij;
	double?h=0;

	if（topflag==0）{
		if（（fp=fopen（file“rb“））==NULL）{
		fprintf（stderr“can‘t?open?%s\n“file）;
		topflag=-1;
		return?0;?}
		fread（topdatasizeof（short）topnnfp）;
		fclose（fp）;
		topflag=1;
	}
	if（topflag==-1）?return?0;
	j=（lon-topulon）/topdx+0.5;
	i=（topulat-lat）/topdy+0.5;
	if（i>=0?&&?i<8400?&&?j>=0?&&?j<8400）{
		if（topdata[i*8400+j]>0?&&?topdata[i*8400+j]<10000）?h=topdata[i*8400+j]/1000.0;
	}
	return?h;

}

string?getFilePath（char*?fullPath）;
//#define?DEBUG?1

/*****************************************************
?*頭信息?普通版本128字節，擴展版本?128+1024?字節。
******************************************************/
typedef?struct?swapldf?{
????/*0*/short?wFileID?;//?2字節。文件標識?局地文件為“LA“?植被?NV?--?NDVI??植被?DV?--?DVI?植被??RV?--?RVI?LST?--?LS??SSY?--?SS
????/*2*/short?wSatelliteID;??//?2字節。?衛星標識?NOAA-14:?14?NOAA-12:?12?FY-1C:103?TERRA:1001
????/*4*/short?wOrbitNo;????????????/*?2字節。?軌道號?*/
????/*6*/short??wUporDown;????????????/*??2字節。升降軌標記?1:?升軌?0:?降軌?*/
????