Филтриране и индексиране на GIS векторни слоеве в C#

С Aspose.GIS за .NET можете да филтрирате слоеве по стойности на атрибути или пространствени граници. Можете също така да използвате индекси, за да ускорите филтрирането и пространствените заявки.

Атрибутен индекс

Филтриране без индекс

Ето как да филтрирате слой по стойности на атрибут:

// For complete examples and data files, please go to https://github.com/aspose-gis/Aspose.GIS-for-.NET
using (var layer = VectorLayer.Open(citiesPath, Drivers.GeoJson))
{
    // Select features based on values of the attribute. This code enumerates all features in the layer
    // and selects all features that match a condition.
    var features = layer.WhereGreaterOrEqual("population", 2000).WhereSmaller("population", 5000);

    // Print results.
    Console.WriteLine("Cities where population >= 2000 and population < 5000:");
    foreach (var feature in features)
    {
        var name = feature.GetValue<string>("name");
        var population = feature.GetValue<int>("population");
        Console.WriteLine(name + ", " + population);
    }
    Console.WriteLine();
}

Филтриране с индекс

Горният код е добър, докато слойът се филтрира само веднъж. Но ако слоят вероятно ще бъде заявен няколко пъти, можем да се възползваме от атрибутните индекси. Отнема известно време за изграждане на атрибутен индекс, но той може да бъде използван многократно за ускоряване на филтрирането.

Следният пример използва файлов атрибутен индекс, за да ускори филтрирането на слой по стойности на атрибута:

// For complete examples and data files, please go to https://github.com/aspose-gis/Aspose.GIS-for-.NET
using (var layer = VectorLayer.Open(citiesPath, Drivers.GeoJson))
{
    // Use attribute index to speed up search by 'population' attribute.
    // Aspose.GIS builds a new index if it doesn't exist, otherwise existing index is reused.
    // Any path can be used.
    var attributeIndexPath = Path.ChangeExtension(citiesPath, "population_out.ix");
    layer.UseAttributesIndex(attributeIndexPath, "population");

    // Select features based on values of the attribute. Since we use attribute index it is not necessary to
    // test all features of the layer and filtering time is reduced significantly.
    var towns = layer.WhereGreaterOrEqual("population", 2000).WhereSmaller("population", 5000);

    // Print results.
    Console.WriteLine("Cities where population >= 2000 and population < 5000:");
    foreach (var town in towns)
    {
        var name = town.GetValue<string>("name");
        var population = town.GetValue<int>("population");
        Console.WriteLine(name + ", " + population);
    }
    Console.WriteLine();
}

Запазване на филтрирани характеристики

Филтрираните характеристики могат да бъдат запазени в слоеве:

// For complete examples and data files, please go to https://github.com/aspose-gis/Aspose.GIS-for-.NET
using (var layer = VectorLayer.Open(citiesPath, Drivers.GeoJson))
{
    // Use attribute index to speed up search by 'population' attribute.
    var attributeIndexPath = Path.ChangeExtension(citiesPath, "population_out.ix");
    layer.UseAttributesIndex(attributeIndexPath, "population");

    // Save all features with population between 2000 and 5000 to the output file.
    layer.WhereGreaterOrEqual("population", 2000)
        .WhereSmaller("population", 5000)
        .SaveTo(DataDir + "towns_out.shp", Drivers.Shapefile);
}

Рендиране на филтрирани характеристики

Възможно е също така да се рендират филтрирани характеристики. Следният пример използва атрибутен индекс, за да избере бързо всички характеристики с население над 2000 и да ги добави към картата:

// For complete examples and data files, please go to https://github.com/aspose-gis/Aspose.GIS-for-.NET
using (var map = new Map(600, 400))
using (var cities = VectorLayer.Open(citiesPath, Drivers.GeoJson))
{
    map.Padding = 20;

    // Use attribute index to speed up search by 'population' attribute.
    var attributeIndexPath = Path.ChangeExtension(citiesPath, "population_out.ix");
    cities.UseAttributesIndex(attributeIndexPath, "population");

    // Render all cities with population greater than 2000.
    map.Add(cities.WhereGreater("population", 2000), new SimpleMarker { FillColor = Color.Red });
    map.Render(outputPath, Renderers.Svg);
}

Пространствен индекс

Пространствените индекси се използват за ускоряване на пространствените заявки. Точно както атрибутните индекси, пространствените индекси се преизползват след създаването си.

Намиране на характеристики най-близки до точка

Ето как да използвате пространствен индекс, за да ускорите търсенето на характеристиката, най-близка до определена точка:

// For complete examples and data files, please go to https://github.com/aspose-gis/Aspose.GIS-for-.NET
using (var layer = VectorLayer.Open(path, Drivers.GeoJson))
{
    // Use spatial index to speed up spatial queries.
    // Aspose.GIS builds a new index if it doesn't exist, otherwise existing index is reused.
    // Any path can be used.
    var spatialIndexPath = Path.ChangeExtension(path, ".spatial_out.ix");
    layer.UseSpatialIndex(spatialIndexPath);

    var point = new Point(12.30, 50.33);

    // Since we use spatial index, nearest-to finds the closest feature much faster.
    var nearest = layer.NearestTo(point);
    Console.WriteLine("City nearest to (12.30 50.33) is " + nearest.GetValue<string>("name"));
    Console.WriteLine();
}

Избор на характеристики, пресичащи се с геометрия

Следният пример използва пространствен индекс, за да ускори избора на характеристики, които се пресичат с геометрия:

// For complete examples and data files, please go to https://github.com/aspose-gis/Aspose.GIS-for-.NET
using (var layer = VectorLayer.Open(path, Drivers.GeoJson))
{
    // Use spatial index to speed up 'WhereIntersects'.
    var spatialIndexPath = Path.ChangeExtension(path, ".spatial_out.ix");
    layer.UseSpatialIndex(spatialIndexPath);

    var polygon = Geometry.FromText("Polygon((12.30 50.33, 22.49 54.87, 21.92 42.53, 12.30 50.33))");
    var intersecting = layer.WhereIntersects(polygon);

    Console.WriteLine("Cities within " + polygon.AsText() + ":");
    foreach (var feature in intersecting)
    {
        var name = feature.GetValue<string>("name");
        var location = (IPoint) feature.Geometry;
        Console.WriteLine($"{name} at ({location.X}, {location.Y})");
    }
    Console.WriteLine();
}