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; } } }

Find and fix images ith Spatial Reference Undefined.

PaulHuffman

I have a directory with 370 tifs. I just found out that a few, an unknown number, have the projection undefined, but they draw aligned with the know projection tifs, so no mystery there. I thought I could loop through a list of tifs in the folder and set the definition on the files missing, similar to https://community.esri.com/thread/13078 but wondered if I should add an "if" to check if the image is one of the ones with projection undefined. But what would be that test? Do you think the test would speed up execution or not?

Najděte více příspěvků označených jako

Accepted answers

PaulHuffman

Yes, the sample code was very helpful. Spent too much time looking through the Syntax and Properties section, when the sample had a clear example of the test. And yes lower case w on workspace was also needed. I couldn't figure out why my raster list seemed to be empty when it hit line 20.

This code worked on a sample data folder. Now I'll run it on the rest. I'm not sure if this testing if the coordinate system is undefined decreased run time or if I could even notice a difference with 370 files. Maybe 3000 files.

import sys, os, arcpy
from arcpy import env
# Load required toolboxes...  
#gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx")  
         
#workspace  
env.workspace = "G:\\usgs_quads\\test"  
#coordsys = "C:\\Program Files\\ArcGISCoordinate Systems\\Projected Coordinate Systems\\National Grids\\British National Grid.prj"  
      
#List all features in  gp.workspace  
rasters = arcpy.ListRasters("*", "TIF")  
#raster = rasters.Next()


# get the coordinate system by describing a feature class
dsc = arcpy.Describe("YESMOWIT CANYON.tif")
coord_sys = dsc.spatialReference

  
for raster in rasters:
    # Create the spatial reference object
    spatial_ref = arcpy.Describe(raster).spatialReference

    # If the spatial reference is unknown
    if spatial_ref.name == "Unknown":
     arcpy.DefineProjection_management(raster, coord_sys)  
     print raster + " is defined."

All comments

PaulHuffman

import sys, os, arcpy
from arcpy import env
# Load required toolboxes...  
#gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx")  
         
#workspace  
env.workspace = "G:\\usgs_quads\\test"  
#coordsys = "C:\\Program Files\\ArcGISCoordinate Systems\\Projected Coordinate Systems\\National Grids\\British National Grid.prj"  
      
#List all features in  gp.workspace  
rasters = arcpy.ListRasters("*", "TIF")  
#raster = rasters.Next()


# get the coordinate system by describing a feature class
dsc = arcpy.Describe("YESMOWIT CANYON.tif")
coord_sys = dsc.spatialReference

  
for raster in rasters:
    # Create the spatial reference object
    spatial_ref = arcpy.Describe(raster).spatialReference

    # If the spatial reference is unknown
    if spatial_ref.name == "Unknown":
     arcpy.DefineProjection_management(raster, coord_sys)  
     print raster + " is defined."

DanPatterson_Retired

you did check the code sample in my previous post from the help file?. You want to get a list of them first ... before ... you define them. Just to make sure we are on the same page

ps you don't need arcgisscripting and check to see whether workspace is capitalized in env.w...

PaulHuffman

I figure the code would be something like this, revised slightly from something I had laying around for fixing shapefiles, trying to remove the old style gp stuff. I'm still trying to figure out how I might test for projection undefined. but it might run just as fast without the test.

import arcgisscripting, sys, os, arcpy
from arcpy import env
# Load required toolboxes...  
#gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx")  
         
#workspace  
env.Workspace = "G:\\usgs_quads\\test"  
#coordsys = "C:\\Program Files\\ArcGISCoordinate Systems\\Projected Coordinate Systems\\National Grids\\British National Grid.prj"  
      
#List all features in  gp.workspace  
rasters = arcpy.ListRasters("*", "TIF")  
#raster = rasters.Next()


# get the coordinate system by describing a feature class
dsc = arcpy.Describe("YESMOWIT CANYON.tif")
coord_sys = dsc.spatialReference

  
for raster in rasters:  
     arcpy.defineprojection(raster, coord_sys)  
     print fc + " is defined."

DanPatterson_Retired

since you don't appear to be afraid of scripting, I will refer you to SpatialReference—Help | ArcGIS for Desktop

which contains all the appropriate information of coordinate systems and their associated properties. The code sample at the bottom will be of use to you

Připojte se k diskusi

Přihlášení členové mohou přispívat, sledovat aktualizace a další. Jste tu noví? Zaregistrujte si bezplatný účet.

Rychlé odkazy

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