I am working on a pilot for a system that will eventually loop over a number of raster files to produce a single output one. If I have two input rasters, the routine works reliably. If I have more than three, it always fails. If I try three after a failure, it fails, but it can be made to work by first resorting to just two inputs and then adding the third. This is my (anonymised) code: import arcpy from arcpy . sa import * arcpy . CheckOutExtension ( "Spatial" ) arcpy . env . extent = "MAXOF" arcpy . env . overwriteOutput = True arcpy . env . qualifiedFieldNames = False arcpy . env . workspace = "C:\aaaa" inRasDir = bbb\\" outRasDir = ccc\\" terms = { } terms = { 'first' : [ '0.025' , '*' , '+' , ] , 'second' : [ '0.025' , '*' , '+' , ] , 'third' : [ '0.025' , '*' , '+' , ] , 'fourth' : [ '0.15' , '*' , '-' ] } def coeff ( lyr ) : if terms [ lyr ] [ 2 ] == '+' : c = float ( terms [ lyr ] [ 0 ] ) else : c = - float ( terms [ lyr ] [ 0 ] ) return c ; outRas = ( ( ( coeff ( 'first' ) * Raster ( inRasDir + "first" ) ) + ( coeff ( 'second' ) * Raster ( inRasDir + "second" ) ) + ( coeff ( 'third' ) * Raster ( inRasDir + "third" ) ) # + (coeff('fourth') * Raster(inRasDir + "fourth")) ) # * 'fifth' ) outRas . save ( outRasDir + "finalRas" ) The error on the last line is: outRas . save ( outRasDir + "finalRas" ) RuntimeError : ERROR 999998 : Unexpected Error . I also tried using CellStatistics but that also failed on the save statement. Any comments would be gratefully received. Regards, Ian

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

arcpy.sa save random ERROR 999998: Unexpected Error

IanMCDONALD

I am working on a pilot for a system that will eventually loop over a number of raster files to produce a single output one. If I have two input rasters, the routine works reliably. If I have more than three, it always fails. If I try three after a failure, it fails, but it can be made to work by first resorting to just two inputs and then adding the third. This is my (anonymised) code:

<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">from</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa <SPAN class="keyword token">import</SPAN> <SPAN class="operator token">*</SPAN>

arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Spatial"</SPAN><SPAN class="punctuation token">)</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>extent <SPAN class="operator token">=</SPAN> <SPAN class="string token">"MAXOF"</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>qualifiedFieldNames <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">False</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> <SPAN class="string token">"C:\aaaa"</SPAN>

inRasDir <SPAN class="operator token">=</SPAN> bbb\\"
outRasDir <SPAN class="operator token">=</SPAN> ccc\\"

terms <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">{</SPAN><SPAN class="punctuation token">}</SPAN>
terms <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">{</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="string token">'first'</SPAN>&nbsp; <SPAN class="punctuation token">:</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="string token">'0.025'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'*'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'+'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="string token">'second'</SPAN> <SPAN class="punctuation token">:</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="string token">'0.025'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'*'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'+'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="string token">'third'</SPAN>&nbsp; <SPAN class="punctuation token">:</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="string token">'0.025'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'*'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'+'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="string token">'fourth'</SPAN>&nbsp; <SPAN class="punctuation token">:</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="string token">'0.15'</SPAN><SPAN class="punctuation token">,</SPAN>&nbsp; <SPAN class="string token">'*'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'-'</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="punctuation token">}</SPAN>

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">coeff</SPAN><SPAN class="punctuation token">(</SPAN>lyr<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> terms<SPAN class="punctuation token">[</SPAN>lyr<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">2</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">==</SPAN> <SPAN class="string token">'+'</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; c <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN>terms<SPAN class="punctuation token">[</SPAN>lyr<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; c <SPAN class="operator token">=</SPAN> <SPAN class="operator token">-</SPAN>float<SPAN class="punctuation token">(</SPAN>terms<SPAN class="punctuation token">[</SPAN>lyr<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">return</SPAN> c<SPAN class="punctuation token">;</SPAN>

outRas <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">(</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="punctuation token">(</SPAN>coeff<SPAN class="punctuation token">(</SPAN><SPAN class="string token">'first'</SPAN><SPAN class="punctuation token">)</SPAN>&nbsp; <SPAN class="operator token">*</SPAN> Raster<SPAN class="punctuation token">(</SPAN>inRasDir <SPAN class="operator token">+</SPAN> <SPAN class="string token">"first"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="operator token">+</SPAN> <SPAN class="punctuation token">(</SPAN>coeff<SPAN class="punctuation token">(</SPAN><SPAN class="string token">'second'</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">*</SPAN> Raster<SPAN class="punctuation token">(</SPAN>inRasDir <SPAN class="operator token">+</SPAN> <SPAN class="string token">"second"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="operator token">+</SPAN> <SPAN class="punctuation token">(</SPAN>coeff<SPAN class="punctuation token">(</SPAN><SPAN class="string token">'third'</SPAN><SPAN class="punctuation token">)</SPAN>&nbsp; <SPAN class="operator token">*</SPAN> Raster<SPAN class="punctuation token">(</SPAN>inRasDir <SPAN class="operator token">+</SPAN> <SPAN class="string token">"third"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="comment token">#&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; + (coeff('fourth')&nbsp; * Raster(inRasDir + "fourth"))</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="punctuation token">)</SPAN>
<SPAN class="comment token">#&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; * 'fifth'</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="punctuation token">)</SPAN>

outRas<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>outRasDir <SPAN class="operator token">+</SPAN> <SPAN class="string token">"finalRas"</SPAN><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>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

The error on the last line is:

&nbsp;&nbsp;&nbsp; outRas<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>outRasDir <SPAN class="operator token">+</SPAN> <SPAN class="string token">"finalRas"</SPAN><SPAN class="punctuation token">)</SPAN>
RuntimeError<SPAN class="punctuation token">:</SPAN> ERROR <SPAN class="number token">999998</SPAN><SPAN class="punctuation token">:</SPAN> Unexpected Error<SPAN class="punctuation token">.</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

I also tried using CellStatistics but that also failed on the save statement.

Any comments would be gratefully received.

Regards,

Ian

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

Accepted answers

IanMCDONALD

Yes, I was beginning to think memory management might be a problem. I have taken your advice to simplify my test system. It is now just a simple loop summing all the rasters. The erratic behaviour appears to have been resolved by using a 64-bit version of Python.

All comments

IanMCDONALD

DanPatterson_Retired

I suspect system resources then...

but simplify things

<SPAN class="keyword token">import</SPAN> arcpy

a <SPAN class="operator token">=</SPAN> np<SPAN class="punctuation token">.</SPAN>random<SPAN class="punctuation token">.</SPAN>randint<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">10</SPAN><SPAN class="punctuation token">,</SPAN> size<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="number token">3</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">3</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN><SPAN class="number token">0.1</SPAN>

b <SPAN class="operator token">=</SPAN> a<SPAN class="operator token">*</SPAN><SPAN class="number token">2</SPAN>

c <SPAN class="operator token">=</SPAN> a<SPAN class="operator token">*</SPAN>np<SPAN class="punctuation token">.</SPAN>random<SPAN class="punctuation token">.</SPAN>rand<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>

d <SPAN class="operator token">=</SPAN> a<SPAN class="operator token">*</SPAN>np<SPAN class="punctuation token">.</SPAN>random<SPAN class="punctuation token">.</SPAN>rand<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>

func <SPAN class="operator token">=</SPAN> a <SPAN class="operator token">+</SPAN> <SPAN class="number token">1.0</SPAN> <SPAN class="operator token">+</SPAN> b <SPAN class="operator token">*</SPAN> <SPAN class="number token">2</SPAN> <SPAN class="operator token">-</SPAN> c <SPAN class="operator token">+</SPAN> d<SPAN class="operator token">/</SPAN><SPAN class="number token">2</SPAN>

func
Out<SPAN class="punctuation token">[</SPAN><SPAN class="number token">14</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">:</SPAN> 
array<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">3.24388419</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">5.03899154</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">2.34633051</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="punctuation token">[</SPAN><SPAN class="number token">3.24388419</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">5.03899154</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">4.5902147</SPAN> <SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="punctuation token">[</SPAN><SPAN class="number token">4.5902147</SPAN> <SPAN class="punctuation token">,</SPAN> <SPAN class="number token">3.69266102</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">1.44877684</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">]</SPAN><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>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Now, replace a, b, c, d with

RasterToNumPyArray—Help | ArcGIS Desktop

RasterToNumPyArray (in_raster, {lower_left_corner}, {ncols}, {nrows}, {nodata_to_value})

and func with

NumPyArrayToRaster—Help | ArcGIS Desktop

NumPyArrayToRaster (in_array, {lower_left_corner}, {x_cell_size}, {y_cell_size}, {value_to_nodata})

Equivalent functions exist in PRO.

I have tons of NumPy examples on my blog

IanMCDONALD

Hi Dan. Thanks for your prompt response.

It was working, as far as the path definition was concerned (see the '\\' at the end of the inRasDir and outRasDir definitions). It was generating an output raster in the directory I expected and I could view it in ArcMap.

I have taken your advice and tried both the raw formatting and forward slashes but the outcome is the same. I can produce an output raster from up to three inputs but not four. If it fails with four, and I immediately try with three, it fails. If I wait for about a minute, without altering the code, it succeeds with three.

DanPatterson_Retired

It is hard to tell with your anon-stuff but this certainly doesn't work...

You need to use raw formatting for paths or os.join

<SPAN class="string token">"C:\aaaa"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"finalRas"</SPAN>
<SPAN class="string token">'C:\x07aaafinalRas'</SPAN>&nbsp; <SPAN class="comment token"># ---- not good</SPAN>

<SPAN class="string token">"C:\aaaa\" + "</SPAN>finalRas"
<SPAN class="string token">'C:\x07aaafinalRas'</SPAN>
&nbsp; File <SPAN class="string token">"<ipython-input-11-8b220ffe2f0d>"</SPAN><SPAN class="punctuation token">,</SPAN> line <SPAN class="number token">1</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="string token">"C:\aaaa\" + "</SPAN>finalRas"
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="operator token">^</SPAN>
SyntaxError<SPAN class="punctuation token">:</SPAN> invalid syntax

<SPAN class="comment token"># This works</SPAN>
<SPAN class="comment token"># ---- you need raw formatting or use forward slashes in the pathe</SPAN>
r<SPAN class="string token">"C:\aaaa"</SPAN> <SPAN class="operator token">+</SPAN> r<SPAN class="string token">"\finalRas"</SPAN>
Out<SPAN class="punctuation token">[</SPAN><SPAN class="number token">12</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="string token">'C:\\aaaa\\finalRas'</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>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

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.