FarmMaps/FarmMapsApiClient
10
0
Fork 1
Code Issues Pull Requests Releases Activity

update blight

Browse Source
This commit is contained in:
Fedde Sijbrandij
2020-12-15 17:10:42 +01:00
parent f3a2b33810 02708ba01d
commit ad49bafa07
7 changed files with 367 additions and 94 deletions
Download Patch File Download Diff File Expand all files Collapse all files

3
FarmmapsNbs/Models/NitrogenInput.cs
View File

@@ -13,5 +13,8 @@ namespace FarmmapsNbs.Models
public string PotatoPurposeType { get; set; }
public int TargetYield { get; set; }
public JObject GeometryJson { get; set; }
public string InputLayerName { get; set; }
}
}

127
FarmmapsNbs/NbsApplication.cs
View File

@@ -9,6 +9,7 @@ using FarmmapsApi.Services;
using FarmmapsNbs.Models;
using Microsoft.Extensions.Logging;
using Newtonsoft.Json;
using Newtonsoft.Json.Linq;
using static FarmmapsApiSamples.Constants;
namespace FarmmapsNbs
@@ -16,12 +17,15 @@ namespace FarmmapsNbs
public class NbsApplication : IApplication
{
private const string DownloadFolder = "Downloads";
private const string SettingsFile = "settings.json";
private readonly ILogger<NbsApplication> _logger;
private readonly FarmmapsApiService _farmmapsApiService;
private readonly NitrogenService _nitrogenService;
private readonly GeneralService _generalService;
private Settings _settings;
public NbsApplication(ILogger<NbsApplication> logger, FarmmapsApiService farmmapsApiService,
GeneralService generalService, NitrogenService nitrogenService)
{
@@ -40,6 +44,10 @@ namespace FarmmapsNbs
if (!Directory.Exists(DownloadFolder))
Directory.CreateDirectory(DownloadFolder);
LoadSettings();
// !! this call is needed the first time an api is called with a fresh clientid and secret !!
await _farmmapsApiService.GetCurrentUserCodeAsync();
var roots = await _farmmapsApiService.GetCurrentUserRootsAsync();
@@ -59,6 +67,8 @@ namespace FarmmapsNbs
private async Task Process(List<UserRoot> roots, NitrogenInput input)
{
// !!specify if you are using an already created cropfield:
bool useCreatedCropfield = false;
var plantingDate = input.PlantingDate;
var measurementDate = input.MeasurementDate;
@@ -75,24 +85,83 @@ namespace FarmmapsNbs
_logger.LogError("Could not find a needed root item");
return;
}
var cropfieldItem = await _generalService.CreateCropfieldItemAsync(myDriveRoot.Code,
$"VRA NBS cropfield {input.OutputFileName}", plantingDate.Year, input.GeometryJson.ToString(Formatting.None));
// Use already created cropfield or create new one
Item cropfieldItem;
if (useCreatedCropfield == false || string.IsNullOrEmpty(_settings.CropfieldItemCode))
{
_logger.LogInformation("Creating cropfield");
cropfieldItem = await _generalService.CreateCropfieldItemAsync(myDriveRoot.Code,
$"VRA NBS cropfield {input.OutputFileName}", plantingDate.Year, input.GeometryJson.ToString(Formatting.None));
_settings.CropfieldItemCode = cropfieldItem.Code;
SaveSettings();
}
else
{
_logger.LogInformation("Cropfield already exists, trying to get it");
cropfieldItem = await _farmmapsApiService.GetItemAsync(_settings.CropfieldItemCode);
}
var geotiffItem = (Item)null;
if (input.File.Contains(".tif") || input.File.Contains(".geotiff")) {
_logger.LogInformation("input = tiff data");
var dataPath = Path.Combine("Data", input.File);
geotiffItem = await _generalService.UploadDataAsync(uploadedRoot, GEOTIFF_PROCESSED_ITEMTYPE, dataPath,
Path.GetFileNameWithoutExtension(input.File));
// No file input, use most recent satellite image
if (string.IsNullOrEmpty(input.File))
{
_logger.LogInformation("No specific data given, retrieving most recent satellite image");
if (geotiffItem == null) {
_logger.LogError("Could not find item for uploaded data");
return;
}
// check if satellite task not yet done, do here and save taskcode
if (useCreatedCropfield == false || string.IsNullOrEmpty(_settings.SatelliteTaskCode))
{
var satelliteTaskCode = await _generalService.RunSatelliteTask(cropfieldItem);
_settings.SatelliteTaskCode = satelliteTaskCode;
SaveSettings();
}
// Select a particular satellite item from satelliteTask
Item satalliteItem = await _generalService.FindSatelliteItem(cropfieldItem, _settings.SatelliteTaskCode);
// must be wdvi[1]
var inputType = (satalliteItem.Data["layers"] as JArray)?[1]["name"].ToString();
if (string.IsNullOrEmpty(inputType))
{
_logger.LogError("Could not get the input type name from the satellite item");
return;
}
// download the geotiff
_logger.LogInformation("Downloading geotiff file");
await _farmmapsApiService.DownloadItemAsync(satalliteItem.Code,
Path.Combine(DownloadFolder, $"nbs_inputSatelliteGeotiff_{input.OutputFileName}.zip"));
// overwrite measurement date by date of satellite item
measurementDate = satalliteItem.DataDate.Value;
geotiffItem = satalliteItem;
}
// (geo)tiff input:
else if (input.File.Contains(".tif") || input.File.Contains(".geotiff")) {
_logger.LogInformation("input = tiff data");
var dataPath = Path.Combine("Data", input.File);
geotiffItem = await _generalService.UploadDataAsync(uploadedRoot, GEOTIFF_PROCESSED_ITEMTYPE, dataPath,
Path.GetFileNameWithoutExtension(input.File));
if (geotiffItem == null) {
_logger.LogError("Could not find item for uploaded data");
return;
}
}
// json/shape input
else {
var isGeoJson = input.File.Contains("json");
var dataPath = Path.Combine("Data", input.File);
@@ -117,6 +186,12 @@ namespace FarmmapsNbs
Path.Combine(DownloadFolder, $"{input.OutputFileName}.input_geotiff.zip"));
}
_logger.LogInformation($"Calculating targetN with targetYield: {input.TargetYield}");
var targetNItem = await _nitrogenService.CreateTargetNItem(cropfieldItem);
var targetNData = await _nitrogenService.CalculateTargetN(cropfieldItem, targetNItem, plantingDate,
@@ -139,7 +214,7 @@ namespace FarmmapsNbs
_logger.LogInformation("Calculating uptake map");
var uptakeMapItem =
await _nitrogenService.CalculateUptakeMap(cropfieldItem, geotiffItem, plantingDate,
measurementDate, input.InputVariable);
measurementDate, input.InputVariable, input.InputLayerName);
if (uptakeMapItem == null)
{
_logger.LogError("Something went wrong with creating the uptakeMap");
@@ -154,8 +229,7 @@ namespace FarmmapsNbs
_logger.LogInformation("Calculating application map");
var applicationMapItem =
await _nitrogenService.CalculateApplicationMap(cropfieldItem, geotiffItem, plantingDate,
measurementDate,
input.InputVariable, targetNData.TargetN);
measurementDate, input.InputVariable, targetNData.TargetN, input.InputLayerName);
if (applicationMapItem == null)
{
@@ -189,5 +263,28 @@ namespace FarmmapsNbs
}
// Functions to save previously created cropfields
private void LoadSettings()
{
if (File.Exists(SettingsFile))
{
var jsonText = File.ReadAllText(SettingsFile);
_settings = JsonConvert.DeserializeObject<Settings>(jsonText);
}
else
{
_settings = new Settings();
}
}
private void SaveSettings()
{
if (_settings == null)
return;
var json = JsonConvert.SerializeObject(_settings);
File.WriteAllText(SettingsFile, json);
}
}
}

