Filtern und Indizieren von GIS-Vektorlayern in C#

Mit Aspose.GIS für .NET können Sie Layer nach Attributwerten oder räumlichen Grenzen filtern. Sie können auch Indizes verwenden, um das Filtern und die räumlichen Abfragen zu beschleunigen.

Attributindex

Filtern ohne Index

Hier erfahren Sie, wie Sie einen Layer nach Werten eines Attributs filtern:

// 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();
}

Filtern mit Index

Der obige Code ist in Ordnung, solange der Layer nur einmal gefiltert wird. Wenn der Layer jedoch wahrscheinlich mehrmals abgefragt wird, können wir von Attributindizes profitieren. Der Aufbau eines Attributindex dauert einige Zeit, kann aber mehrfach wiederverwendet werden, um das Filtern zu beschleunigen.

Das folgende Beispiel verwendet eine Attributindexdatei, um das Filtern von Layern nach Werten des Attributs zu beschleunigen:

// 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();
}

Gefilterte Features speichern

Gefilterte Features können in Layern gespeichert werden:

// 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);
}

Gefilterte Features rendern

Es ist auch möglich, gefilterte Features zu rendern. Das folgende Beispiel verwendet einen Attributindex, um schnell alle Features mit einer Bevölkerung von mehr als 2000 auszuwählen und sie der Karte hinzuzufügen:

// 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);
}

Räumlicher Index

Räumliche Indizes werden verwendet, um räumliche Abfragen zu beschleunigen. Genau wie Attributindizes werden räumliche Indizes nach der Erstellung wiederverwendet.

Features in der Nähe eines Punktes finden

Hier erfahren Sie, wie Sie einen räumlichen Index verwenden, um die Suche nach dem Feature in der Nähe eines bestimmten Punktes zu beschleunigen:

// 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();
}

Features auswählen, die sich mit einer Geometrie schneiden

Das folgende Beispiel verwendet einen räumlichen Index, um die Auswahl von Features zu beschleunigen, die sich mit einer Geometrie schneiden:

// 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();
}