private void QueryIDs(string[] dataSets)
{
long setID, prdsID;
double selfFillValue, dayFillValue, dayInvalidValue;
string set;
foreach (string setName in dataSets)
{
set = setName.Replace("_", "");
if (_dbCon.QueryDatasetsID("AIRS", set, out prdsID, out setID))
{
_pro._datasetsIDs.Add(set, setID);
_pro._datasetsPrdsIDs.Add(set, prdsID);
}
if (_dbCon.QueryDatasetsInvalidValue("AIRS", set, out selfFillValue, out dayFillValue, out dayInvalidValue))
{
_pro._selfFillValue.Add(set, selfFillValue);
_pro._dayFillValue.Add(set, dayFillValue);
_pro._dayInvalidValue.Add(set, dayInvalidValue);
}
}
long regionID;
foreach (PrjEnvelopeItem prjItem in _prjEnvelopes)
{
if (!_dbCon.IshasDataRegionRecord(prjItem, out regionID)) //!_dbCon.IshasRecord(tableName, "regionName", prjItem.Name.ToLower()))
{
regionID = _dbCon.InsertRegionTable(prjItem.Name.ToLower(), prjItem.PrjEnvelope.MinX, prjItem.PrjEnvelope.MaxX, prjItem.PrjEnvelope.MinY, prjItem.PrjEnvelope.MaxY);
}
_pro._regionIDs.Add(prjItem.Name, regionID);
}
}
public void TryCreatePeriodRecord(string file, Action <int, string> progressCallback = null)
{
string fname, localfname, region, sensor = "", statName, datasetname, period;
string[] parts;
long prdID = 0, peridID, staticsID, datasetID;
int year = 0, month = 0, xun = 0, regionID;
//double validPercent;
//CloudTopTemperature_MOD06_china_Ten_2011-11-2_AVG_day_0.05.LDF
//CloudTopTemperature_MOD06_china_Month_2011-11_AVG_day_0.05.LDF
//CloudTopTemperature_MOD06_china_Year_2011_AVG_day_0.05.LDF
if (progressCallback != null)
{
progressCallback(-1, "\t\t\t\t入库开始...," + file);
}
fname = Path.GetFileNameWithoutExtension(file);
localfname = file.Replace(_inputDIR, "");
parts = fname.Split('_');
period = parts[3].ToLower();
if (period != "ten" && period != "month" && period != "year")
{
return;
}
if (!_periodTypeID.ContainsKey(period))
{
return;
}
peridID = _periodTypeID[period];//周期类型
string sta = parts[5].ToUpper();
if (sta != "MIN" && sta != "MAX" && sta != "AVG")
{
return;
}
statName = sta;
staticsID = _statName2ID[sta];//统计类型
string label = parts[1].ToUpper();
if (label == "MOD06" || label == "MYD06" || label == "AIRS")
{
sensor = label;
}
else
{
return;
}
datasetname = parts[0];
if (_dbConnect.QueryDatasetsID(sensor, datasetname, out prdID, out datasetID))
{
region = parts[2];
if (!_regionID.Keys.Contains(region))
{
if (!_dbConnect.QueryRegionID(region.ToLower(), out regionID))
{
return;
}
_regionID.Add(region, regionID);
}
else
{
regionID = _regionID[region];
}
if (regionID == 0)
{
regionID = 1;
}
if (period == "ten")
{
Regex XunDateReg = new Regex(@"(?<year>\d{4})-(?<month>\d{2})-(?<Xun>\d{1})", RegexOptions.Compiled);
Match match = XunDateReg.Match(parts[4]);
if (!match.Success)
{
return;
}
year = Int16.Parse(match.Groups["year"].Value);
month = Int16.Parse(match.Groups["month"].Value);
xun = Int16.Parse(match.Groups["Xun"].Value);
}
else if (period == "month")
{
Regex MonthDateReg = new Regex(@"(?<year>\d{4})-(?<month>\d{2})", RegexOptions.Compiled);
Match match = MonthDateReg.Match(parts[4]);
if (!match.Success)
{
return;
}
year = Int16.Parse(match.Groups["year"].Value);
month = Int16.Parse(match.Groups["month"].Value);
}
else//period == "year"
{
year = Int16.Parse(parts[4]);
}
//GetRsEnvolop(file, out resol, out validPercent);
float resol = float.Parse(parts[parts.Length - 1]);
string datasource = "";
if (parts[parts.Length - 2].ToLower() == "day")
{
datasource = "D";
}
else
{
datasource = "N";
}
if (!_dbConnect.IshasRecord(_tableName, "ImageName", Path.GetFileName(file)))
{
_dbConnect.InsertPeriodicSynProductsTable(prdID, datasetID, localfname, regionID, region, sensor, resol, peridID, year, month, xun, _statName2ID[statName], datasource);
}
else
{
_dbConnect.DeleteCLDParatableRecord(_tableName, "ImageName", Path.GetFileName(file));
_dbConnect.InsertPeriodicSynProductsTable(prdID, datasetID, localfname, regionID, region, sensor, resol, peridID, year, month, xun, _statName2ID[statName], datasource);
}
}
}
private void DoProcess()
{
if (_originFiles2Base.Count != 0)
{
_originFiles2Base.Clear();
}
if (cbxPrdsLevl.SelectedIndex == 0)//原始数据
{
try
{
string subdataDocDir = Path.Combine(_dataDocDir, "原始数据");
#region 原始数据入库
//_logfName = "RawData2Database";
if (radiMODIS.Checked)
{
GetMODOriginDataFiles(_state);
if (!_originFiles2Base.Keys.Contains("MOD06") || _originFiles2Base["MOD06"].Count < 1)
{
txtErrorLog.Text += "当前路径不存在可入库的MODIS数据！\r\n";
}
}
else if (radiAIRS.Checked)
{
GetAIRSOriginDataFiles(_state);
if (!_originFiles2Base.Keys.Contains("AIRS") || _originFiles2Base["AIRS"].Count < 1)
{
txtErrorLog.Text += "当前路径不存在可入库的AIRS数据！\r\n";
}
}
else if (radiISCCP.Checked)
{
GetISCCPOriginDataFiles(_state);
if (_originFiles2Base.Keys.Contains("ISCCP"))
{
if (_state != null)
{
_state("共计" + _originFiles2Base["ISCCP"].Count + "个待入库ISCCP文件！");
}
else
{
txtErrorLog.Text += "共计" + _originFiles2Base["ISCCP"].Count + "个待入库ISCCP文件！\r\n";
}
}
}
else if (radiCloudSAT.Checked)
{
GetCloudSATOriginDataFiles(_state);
if (!_originFiles2Base.Keys.Contains("CloudSAT") || _originFiles2Base["CloudSAT"].Count < 1)
{
txtErrorLog.Text += "当前路径不存在可入库的CloudSAT数据！\r\n";
}
}
if (_originFiles2Base.Count == 0)
{
MessageBox.Show("当前目录不存在可入库的文件！请重新选择");
return;
}
#region 数据归档
//Dictionary<string, List<string>> uniformOriginFiles = new Dictionary<string, List<string>>();
//foreach (string mode in _originFiles2Base.Keys)
//{
// if (_state != null)
// _state(string.Format("正在归档{0}数据，请稍候...", mode));
// List<string> unifiles = FileToDatabase.Files2UniformDir(_originFiles2Base[mode], _dataDocDir, mode, _state);
// uniformOriginFiles.Add(mode, unifiles);
//}
#endregion
try
{
OriginData2Database OriDbase = new OriginData2Database(subdataDocDir, _originFiles2Base, cbxOverrideRecord.Checked, _state);
OriDbase.IsFiles2UniformDir = cbxData2DocDir.Checked;
OriDbase._DocDir = _dataDocDir;
runTaskThread1 = new Thread(new ThreadStart(OriDbase.CheckFile2Table));
runTaskThread1.IsBackground = true;
if (_state != null)
{
_state(string.Format("开始将数据归档并入库，请稍候..."));
}
runTaskThread1.Start();
}
catch (System.Exception ex)
{
_state(ex.Message);
return;
}
#endregion
}
catch (System.Exception ex)
{
_state(ex.Message);
return;
}
}
else if (cbxPrdsLevl.SelectedIndex == 1)//历史日产品数据
{
#region 日产品数据入库
string fileFilter, sensor;
if (radiAIRS.Checked)
{
fileFilter = "*AIRS*day*.ldf";//文件名的格式有待确认
sensor = "AIRS";
}
else
{
fileFilter = "mod06_*_*.dat";//文件名的格式有待确认
sensor = "MODIS";
}
List <string> dayprdsfiles = Directory.GetFiles(_inputDir, fileFilter, SearchOption.AllDirectories).ToList();
//归档；
//从文件名中解析数据的时间信息、数据集信息、区域信息、分辨率信息、日夜标识信息
//拼接生成数据的归档路径，对数据按照现有的规则进行重命名；//日拼接产品\CloudTopTemperature\TERRA\MODIS\2011\1\day\0.05
List <string> unidayprdsfiles = new List <string>();
//string newfnameformat = "{0}_MOD06_china_day_{1}_0.01.dat";
string newfname = "";
foreach (string oldfile in dayprdsfiles)
{
//mod06_20131201_CMLF.dat
//CloudOpticalThickness_MOD06_china_day_20110101_0.01.LDF
System.IO.File.Move(oldfile, newfname);
}
//入库---注意数据中的无效值问题
//需要查询数据集对应的产品ID及数据集ID，regionID，
//double selfFillValue = _selfFillValue[setName];
//double dayFillValue = _dayFillValue[setName];
//double dayInvalidValue = _dayInvalidValue[setName];
//重新计算有效百分比，granulesCount/granulesTimes、日夜标识DataSource
#region 入库
//归档路径
string tableName = "CP_DayMergeProducts_TB";
ConnectMySqlCloud dbcon = new ConnectMySqlCloud(_dataBaseXml);
string setName, imagedata, regionName, fname;
float resl;
double invalidValue = 32767;
int validPercent, regionID = 1;;
long prdID, datasetID = 0;
string datasource = "";
string[] parts;
DateTime fdate;
DataProcesser datapro = new DataProcesser();
foreach (string newFileName in unidayprdsfiles)
{
imagedata = newFileName.Replace(_dataDocDir, "");
fname = Path.GetFileNameWithoutExtension(newFileName);
parts = fname.Split('_');
if (parts.Length != 6)
{
continue;
}
setName = parts[0];
if (!dbcon.QueryDatasetsID(sensor, setName, out prdID, out datasetID))
{
continue;
}
regionName = parts[2];
if (!dbcon.QueryRegionID(regionName, out regionID))
{
}
if (parts[3].ToLower() == "day")
{
datasource = "D";
}
else
{
datasource = "N";
}
//fdate =DataProcesser.getDayTime(newFileName);
fdate = DateTime.TryParseExact(Path.GetFileNameWithoutExtension(newFileName).Split('_')[4], "yyyyMMdd", System.Globalization.CultureInfo.CurrentCulture, System.Globalization.DateTimeStyles.None, out fdate) ? fdate : DateTime.MinValue;
validPercent = ComputeValidPercent(datapro, newFileName, invalidValue, out resl);
if (!dbcon.IshasRecord(tableName, "ImageName", Path.GetFileName(newFileName)))
{
dbcon.InsertNDayMergeProductsTable(fdate, prdID, datasetID, imagedata, regionID, regionName, sensor, resl, validPercent, 0, "", datasource);
_state(newFileName + "入库完成！");
}
else
{
dbcon.DeleteCLDParatableRecord(tableName, "ImageName", Path.GetFileName(newFileName));
dbcon.InsertNDayMergeProductsTable(fdate, prdID, datasetID, imagedata, regionID, regionName, sensor, resl, validPercent, 0, "", datasource);
_state(newFileName + "入库更新完成！");
}
}
#endregion
#endregion
}
else if (cbxPrdsLevl.SelectedIndex == 2)//周期合成产品
{
//PeriodicSynPrds2Base psp = new PeriodicSynPrds2Base(_inputDir, _logfName, _dataBaseXml);
//psp.StartComp();
}
}
public void Compute(Action <int, string> progressCallback)
{
if (File.Exists(_dataBaseXml))
{
_con = new ConnectMySqlCloud();
}
else
{
if (progressCallback != null)
{
progressCallback(-1, "数据库配置文件" + @"\SystemData\ProductArgs\CLD\CPDataBaseArgs.xml" + "不存在！");
}
return;
}
_periodicsyn = new PeriodicSynPrds2Base(_args.InputDir, "周期合成入库", _con);
_statics = _args.StatisticsTypes;//统计类型
_overwritefiles = _args.OverWriteHistoryFiles;
string[] periods = null;
if (_periodPrdsArgs == "YEAR")
{
periods = new string[3] {
"YEAR", "MONTH", "TEN"
}
}
;
else if (_periodPrdsArgs == "MONTH")
{
periods = new string[2] {
"MONTH", "TEN"
}
}
;
//string[] periods = _args.PeriodTypes;
_inputDataPath = Path.Combine(_args.InputDir, "日拼接产品");
if (!Directory.Exists(_inputDataPath))
{
if (progressCallback != null)
{
progressCallback(-1, "输入文件目录" + _inputDataPath + "不存在！请重试！");
}
return;
}
if (progressCallback != null)
{
progressCallback(-1, "输入文件夹：" + _inputDataPath);
}
string outputdir = Path.Combine(_args.OutputDir, "周期合成产品");
if (progressCallback != null)
{
progressCallback(-1, "输出文件夹：" + outputdir);
}
if (!Directory.Exists(outputdir))
{
Directory.CreateDirectory(outputdir);
}
//日拼接产品\CloudEffectiveEmissivity\TERRA\MODIS\2011\1\day或night\resl
int al = _validsets.Length;
if (al == 0)
{
return;
}
float interval = 99 / al / 2;
int i = -1, pp;
string childpath;
List <string> childPath = new List <string>();
foreach (string set in _validsets)
{
i++;
pp = i + 1;
string setName = set.Replace("_", "");
double filvalue, invalidvalue, setfill;
//childpath = Path.Combine(_inputDataPath, set.Replace("_", ""));
childpath = Path.Combine(new string[] { _inputDataPath, set.Replace("_", ""), _satellite, _sensor, _singleyear.ToString() });
if (progressCallback != null)
{
progressCallback((int)(i * interval * 2) + 1, "共" + _validsets.Length + "个数据集,开始处理第" + pp + "个，" + set + ",文件夹：" + childpath);
}
if (progressCallback != null)
{
progressCallback(-1, "\t链接数据库,查询数据集ID及无效值...");
}
if (!_con.QueryDatasetsID(_dataOrigin, setName, out _prdID, out _datasetID))//_sensor
{
if (progressCallback != null)
{
progressCallback(-1, "\t链接数据库,查询数据集ID失败！");
}
continue;
}
if (!_con.QueryDatasetsInvalidValue(_dataOrigin, setName, out setfill, out filvalue, out invalidvalue))//_sensor
{
if (progressCallback != null)
{
progressCallback(-1, "\t链接数据库,查询数据集无效值失败！");
}
continue;
}
RasterStatics sta = new RasterStatics();
sta.SetDstFillValues(filvalue.ToString());
sta.SetDstInvalidValues(invalidvalue.ToString());
//sta.SetDstInvalidValues("0");
#region 原有算法
string dayfileformat = null;
//CloudFractionDay_MOD06_china_day_20060101_0.05.LDF
string[] danightlabels = new string[2] {
"day", "night"
};
int dnc = 0, kk;
foreach (string label in danightlabels)
{
if (progressCallback != null)
{
progressCallback((int)(interval * (i + dnc++ / 2.0f)) + 1, "\t开始扫描数据集" + set + "的" + _singleyear + "年的" + label + "文件...");
}
kk = dnc;
//CloudFractionDay_MOD06_china_day_20060101_0.05.LDF
DirectoryInfo yearDirInfo = new DirectoryInfo(childpath);
int monthCount = yearDirInfo.GetDirectories().Length;
if (progressCallback != null)
{
progressCallback(-1, "\t共" + monthCount + "个月待处理...");
}
if (monthFiles != null)
{
monthFiles.Clear();
}
dayfileformat = string.Format("{0}*{1}*{2}_{3}*.ldf", setName, _dataOrigin, label, _singleyear);
foreach (DirectoryInfo d in yearDirInfo.GetDirectories()) //查找子目录，月
{
//日拼接产品\CloudEffectiveRadius\TERRA\MODIS\2006\1\day\0.01
if (progressCallback != null)
{
progressCallback(-1, "\t开始扫描" + _singleyear + "年" + d.Name + "月的旬文件...");
}
TryComputeXunFiles(Path.Combine(yearDirInfo.FullName, d.ToString()), dayfileformat, progressCallback);
}
//string[] parafiles = Directory.GetFiles(childpath, string.Format(dayfileformat, setName,_dataOrigin, label,_singleyear), SearchOption.AllDirectories);
//if (parafiles.Length == 0)
//{
// if (progressCallback != null)
// progressCallback(-1, "\t没有符合数据集" + set + "的" + _singleyear + "年的" + label + "的文件！");
// continue;
//}
////按月进行分组
//Dictionary<string, List<string>> monthFiles = new Dictionary<string, List<string>>();
//foreach (string file in parafiles.ToArray())
//{
// string path = Path.GetDirectoryName(file);
// if (!path.Replace(childpath, "").Contains("\\" + _singleyear.ToString() + "\\"))
// continue;
// if (!monthFiles.ContainsKey(path))
// monthFiles.Add(path, new List<string>(new string[1] { file }));
// else
// monthFiles[path].Add(file);
//}
//if (progressCallback != null)
// progressCallback(-1, "\t共" + monthFiles.Count + "个月的" + parafiles.Length + "个日拼接文件...");
//旬文件名，产品类型_MOD06_region_旬产品_yyyy-mm-N_MAX/MIN/AVG_day或night_0.05.LDF
Dictionary <string, string> dirxunlist = new Dictionary <string, string>();
//将数据按区域进行分组；
float m = -1;
int ml = monthFiles.Keys.Count;
if (ml == 0)
{
if (progressCallback != null)
{
progressCallback(-1, "\t没有符合数据集" + set + "的" + _singleyear + "年的" + label + "文件！");
}
continue;
}
string[] monthdirparts = null;
float reslf = 0;
foreach (string monthdir in monthFiles.Keys)
{
int start = monthdir.LastIndexOf("\\") + 1;
string resl = monthdir.Substring(start, monthdir.Length - start);
monthdirparts = monthdir.Split('\\');
if (!float.TryParse(monthdirparts[monthdirparts.Length - 1], out reslf) || !danightlabels.Contains(monthdirparts[monthdirparts.Length - 2].ToLower()))
{
continue;
}
string dirXun = "{0}\\{1}\\{2}\\{3}\\{4}";
dirXun = string.Format(dirXun, outputdir, set.Replace("_", ""), _satellite, _sensor, label + "\\" + resl + "\\Ten");
if (!dirxunlist.ContainsKey(dirXun))
{
dirxunlist.Add(dirXun, resl);
}
m++;
Dictionary <string, List <string> > regionFiles = new Dictionary <string, List <string> >();
foreach (string file in monthFiles[monthdir])
{
string fileName = Path.GetFileNameWithoutExtension(file);
string[] parts = fileName.Split('_');
//CloudMultiLayerFlag_MOD06_china_day_20110101_0.01.LDF
if (parts.Length < 6)
{
continue;
}
string region = parts[2];
if (!regionFiles.ContainsKey(region))
{
regionFiles.Add(region, new List <string>(new string[1] {
file
}));
}
else
{
regionFiles[region].Add(file);
}
}
if (progressCallback != null)
{
progressCallback((int)(interval * (i + kk / 2.0f + m / ml * 0.8)) + 1, "\t\t正在合成" + monthdir + "的旬产品...");
}
ComputeTenPeriodSyn(setName, dirXun, regionFiles, sta, resl, progressCallback, label);
}
string dirnext;
if (periods.Contains("MONTH"))
{
if (dirxunlist.Count < 1)
{
continue;
}
if (progressCallback != null)
{
progressCallback(-1, "\t开始合成" + setName + "的" + _singleyear + "年" + label + "的月产品...");
}
int xlength = dirxunlist.Count;
int xcur = 0;
foreach (string dirXun in dirxunlist.Keys)
{
string resl = dirxunlist[dirXun];
string flabel = string.Format("_{0}_{1}.LDF", label, resl);
if (progressCallback != null)
{
progressCallback((int)(interval * (i + kk / 2.0f + xcur++ / (1.0f * xlength) * 0.2 + 0.8)) + 1, "\t\t正在合成" + dirXun + "的月产品...");
}
UpdateDataNext(dirXun, "Ten", "Month", out dirnext, sta, progressCallback, flabel);
if (periods.Contains("YEAR"))
{
if (progressCallback != null)
{
progressCallback(-1, "\t\t正在合成" + dirnext + "的年产品...");
}
UpdateDataNext(dirnext, "Month", "Year", out dirnext, sta, progressCallback, flabel);
}
}
}
}
#endregion
}
if (progressCallback != null)
{
progressCallback(100, "周期合成完成！");
}
return;
}
