gp = GeoProcess()
#########################################################
# Fix geometries
#########################################################
print('fixing geometries')
fix_geo = gp.fix_geometry(wldsin)
#########################################################
# Add autoincremental field
#########################################################
print('adding autoincremental id-field')
autoinc_id = gp.add_autoincremental_id(fix_geo, 'GID')
#########################################################
# Copy attribute
#########################################################
print('copying language name into a field with shorter attribute name')
language_name_copied = gp.copy_attribute(autoinc_id, 'NAME_PROP', 'lnm')
#########################################################
# Drop field(s)
#########################################################
print('dropping fields except GID, ID, lnm')
keepfields = ['GID', 'ID', 'lnm']
gp.drop_fields(language_name_copied,
keep_fields=keepfields,
output_object=wldsout)
print('DONE!')
##################################################################
print('extracting vertices (get endpoints of each line)')
extract_vertices = gp.extract_vertices(centroids_nearest_coast_distance_joined)
##################################################################
# Extract by attribute
##################################################################
print('keeping only vertices on coast')
extract_by_attribute = gp.extract_by_attribute(extract_vertices, 'distance',
'>', 0)
#########################################################
# Copy attribute
#########################################################
print('creating new field: centroid latitude (keep field names straight)')
added_field_cent_lat = gp.copy_attribute(extract_by_attribute, 'ycoord',
'cent_lat')
#########################################################
# Copy attribute
#########################################################
print('creating new field: centroid longitude (keep field names straight)')
added_field_cent_lon = gp.copy_attribute(added_field_cent_lat, 'xcoord',
'cent_lon')
#########################################################
# Drop fields
#########################################################
print('dropping unnecessary fields')
keepfields = ['ne_10m_adm', 'ADMIN', 'ISO_A3', 'cent_lat', 'cent_lon']
centroids_lat_lon_drop_fields = gp.drop_fields(added_field_cent_lon,
keep_fields=keepfields)