static void Main(string[] args)
{
///////////////////////////////////////////////////////////////////////////////
//介绍性开头
Console.WriteLine("++++++++++++++++++++++++++++++++++++++++++++++++++++++");
Console.WriteLine("+++ Hour QPE Analysis V1.0 +++");
Console.WriteLine("+++++ Supproted By CaoYong 2018.11.12 +++++++++");
Console.WriteLine("+++++ QQ: 403637605 +++++++++");
Console.WriteLine("++++++++++++++++++++++++++++++++++++++++++++++++++++++");
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
//打开计时器
Stopwatch sw = new Stopwatch(); //创建计时器
sw.Start(); //开启计数器
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
string appDir = System.AppDomain.CurrentDomain.SetupInformation.ApplicationBase; //程序启动文件夹
//string appDir = @"E:/HourQPE_V1/";
System.Environment.CurrentDirectory = appDir; //设置shell所在文件夹
string logPath = appDir + @"log/" + DateTime.Now.Year.ToString("d4") + DateTime.Now.Month.ToString("d2") + DateTime.Now.Day.ToString("d2") + @".txt"; //日志文件夹地址
Log simpleLog = new Log(logPath); //建立log对象，用于日志的记录 //输出站点ID计算信息
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
try
{
///////////////////////////////////////////////////////////////////////////
//时间处理(北京时)
DateTime dtNow = DateTime.Now; //程序启动时间（北京时）
if (args.Length == 0) //实时运算处理
{
dtNow = DateTime.Now;
//dtNow = new DateTime(2018, 11, 12, 14, 00, 00);
}
else if (args.Length == 1 && args[0].Length == 12) //指定日期运算处理
{
try
{
int argYr = int.Parse(args[0].Substring(0, 4));
int argMo = int.Parse(args[0].Substring(4, 2));
int argDy = int.Parse(args[0].Substring(6, 2));
int argHr = int.Parse(args[0].Substring(8, 2));
int argMn = int.Parse(args[0].Substring(10, 2));
dtNow = new DateTime(argYr, argMo, argDy, argHr, argMn, 0);
}
catch (Exception ex)
{
simpleLog.WriteError(ex.Message, 1);
simpleLog.WriteError("date args content is not right!", 1);
return;
}
}
else
{
simpleLog.WriteError("date args is not right!", 1);
return;
}
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
//读取控制文件
string paraFilePath = appDir + @"para/para.ini"; //控制文件地址
string inputSamplePath = null; //自动站站点数据输入地址
string outputSamplePath = null; //QPE分析格点输出地址
if (!File.Exists(paraFilePath))
{
simpleLog.WriteError("para file is not exist!", 1);
return;
}
else
{
FileStream paraFS = new FileStream(paraFilePath, FileMode.Open, FileAccess.Read);
StreamReader paraSR = new StreamReader(paraFS, Encoding.GetEncoding("gb2312"));
{
try
{
string strTmp = paraSR.ReadLine();
string[] strArrayTmp = strTmp.Split(new char[] { '=' });
inputSamplePath = strArrayTmp[1].Trim(); //获取输入地址
strTmp = paraSR.ReadLine();
strArrayTmp = strTmp.Split(new char[] { '=' });
outputSamplePath = strArrayTmp[1].Trim(); //获取输出地址
}
catch
{
simpleLog.WriteError("para content is not right!", 1);
return;
}
}
paraSR.Close();
paraFS.Close();
}
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
string stainfoPath = appDir + @"info/sta.info"; //此处以完整第3类数据站点数据为站点信息
DateTime dtStart = dtNow.AddHours(-6);
DateTime dtEnd = dtNow;
DateTime dtCurrent = dtStart;
while (DateTime.Compare(dtCurrent, dtEnd) <= 0) //冗余设置，分析之前6小时至当前的结果
{
DateTime dtInput = dtCurrent.AddHours(-8); //原始数据为世界时
DateTime dtOutput = dtCurrent.AddHours(-8); //输出结果为世界时
string inputFilePath = StringProcess.DateReplace(inputSamplePath, dtInput, 000);
string outputFilePath = StringProcess.DateReplace(outputSamplePath, dtOutput, 000);
if (File.Exists(outputFilePath + ".m4")) //存在输出文件则不进行分析
{
dtCurrent = dtCurrent.AddHours(1);
continue;
}
Console.WriteLine("inputPath=" + inputFilePath);
Console.WriteLine("OutputPath=" + outputFilePath);
ScatterData sdInputData = new ScatterData(stainfoPath);
sdInputData.ClearToNum(0.0);
sdInputData.ReadValFromMicaps3(inputFilePath);
sdInputData.ClearToNumGreaterThan(0.0, 100.0); //简单质量控制
sdInputData.ClearToNumLessThan(0.0, 0.10); //简单质量控制，注意冬季情况的适用性
GridData gdOutputData = new GridData(70.0, 140.0, 0.0, 60.0, 0.05, 0.05);
gdOutputData.ClearToNum(0.0);
gdOutputData = SpatialAnalisis.GressManInterpolationForRain(sdInputData, gdOutputData, 0.2, 2.0, new double[] { 1.0, 0.8, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1 }); //分析25-5公里尺度
gdOutputData.WriteFloatValToBin(outputFilePath + ".bin");
string strHeader = StringProcess.DateReplace(@"diamond 4 YYYYMMDDHH00_01hour_QPE YYYY MM DD HH NN 000 0.05 0.05 70.0 140.0 0.0 60.0 1401 1201 2.0 -2.0 20.0 1 00", dtOutput, 000);
gdOutputData.WriteValToMicaps4(outputFilePath + ".m4", strHeader);
///////////////////////////////////////////////////////////////////////////////
//检验程序
ScatterData sdFromQPE = new ScatterData(stainfoPath);
sdFromQPE.ClearToNum(0.0);
sdFromQPE.BilinearInterpolationFromGridData(gdOutputData);
double varError = 0.0;
for (int n = 0; n < sdFromQPE.StaData.Length; n++)
{
varError = varError + (sdFromQPE.StaData[n].Value - sdInputData.StaData[n].Value) * (sdFromQPE.StaData[n].Value - sdInputData.StaData[n].Value);
}
varError = varError / sdFromQPE.StaData.Length;
varError = Math.Sqrt(varError);
simpleLog.WriteInfo("Error is " + varError, 1);
///////////////////////////////////////////////////////////////////////////////
dtCurrent = dtCurrent.AddHours(1);
}
///////////////////////////////////////////////////////////////////////////////
}
catch (Exception ex)
{
simpleLog.WriteError(ex.Message, 1);
}
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
sw.Stop();
simpleLog.WriteInfo("time elasped: " + sw.Elapsed, 1);
///////////////////////////////////////////////////////////////////////////////
}
