I have quite a few hundred rasters that I am looking to generate contour lines off. I have successfully run the Contour tool on individual rasters located in the 3D Analyst Tools -> Raster Surface -> toolkit but can anyone advise if there is a script or other tool that I can point to a series of directories to process quite a few hundred sets of contour lines off? I've seen the script in the help file which is below but is there a fileobject algorithm that can be run over a directory with to extract filenames with many files in it? Contour (3D Analyst) # Name: Contour_3d_Ex_02.py # Description: Creates contours or isolines from a raster surface. # Requirements: 3D Analyst Extension # Import system modules import arcpy from arcpy import env # Set environment settings env . workspace = "C:/data" # Set local variables inRaster = "elevation" contourInterval = 200 baseContour = 0 outContours = "C:/sapyexamples/output/outcontours02.shp" # Check out the ArcGIS 3D Analyst extension license arcpy . CheckOutExtension ( "3D" ) # Execute Contour arcpy . Contour_3d ( inRaster , outContours , contourInterval , baseContour )

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

Running Contour Tool on Multiple Rasters

PhillipCouch

I have quite a few hundred rasters that I am looking to generate contour lines off. I have successfully run the Contour tool on individual rasters located in the 3D Analyst Tools -> Raster Surface -> toolkit but can anyone advise if there is a script or other tool that I can point to a series of directories to process quite a few hundred sets of contour lines off?

I've seen the script in the help file which is below but is there a fileobject algorithm that can be run over a directory with to extract filenames with many files in it?

Contour (3D Analyst)

<SPAN class=""># Name: Contour_3d_Ex_02.py</SPAN><SPAN class=""># Description: Creates contours or isolines from a raster surface.</SPAN><SPAN class=""># Requirements: 3D Analyst Extension</SPAN><SPAN class=""># Import system modules</SPAN><SPAN class="">import</SPAN> <SPAN class="">arcpy</SPAN><SPAN class="">from</SPAN> <SPAN class="">arcpy</SPAN> <SPAN class="">import</SPAN> <SPAN class="">env</SPAN><SPAN class=""># Set environment settings</SPAN><SPAN class="">env</SPAN><SPAN class="">.</SPAN><SPAN class="">workspace</SPAN> <SPAN class="">=</SPAN> <SPAN class="">"C:/data"</SPAN><SPAN class=""># Set local variables</SPAN><SPAN class="">inRaster</SPAN> <SPAN class="">=</SPAN> <SPAN class="">"elevation"</SPAN><SPAN class="">contourInterval</SPAN> <SPAN class="">=</SPAN> <SPAN class="">200</SPAN><SPAN class="">baseContour</SPAN> <SPAN class="">=</SPAN> <SPAN class="">0</SPAN><SPAN class="">outContours</SPAN> <SPAN class="">=</SPAN> <SPAN class="">"C:/sapyexamples/output/outcontours02.shp"</SPAN><SPAN class=""># Check out the ArcGIS 3D Analyst extension license</SPAN><SPAN class="">arcpy</SPAN><SPAN class="">.</SPAN><SPAN class="">CheckOutExtension</SPAN><SPAN class="">(</SPAN><SPAN class="">"3D"</SPAN><SPAN class="">)</SPAN><SPAN class=""># Execute Contour</SPAN><SPAN class="">arcpy</SPAN><SPAN class="">.</SPAN><SPAN class="">Contour_3d</SPAN><SPAN class="">(</SPAN><SPAN class="">inRaster</SPAN><SPAN class="">,</SPAN> <SPAN class="">outContours</SPAN><SPAN class="">,</SPAN> <SPAN class="">contourInterval</SPAN><SPAN class="">,</SPAN> <SPAN class="">baseContour</SPAN><SPAN class="">)</SPAN>

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

Accepted answers

All comments

DanPatterson_Retired

And this missive on paths.... /blogs/dan_patterson/2016/08/14/filenames-and-file-paths-in-python

spaces are the killer when working with the spatial analyst... simply moving folder contents goes a long way.