132
FarmmapsNbs/NitrogenInput.json
View File

@@ -35,7 +35,7 @@
[ 3.40843828875524, 50.638966444680605 ]
]
]
//},
},
//{
//"file": "[...].json",
//"inputVariable": "wdvi",
@@ -62,7 +62,7 @@
// ]
// ]
//}
}
//}
//},
//{
@@ -106,6 +106,124 @@
//{
// "file": "Scan_1_20190605.zip",
// "inputVariable": "irmi",
{
"file": "", // keep emptpy to use satellite image
"inputVariable": "wdvi",
"InputLayerName": "wdvi",
"outputFileName": "rtest1",
"plantingDate": "2020-05-01",
"measurementDate": "2020-06-14",
"potatoPurposeType": "consumption",
"targetYield": 60,
"geometryJson": {
"type": "Polygon",
"coordinates": [
[
[ 4.960707146896585, 52.800583669708487 ],
[ 4.960645975538824, 52.800470217610922 ],
[ 4.962140695752897, 52.799177147194797 ],
[ 4.967523821195745, 52.801502400041208 ],
[ 4.966336768950911, 52.802543735879809 ],
[ 4.961711880764330, 52.801009996856429 ],
[ 4.960707146896585, 52.800583669708487 ]
]
]
}
}
//,
//{
// "file": "", // keep emptpy to use satellite image
// "inputVariable": "wdvi",
// "InputLayerName": "wdvi",
// "outputFileName": "mtest1",
// "plantingDate": "2020-04-01",
// "measurementDate": "2020-06-24",
// "potatoPurposeType": "consumption",
// "targetYield": 60,
// "geometryJson": {
// "type": "Polygon",
// "coordinates": [
// [
// [ 3.575458869234128, 51.308707885669762 ],
// [ 3.575508957999423, 51.30878478562019 ],
// [ 3.576188404403633, 51.309372997559777 ],
// [ 3.576188872410267, 51.309374219701091 ],
// [ 3.576210290749152, 51.309430091473608 ],
// [ 3.57621266537704, 51.309483685674898 ],
// [ 3.576213477455834, 51.309502027981374 ],
// [ 3.577543136860447, 51.310682367527122 ],
// [ 3.57796629328866, 51.31104321405175 ],
// [ 3.578442479292087, 51.311449273042747 ],
// [ 3.57866702353106, 51.311636072786726 ],
// [ 3.57880446997978, 51.31157117599529 ],
// [ 3.579155910205885, 51.311863542729718 ],
// [ 3.579175814007489, 51.311875435159394 ],
// [ 3.579293885246395, 51.311936532835396 ],
// [ 3.579413896180069, 51.311998649478575 ],
// [ 3.579514543462617, 51.312041110734917 ],
// [ 3.579611760655688, 51.312082118352606 ],
// [ 3.579635115371588, 51.312093949652223 ],
// [ 3.579793143414486, 51.312189437140432 ],
// [ 3.579966991648108, 51.312286148850511 ],
// [ 3.580079704980967, 51.312332458349751 ],
// [ 3.580203717638148, 51.312336471368539 ],
// [ 3.580307101018293, 51.312330239539847 ],
// [ 3.580383836270609, 51.312317097185243 ],
// [ 3.580505207977176, 51.312279163554869 ],
// [ 3.580610387713855, 51.312233723091026 ],
// [ 3.5806309754483, 51.312226093729677 ],
// [ 3.580638516049738, 51.312223727082049 ],
// [ 3.58075536599681, 51.312186990706344 ],
// [ 3.580787745633303, 51.312176820129551 ],
// [ 3.580829682241423, 51.312167665428959 ],
// [ 3.58086828456562, 51.312162614898625 ],
// [ 3.580980493636721, 51.312147935609723 ],
// [ 3.581014352632766, 51.312145662592656 ],
// [ 3.581028980583245, 51.312145592849248 ],
// [ 3.581119189823368, 51.312145103215911 ],
// [ 3.581195330198145, 51.312144693075908 ],
// [ 3.581243537809229, 51.312148741603245 ],
// [ 3.58132480221972, 51.312163110548383 ],
// [ 3.581426517039001, 51.312181089466016 ],
// [ 3.581448095953263, 51.312184910295024 ],
// [ 3.581474337475052, 51.312191038736145 ],
// [ 3.581709405982819, 51.312260068944951 ],
// [ 3.581727319558337, 51.312266163697757 ],
// [ 3.581753583356718, 51.312276407881612 ],
// [ 3.581841655772683, 51.312310743468075 ],
// [ 3.581878851795624, 51.312325234405343 ],
// [ 3.581905889860924, 51.312338260798654 ],
// [ 3.581906254501594, 51.312338472113815 ],
// [ 3.582048236499295, 51.312422439022804 ],
// [ 3.58189849928915, 51.312481944577037 ],
// [ 3.583044298383354, 51.312095780444281 ],
// [ 3.582984006671231, 51.312017283925279 ],
// [ 3.582743535862999, 51.311699343434064 ],
// [ 3.582628599916243, 51.311582190756774 ],
// [ 3.581446834435509, 51.310511259982569 ],
// [ 3.580621864908701, 51.309767270412806 ],
// [ 3.579610575760466, 51.308860440593946 ],
// [ 3.579112608916012, 51.308394999612226 ],
// [ 3.578688808506157, 51.307968441218165 ],
// [ 3.578394256007207, 51.307644098092617 ],
// [ 3.578355980318371, 51.307607614964702 ],
// [ 3.578217977585775, 51.307654179547846 ],
// [ 3.577480636332469, 51.307921035607997 ],
// [ 3.575695560441006, 51.308577167973212 ],
// [ 3.575668643609169, 51.30855384769157 ],
// [ 3.575666204524265, 51.308551734020703 ],
// [ 3.575506397192348, 51.308609906947261 ],
// [ 3.575459139533024, 51.308653178431456 ],
// [ 3.575458869234128, 51.308707885669762 ]
// ]
// ]
// }
//}
//,
//{
// "file": "Scan_1_20190605.zip",
// "inputVariable": "irmi",
// "InputLayerName": "IRMI",
// "outputFileName": "vranbs3",
// "plantingDate": "2019-04-18",
// "measurementDate": "2019-06-20",
@@ -144,6 +262,12 @@
//{
// "file": "Scan_1_20190605.zip",
// "inputVariable": "irmi",
//}
//,
//{
// "file": "Scan_1_20190605.zip",
// "inputVariable": "irmi",
// "InputLayerName": "",
// "outputFileName": "vranbs4",
// "plantingDate": "2019-04-18",
// "measurementDate": "2019-07-03",
@@ -178,5 +302,9 @@
// ]
// ]
// }
}
//}
]

