Hfrhyu

Is this a crack on the carbon frame?

Is it tax fraud for an individual to declare non-taxable revenue as taxable income? (US tax laws)

Risk of getting Chronic Wasting Disease (CWD) in the United States?

What defenses are there against being summoned by the Gate spell?

What typically incentivizes a professor to change jobs to a lower ranking university?

Test whether all array elements are factors of a number

Prove that NP is closed under karp reduction?

How to format long polynomial?

Do VLANs within a subnet need to have their own subnet for router on a stick?

How does strength of boric acid solution increase in presence of salicylic acid?

How is the claim "I am in New York only if I am in America" the same as "If I am in New York, then I am in America?

Email Account under attack (really) - anything I can do?

Writing rule stating superpower from different root cause is bad writing

Why was the small council so happy for Tyrion to become the Master of Coin?

Have astronauts in space suits ever taken selfies? If so, how?

How does one intimidate enemies without having the capacity for violence?

Why does Kotter return in Welcome Back Kotter?

What do the dots in this tr command do: tr .............A-Z A-ZA-Z <<< "JVPQBOV" (with 13 dots)

If I cast Expeditious Retreat, can I Dash as a bonus action on the same turn?

Is it important to consider tone, melody, and musical form while writing a song?

How do we improve the relationship with a client software team that performs poorly and is becoming less collaborative?

Why not use SQL instead of GraphQL?

What are the differences between the usage of 'it' and 'they'?

