nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Reprojecting Rasters Issue

EricaNewton

Hi there,

I am new to python and attempting to set up a batch process to reproject rasters in a folder.

This is the script and error message I am getting. I defined the projection of the rasters in the folder to WGS 1984 UTM Zone 16N before starting.

>>> import arcpy
... arcpy.env.workspace = r"C:/GIS/NDVI/Georef_27_21_TIFFs/test"
... rasterList = arcpy.ListRasters("*", "All")
... for raster in rasterList:
... print raster
... outraster = raster+"_ontlambproj"
... arcpy.ProjectRaster_management("raster", "outraster", "MNR_LAMBERT_PROJECTION_3161.prj", "NEAREST", "30", "WGS_1984_(ITRF00)_To_NAD_1983", "#", r"C:/Program Files/ArcGIS/Desktop10.0/Coordinate Systems/Projected Coordinate Systems/UTM/WGS 1984/Northern Hemisphere/WGS 1984 UTM Zone 16N.prj")
...
L5027021_02119850802_B30.TIF

LT5027021008516650_B3.tif

Runtime error <class 'arcgisscripting.ExecuteError'>: Undefined coordinate system for input dataset.

Can anyone help me determine why I am getting this error?

Thank you so much!

Erica

Finde mehr Beiträge mit dem Tag

Accepted answers

All comments

EricaNewton


[HTML]

Apparently I just wasn't waiting long enough.  The second code (below) works for both. 

import arcpy
arcpy.env.workspace = r"C:/GIS/NDVI/Georef_27_21_TIFFs/test"
rasters = arcpy.ListRasters("*", "TIF")
for raster in rasters:
    print raster
    outraster = raster.rpartition(".")[0]+"_olp.tif"
    arcpy.ProjectRaster_management(raster, outraster, "MNR_LAMBERT_PROJECTION_3161.prj", "NEAREST", "30", "WGS_1984_(ITRF00)_To_NAD_1983", "#", r"C:/Program Files/ArcGIS/Desktop10.0/Coordinate Systems/Projected Coordinate Systems/UTM/WGS 1984/Northern Hemisphere/WGS 1984 UTM Zone 16N.prj")

[/HTML]

MathewCoyle

Wrapping your code in [ CODE][/CODE ] tags lets you preserve your formatting.

As to the question, do both Rasters have unique names? Same coordinate systems? When you list the rasters do both come up? The syntax for limiting the list by .tif images is "TIFF" I believe.

EricaNewton

Sorry, the first script only has one indent for "print raster"

The second script is indented at "print raster" and all following lines.

EricaNewton

There is one final tweak that I can't seem to figure out with this script. Both of the following only return one reprojected file, not both that are in the "test" folder. Any ideas on how to get it to process both .tif files?

import arcpy
arcpy.env.workspace = r"C:/GIS/NDVI/Georef_27_21_TIFFs/test"
rasters = arcpy.ListRasters("*", "All")
for raster in rasters:
    print raster
outraster = raster.rpartition(".")[0]+"_olp.tif"
arcpy.ProjectRaster_management(raster, outraster, "MNR_LAMBERT_PROJECTION_3161.prj", "NEAREST", "30", "WGS_1984_(ITRF00)_To_NAD_1983", "#", r"C:/Program Files/ArcGIS/Desktop10.0/Coordinate Systems/Projected Coordinate Systems/UTM/WGS 1984/Northern Hemisphere/WGS 1984 UTM Zone 16N.prj")

import arcpy
arcpy.env.workspace = r"C:/GIS/NDVI/Georef_27_21_TIFFs/test"
rasters = arcpy.ListRasters("*", "TIF")
for raster in rasters:
    print raster
    outraster = raster.rpartition(".")[0]+"_olp.tif"
    arcpy.ProjectRaster_management(raster, outraster, "MNR_LAMBERT_PROJECTION_3161.prj", "NEAREST", "30", "WGS_1984_(ITRF00)_To_NAD_1983", "#", r"C:/Program Files/ArcGIS/Desktop10.0/Coordinate Systems/Projected Coordinate Systems/UTM/WGS 1984/Northern Hemisphere/WGS 1984 UTM Zone 16N.prj")

EricaNewton

Perfect. Thanks to both of you very much. Now I can keep going !

Erica

JamieDrisdelle

You can also use this to get the base name of a raster without its extension like this. Both options will work.

desc = arcpy.Describe(raster)
name = desc.baseName

Jamie

MathewCoyle

Yes that is just simple string manipulation, since raster = "originalfilename.tif" when you do outraster = raster+"_olp.tif" you get the logical originalfilename.tif_olp.tif

There are a couple ways to solve that issue.

raster = raster.rpartition(".")[0]

EricaNewton

THANK YOU!!! I thought it must have been a simple solution, but being very new to this, I have been working on it for days. I'm still having a problem with the file name (it renames the files as originalfilename.tif_olp.tif - need to fix that to originalfilename_olp.tif but I think I can figure that out - here is the final script that properly reprojects the tif.

... import arcpy
... arcpy.env.workspace = r"C:/GIS/NDVI/Georef_27_21_TIFFs/test"
... rasterList = arcpy.ListRasters("*", "All")
... for raster in rasterList:
... print raster
... outraster = raster+"_olp.tif"
... arcpy.ProjectRaster_management(raster, outraster, "MNR_LAMBERT_PROJECTION_3161.prj", "NEAREST", "30", "WGS_1984_(ITRF00)_To_NAD_1983", "#", r"C:/Program Files/ArcGIS/Desktop10.0/Coordinate Systems/Projected Coordinate Systems/UTM/WGS 1984/Northern Hemisphere/WGS 1984 UTM Zone 16N.prj")
...
L5027021_02119850802_B30.TIF

LT5027021008516650_B3.tif

MathewCoyle

"raster" and "outraster" are variables not string names, so remove the quotes and try that.

Nimm an der Unterhaltung teil

Angemeldete Mitglieder können Beiträge verfassen, Updates folgen und mehr. Neu hier? Registriere ein kostenloses Konto.

Schnellzugriffe

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.