if not os.path.exists(junkpath):
os.mkdir(junkpath)
gp = GeoProcess()
###########################################################################
###########################################################################
###########################################################################
######################################################################
# 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)
crs_wcea = QgsCoordinateReferenceSystem('ESRI:54034')
gp = GeoProcess()
##################################################################
# Drop field(s)
##################################################################
print('dropping unnecessary fields')
keepfields = ['ne_10m_adm', 'ADMIN', 'ISO_A3']
countries_drop_fields = gp.drop_fields(admin_in, keep_fields=keepfields)
##################################################################
# Reproject layer
##################################################################
print('projecting to world cylindical equal area')
countries_reprojected = gp.reproject_layer(countries_drop_fields, 'ESRI:54034')
##################################################################
# Fix geometries
##################################################################
print('fixing geometries')
countries_fix_geo = gp.fix_geometry(countries_reprojected)
##################################################################
# Field calculator, output to csv
##################################################################
print('calculating areas, outputting to csv')
gp.add_area_attribute(countries_fix_geo, 'area', output_object=areas_out)
print('DONE!')