topic Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240 in Python Questions

Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

GurminderBharani1 — Sat, 11 Dec 2021 06:13:58 GMT

Unable to save a 3 dimensional array to stacked raster .img file.

I have this array:

arr.shape (10, 148, 160)

Now when I try to save this into raster .img format:

out_ras = arcpy.NumPyArrayToRaster(arr) 
out_ras.save(r"C:\Users\Gurminder\Documents\study\IWMI\HomeTools\Dump\Converted_arr\stack.img")

I get the following error:

Runtime error Traceback (most recent call last): File "<string>", line 1, in <module> RuntimeError: 
ERROR 010240: Could not save raster dataset to C:\Users\Gurminder\Documents\study\IWMI\HomeTools\Dump\
Converted_arr\stack.img with output format IMAGINE Image.

Although a file named stack.img is created in the given directory but all it contains is garbage values.

I have been facing this issue since I have upgraded from ArcMap 10.3 to 10.4, before upgrade everything was working just fine.

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

DanPatterson_Retired — Sat, 11 Dec 2021 06:14:01 GMT

That error message is a generic one normally associated with memory issues or file path issues when working with raster data. (search on error 010240 to see what I mean) But your array is only 10 bands with 148x160 in size, so that shouldn't be an issue...

I would suggest saving the file to a shorter path as a first step.

Secondly, try saving one of the bands to a new shorter path.

Third, you could definitely have something wrong with one of the bands that is producing garbage in the final array, so you might want to throw in some print statements

Finally create this 10 band raster and see if you can save it

>>> import numpy as np
>>> np.set_printoptions(edgeitems=2,linewidth=80,precision=2,
...                     suppress=True,threshold=20,
...                     formatter={'float': '{: 0.3f}'.format})
>>> 
>>> a = np.arange(10*5*6).reshape(10,5,6)
>>> a
array([[[  0,   1, ...,   4,   5],
        [  6,   7, ...,  10,  11],
        ..., 
        [ 18,  19, ...,  22,  23],
        [ 24,  25, ...,  28,  29]],


       [[ 30,  31, ...,  34,  35],
        [ 36,  37, ...,  40,  41],
        ..., 
        [ 48,  49, ...,  52,  53],
        [ 54,  55, ...,  58,  59]],


       ..., 
       [[240, 241, ..., 244, 245],
        [246, 247, ..., 250, 251],
        ..., 
        [258, 259, ..., 262, 263],
        [264, 265, ..., 268, 269]],


       [[270, 271, ..., 274, 275],
        [276, 277, ..., 280, 281],
        ..., 
        [288, 289, ..., 292, 293],
        [294, 295, ..., 298, 299]]])
‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Some test results that you can also run on yours

>>> a = np.arange(10*5*6).reshape(10,5,6)
>>> a[0]
array([[ 0,  1,  2,  3,  4,  5],
       [ 6,  7,  8,  9, 10, 11],
       [12, 13, 14, 15, 16, 17],
       [18, 19, 20, 21, 22, 23],
       [24, 25, 26, 27, 28, 29]])
>>> a[-1]
array([[270, 271, 272, 273, 274, 275],
       [276, 277, 278, 279, 280, 281],
       [282, 283, 284, 285, 286, 287],
       [288, 289, 290, 291, 292, 293],
       [294, 295, 296, 297, 298, 299]])
>>> np.mean(a,axis=0)
array([[ 135.,  136.,  137.,  138.,  139.,  140.],
       [ 141.,  142.,  143.,  144.,  145.,  146.],
       [ 147.,  148.,  149.,  150.,  151.,  152.],
       [ 153.,  154.,  155.,  156.,  157.,  158.],
       [ 159.,  160.,  161.,  162.,  163.,  164.]])
>>> np.ptp(a,axis=0)
array([[270, 270, 270, 270, 270, 270],
       [270, 270, 270, 270, 270, 270],
       [270, 270, 270, 270, 270, 270],
       [270, 270, 270, 270, 270, 270],
       [270, 270, 270, 270, 270, 270]])
‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

GurminderBharani1 — Sat, 11 Dec 2021 06:14:03 GMT

Hello Dan,

Thank you for the prompt reply. I have made the array suggested by you and still I am getting the same error while saving it.

>>> import numpy as np  
np.set_printoptions(edgeitems=2,linewidth=80,precision=2,  
                    suppress=True,threshold=20,  
                    formatter={'float': '{: 0.3f}'.format})  
  
a = np.arange(10*5*6).reshape(10,5,6)  
out_ras = arcpy.NumPyArrayToRaster(a) 
out_ras.save(r"C:\Users\Gurminder\Documents\study\IWMI\HomeTools\Dump\Converted_arr\stack_geonet.img") 
Runtime error 
Traceback (most recent call last):
  File "<string>", line 8, in <module>
RuntimeError: ERROR 010240: Could not save raster dataset to C:\Users\Gurminder\Documents\study\IWMI\HomeTools\Dump\Converted_arr\stack_geonet.img with output format IMAGINE Image.

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

