Error when Projecting Rasters using ArcPy?

Question

Using ArcGIS 10.3 I'm receiving the dreaded 'ERROR 999999' when attempting to project a geodatabase full of rasters from WGS84_UTM_zone_10N to Nad83_BC_Environment_Albers

import arcpy

InRasterGDB = r"C:\Users\PoirierM\Documents\GISData\LandsatSSTImagery.gdb"
OutRasterGDB = r"C:\Users\PoirierM\Documents\GISData\LandsatSSTImageryAlbers.gdb"
OutCoordSystem = "3005"
ResampleType = "BILINEAR"
CellSize = "30 30"
GeoTransform = "WGS_1984_(ITRF00)_To_NAD_1983"
InCoordSystem = "32610"

arcpy.env.workspace = InRasterGDB

for fname in arcpy.ListRasters():
    print fname
    print OutRasterGDB + "\\"  + fname
    arcpy.ProjectRaster_management(InRasterGDB + "\\" + fname, OutRasterGDB + "\\" + fname, OutCoordSystem, ResampleType, CellSize, GeoTransform, "", InCoordSystem)

This should work, but it returns the following error when run in the Python Window in ArcCatalog:

Runtime error
  Traceback (most recent call last):
   File "<string>", line 16, in <module>
   File "c:\program files (x86)\arcgis\desktop10.3\arcpy\arcpy\management.py", line 9009, in ProjectRaster
     raise e ExecuteError: ERROR 999999: Error executing function. Failed to execute (ProjectRaster). 

I can do a single raster manually through either the tool interface or within modelbuilder. Can anyone suggest why it won't work as a Python script?

---

I've got a working script now - see my (MPoirier) answer below.

Answer 1

I believe the problem might be in your coordinate system definitions -- ref. the Project Raster help page.

The coordinate system to which the input raster will be projected. The default value is set based on the Output Coordinate System environment setting.

Valid values for this parameter are

• A file with the ".prj" extension (the prj files which ship with ArcGIS can be found in "C:\Program Files\ArcGIS\Coordinate Systems").
• An existing feature class, feature dataset, raster catalog (basically anything with a coordinate system).
• The string representation of a coordinate system. These lengthy strings can be generated by adding a coordinate system variable to ModelBuilder, setting the variable's value as desired, then exporting the model to a Python script.

You're giving it WKIDs directly, which isn't listed as an option.

But, you can make a spatial reference via the WKID using the SpatialReference property.

OutCoordSystem = arcpy.SpatialReference(3005)
InCoordSystem = arcpy.SpatialReference(32610)

Answer 2

Here's a working script. I started from the modelbuilder code for reprojecting a single raster, and incrementally changed inputs to variables testing at each step. Seemingly the only difference between the working script below and my previous non-functioning script is that the variables are defined within the body of the for loop.

I've also added some comments to make things clear.

##Script for looping through a file geodatabase to find rasters and reproject them from WGS_1984_UTM_zone_10N to NAD_1983_BC_Environment_Albers
    ## Matthew Poirier, Last modified February 18th 2015

    # Import arcpy module
    import arcpy

    #Define workspace (geodatabase location)
    InRasterGDB = r"C:\Users\PoirierM\Documents\GISData\LandsatSSTImagery.gdb"
    arcpy.env.workspace = InRasterGDB

    #Loop through list of all rasters in the workspace
    for Ras in arcpy.ListRasters():

    # Local variables for input (optional) and output:
        InRas = "C:\\Users\\PoirierM\\Documents\\GISData\\LandsatSSTImagery.gdb\\" + Ras 
        OutRas = "C:\\Users\\PoirierM\\Documents\\GISData\\LandsatSSTImageryAlbers.gdb\\" + Ras

    # Define the projection strings - could also use string longform copied from modelbuilder script (tested) or an existing .prj file
        InProj = "PROJCS['WGS_1984_UTM_zone_10N',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['false_easting',500000.0],PARAMETER['false_northing',0.0],PARAMETER['central_meridian',-123.0],PARAMETER['scale_factor',0.9996],PARAMETER['latitude_of_origin',0.0],UNIT['Meter',1.0]]"
        OutProj = "PROJCS['NAD_1983_BC_Environment_Albers',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',1000000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',-126.0],PARAMETER['Standard_Parallel_1',50.0],PARAMETER['Standard_Parallel_2',58.5],PARAMETER['Latitude_Of_Origin',45.0],UNIT['Meter',1.0]]"

    # Specify the geographic datum transformation method to be used
        GeoTran = "WGS_1984_(ITRF00)_To_NAD_1983"

    # Specify the resampling type (bilinear or cubic for continuous data, defaults to nearest)
        ReSample = "BILINEAR"

    # Define the output cell size (defaults to that of the input raster
        CellSize = "30 30"

    # Process: Project Raster - note empty quotes are placeholder for optional 'registration_point' parameter
        arcpy.ProjectRaster_management(InRas, OutRas, OutProj, ReSample, CellSize, GeoTran, "", InProj)

Answer 3

Try passing cell size as one dimension.

like:

CellSize = "30"

Otherwise try testing the script without any of the optional parameters and see what happens.