###########################################################################
######################################################################
# reproject rivers, rail, and counties to WGS84
######################################################################
print('reprojecting rivers, rail, and counties to WGS84')
rail70wgs84 = gp.reproject_layer(rail70, 'EPSG:4326')
rail90wgs84 = gp.reproject_layer(rail90, 'EPSG:4326')
riverswgs84 = gp.reproject_layer(rivers, 'EPSG:4326')
countieswgs84 = gp.reproject_layer(counties, 'EPSG:4326')
######################################################################
# county centroids
######################################################################
print('finding county centroids')
centroids = gp.centroids(countieswgs84)
######################################################################
# adding indicators and dropping unnecessary fields
######################################################################
print('adding indicators')
rail70wgs84 = gp.add_constant_field(rail70wgs84, 'rail', 1)
rail90wgs84 = gp.add_constant_field(rail90wgs84, 'rail', 1)
riverswgs84 = gp.add_constant_field(riverswgs84, 'river', 1)
centroids = gp.add_constant_field(centroids, 'centroids', 1)
print('dropping unnecessary fields')
rail70wgs84 = gp.drop_fields(rail70wgs84, keep_fields=['rail'])
rail90wgs84 = gp.drop_fields(rail90wgs84, keep_fields=['rail'])
riverswgs84 = gp.drop_fields(riverswgs84, keep_fields=['river'])
centroids = gp.drop_fields(centroids, keep_fields=['centroids'])
# Fix geometries
#########################################################
print('fixing geometries, coast')
fixgeo_coast = gp.fix_geometry(coastin)
#########################################################
# Fix geometries
#########################################################
print('fixing geometries, countries')
fixgeo_countries = gp.fix_geometry(adminin)
#########################################################
# Centroids
#########################################################
print('finding country centroids')
country_centroids = gp.centroids(fixgeo_countries)
#########################################################
# Add geometry attributes
#########################################################
print('adding co-ordinates to centroids')
centroids_with_coordinates = gp.add_xycoordinates(country_centroids)
#########################################################
# Drop fields
#########################################################
print('dropping unnecessary fields, coast')
keepfields = ['featurecla']
gp.drop_fields(fixgeo_coast, keep_fields=keepfields, output_object=coastout)
##################################################################