DanPatterson_Retired — Mon, 27 Jun 2016 09:57:18 GMT

what about my other suggestions? Use the array I gave you to try...also, try something other than *.img

I would suggest saving the file to a shorter path as a first step.

Secondly, try saving one of the bands to a new shorter path.

Third, you could definitely have something wrong with one of the bands that is producing garbage in the final array, so you might want to throw in some print statements

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

NeilAyres — Mon, 27 Jun 2016 09:58:18 GMT

As Dan has said, try a much shorter path to save.

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

GurminderBharani1 — Tue, 28 Jun 2016 04:05:50 GMT

I tried shorter path with img as well as tif but still I get the same error

>>> import numpy as np
np.set_printoptions(edgeitems=2,linewidth=80,precision=2,
                    suppress=True,threshold=20,
                    formatter={'float': '{: 0.3f}'.format})

a = np.arange(10*5*6).reshape(10,5,6)
out_ras = arcpy.NumPyArrayToRaster(a)
out_ras.save("C:\\Users\\Gurminder\\Desktop\\stack_geonet.img")
Runtime error
Traceback (most recent call last):
File "<string>", line 8, in <module>
RuntimeError: ERROR 010240: Could not save raster dataset to C:\Users\Gurminder\Desktop\stack_geonet.img with output format IMAGINE Image.
>>> import numpy as np
np.set_printoptions(edgeitems=2,linewidth=80,precision=2,
                    suppress=True,threshold=20,
                    formatter={'float': '{: 0.3f}'.format})

a = np.arange(10*5*6).reshape(10,5,6)
out_ras = arcpy.NumPyArrayToRaster(a)
out_ras.save("C:\\Users\\Gurminder\\Desktop\\stack_geonet.tif")
Runtime error
Traceback (most recent call last):
File "<string>", line 8, in <module>
RuntimeError: ERROR 010240: Could not save raster dataset to C:\Users\Gurminder\Desktop\stack_geonet.tif with output format TIFF.

I have tried converting single band raster to array and back to raster. Everything works fine when it comes to single band.

Taking the array as you have mentioned cannot contain garbage values still it gives the same error so garbage values does not exist in the array "a".

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

GurminderBharani1 — Tue, 28 Jun 2016 04:08:10 GMT

The small path that tried was up to three directories, still I get the error. Is it short enough. I have printed the output to Dan's reply you have check the code there.

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

DanPatterson_Retired — Tue, 28 Jun 2016 04:14:38 GMT

try saving to

r"c:\test\x.img"

r"c:\test\x.tif"

shorter and not in the same path and your file may be barking at the length of its name or the underscore or anything else... let's stretch a little in the things to test

In short, it doesn't like your filename or the saving location, or you don't have permissions OR are you saving to the DESKTOP????

If so .... no ... Bill Gates may say it is ok, but it isn't

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

GurminderBharani1 — Tue, 28 Jun 2016 06:32:14 GMT

Still having the same issue.

>>> import numpy as np
... np.set_printoptions(edgeitems=2,linewidth=80,precision=2,
...                     suppress=True,threshold=20,
...                     formatter={'float': '{: 0.3f}'.format})
...
... a = np.arange(10*5*6).reshape(10,5,6)
... out_ras = arcpy.NumPyArrayToRaster(a)
... out_ras.save(r"C:\test\x.img")
...
Runtime error
Traceback (most recent call last):
File "<string>", line 8, in <module>
RuntimeError: ERROR 010240: Could not save raster dataset to C:\test\x.img with output format IMAGINE Image.
>>> import numpy as np
... np.set_printoptions(edgeitems=2,linewidth=80,precision=2,
...                     suppress=True,threshold=20,
...                     formatter={'float': '{: 0.3f}'.format})
...
... a = np.arange(10*5*6).reshape(10,5,6)
... out_ras = arcpy.NumPyArrayToRaster(a)
... out_ras.save(r"C:\test\x.tif")
...
Runtime error
Traceback (most recent call last):
File "<string>", line 8, in <module>
RuntimeError: ERROR 010240: Could not save raster dataset to C:\test\x.tif with output format TIFF.

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

DanPatterson_Retired — Tue, 28 Jun 2016 10:32:28 GMT

Ahhhh

I will let you figure out what line you left out

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

GurminderBharani1 — Tue, 28 Jun 2016 10:51:19 GMT

Hello Dan,

Are you suggesting that import arcpy was missing.

I have tried with import statement still its giving me the same error.

One question Dan, what version of ArcMap are you using?

>>> import numpy as np
... import arcpy
... np.set_printoptions(edgeitems=2,linewidth=80,precision=2,
...                     suppress=True,threshold=20,
...                     formatter={'float': '{: 0.3f}'.format})
...
... a = np.arange(10*5*6).reshape(10,5,6)
... out_ras = arcpy.NumPyArrayToRaster(a)
... out_ras.save(r"C:\test\x.img")
...
Runtime error
Traceback (most recent call last):
File "<string>", line 9, in <module>
RuntimeError: ERROR 010240: Could not save raster dataset to C:\test\x.img with output format IMAGINE Image.
>>> import numpy as np
... import arcpy
... a = np.arange(10*5*6).reshape(10,5,6)
... out_ras = arcpy.NumPyArrayToRaster(a)
... out_ras.save(r"C:\test\x.tif")
...
Runtime error
Traceback (most recent call last):
File "<string>", line 5, in <module>
RuntimeError: ERROR 010240: Could not save raster dataset to C:\test\x.tif with output format TIFF.
>>>

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

