Vector Layer Editing

import os
import globalmapper as gm

layer_data = "Z:/GM-20069/edit_vector_features_by_attribute_sample_data/gfveg.shp"

# This script iterates through all of the area features in a shapefile. It checks for names of vector areas with a specific
# attribute value (GREAT BASIN CONIFER WOODLAND) and display name (COMMUNITY). It changes the area's classification to wood
# and adds a label to that area.

# Closes all layers that may be open.
layers_are_open = gm.GetLoadedLayerList()
if layers_are_open:
    arr_ptr, arr_size = layers_are_open
    initial_layers = gm.GM_LayerHandle_array_frompointer(arr_ptr)
    for i in range(arr_size):
        gm.CloseLayer(initial_layers[i])

err_LoadLayerList, arr_ptr, arr_size = gm.LoadLayerList(
    layer_data, gm.GM_LoadFlags_UseDefaultLoadOpts
)
layer_arr = gm.GM_LayerHandle_array_frompointer(arr_ptr)

### START SAMPLE
i = 0
modified_ct = 0
feature = gm.GetAreaFeature(layer_arr[0], i)

# Iterates over area features
while feature:
    # Converts from C++ pointer to an array
    attr_array = gm.GM_AttrValue_array_frompointer(feature.mFeatureInfo.mAttrList)
    # Converts array to a python list
    attribute_values = gm.carray_to_list(attr_array, feature.mFeatureInfo.mNumAttrs)

    # Iterates through python list of attribute names and attribute values
    for attribute in attribute_values:
        attr_name = attribute.mName
        attr_val = attribute.mVal

        # If conditions are met, area feature is assigned to feature class AFC_WOODS
        if attr_name == "COMMUNITY" and attr_val == "GREAT BASIN CONIFER WOODLAND":
            err_SetFeatureClass = gm.SetFeatureClass(
                layer_arr[0], gm.GM_FeatureClass_Area, i, gm.AFC_WOODS
            )
            modified_ct += 1

        # If conditions are met, area feature is given a label
        if attr_name == "COMMUNITY" and attr_val == "GREAT BASIN CONIFER WOODLAND":
            err_SetFeatureLabel = gm.SetFeatureLabel(
                layer_arr[0],
                gm.GM_FeatureClass_Area,
                i,
                "Great Basin Conifer Woodland",
            )

    # Frees deep copies from memory
    # Without this step the script may freeze
    gm.FreeAreaFeature(feature)
    i += 1
    feature = gm.GetAreaFeature(layer_arr[0], i)

# Displays the number of area features that were set to the Feature Class "Area Woods".
print("Set", modified_ct, "area features to the type WOODS")

### END SAMPLE

assert (
    gm.GM_Error_None == err_LoadLayerList
), f"LoadLayerList: Expected {gm.GM_Error_None}, got {err_LoadLayerList}"
assert (
    gm.GM_Error_None == err_SetFeatureClass
), f"SetFeatureClass: Expected {gm.GM_Error_None}, got {err_SetFeatureClass}"
assert (
    gm.GM_Error_None == err_SetFeatureLabel
), f"SetFeatureLabel: Expected {gm.GM_Error_None}, got {err_SetFeatureLabel}"

# The original amount of features modified when this sample was created
if modified_ct != 71:
    assert (
        False
    ), f"Expected 71 area features to be modified, but {modified_ct} were instead"

assert os.path.exists(layer_data), f"Sample data doesn't exist: {layer_data}"

import os
import globalmapper as gm

layer_data = "Z:/GM-20069/vector_layer_editing_sample_data/gfveg.shp"

# A sample script that deletes all holes from a shapefile by iterating through each area feature. For 
# each feature that hasn't already been deleted, it deletes that hole. Afterward, the total number of 
# holes deleted and the layer from which they were removed are output.

# Closes all layers that may be open.
layers_are_open = gm.GetLoadedLayerList()
if layers_are_open:
    arr_ptr, arr_size = layers_are_open
    initial_layers = gm.GM_LayerHandle_array_frompointer(arr_ptr)
    for i in range(arr_size):
        gm.CloseLayer(initial_layers[i])

