private void button6_Click(object sender, EventArgs e)
{
FileCalc lfilecalc = new FileCalc();
lfilecalc.calculate(calcRootPath);
MessageBox.Show("Completed!!");
}
public void FileCalcTestDivide()
{
var target = new FileCalc("..\\..\\..\\Calculator\\inputDivide.txt", new OperationRegistryImpl());
Assert.AreEqual(target.Compute(), 1);
}
public void FileCalcTestMultiply()
{
var target = new FileCalc("..\\..\\..\\Calculator\\inputMultiply.txt", new OperationRegistryImpl());
Assert.AreEqual(target.Compute(), 325);
}
private static void Anomaly(string Folder)
{
foreach (string FileCalc in Directory.EnumerateFiles(GeoServerModisDataDir, "*.tif", SearchOption.TopDirectoryOnly))
{
if (FileCalc.Contains("Anomaly") || FileCalc.Contains("BASE"))
{
continue;
}
// calculate anomaly
string fileNameWithoutExtension = Path.GetFileNameWithoutExtension(FileCalc),
anomalyFile = Path.Combine(GeoServerModisDataDir, Path.ChangeExtension(fileNameWithoutExtension + "_Anomaly", "tif")),
baseFile = Path.Combine(GeoServerModisDataDir, Path.ChangeExtension($"ABASE{fileNameWithoutExtension.Substring(5)}", "tif")),
letters = "ABCDEFGHIJKLMNOPQRSTUVWXY";
// check if anomaly already exists
if (File.Exists(anomalyFile))
{
continue;
}
// check if base layer for anomaly calculation already exists, if no then try to create it
string arguments = "--co COMPRESS=LZW";
if (!File.Exists(baseFile))
{
// check if base layers for base calculation already exist
bool baseExists = true;
for (int year = AnomalyStartYear; year <= AnomalyFinishYear; year++)
{
string baseYearFile = Path.Combine(GeoServerModisDataDir, Path.ChangeExtension(fileNameWithoutExtension.Remove(1, 4).Insert(1, year.ToString()), "tif"));
if (!File.Exists(baseYearFile))
{
baseExists = false;
break;
}
}
if (!baseExists)
{
continue;
}
// create base file to day
for (int year = AnomalyStartYear; year <= AnomalyFinishYear; year++)
{
int letterIndex = year - AnomalyStartYear;
string baseYearFile = Path.ChangeExtension(fileNameWithoutExtension.Remove(1, 4).Insert(1, year.ToString()), "tif");
arguments += $" -{letters[letterIndex]} {Path.GetFileName(baseYearFile)}";
}
arguments += $" --outfile={Path.GetFileName(baseFile)}";
arguments += $" --calc=\"((";
for (int year = AnomalyStartYear; year <= AnomalyFinishYear; year++)
{
int letterIndex = year - AnomalyStartYear;
arguments += $"{letters[letterIndex]}+";
}
arguments = arguments.Remove(arguments.Length - 1);
arguments += $")/{(AnomalyFinishYear - AnomalyStartYear + 1).ToString()})\"";
GDALExecute(
CMDPath,
"gdal_calc.py",
Folder,
arguments);
}
// calculate
arguments = "--co COMPRESS=LZW";
arguments += $" -{letters[0]} {Path.GetFileName(baseFile)}";
arguments += $" -{letters[1]} {Path.GetFileName(FileCalc)} ";
arguments += $"--outfile={Path.GetFileName(anomalyFile)} ";
arguments += $"--calc=\"(B-A)*0.01\"";
GDALExecute(
CMDPath,
"gdal_calc.py",
Folder,
arguments);
// publish
string layerName = Path.GetFileNameWithoutExtension(anomalyFile);
// store
string publishParameters = $" -v -u" +
$" {GeoServerUser}:{GeoServerPassword}" +
$" -POST -H \"Content-type: text/xml\"" +
$" -d \"<coverageStore><name>{layerName}</name><type>GeoTIFF</type><enabled>true</enabled><workspace>{GeoServerWorkspace}</workspace><url>" +
$"/data/{GeoServerWorkspace}/{layerName}.tif</url></coverageStore>\"" +
$" {GeoServerURL}rest/workspaces/{GeoServerWorkspace}/coveragestores?configure=all";
CurlExecute(
CMDPath,
publishParameters);
// layer
publishParameters = $" -v -u" +
$" {GeoServerUser}:{GeoServerPassword}" +
$" -PUT -H \"Content-type: text/xml\"" +
$" -d \"<coverage><name>{layerName}</name><title>{layerName}</title><defaultInterpolationMethod><name>nearest neighbor</name></defaultInterpolationMethod></coverage>\"" +
$" \"{GeoServerURL}rest/workspaces/{GeoServerWorkspace}/coveragestores/{layerName}/coverages?recalculate=nativebbox\"";
CurlExecute(
CMDPath,
publishParameters);
// style
string index = ModisDataSetIndex.ToString().PadLeft(2, '0'),
style = $"{GeoServerWorkspace}:{ModisSource}_{ModisProduct.Replace(".", "")}_B{index}_{ModisDataSet}_Anomaly";
publishParameters = $" -v -u" +
$" {GeoServerUser}:{GeoServerPassword}" +
$" -X PUT -H \"Content-type: text/xml\"" +
$" -d \"<layer><defaultStyle><name>{style}</name></defaultStyle></layer>\"" +
$" {GeoServerURL}rest/layers/{GeoServerWorkspace}:{layerName}.xml";
CurlExecute(
CMDPath,
publishParameters);
}
}