DanPatterson_Retired — Tue, 28 Jun 2016 10:57:02 GMT

If you are doing it inside of arcmap's python IDE, you shouldn't need it... but you didn't specify where you were running it.

So... it should work without the import arcpy.

I am using ArcMap 10.4.1 and I did the calculation in pythonwin running python 3.4, but it worked in 2.7 as well

I basically think your setup is ....

So the only thing I can suggest is go try it on someone else's machine. The code works, an image was produced and you will just keep spinning your wheels until you can prove that you can save that basic stacked array.

Good luck

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

DanPatterson_Retired — Tue, 28 Jun 2016 11:04:55 GMT

Xander Bakker can you try the code? I can't flag or redirect elsewhere since you are owner of the group

Neil Ayres or Peter Wilson can one of you guys on older Arc* installs, give it a shot, since I know you have worked with arrays

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

XanderBakker — Sat, 11 Dec 2021 06:14:06 GMT

I have moved the thread to the Python space and tested the code in PyScripter with Python 2.7.10 and ArcGIS 10.4.1. Writing to TIF, IMG and FGDB raster all resulted in the same error message. In case of the IMG format the raster was created with all bands filled with NoData, in case of the TIFF all bands had value 0 and in case of the FGDB raster the output was not created.

Reducing the number of bands to 3 did not yield a different result (same error), but when using a single band it will be created without problem. You could write each band separately and use arcpy.CompositeBands_management to create the multiband output, but I'm pretty sure that's not what you're after.

When you try to create a composite bands using the bands of the multi band temp output, it will create the raster without error message, but it will contain all NoData values:

import numpy as np
import arcpy
import os
np.set_printoptions(edgeitems=2,linewidth=80,precision=2,
    suppress=True,threshold=20,
    formatter={'float': '{: 0.3f}'.format})
a = np.arange(10*5*6).reshape(10,5,6)
out_ras = arcpy.NumPyArrayToRaster(a)

desc = arcpy.Describe(out_ras)
lst_ras = [child.catalogPath for child in desc.children]
in_ras = ';'.join(lst_ras)
arcpy.CompositeBands_management(in_rasters=in_ras, out_raster=r"D:\Temp\ras10bands.tif")

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

DanPatterson_Retired — Tue, 28 Jun 2016 12:49:45 GMT

Bizarre... so the key seems to be having Unicode (with 2.7) or just use python 3.4.x since it sure has worked for me, regardless of the output format (*.img or *.tiff. I get a 10 band raster. I suppose another reason to upgrade python.

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

NeilAyres — Tue, 28 Jun 2016 14:12:50 GMT

Just moved to 10.4.1 last week.

Seems that Xander's work around will do the job. But perhaps this is an undocumented feature in 10.4.1

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

DanPatterson_Retired — Tue, 28 Jun 2016 14:27:04 GMT

Neil... Xander is using 10.4.1 but python 2.7.x

It is either the python version (3.4.x ergo Unicode enabled) or the user

I can't narrow it down further... I normally use *.npy files to store multiband rasters/arrays, they are way smaller and I rarely need to map something.

These three are the same array in npy, img and tif formats

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

DanPatterson_Retired — Sat, 11 Dec 2021 06:14:09 GMT

WELLLL

Python version....

>>> import numpy as np
>>> import arcpy
>>> a = np.arange(10*5*6).reshape(10,5,6)
>>> import sys; print(sys.version)
2.7.10 (default, May 23 2015, 09:40:32) [MSC v.1500 32 bit (Intel)]
>>> 
>>> out_ras = arcpy.NumPyArrayToRaster(a)
RuntimeError: ERROR 010240: Could not save raster dataset to C:\Users\Dan\Documents\ArcGIS\Default.gdb\numpy with output format FGDBR.
>>>
‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Install ArcGIS PRO and you will get Python 3.4.3 (or you can install it yourself and fiddle with associations)

Anaconda is coming soon.

If you don't need to 'see your stuff' use the *.npy format

EDIT

python 2.7 still being distributed with arc 10.5, so the PRO install will do until 3.x is distributed with arcmap

Re: Unable to save 3 dimensional array as stacked raster .img file in arcpy, RuntimeError: ERROR 010240

DanPatterson_Retired — Sat, 17 Feb 2018 01:37:50 GMT

Some time later... I new I saw it somewhere!

This solves the dilemma of where to write temporary tif files... The describe object was the key

Nice... Thanks Xander Bakker

E‌xcept when you bring the stack back in with RasterToNumPyArray, only the first three bands have data, the rest are zero...

Very strange