private bool AddResultToUserSession(VisibilityCalcResults visibilityResults)
{
try
{
if (!context.MilSp_VisibilityUserSessions.Any(
r => r.userName == visibilityResults.UserName && r.visibilityResultId == visibilityResults.Id))
{
var sessionValue = new MilSp_VisibilityUserSession
{
userName = visibilityResults.UserName,
visibilityResultId = visibilityResults.Id
};
context.MilSp_VisibilityUserSessions.InsertOnSubmit(sessionValue);
Submit();
}
return(true);
}
catch (Exception ex)
{
log.WarnEx($"Unexpected exception:{ex.Message}");
}
return(false);
}
public VisibilityCalcResults UpdateVisibilityResults(VisibilityCalcResults visibilityResults)
{
try
{
var resultsEntity = context.MilSp_VisiblityResults.FirstOrDefault(res => res.Id.Trim() == visibilityResults.Id);
if (resultsEntity != null)
{
resultsEntity.Update(visibilityResults);
Submit();
log.InfoEx($"Visibility results {visibilityResults.Id} was successfully updated");
return(context.MilSp_VisiblityResults.First(session => session.Id.Trim() == visibilityResults.Id).Get());
}
log.WarnEx($"Visibility results {visibilityResults.Id} not found");
}
catch (MilSpaceDataException ex)
{
log.WarnEx(ex.Message);
if (ex.InnerException != null)
{
log.WarnEx(ex.InnerException.Message);
}
}
catch (Exception ex)
{
log.WarnEx($"Unexpected exception:{ex.Message}");
}
return(null);;
}
public static VisibilityCalcResults UpdateVisibilityResults(VisibilityCalcResults visibilityResults)
{
using (var accessor = new VisibilityDataAccess())
{
var res = accessor.UpdateVisibilityResults(visibilityResults);
return(res);
}
}
public static bool AddSharedVisibilityResultsToUserSession(VisibilityCalcResults visibilityResults)
{
using (var accessor = new VisibilityDataAccess())
{
var res = accessor.AddShareddResultsToUserSession(visibilityResults);
return(res);
}
}
private static KeyValuePair <ILayer, LayerTypeEnum> GetVisibilityLayer(VisibilityResultInfo info, IObjectFactory objFactory, IFunctionalSurface baseSurface)
{
KeyValuePair <ILayer, LayerTypeEnum> layerKeyValuePair = new KeyValuePair <ILayer, LayerTypeEnum>();
Type factoryType = Type.GetTypeFromProgID("esriDataSourcesGDB.FileGDBWorkspaceFactory");
string typeFactoryID = factoryType.GUID.ToString("B");
IWorkspaceFactory workspaceFactory = (IWorkspaceFactory)objFactory.Create(typeFactoryID);
IWorkspace2 workspace = (IWorkspace2)workspaceFactory.OpenFromFile(info.GdbPath, 0);
var rastersTypes = VisibilityCalcResults.GetRasterResults();
if (rastersTypes.Any(type => type == info.RessutType))
{
var rasterLayer = CreateRasterLayer(info.ResultName, workspace, objFactory, info.GdbPath);
if (rasterLayer != null)
{
SetVisibilitySessionRaster3DProperties(rasterLayer, objFactory, baseSurface);
layerKeyValuePair = new KeyValuePair <ILayer, LayerTypeEnum>(rasterLayer, LayerTypeEnum.Raster);
}
}
if (info.RessutType == VisibilityCalculationResultsEnum.ObservationPoints || info.RessutType == VisibilityCalculationResultsEnum.ObservationPointSingle)
{
var pointFeatureLayer = CreateFeatureLayer(info.ResultName, workspace, objFactory);
if (pointFeatureLayer != null)
{
SetFeatures3DProperties(pointFeatureLayer, objFactory, baseSurface);
layerKeyValuePair = new KeyValuePair <ILayer, LayerTypeEnum>(pointFeatureLayer, LayerTypeEnum.PointFeature);
}
}
if (info.RessutType == VisibilityCalculationResultsEnum.VisibilityAreaPolygons || info.RessutType == VisibilityCalculationResultsEnum.ObservationObjects)
{
var polygonFeatureLayer = CreateFeatureLayer(info.ResultName, workspace, objFactory);
if (polygonFeatureLayer != null)
{
SetFeatures3DProperties(polygonFeatureLayer, objFactory, baseSurface);
layerKeyValuePair = new KeyValuePair <ILayer, LayerTypeEnum>(polygonFeatureLayer, LayerTypeEnum.PolygonFeature);
}
}
Marshal.ReleaseComObject(workspaceFactory);
return(layerKeyValuePair);
}
public VisibilityCalcResults AddVisibilityResults(VisibilityCalcResults visibilityResults)
{
try
{
if (!context.MilSp_VisiblityResults.Any(res => res.Id == visibilityResults.Id))
{
var resultsEntity = visibilityResults.Get();
context.MilSp_VisiblityResults.InsertOnSubmit(resultsEntity);
if (!AddResultToUserSession(visibilityResults))
{
throw new MilSpaceDataException("VisibilityUserSession", DataOperationsEnum.Insert);
}
Submit();
log.InfoEx($"Session {visibilityResults.Id} was successfully added");
}
else
{
log.WarnEx($"Session {visibilityResults.Id} cannot added because of it already exists");
}
return(context.MilSp_VisiblityResults.First(session => session.Id == visibilityResults.Id).Get());
}
catch (MilSpaceDataException ex)
{
log.WarnEx(ex.Message);
if (ex.InnerException != null)
{
log.WarnEx(ex.InnerException.Message);
}
}
catch (Exception ex)
{
log.WarnEx($"Unexpected exception:{ex.Message}");
}
return(null);
}
public bool AddShareddResultsToUserSession(VisibilityCalcResults visibilityResults)
{
visibilityResults.UserName = Environment.UserName;
return(AddResultToUserSession(visibilityResults));
}
private IEnumerable <string> ProcessObservationPoint(List <string> results)
{
logger.InfoEx($"> ProcessObservationPoint START with parameters {calcResults}");
Dictionary <int, List <string> > generatedResults = null;
if (!showAllResults)
{
logger.InfoEx("Show all results was not set");
generatedResults = new Dictionary <int, List <string> >();
}
int iStepNum = 1;
IEnumerable <string> messages = null;
string oservStationsFeatureClassName = null;
//Used to clip the base immge to make it smaller and reduce the culculation time
int currentPointIdForBestParamsCalculations = pointsFilteringIds.Last();
//Handle Observation Objects
if (calcResults.HasFlag(VisibilityCalculationResultsEnum.ObservationObjects))
{
oservStationsFeatureClassName = VisibilityTask.GetResultName(VisibilityCalculationResultsEnum.ObservationObjects, outputSourceName);
logger.InfoEx($"Export observation objects to {oservStationsFeatureClassName}");
//Add result To Delete if not to show all results. Add OP result
if (!showAllResults)
{
//pointResults =
generatedResults.Add(-1, new List <string>()
{
oservStationsFeatureClassName
});
}
var exportedFeatureClass = GdbAccess.Instance.ExportObservationFeatureClass(
observStationsfeatureClass as IDataset,
oservStationsFeatureClassName,
stationsFilteringIds);
if (!string.IsNullOrWhiteSpace(exportedFeatureClass))
{
logger.InfoEx($"Observation points feature class was created:{exportedFeatureClass}");
results.Add(iStepNum.ToString() + ". " + "Створено копію ОН для розрахунку: " + exportedFeatureClass);
iStepNum++;
}
else
{
string errorMessage = $"The feature calss {oservStationsFeatureClassName} was not exported";
result.ErrorMessage = errorMessage;
logger.ErrorEx("> ProcessObservationPoint ERROR ExportObservationFeatureClass (1). errorMessage:{0}", errorMessage);
messages.Append(errorMessage);
results.Add("Помилка: " + errorMessage);
return(messages);
}
}
else if (calcResults.HasFlag(VisibilityCalculationResultsEnum.BestParametersTable))
{
oservStationsFeatureClassName = observStationsfeatureClass.AliasName;
logger.InfoEx($"Export observation objects to {oservStationsFeatureClassName}");
}
//Handle Observation Points
List <KeyValuePair <VisibilityCalculationResultsEnum, int[]> > pointsIDs =
new List <KeyValuePair <VisibilityCalculationResultsEnum, int[]> >();
if (calcResults.HasFlag(VisibilityCalculationResultsEnum.VisibilityAreaRaster) &&
!calcResults.HasFlag(VisibilityCalculationResultsEnum.BestParametersTable))
{
pointsIDs.Add(
new KeyValuePair <VisibilityCalculationResultsEnum, int[]>(
VisibilityCalculationResultsEnum.ObservationPoints,
pointsFilteringIds));
}
if (calcResults.HasFlag(VisibilityCalculationResultsEnum.ObservationPointSingle))
{
if (pointsFilteringIds.Length > 1)
{
pointsIDs.AddRange(
pointsFilteringIds.Select(
id => new KeyValuePair <VisibilityCalculationResultsEnum, int[]>(
VisibilityCalculationResultsEnum.ObservationPointSingle, new int[] { id }))
.ToArray());
}
else
{
calcResults &= ~VisibilityCalculationResultsEnum.ObservationPointSingle;
calcResults &= ~VisibilityCalculationResultsEnum.VisibilityAreaRasterSingle;
}
}
if (calcResults.HasFlag(VisibilityCalculationResultsEnum.BestParametersTable))
{
pointsIDs.AddRange(
pointsFilteringIds.Select(
id => new KeyValuePair <VisibilityCalculationResultsEnum, int[]>(
VisibilityCalculationResultsEnum.ObservationPointSingle, new int[] { id }))
.ToArray());
}
int index = -1;
bool removeFullImageFromresult = false;
int[] objIds = null;
if (calcResults.HasFlag(VisibilityCalculationResultsEnum.ObservationObjects))
{
objIds = stationsFilteringIds;
}
var coverageTableManager = new CoverageTableManager();
BestOPParametersManager bestOPParametersManager = null;
if (calcResults.HasFlag(VisibilityCalculationResultsEnum.BestParametersTable))
{
bestOPParametersManager = new BestOPParametersManager(visibilityPercent, currentPointIdForBestParamsCalculations, showAllResults);
}
// There is no reason to clip image for the BestParameters mode, since the point just shance its Z value
var clipSourceImageForEveryPoint = true;
//Used to define a result layer name;
int pointIndex = -1;
var orifinalResterSource = rasterSource;
foreach (var curPoints in pointsIDs)
{
int pointId;
if (calcResults.HasFlag(VisibilityCalculationResultsEnum.BestParametersTable))
{
pointId = currentPointIdForBestParamsCalculations;
pointIndex = pointId - 1;
}
else
{
pointId = curPoints.Key == VisibilityCalculationResultsEnum.ObservationPoints ? -1 : ++index;
pointIndex = pointId;
}
logger.InfoEx($"Start processing point id {pointId} with index {pointIndex}");
if (!showAllResults)// Cretae the list to add datasets to delete
{
generatedResults.Add(pointId, new List <string>());
}
int[] curPointsValue = (calcResults.HasFlag(VisibilityCalculationResultsEnum.BestParametersTable)) ?
new int[] { pointId } : curPoints.Value;
var observPointFeatureClassName = VisibilityTask.GetResultName(curPoints.Key, outputSourceName, pointIndex);
logger.InfoEx($"Processing point(s) into {observPointFeatureClassName}");
//Add result To Delete if not to show all results. Add OP result
if (!showAllResults && generatedResults.ContainsKey(pointId))
{
generatedResults[pointId].Add(observPointFeatureClassName);
}
var exportedFeatureClass = GdbAccess.Instance.ExportObservationFeatureClass(
observPointsfeatureClass as IDataset,
observPointFeatureClassName,
curPointsValue);
results.Add(iStepNum.ToString() + ". " + "Створено копію ПС для розрахунку: " + exportedFeatureClass);
logger.InfoEx($"Point(s) was exported into {observPointFeatureClassName} to process visibility");
iStepNum++;
if (string.IsNullOrWhiteSpace(exportedFeatureClass))
{
string errorMessage = $"The feature calss {observPointFeatureClassName} was not exported";
result.Exception = new MilSpaceVisibilityCalcFailedException(errorMessage);
result.ErrorMessage = errorMessage;
logger.ErrorEx("> ProcessObservationPoint ERROR ExportObservationFeatureClass. errorMessage:{0}", errorMessage);
results.Add("Помилка: " + errorMessage + " ПС: " + pointId.ToString());
return(messages);
}
bool clipSoutceImageByPotentialArea = clipSourceImageForEveryPoint && (curPoints.Key == VisibilityCalculationResultsEnum.ObservationPoints || curPoints.Key == VisibilityCalculationResultsEnum.ObservationPointSingle);
if (clipSoutceImageByPotentialArea)
{
logger.InfoEx($"Clipping calc image to potential area");
string featureClassName = exportedFeatureClass;
if (calcResults.HasFlag(VisibilityCalculationResultsEnum.BestParametersTable) && curPoints.Value.First() == 1)
{
featureClassName = observPointsfeatureClass.AliasName;
}
logger.InfoEx($"Calculate potential area for observation point in {featureClassName} and objects from {oservStationsFeatureClassName}");
coverageTableManager.SetCalculateAreas(featureClassName, oservStationsFeatureClassName, !calcResults.HasFlag(VisibilityCalculationResultsEnum.BestParametersTable));
clipSourceImageForEveryPoint = false;// !calcResults.HasFlag(VisibilityCalculationResultsEnum.BestParametersTable);
var visibilityPotentialAreaFCName = VisibilityCalcResults.GetResultName(VisibilityCalculationResultsEnum.VisibilityAreasPotential, outputSourceName);
coverageTableManager.AddPotentialArea(
visibilityPotentialAreaFCName,
(curPoints.Key == VisibilityCalculationResultsEnum.ObservationPoints ||
curPoints.Key == VisibilityCalculationResultsEnum.ObservationPointSingle), pointIndex + 1);
results.Add(iStepNum.ToString() + ". " + "Розраховано потенційне покриття: " + visibilityPotentialAreaFCName + " ПС: " + pointId.ToString());
//Calc protentilly visiblw area as Image
var fc = GenerateWorknArea(visibilityPotentialAreaFCName, observPointFeatureClassName, outputSourceName);
results.Add(iStepNum.ToString() + ". " + "Розраховано потенційне покриття для розрахунку: " + fc.AliasName + " ПС: " + pointId.ToString());
//Try to clip the raster source by visibilityPotentialAreaFCName
var visibilityPotentialAreaImgName =
VisibilityCalcResults.GetResultName(pointIndex > -1 ?
VisibilityCalculationResultsEnum.VisibilityRastertPotentialAreaSingle :
VisibilityCalculationResultsEnum.VisibilityRastertPotentialArea,
outputSourceName, pointIndex);
//Add result To Delete if not to show all results. Add OP result
if (!showAllResults && generatedResults.ContainsKey(pointId))
{
generatedResults[pointId].Add(visibilityPotentialAreaImgName);
}
if (!CalculationLibrary.ClipVisibilityZonesByAreas(
rasterSource,
visibilityPotentialAreaImgName,
fc.AliasName,
out messages, "NONE"))
{
string errorMessage = $"The result {visibilityPotentialAreaImgName} was not generated";
result.Exception = new MilSpaceVisibilityCalcFailedException(errorMessage);
result.ErrorMessage = errorMessage;
logger.ErrorEx("> ProcessObservationPoint ERROR ClipVisibilityZonesByAreas. errorMessage:{0}", errorMessage);
results.Add("Помилка: " + errorMessage + " ПС: " + pointId.ToString());
return(messages);
}
else
{
results.Add(iStepNum.ToString() + ". " + "Розраховано покриття для розрахунку на основі потенційноі видимості: " + visibilityPotentialAreaImgName + " ПС: " + pointId.ToString());
rasterSource = visibilityPotentialAreaImgName;
//Delete temporary Featureclass usewd for clipping the base image by potential area
GdbAccess.Instance.RemoveFeatureClass(fc.AliasName);
}
}
//Generate Visibility Raster
string featureClass = observPointFeatureClassName;
string outImageName = VisibilityTask.GetResultName(
curPoints.Key == VisibilityCalculationResultsEnum.ObservationPoints ?
VisibilityCalculationResultsEnum.VisibilityAreaRaster :
VisibilityCalculationResultsEnum.VisibilityAreaRasterSingle,
outputSourceName,
pointIndex);
//Add result To Delete if not to show all results. Add Visibility raster result
if (!showAllResults && generatedResults.ContainsKey(pointId))
{
generatedResults[pointId].Add(outImageName);
}
if (!CalculationLibrary.GenerateVisibilityData(
rasterSource,
featureClass,
VisibilityAnalysisTypesEnum.Frequency,
outImageName,
out messages))
{
string errorMessage = $"The result {outImageName} was not generated";
result.Exception = new MilSpaceVisibilityCalcFailedException(errorMessage);
result.ErrorMessage = errorMessage;
logger.ErrorEx("> ProcessObservationPoint ERROR ConvertRasterToPolygon. errorMessage:{0}", errorMessage);
logger.ErrorEx("> ProcessObservationPoint ERROR ConvertRasterToPolygon. errorMessage:{0}", errorMessage);
results.Add("Помилка: " + errorMessage + " ПС: " + pointId.ToString());
return(messages);
}
else
{
results.Add(iStepNum.ToString() + ". " + "Розраховано видимість: " + outImageName + " ПС: " + pointId.ToString());
iStepNum++;
logger.InfoEx($"Visibility image {outImageName} for point {pointId} was generated");
string visibilityArePolyFCName = null;
//ConvertToPolygon full visibility area
if (calcResults.HasFlag(VisibilityCalculationResultsEnum.VisibilityAreaPolygons) &&
!calcResults.HasFlag(VisibilityCalculationResultsEnum.ObservationObjects))
{
logger.InfoEx($"Calculation of {VisibilityCalculationResultsEnum.VisibilityAreaPolygons} (Convert visible areas to polygon) was swithced on");
visibilityArePolyFCName =
VisibilityTask.GetResultName(pointIndex > -1 ?
VisibilityCalculationResultsEnum.VisibilityAreaPolygonSingle :
VisibilityCalculationResultsEnum.VisibilityAreaPolygons, outputSourceName, pointIndex);
//Add result To Delete if not to show all results. Add Visibility polygon result
if (!showAllResults && generatedResults.ContainsKey(pointId))
{
generatedResults[pointId].Add(visibilityArePolyFCName);
}
if (!CalculationLibrary.ConvertRasterToPolygon(outImageName, visibilityArePolyFCName, out messages))
{
if (!messages.Any(m => m.StartsWith("ERROR 010151"))) // Observatioj areas dont intersect Visibility aresa
{
string errorMessage = $"The result {visibilityArePolyFCName} was not generated";
result.Exception = new MilSpaceVisibilityCalcFailedException(errorMessage);
result.ErrorMessage = errorMessage;
logger.ErrorEx("> ProcessObservationPoint ERROR ConvertRasterToPolygon. errorMessage:{0}", errorMessage);
results.Add("Помилка: " + errorMessage + " ПС: " + pointId.ToString());
return(messages);
}
results.Add(iStepNum.ToString() + ". " + "Видимисть відсутня. Полігони не були конвертовані: " + visibilityArePolyFCName + " ПС: " + pointId.ToString());
visibilityArePolyFCName = string.Empty;
logger.InfoEx($"There is no visibile areas for OP {pointId}");
}
else
{
results.Add(iStepNum.ToString() + ". " + "Конвертовано у полігони: " + visibilityArePolyFCName + " ПС: " + pointId.ToString());
logger.InfoEx($"Visibile areas for OP {pointId} was converted into {visibilityArePolyFCName} ");
}
iStepNum++;
}
//Clip
if (!string.IsNullOrEmpty(oservStationsFeatureClassName))
{
var inClipName = outImageName;
var resultLype =
pointId > -1 ?
VisibilityCalculationResultsEnum.VisibilityObservStationClipSingle :
VisibilityCalculationResultsEnum.VisibilityObservStationClip;
var outClipName = VisibilityTask.GetResultName(resultLype,
outputSourceName, pointIndex);
//Add result To Delete if not to show all results. Add Visibility polygon result
if (!showAllResults && generatedResults.ContainsKey(pointId))
{
generatedResults[pointId].Add(outClipName);
}
logger.InfoEx($"Try to clip the visible areas {inClipName} by orservation object {oservStationsFeatureClassName}");
if (!CalculationLibrary.ClipVisibilityZonesByAreas(
inClipName,
outClipName,
oservStationsFeatureClassName,
out messages))
{
string errorMessage = $"The result {outClipName} was not generated";
result.Exception = new MilSpaceVisibilityCalcFailedException(errorMessage);
result.ErrorMessage = errorMessage;
logger.ErrorEx("> ProcessObservationPoint ERROR ClipVisibilityZonesByAreas. errorMessage:{0}", errorMessage);
results.Add("Помилка: " + errorMessage + " ПС: " + pointId.ToString());
return(messages);
}
else
{
results.Add(iStepNum.ToString() + ". " + "Зона видимості зведена до дійсного розміру: " + outClipName + " ПС: " + pointId.ToString());
logger.InfoEx($"Visible area {inClipName} was clipped by orservation objects {oservStationsFeatureClassName}");
iStepNum++;
if (!calcResults.HasFlag(resultLype))
{
calcResults |= resultLype;
}
//Change to VisibilityAreaPolygonForObjects
var curCulcRResult =
pointId > -1 ?
VisibilityCalculationResultsEnum.VisibilityAreaPolygonSingle :
VisibilityCalculationResultsEnum.VisibilityAreaPolygons;
visibilityArePolyFCName = VisibilityTask.GetResultName(curCulcRResult, outputSourceName, pointIndex);
var rasterDataset = GdbAccess.Instance.GetDatasetFromCalcWorkspace(
outClipName, VisibilityCalculationResultsEnum.VisibilityAreaRaster);
bool isEmpty = EsriTools.IsRasterEmpty((IRasterDataset2)rasterDataset);
if (isEmpty)
{
if (calcResults.HasFlag(curCulcRResult))
{
calcResults &= ~curCulcRResult;
}
results.Add($"{iStepNum.ToString()}. Видимість відсутня. Полігони {visibilityArePolyFCName} не було сформовано. ПС: {pointId.ToString()}");
logger.InfoEx($"There is no visible areas in {visibilityArePolyFCName}. Point {pointId}");
}
else
{
if (!CalculationLibrary.ConvertRasterToPolygon(outClipName, visibilityArePolyFCName, out messages))
{
if (!messages.Any(m => m.StartsWith("ERROR 010151"))) // Observatioj areas dont intersect Visibility area
{
string errorMessage = $"The result {visibilityArePolyFCName} was not generated";
result.Exception = new MilSpaceVisibilityCalcFailedException(errorMessage);
result.ErrorMessage = errorMessage;
logger.ErrorEx("> ProcessObservationPoint ERROR ConvertRasterToPolygon. errorMessage:{0}", errorMessage);
results.Add("Помилка: " + errorMessage + " ПС: " + pointId.ToString());
return(messages);
}
results.Add(iStepNum.ToString() + ". " + "Видимість відсутня. Полігони не було сформовано: " + visibilityArePolyFCName + " ПС: " + pointId.ToString());
logger.InfoEx($"There is no visible areas. {visibilityArePolyFCName} was not generated. Point {pointId}");
}
else
{
//Add result To Delete if not to show all results. Add Visibility polygon result
if (!showAllResults && generatedResults.ContainsKey(pointId))
{
generatedResults[pointId].Add(visibilityArePolyFCName);
}
results.Add(iStepNum.ToString() + ". " + "Конвертовано у полігони: " + visibilityArePolyFCName + " ПС: " + pointId.ToString());
}
}
iStepNum++;
}
}
else if (calcResults.HasFlag(VisibilityCalculationResultsEnum.VisibilityAreasTrimmedByPoly) &&
!string.IsNullOrEmpty(visibilityArePolyFCName))
{
//Clip visibility images to valuable extent
var inClipName = outImageName;
var outClipName = VisibilityTask.GetResultName(pointIndex > -1 ?
VisibilityCalculationResultsEnum.VisibilityAreaTrimmedByPolySingle :
VisibilityCalculationResultsEnum.VisibilityAreasTrimmedByPoly, outputSourceName, pointIndex);
if (!CalculationLibrary.ClipVisibilityZonesByAreas(
inClipName,
outClipName,
visibilityArePolyFCName,
out messages,
"NONE"))
{
string errorMessage = $"The result {outClipName} was not generated";
result.Exception = new MilSpaceVisibilityCalcFailedException(errorMessage);
result.ErrorMessage = errorMessage;
logger.ErrorEx("> ProcessObservationPoint ERROR ClipVisibilityZonesByAreas. errorMessage:{0}", errorMessage);
results.Add("Помилка: " + errorMessage + " ПС: " + pointId.ToString());
return(messages);
}
else
{
//Add result To Delete if not to show all results. Add Visibility polygon result
if (!showAllResults && generatedResults.ContainsKey(pointId))
{
generatedResults[pointId].Add(outClipName);
}
results.Add(iStepNum.ToString() + ". " + "Зона видимості зведена до дійсного розміру: " + outClipName + " ПС: " + pointId.ToString());
iStepNum++;
removeFullImageFromresult = true;
}
}
if (calcResults.HasFlag(VisibilityCalculationResultsEnum.CoverageTable))
{
var pointsCount = pointsFilteringIds.Where(id => id > -1).Count();
coverageTableManager.CalculateCoverageTableDataForPoint(
(pointId == -1),
visibilityArePolyFCName,
pointsCount,
exportedFeatureClass,
pointIndex + 1);
results.Add(iStepNum.ToString() + ". " + "Сформовані записи таблиці покриття. для ПС: " + pointId.ToString());
iStepNum++;
}
if (calcResults.HasFlag(VisibilityCalculationResultsEnum.BestParametersTable))
{
if (pointId != -1)
{
currentPointIdForBestParamsCalculations = bestOPParametersManager.FindVisibilityPercent(visibilityArePolyFCName, observStationsfeatureClass,
stationsFilteringIds, pointId);
results.Add(iStepNum.ToString() + ". " + "Знайдено відсоток покриття для кроку " + pointId.ToString());
iStepNum++;
if (currentPointIdForBestParamsCalculations < 0)
{
break;
}
}
}
}
}
if (calcResults.HasFlag(VisibilityCalculationResultsEnum.CoverageTable))
{
var coverageTable = VisibilityTask.GetResultName(VisibilityCalculationResultsEnum.CoverageTable, outputSourceName);
coverageTableManager.SaveDataToCoverageTable(coverageTable);
results.Add(iStepNum.ToString() + ". " + "Збережена загальна таблиця покриття: " + coverageTable);
iStepNum++;
}
if (calcResults.HasFlag(VisibilityCalculationResultsEnum.BestParametersTable))
{
var bestParamsTableName = VisibilityTask.GetResultName(VisibilityCalculationResultsEnum.BestParametersTable, outputSourceName);
if (bestOPParametersManager.CreateVOTable(observPointsfeatureClass,
visibilityPercent, bestParamsTableName))
{
//Delete unnecessary datasets
if (!showAllResults && bestOPParametersManager.AppropriateedParams != null)
{
foreach (var key in generatedResults.Keys)
{
if (!bestOPParametersManager.AppropriateedParams.Any(k => k.Key == key))
{
foreach (var fc in generatedResults[key])
{
//try to remove feature class
if (!GdbAccess.Instance.RemoveFeatureClass(fc))
{
// it is not FC try to remove raster
GdbAccess.Instance.RemoveRasterDataset(fc);
}
}
}
}
}
results.Add(iStepNum.ToString() + ". " + "Збережена таблиця накращих параметрів ПН для мінімальної видимості " + visibilityPercent + "%: " + bestParamsTableName);
iStepNum++;
}
else
{
results.Add(iStepNum.ToString() + ". " + "Відсутні результати для видимості " + visibilityPercent +
"%. Збережена таблиця з накращими можливими параметрами ПН: " + bestParamsTableName);
iStepNum++;
}
}
//Set real results to show
if (removeFullImageFromresult)
{
if (calcResults.HasFlag(VisibilityCalculationResultsEnum.VisibilityAreaRaster))
{
calcResults &= ~VisibilityCalculationResultsEnum.VisibilityAreaRaster;
}
if (calcResults.HasFlag(VisibilityCalculationResultsEnum.VisibilityAreaRasterSingle))
{
calcResults &= ~VisibilityCalculationResultsEnum.VisibilityAreaRasterSingle;
}
}
session.CalculatedResults = (int)calcResults;
logger.InfoEx("> ProcessObservationPoint END");
return(messages);
}
public static IEnumerable <string> GetCalculatedObserObjectsName(string resultId)
{
var calcObservObjectsFeatureClass = VisibilityCalcResults.GetResultName(VisibilityCalculationResultsEnum.ObservationObjects, resultId);
return(GdbAccess.Instance.GetCalcEntityNamesFromFeatureClass(calcObservObjectsFeatureClass, "sTitleOO"));
}
public static IFeatureClass CreateOPFeatureClass(WizardResult calcResult, IFeatureClass observatioPointsFeatureClass,
IActiveView activeView, IRaster raster)
{
var calcObservPointsFeatureClass = VisibilityCalcResults.GetResultName(VisibilityCalculationResultsEnum.ObservationPoints, calcResult.TaskName);
var observPointTemporaryFeatureClass =
GdbAccess.Instance.GenerateTemporaryFeatureClassWithRequitedFields(observatioPointsFeatureClass.Fields, calcObservPointsFeatureClass);
bool isCircle;
double maxAzimuth = 0;
double minAzimuth = 360;
double maxDistance = 0;
calcResult.ObservationStation.Project(activeView.FocusMap.SpatialReference);
var observationStationPolygon = calcResult.ObservationStation as IPolygon;
var observStationEnvelope = observationStationPolygon.Envelope;
var observStationEnvelopePoints = new IPoint[] { observStationEnvelope.LowerLeft, observStationEnvelope.LowerRight,
observStationEnvelope.UpperRight, observStationEnvelope.UpperLeft };
var observerPointGeometry = new PointClass
{
X = calcResult.ObservationPoint.X.Value,
Y = calcResult.ObservationPoint.Y.Value,
SpatialReference = EsriTools.Wgs84Spatialreference
};
observerPointGeometry.Project(activeView.FocusMap.SpatialReference);
observerPointGeometry.AddZCoordinate(raster);
observerPointGeometry.ZAware = true;
if (double.IsNaN(observerPointGeometry.Z))
{
throw new MilSpacePointOutOfRatserException(calcResult.ObservationPoint.Objectid, calcResult.RasterLayerName);
}
var minDistance = FindMinDistance(observStationEnvelopePoints, observerPointGeometry);
// ---- Get min and max azimuths ----
var points = EsriTools.GetPointsFromGeometries(new IGeometry[] { calcResult.ObservationStation },
observerPointGeometry.SpatialReference,
out isCircle).ToArray();
bool isPointInside = EsriTools.IsPointOnExtent(observStationEnvelope, observerPointGeometry);
// Set azimuth for circle polygon
if (isCircle && !isPointInside)
{
for (int i = 0; i < observStationEnvelopePoints.Length; i++)
{
var line = new Line()
{
FromPoint = observerPointGeometry,
ToPoint = observStationEnvelopePoints[i],
SpatialReference = observerPointGeometry.SpatialReference
};
if (i == 0)
{
maxDistance = line.Length;
}
else if (maxDistance < line.Length)
{
maxDistance = line.Length;
}
if (minAzimuth > line.PosAzimuth())
{
minAzimuth = line.PosAzimuth();
}
if (maxAzimuth < line.PosAzimuth())
{
maxAzimuth = line.PosAzimuth();
}
}
}
else
{
EsriTools.CreateDefaultPolylinesForFun(observerPointGeometry, points, new IGeometry[] { calcResult.ObservationStation },
isCircle, isPointInside, -1, out minAzimuth, out maxAzimuth, out maxDistance).ToList();
}
// ---- End. Get min and max azimuths ----
// observerPointGeometry.Project(activeView.FocusMap.SpatialReference);
var pointCopy = observerPointGeometry.CloneWithProjecting();
(pointCopy as PointClass).ZAware = true;
int idObserPooint = 0;
for (var currentHeight = calcResult.FromHeight; currentHeight <= calcResult.ToHeight; currentHeight += calcResult.Step)
{
double maxTiltAngle = -90;
double minTiltAngle = 90;
var height = currentHeight;
// Set parameters for case if point located into the polygon
if (isCircle && isPointInside)
{
minTiltAngle = -90;
minDistance = 0;
}
for (int i = 0; i < observStationEnvelopePoints.Length; i++)
{
var currentTitleAngle = EsriTools.FindAngleByDistanceAndHeight(height, observerPointGeometry, observStationEnvelopePoints[i], raster);
if (minTiltAngle > currentTitleAngle)
{
minTiltAngle = currentTitleAngle;
}
if (maxTiltAngle < currentTitleAngle)
{
maxTiltAngle = currentTitleAngle;
}
}
// Create observation point copy with changing height, distance and angles values
ObservationPoint currentObservationPoint = calcResult.ObservationPoint.ShallowCopy();
currentObservationPoint.RelativeHeight = currentHeight;
currentObservationPoint.AngelMinH = minTiltAngle;
currentObservationPoint.AngelMaxH = maxTiltAngle;
currentObservationPoint.AzimuthStart = minAzimuth;
currentObservationPoint.AzimuthEnd = maxAzimuth;
currentObservationPoint.InnerRadius = minDistance;
currentObservationPoint.OuterRadius = maxDistance;
currentObservationPoint.AzimuthMainAxis = calcResult.ObservationPoint.AzimuthMainAxis;
currentObservationPoint.Id = idObserPooint.ToString();
GdbAccess.Instance.AddObservPoint(pointCopy, currentObservationPoint, observPointTemporaryFeatureClass);
idObserPooint++;
}
return(observPointTemporaryFeatureClass);
}