err_LoadLayerList, arr_ptr, arr_size = gm.LoadLayerList(
    layer_data, gm.GM_LoadFlags_UseDefaultLoadOpts
)

arr_ptr, arr_size = gm.GetLoadedLayerList()
layers_array = gm.GM_LayerHandle_array_frompointer(arr_ptr)
layers_list = gm.carray_to_list(layers_array, arr_size)

### START SAMPLE
# Iterates over area features
i = 0
holes_deleted = 0
feature = gm.GetAreaFeature(layers_list[0], i)
# Goes through all available area features in the provided layer
while feature:
    # For each hole in that feature, if it's not deleted, delete it
    for j in range(feature.mNumHoles):
        err_IsHoleInAreaDeleted, val = gm.IsHoleInAreaDeleted(layers_list[0], i, j)
        if not val:
            err_DeleteHoleInArea = gm.DeleteHoleInArea(layers_list[0], i, j, True)
            holes_deleted += 1

    # Frees deep copies from memory
    # Without this step the script causes the program to freeze
    gm.FreeAreaFeature(feature)
    i += 1
    feature = gm.GetAreaFeature(layers_list[0], i)

layer_name = gm.GetLayerInfo(layers_list[0]).mDescription
print("\n" + str(holes_deleted), "holes deleted from", layer_name, "\n")

### END SAMPLE

assert (
    gm.GM_Error_None == err_LoadLayerList
), f"LoadLayerList: Expected {gm.GM_Error_None}, got {err_LoadLayerList}"
assert (
    gm.GM_Error_None == err_DeleteHoleInArea
), f"DeleteHoleInArea: Expected {gm.GM_Error_None}, got {err_DeleteHoleInArea}"

assert os.path.exists(layer_data), f"Sample data doesn't exist: {layer_data}"

assert holes_deleted == 191, "Incorrect number of holes deleted"

import os
import globalmapper as gm

import_dir = "Z:/GM-20069/vector_layer_editing_sample_data/"

# A list of attributes is created for a given feature and then assigned to it.

# Close all layers that may be open
layers_are_open = gm.GetLoadedLayerList()
if layers_are_open:
    arr_ptr, arr_size = layers_are_open
    initial_layers = gm.GM_LayerHandle_array_frompointer(arr_ptr)
    for i in range(arr_size):
        gm.CloseLayer(initial_layers[i])

err_LoadLayerList, array_pointer, array_size_param = gm.LoadLayerList(
    os.path.join(import_dir, "area.shp"),
    gm.GM_LoadFlags_HideAllPrompts,
)

arr_ptr, arr_size = gm.GetLoadedLayerList()
layers_array = gm.GM_LayerHandle_array_frompointer(arr_ptr)
layers_list = gm.carray_to_list(layers_array, arr_size)

### START SAMPLE
# Adds two attributes to an area feature
attr_1 = gm.GM_AttrValue_t("GM_TYPE", "Building")
attr_2 = gm.GM_AttrValue_t("AREA_TYPE", "Residential")
attr_list = [attr_1, attr_2]
size = len(attr_list)
attr_array = gm.GM_AttrValue_array(size)
for i in range(size):
    attr_array[i] = attr_list[i]
err_SetFeatureAttrList = gm.SetFeatureAttrList(
    layers_list[0], gm.GM_FeatureClass_Area, 0, attr_array, size
)

### END SAMPLE

assert (
    gm.GM_Error_None == err_LoadLayerList
), f"LoadLayerList: Expected {gm.GM_Error_None}, got {err_LoadLayerList}"
assert (
    gm.GM_Error_None == err_SetFeatureAttrList
), f"SetFeatureAttrList: Expected {gm.GM_Error_None}, got {err_SetFeatureAttrList}"

assert os.path.isdir(import_dir), f"Import directory doesn't exist: {import_dir}"

assert os.path.exists(
    os.path.join(import_dir, "area.shp")
), "Sample data doesn't exist"