23
FarmmapsNbs/NitrogenService.cs
View File

@@ -74,7 +74,7 @@ namespace FarmmapsNbs
var item = await _farmmapsApiService.GetItemAsync(targetNItem.Code);
return item.Data.ToObject<TargetNData>();
}
/// <summary>
/// Calculates the uptake map based on the given inputs
/// </summary>
@@ -82,9 +82,11 @@ namespace FarmmapsNbs
/// <param name="inputItem"></param>
/// <param name="plantingDate">The date the crop is planted</param>
/// <param name="measurementDate">The date the measurements are taken</param>
/// <param name="inputType">Data type, could be yara, ci, irmi or wdvi</param>
/// <param name="inputLayerName">Column name in which the sensor value is stored</param>
/// <returns></returns>
public async Task<Item> CalculateUptakeMap(Item cropfieldItem, Item inputItem, DateTime plantingDate,
DateTime measurementDate, string inputType)
DateTime measurementDate, string inputType, string inputLayerName)
{
var nbsUptakeMapRequest = new TaskRequest {TaskType = VRANBS_TASK};
nbsUptakeMapRequest.attributes["operation"] = "uptake";
@@ -92,10 +94,9 @@ namespace FarmmapsNbs
nbsUptakeMapRequest.attributes["plantingDate"] = plantingDate.ToString("o");
nbsUptakeMapRequest.attributes["measurementDate"] = measurementDate.ToString("o");
nbsUptakeMapRequest.attributes["inputType"] = inputType.ToLower();
nbsUptakeMapRequest.attributes["inputLayerName"] = "IRMI"; //toevoeging FS. Kolom IRMI hernoemd als IMI. Deze wordt niet automatisch herkend. En moet dus gespecificeerd worden.
if (!(string.IsNullOrEmpty(inputLayerName))) nbsUptakeMapRequest.attributes["inputLayerName"] = inputLayerName;
//toevoeging FS. Kolom IRMI hernoemd als IMI. Deze wordt niet automatisch herkend. En moet dus gespecificeerd worden.
//var layers = inputItem.Data["layers"]; //toevoeging FS, check welke data lagen worden omgezet
//_logger.LogInformation($"DataLayers: {layers}"); //toevoeging FS check welke data lagen worden omgezet
string itemTaskCode = await _farmmapsApiService.QueueTaskAsync(cropfieldItem.Code, nbsUptakeMapRequest);
@@ -116,9 +117,10 @@ namespace FarmmapsNbs
var itemName = "VRANbs uptake";
var uptakeMapItem = await _generalService.FindChildItemAsync(cropfieldItem.Code,
GEOTIFF_PROCESSED_ITEMTYPE, itemName,
GEOTIFF_PROCESSED_ITEMTYPE, itemName,
i => i.Updated >= itemTask.Finished.GetValueOrDefault(DateTime.UtcNow) &&
i.Name.ToLower().Contains(itemName.ToLower()));
i.Name.ToLower().Contains(itemName.ToLower()));
if (uptakeMapItem == null)
{
_logger.LogError("Could not find the uptake geotiff child item under cropfield");
@@ -136,11 +138,11 @@ namespace FarmmapsNbs
/// <param name="inputItem">The farmmaps item containing the geotiff data</param>
/// <param name="plantingDate">The date the crop is planted</param>
/// <param name="measurementDate">The date the measurements are taken</param>
/// <param name="inputType">The inputtype to use</param>
/// <param name="inputType">The inputtype to use, could be yara, ci, irmi or wdvi</param>
/// <param name="targetN">The target nitrogen to use for the calculations</param>
/// <returns></returns>
public async Task<Item> CalculateApplicationMap(Item cropfieldItem, Item inputItem, DateTime plantingDate,
DateTime measurementDate, string inputType, double targetN)
DateTime measurementDate, string inputType, double targetN, string inputLayerName)
{
var nbsApplicationMapRequest = new TaskRequest {TaskType = VRANBS_TASK};
nbsApplicationMapRequest.attributes["operation"] = "application";
@@ -150,7 +152,8 @@ namespace FarmmapsNbs
nbsApplicationMapRequest.attributes["inputCode"] = inputItem.Code;
nbsApplicationMapRequest.attributes["inputType"] = inputType.ToLower();
nbsApplicationMapRequest.attributes["targetN"] = targetN.ToString(CultureInfo.InvariantCulture);
if (!(string.IsNullOrEmpty(inputLayerName))) nbsApplicationMapRequest.attributes["inputLayerName"] = inputLayerName;
string itemTaskCode = await _farmmapsApiService.QueueTaskAsync(cropfieldItem.Code, nbsApplicationMapRequest);
await PollTask(TimeSpan.FromSeconds(5), async (tokenSource) =>