Accidentally leaked the solution to an assignment, what to do now? (I'm the prof)

Projecting many raster files to coordinate system using ArcPy gives ERROR 999999?

Error when Projecting Rasters using ArcPy?Project a set of coordinates - points from Excel to a shapefile with unknown coordinate system (map-match)Projected raster data has inconsistent extent but not shapefileDo I need to reproject a NAD83 raster (with no projection) before using it with a shapefile that has a specific NAD83 projection?Why does the difference between two overlapping rasters exist?Error 999999 with Project RasterRe-projecting Multiple Shapefiles using ArcPy?ArcPy: 999999 ERROR when using AddJoin_managementArcGIS Desktop enforcing Vertical Coordinate System properties for rasters?Executing CostDistance Function in ArcPy gives Error 999999?

.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;

-1

I am very new to Python and trying to project many raster files using arcpy. Here is the screenshot of my raster files.

These raster files have an unknown coordinate system but I can project them with "ProjectRaster" in ArcGIS pro using 'GCS_WGS_1984'. I looked several online examples and revise one of them but could not work. Here is my scripts.

import acrpy
import os 
arcpy.env.workspace = "D:/STAR/assignment/assignment7/info_data_test" #The 
file where your files are
arcpy.env.overwriteOutput = True
outWorkspace = "D:/STAR/assignment/assignment7/info_data_test_results" #where 
you want the new rasters to go
arcpy.env.parallelProcessingFactor = "50"
rasters = arcpy.ListRasters()

#Begin loop

for inrl in rasters:

# Determine if the input has a defined coordinate system, can't project it if 
it does not
dsc = arcpy.Describe(inrl)

if dsc.spatialReference.Name == "Unknown":
 # Determine the new output feature class path and name
 outrl = os.path.join(outWorkspace, inrl)

 # Set output coordinate system
 outCS = arcpy.SpatialReference('GCS_WGS_1984')

 # run project tool
 arcpy.AddMessage("Reprojecting" + inrl)
 arcpy.Project_management(inrl, outrl, outCS)

 # check messages
 print(arcpy.GetMessages())

else:

 print ('skipped this raster file: ' + inrl)

 print(arcpy.AddMessage("Reprojection complete"))

 print "Script run ok"

When I ran this script, there appear mistakes. Here is the screenshot of error information.

arcpy coordinate-system error-999999

share|improve this question

edited Apr 3 at 2:34

PolyGeo♦

53.9k1781245

asked Apr 3 at 1:08

Zhenyu YaoZhenyu Yao

111

New contributor

Zhenyu Yao is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  
  Please remember to post error messages as ASCII text in the body of the question. Images are not searchable or legible.
  
  – Vince
  Apr 3 at 1:19
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  If you want to define a projection to the raster data, you need to use arcpy.DefineProjection_management(). Look at this help: pro.arcgis.com/en/pro-app/tool-reference/data-management/…. If you want to change the projection of raster data from one projection to another, you need to use arcpy.ProjectRaster_management. pro.arcgis.com/en/pro-app/tool-reference/data-management/…
  
  – ahmadhanb
  Apr 3 at 1:20
  
  

  
  
  
  

add a comment | 

-1

I am very new to Python and trying to project many raster files using arcpy. Here is the screenshot of my raster files.

These raster files have an unknown coordinate system but I can project them with "ProjectRaster" in ArcGIS pro using 'GCS_WGS_1984'. I looked several online examples and revise one of them but could not work. Here is my scripts.

import acrpy
import os 
arcpy.env.workspace = "D:/STAR/assignment/assignment7/info_data_test" #The 
file where your files are
arcpy.env.overwriteOutput = True
outWorkspace = "D:/STAR/assignment/assignment7/info_data_test_results" #where 
you want the new rasters to go
arcpy.env.parallelProcessingFactor = "50"
rasters = arcpy.ListRasters()

#Begin loop

for inrl in rasters:

# Determine if the input has a defined coordinate system, can't project it if 
it does not
dsc = arcpy.Describe(inrl)

if dsc.spatialReference.Name == "Unknown":
 # Determine the new output feature class path and name
 outrl = os.path.join(outWorkspace, inrl)

 # Set output coordinate system
 outCS = arcpy.SpatialReference('GCS_WGS_1984')

 # run project tool
 arcpy.AddMessage("Reprojecting" + inrl)
 arcpy.Project_management(inrl, outrl, outCS)

 # check messages
 print(arcpy.GetMessages())

else:

 print ('skipped this raster file: ' + inrl)

 print(arcpy.AddMessage("Reprojection complete"))

 print "Script run ok"

When I ran this script, there appear mistakes. Here is the screenshot of error information.

arcpy coordinate-system error-999999

share|improve this question

edited Apr 3 at 2:34

PolyGeo♦

53.9k1781245

asked Apr 3 at 1:08

Zhenyu YaoZhenyu Yao

111

New contributor

Zhenyu Yao is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  
  Please remember to post error messages as ASCII text in the body of the question. Images are not searchable or legible.
  
  – Vince
  Apr 3 at 1:19
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  If you want to define a projection to the raster data, you need to use arcpy.DefineProjection_management(). Look at this help: pro.arcgis.com/en/pro-app/tool-reference/data-management/…. If you want to change the projection of raster data from one projection to another, you need to use arcpy.ProjectRaster_management. pro.arcgis.com/en/pro-app/tool-reference/data-management/…
  
  – ahmadhanb
  Apr 3 at 1:20
  
  

  
  
  
  

add a comment | 

I am very new to Python and trying to project many raster files using arcpy. Here is the screenshot of my raster files.

These raster files have an unknown coordinate system but I can project them with "ProjectRaster" in ArcGIS pro using 'GCS_WGS_1984'. I looked several online examples and revise one of them but could not work. Here is my scripts.

import acrpy
import os 
arcpy.env.workspace = "D:/STAR/assignment/assignment7/info_data_test" #The 
file where your files are
arcpy.env.overwriteOutput = True
outWorkspace = "D:/STAR/assignment/assignment7/info_data_test_results" #where 
you want the new rasters to go
arcpy.env.parallelProcessingFactor = "50"
rasters = arcpy.ListRasters()

#Begin loop

for inrl in rasters:

# Determine if the input has a defined coordinate system, can't project it if 
it does not
dsc = arcpy.Describe(inrl)

if dsc.spatialReference.Name == "Unknown":
 # Determine the new output feature class path and name
 outrl = os.path.join(outWorkspace, inrl)

 # Set output coordinate system
 outCS = arcpy.SpatialReference('GCS_WGS_1984')

 # run project tool
 arcpy.AddMessage("Reprojecting" + inrl)
 arcpy.Project_management(inrl, outrl, outCS)

 # check messages
 print(arcpy.GetMessages())

else:

 print ('skipped this raster file: ' + inrl)

 print(arcpy.AddMessage("Reprojection complete"))

 print "Script run ok"

When I ran this script, there appear mistakes. Here is the screenshot of error information.

arcpy coordinate-system error-999999

share|improve this question

edited Apr 3 at 2:34

PolyGeo♦

53.9k1781245

asked Apr 3 at 1:08

Zhenyu YaoZhenyu Yao

111

New contributor

Zhenyu Yao is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

import acrpy
import os 
arcpy.env.workspace = "D:/STAR/assignment/assignment7/info_data_test" #The 
file where your files are
arcpy.env.overwriteOutput = True
outWorkspace = "D:/STAR/assignment/assignment7/info_data_test_results" #where 
you want the new rasters to go
arcpy.env.parallelProcessingFactor = "50"
rasters = arcpy.ListRasters()

#Begin loop

for inrl in rasters:

# Determine if the input has a defined coordinate system, can't project it if 
it does not
dsc = arcpy.Describe(inrl)

if dsc.spatialReference.Name == "Unknown":
 # Determine the new output feature class path and name
 outrl = os.path.join(outWorkspace, inrl)

 # Set output coordinate system
 outCS = arcpy.SpatialReference('GCS_WGS_1984')

 # run project tool
 arcpy.AddMessage("Reprojecting" + inrl)
 arcpy.Project_management(inrl, outrl, outCS)

 # check messages
 print(arcpy.GetMessages())

else:

 print ('skipped this raster file: ' + inrl)

 print(arcpy.AddMessage("Reprojection complete"))

 print "Script run ok"

share|improve this question

edited Apr 3 at 2:34

PolyGeo♦

53.9k1781245

asked Apr 3 at 1:08

Zhenyu YaoZhenyu Yao

111

New contributor

Zhenyu Yao is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

share|improve this question

edited Apr 3 at 2:34

PolyGeo♦

53.9k1781245

asked Apr 3 at 1:08

Zhenyu YaoZhenyu Yao

111

New contributor

Zhenyu Yao is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

edited Apr 3 at 2:34

PolyGeo♦

53.9k1781245

edited Apr 3 at 2:34

PolyGeo♦

53.9k1781245

asked Apr 3 at 1:08

Zhenyu YaoZhenyu Yao

111

New contributor

Zhenyu Yao is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

asked Apr 3 at 1:08

Zhenyu YaoZhenyu Yao

111