// Copyright (C) Stichting Deltares 2016. All rights reserved. // // This file is part of Ringtoets. // // Ringtoets is free software: you can redistribute it and/or modify // it under the terms of the GNU Lesser General Public License as published by // the Free Software Foundation, either version 3 of the License, or // (at your option) any later version. // // This program is distributed in the hope that it will be useful, // but WITHOUT ANY WARRANTY; without even the implied warranty of // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the // GNU Lesser General Public License for more details. // // You should have received a copy of the GNU Lesser General Public License // along with this program. If not, see . // // All names, logos, and references to "Deltares" are registered trademarks of // Stichting Deltares and remain full property of Stichting Deltares at all times. // All rights reserved. using System; using System.Collections.Generic; using System.Linq; using Core.Common.Base.Geometry; using Core.Common.IO.Exceptions; using Core.Common.Utils.Builders; using Core.Components.Gis.Data; using Core.Components.Gis.Features; using Core.Components.Gis.Geometries; using DotSpatial.Data; using DotSpatial.Topology; using CoreCommonUtilsResources = Core.Common.Utils.Properties.Resources; using GisIOResources = Core.Components.Gis.IO.Properties.Resources; namespace Core.Components.Gis.IO.Readers { /// /// Class to be used to read points from a shapefile. /// public class PointShapeFileReader : ShapeFileReaderBase { private int readIndex; /// /// Creates a new instance of . /// /// The path to the shape file. /// Thrown when is invalid. /// Thrown when /// /// points to a file that doesn't exist. /// The shapefile has non-point geometries in it. /// /// public PointShapeFileReader(string filePath) : base(filePath) { try { ShapeFile = new PointShapefile(filePath); } catch (ApplicationException e) { string message = new FileReaderErrorMessageBuilder(filePath) .Build(GisIOResources.PointShapeFileReader_File_contains_geometries_not_points); throw new CriticalFileReadException(message, e); } } public override FeatureBasedMapData ReadShapeFile(string name = null) { List featureList = new List(); while (readIndex != GetNumberOfLines()) { featureList.Add(ReadFeatureLine()); } string mapFeatureName = !string.IsNullOrWhiteSpace(name) ? name : GisIOResources.PointShapeFileReader_ReadLine_Points; return ConvertPointFeaturesToMapPointData(featureList, mapFeatureName); } public override FeatureBasedMapData ReadLine(string name = null) { if (readIndex == GetNumberOfLines()) { return null; } try { IFeature pointFeature = GetFeature(readIndex); string mapFeatureName = !string.IsNullOrWhiteSpace(name) ? name : GisIOResources.PointShapeFileReader_ReadLine_Points; return ConvertPointFeatureToMapPointData(pointFeature, mapFeatureName); } finally { readIndex++; } } public override IFeature GetFeature(int index) { return ShapeFile.Features[index]; } private IFeature ReadFeatureLine() { try { return GetFeature(readIndex); } finally { readIndex++; } } private FeatureBasedMapData ConvertPointFeatureToMapPointData(IFeature pointFeature, string name) { MapFeature feature = CreateMapFeatureForPointFeature(pointFeature); IEnumerable mapFeatures = new List { feature }; return new MapPointData(mapFeatures, name); } private FeatureBasedMapData ConvertPointFeaturesToMapPointData(IEnumerable featureList, string name) { MapFeature[] mapFeatures = featureList.Select(CreateMapFeatureForPointFeature).ToArray(); return new MapPointData(mapFeatures, name); } private static MapFeature CreateMapFeatureForPointFeature(IFeature pointFeature) { IEnumerable mapGeometries = pointFeature.Coordinates.Select(c => new MapGeometry(GetMapGeometryPointCollections(c))); return new MapFeature(mapGeometries); } private static IEnumerable> GetMapGeometryPointCollections(Coordinate c) { return new[] { new[] { new Point2D(c.X, c.Y) } }; } } }