ChrisDonohue__GISP

Not to get too far off subject, but Python will handle forward and back slashes. But forward slashes are easier to implement.

Also, from an ESRI blog:

Notice that in calling the clip tool, we used forward slashes (/). Back slashes () and forward slashes work equally well in ArcGIS on ALL platforms, so you don’t need to worry about switching back and forth when you change to Windows. You can even use back slashes on Linux and Solaris. The problem is that in Python, back slashes are a special character, so to use them, you have to treat them differently than most other characters.

- You can escape them: “\some\path”

- You can use raw strings: r”somepath”

- You can use unicode strings: u”somepath”

If you stick to forward slashes on all platforms, you improve the look of your code
and reduce the effort it takes to type it.

Source: Tips and tricks – Unix, Python, and Geoprocessing | ArcGIS Blog

Chris Donohue, GISP

PhillipCouch

Thanks Dan, the file names and folder names are huge and renaming 17000 rasters would be a nightmare. Upon a second look I noted the python coding uses forward slashes rather than back slashes for directory paths. I also thought that the "elevation" was a variable call to the contour function but it is actually the name of the raster file you are trying to process. Below is the working code with comments for anyone it can help.

This is basically from the help guide with brief commentary that tripped me up.

<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">from</SPAN> arcpy <SPAN class="keyword token">import</SPAN> env
<SPAN class="keyword token">from</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa <SPAN class="keyword token">import</SPAN> <SPAN class="operator token">*</SPAN>
env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> <SPAN class="string token">"enter the directory path where your files are"</SPAN>
inRaster <SPAN class="operator token">=</SPAN> <SPAN class="string token">"enter your file name with .asc at the end"</SPAN>
contourInterval <SPAN class="operator token">=</SPAN> <SPAN class="number token">5</SPAN>
baseContour <SPAN class="operator token">=</SPAN> <SPAN class="number token">0</SPAN>
outContours <SPAN class="operator token">=</SPAN> <SPAN class="string token">"D:/map/dembula/really long nested folder structure with spaces in the folder name/testcountours.shp"</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Spatial"</SPAN><SPAN class="punctuation token">)</SPAN>
Contour<SPAN class="punctuation token">(</SPAN>inRaster<SPAN class="punctuation token">,</SPAN> outContours<SPAN class="punctuation token">,</SPAN> contourInterval<SPAN class="punctuation token">,</SPAN> baseContour<SPAN class="punctuation token">)</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

DanPatterson_Retired

Yes... there is a limit on path lengths... and if there are any spaces in them, many Spatial Analyst tools will fail.

Move your stuff to a simple folder like...

c:\test\Data

ChrisDonohue__GISP

To expand on Dan's suggestion, here's how to format the code in GeoNet (a process which is not obvious). Since spacing is important for Python, getting the code formatted can aid in troubleshooting by GeoNet responders.

I know you stated you tried the syntax highlighter, but just wanted to be sure the complete instructions were available.

Posting code with Syntax Highlighting on GeoNet

Chris Donohue, GISP

PhillipCouch

Thanks I've used the syntax highlighter for python and keep getting syntax errors on line 4. It's a pretty long directory name and file name but I've used the syntax call from the help guide which is repeated below.

import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "input file directory and filename" Contour("elevation", "D:\test\contour.shp", 5, 0)

I've also tried copying the code from the background process syntax and getting the same error.

import arcpy
from arcpy import env
from arcpy.sa import *
Contour "input file directory and filename" "\test\contour.shp" 5 0 1

Can anyone please tell me if there is a glaring issue in the code?

DanPatterson_Retired

Phillip, it would help if you formatted your code using the syntax highlighter for python ( follow the dots ... in thread tools)

PhillipCouch

Thankyou Abdullah, I'm having trouble creating a raster dataset and using model builder so I'll stick with the coding.

AbdullahAnter

Put all rasters in one (workspace/Raster dataset). Then use ModelBuilder and put Iterate Rasters.

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.