def join_attributes(raster, raster_join_field, attribute_file,
attribute_join_field, attribute_metadata_file):
"""
Join attributes to a raster
Parameters:
-----------
raster : str
raster to join attributes to
raster_join_field : str
name of join field in raster
attribute_file: str
name of file with attributes to join to raster
attribute_join_field: str
name of join field in attribute file
attribute_metadata_file: str
name of file with attribute metadata to decide what variables to
drop from join file
Returns:
--------
None
"""
model_dir = get_path(raster)
gp = geoprocessor.Geoprocessor(model_dir)
# create list of attributes to drop from join file - we only want a
# subset of all variables joined to the NN grids (PROJECT_ATTR = 1), so
# we need to specify all variables that have PROJECT_ATTR = 0 in the
# metadata
mp = xsmp.XMLStandMetadataParser(attribute_metadata_file)
drop_fields = [x.field_name for x in mp.attributes if x.project_attr == 0]
try:
gp.join_attributes(raster,
raster_join_field,
attribute_file,
attribute_join_field,
drop_fields=drop_fields)
except:
print(sys.exc_info())
def build_vats(rasters):
"""
Build value attribute tables
Parameters:
-----------
rasters : list
list of rasters to build vats (or update if vat already built)
Returns:
--------
None
"""
for raster in rasters:
try:
model_dir = get_path(raster)
gp = geoprocessor.Geoprocessor(model_dir)
gp.build_vat(raster)
except:
print(sys.exc_info())
continue
def copy_rasters(rasters):
"""
Copy rasters without the attribute table. The rasters are overwritten
to the same input names
Parameters:
-----------
rasters : list
list of rasters to copy
Returns:
--------
None
"""
for raster in rasters:
try:
model_dir = get_path(raster)
gp = geoprocessor.Geoprocessor(model_dir)
gp.overwrite(gp.copy_raster_no_attributes, raster)
except:
print(sys.exc_info())
continue
def integerize_rasters(rasters):
"""
Convert floating point rasters to integer rasters. The rasters are
overwritten to the same input names
Parameters:
-----------
rasters : list
list of rasters to convert to integers
Returns:
--------
None
"""
for raster in rasters:
try:
model_dir = get_path(raster)
gp = geoprocessor.Geoprocessor(model_dir)
gp.overwrite(gp.convert_to_integer, raster)
except:
print(sys.exc_info())
continue
def define_projections(rasters, projection_file):
"""
Define projections on rasters
Parameters:
-----------
rasters : list
list of rasters to define projections on
projection_file : str
name and full path of file containing projection parameters
Returns:
--------
None
"""
for raster in rasters:
try:
model_dir = get_path(raster)
gp = geoprocessor.Geoprocessor(model_dir)
gp.define_projection(raster, projection_file)
except:
print(sys.exc_info())
continue
def create_masked_raster(in_raster, boundary_raster, mask_raster, out_raster):
"""
Create a masked raster
Parameters:
-----------
in_raster : str
input raster to apply mask
boundary_raster : str
boundary raster used to clip masked output
mask_raster: str
raster to use as a mask
out_raster: str
name of output masked raster
Returns:
--------
None
"""
model_dir = get_path(in_raster)
gp = geoprocessor.Geoprocessor(model_dir)
gp.create_clipped_masked_raster(in_raster, boundary_raster, mask_raster,
out_raster)