Converted polygon shapefiles to a 30-meter raster grid.
An ArcGIS 10.0 model builder script was run to convert each set
of block group boundaries (by state) to polygons. An "exported"
python version of is attached to document the processing steps.
# ---------------------------------------------------------------------------
# script.py
# Created on: 2012-02-23 16:09:53.00000
# (generated by ArcGIS/ModelBuilder)
# Usage: script <Block_group_shapefile> <Block_group_raster>
# Description:
# ---------------------------------------------------------------------------
# Import arcpy module
import arcpy
# Check out any necessary licenses
arcpy.CheckOutExtension("spatial")
# Load required toolboxes
arcpy.ImportToolbox("E:/work/bg2010/Census To Raster Tools.tbx")
# Script arguments
Block_group_shapefile = arcpy.GetParameterAsText(0)
if Block_group_shapefile == '#' or not Block_group_shapefile:
Block_group_shapefile = "tl_2010_11_bg10" # provide a default value if unspecified
Block_group_raster = arcpy.GetParameterAsText(1)
if Block_group_raster == '#' or not Block_group_raster:
Block_group_raster = "E:\\work\\bg10g" # provide a default value if unspecified
# Local variables:
tabblock2010_11_pophu_Projec1_shp = Block_group_shapefile
Raster_table__2_ = Block_group_shapefile
Raster_table = Block_group_raster
Snap_raster = "%scratchworkspace%\\xxsnap"
# Process: Project
arcpy.Project_management(Block_group_shapefile, tabblock2010_11_pophu_Projec1_shp,
"PROJCS['USA_Contiguous_Albers_Equal_Area_Conic_USGS_version',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983'," \
"SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]," \
"PROJECTION['Albers'],PARAMETER['False_Easting',0.0],PARAMETER['False_Northing',0.0]," \
"PARAMETER['Central_Meridian',-96.0],PARAMETER['Standard_Parallel_1',29.5],PARAMETER['Standard_Parallel_2',45.5]," \
"PARAMETER['Latitude_Of_Origin',23.0],UNIT['Meter',1.0]]", "",
"GEOGCS['GCS_North_American_1983'" \
"DATUM['D_North_American_1983',SPHEROID['GRS_1980',6378137.0,298.257222101]]," \
"PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]")
# Process: Create Snap Raster
arcpy.CreateSnapRaster_c2raster("15", "15", "30", Snap_raster,
"PROJCS['USA_Contiguous_Albers_Equal_Area_Conic_USGS_version',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983'," \
"SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]," \
"PROJECTION['Albers'],PARAMETER['False_Easting',0.0],PARAMETER['False_Northing',0.0]," \
"PARAMETER['Central_Meridian',-96.0],PARAMETER['Standard_Parallel_1',29.5],PARAMETER['Standard_Parallel_2',45.5]," \
"PARAMETER['Latitude_Of_Origin',23.0],UNIT['Meter',1.0]]")
# Process: Polygon to Raster
tempEnvironment0 = arcpy.env.snapRaster
arcpy.env.snapRaster = "%scratchworkspace%\\xxsnap"
arcpy.PolygonToRaster_conversion(tabblock2010_11_pophu_Projec1_shp, "FID", Block_group_raster, "MAXIMUM_AREA", "NONE", "30")
arcpy.env.snapRaster = tempEnvironment0
# Process: Make Table View
arcpy.MakeTableView_management(Block_group_raster, Raster_table, "", "", "VALUE VALUE VISIBLE NONE;COUNT COUNT VISIBLE NONE")
# Process: Join Field
arcpy.JoinField_management(Raster_table, "VALUE", Block_group_shapefile, "FID", "GEOID10;ALAND10;AWATER10"