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

RasterToNumPyArray gives an error

GraziaZulian

I'm having problems by using RasterToNumPyArray.

With the smal test data base it works, but if I try with a huge dataset (columns and rows: 67000, 58000, cellsize: 100 m) it gives me this error:

Traceback (most recent call last):
File "C:\Python26\ArcGIS10.0\Lib\site-packages\pythonwin\pywin\framework\scriptutils.py", line 309, in RunScript
    debugger.run(codeObject, __main__.__dict__, start_stepping=0)
File "C:\Python26\ArcGIS10.0\Lib\site-packages\pythonwin\pywin\debugger\__init__.py", line 60, in run
    _GetCurrentDebugger().run(cmd, globals,locals, start_stepping)
File "C:\Python26\ArcGIS10.0\Lib\site-packages\pythonwin\pywin\debugger\debugger.py", line 624, in run
    exec cmd in globals, locals
File "D:\ES_Map\recreationnew\RecNumpy2.py", line 31, in <module>
    in1 = arcpy.RasterToNumPyArray(arcpy.env.workspace + "\\eu_coastna") #arcpy.RasterToNumPyArray(input1)
File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\__init__.py", line 1113, in RasterToNumPyArray
    return _RasterToNumPyArray(*args, **kwargs)
RuntimeError: ERROR 999998: Unexpected Error.

could you help me please.
here is a simple example of my code

RecNumpy4.py.zip

Finde mehr Beiträge mit dem Tag

Accepted answers

All comments

DavidWilton

I ran into this issue so I wrote an iterator class to do this. In my case I wanted to return the X,Y,Z values but you could easily modify to return the numpy array.

It can be used simply like this:

<SPAN class="keyword token">for</SPAN> block <SPAN class="keyword token">in</SPAN> RasterBlockIterator<SPAN class="punctuation token">(</SPAN>raster<SPAN class="punctuation token">,</SPAN> no_data_val<SPAN class="operator token">=</SPAN>self<SPAN class="punctuation token">.</SPAN>no_data_val<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
   <SPAN class="keyword token">for</SPAN> point <SPAN class="keyword token">in</SPAN> block<SPAN class="punctuation token">:</SPAN>
      x <SPAN class="operator token">=</SPAN> point<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN>
      y <SPAN class="operator token">=</SPAN> point<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN> 
      z <SPAN class="operator token">=</SPAN> point<SPAN class="punctuation token">[</SPAN><SPAN class="number token">2</SPAN><SPAN class="punctuation token">]</SPAN>‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

The iterator class. If you want the numpy array just return the arr on line 62. By default it works in blocks of 1000*1000

<SPAN class="keyword token">class</SPAN> <SPAN class="token class-name">RasterBlockIterator</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="string token">"""Iterator that extracts tiles of data from a Esri Raster."""</SPAN>

    <SPAN class="keyword token">def</SPAN> <SPAN class="token function">__init__</SPAN><SPAN class="punctuation token">(</SPAN>self<SPAN class="punctuation token">,</SPAN> raster<SPAN class="punctuation token">:</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">,</SPAN> tile_size<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="number token">1000</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">1000</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> no_data_val<SPAN class="operator token">=</SPAN>np<SPAN class="punctuation token">.</SPAN>nan<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
        <SPAN class="string token">"""
        :param raster: Raster to iterate
        :param tile_size: size of each block/tile to extrac to numpy array
        :param no_data_val: value in the raster which should be treated as no data
        """</SPAN>
        self<SPAN class="punctuation token">.</SPAN>tile_size <SPAN class="operator token">=</SPAN> tile_size
        self<SPAN class="punctuation token">.</SPAN>raster <SPAN class="operator token">=</SPAN> raster
        self<SPAN class="punctuation token">.</SPAN>no_data_val <SPAN class="operator token">=</SPAN> no_data_val
        self<SPAN class="punctuation token">.</SPAN>cell_size <SPAN class="operator token">=</SPAN> raster<SPAN class="punctuation token">.</SPAN>meanCellHeight
        self<SPAN class="punctuation token">.</SPAN>pos <SPAN class="operator token">=</SPAN> <SPAN class="number token">0</SPAN>

        <SPAN class="comment token"># work out the size of all the block to create. They won't be extracted at this point</SPAN>
        t_rows <SPAN class="operator token">=</SPAN> raster<SPAN class="punctuation token">.</SPAN>height
        t_cols <SPAN class="operator token">=</SPAN> raster<SPAN class="punctuation token">.</SPAN>width
        x_min <SPAN class="operator token">=</SPAN> raster<SPAN class="punctuation token">.</SPAN>extent<SPAN class="punctuation token">.</SPAN>XMin
        y_min <SPAN class="operator token">=</SPAN> raster<SPAN class="punctuation token">.</SPAN>extent<SPAN class="punctuation token">.</SPAN>YMin
        self<SPAN class="punctuation token">.</SPAN>tiles <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
        <SPAN class="keyword token">for</SPAN> row <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">,</SPAN> t_rows<SPAN class="punctuation token">,</SPAN> tile_size<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
            <SPAN class="keyword token">for</SPAN> col <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">,</SPAN> t_cols<SPAN class="punctuation token">,</SPAN> tile_size<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
                colend <SPAN class="operator token">=</SPAN> min<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN>t_cols<SPAN class="punctuation token">,</SPAN> col <SPAN class="operator token">+</SPAN> tile_size<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>
                rowend <SPAN class="operator token">=</SPAN> min<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN>t_rows<SPAN class="punctuation token">,</SPAN> row <SPAN class="operator token">+</SPAN> tile_size<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>
                xstart <SPAN class="operator token">=</SPAN> x_min <SPAN class="operator token">+</SPAN> <SPAN class="punctuation token">(</SPAN>col <SPAN class="operator token">*</SPAN> self<SPAN class="punctuation token">.</SPAN>cell_size<SPAN class="punctuation token">)</SPAN>
                ystart <SPAN class="operator token">=</SPAN> y_min <SPAN class="operator token">+</SPAN> <SPAN class="punctuation token">(</SPAN>row <SPAN class="operator token">*</SPAN> self<SPAN class="punctuation token">.</SPAN>cell_size<SPAN class="punctuation token">)</SPAN>
                self<SPAN class="punctuation token">.</SPAN>tiles<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>namedtuple<SPAN class="punctuation token">(</SPAN><SPAN class="string token">'tile'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'colstart colend rowstart rowend xstart ystart'</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">(</SPAN>
                    colstart<SPAN class="operator token">=</SPAN>col<SPAN class="punctuation token">,</SPAN>
                    colend<SPAN class="operator token">=</SPAN>colend<SPAN class="punctuation token">,</SPAN>
                    rowstart<SPAN class="operator token">=</SPAN>row<SPAN class="punctuation token">,</SPAN>
                    rowend<SPAN class="operator token">=</SPAN>rowend<SPAN class="punctuation token">,</SPAN>
                    xstart<SPAN class="operator token">=</SPAN>xstart<SPAN class="punctuation token">,</SPAN>
                    ystart<SPAN class="operator token">=</SPAN>ystart
                <SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="keyword token">def</SPAN> <SPAN class="token function">__iter__</SPAN><SPAN class="punctuation token">(</SPAN>self<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
        <SPAN class="keyword token">return</SPAN> self

    <SPAN class="keyword token">def</SPAN> <SPAN class="token function">__len__</SPAN><SPAN class="punctuation token">(</SPAN>self<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
        <SPAN class="keyword token">return</SPAN> len<SPAN class="punctuation token">(</SPAN>self<SPAN class="punctuation token">.</SPAN>tiles<SPAN class="punctuation token">)</SPAN>

    <SPAN class="keyword token">def</SPAN> <SPAN class="token function">__next__</SPAN><SPAN class="punctuation token">(</SPAN>self<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">-</SPAN><SPAN class="operator token">></SPAN> List<SPAN class="punctuation token">[</SPAN>Tuple<SPAN class="punctuation token">[</SPAN>float<SPAN class="punctuation token">,</SPAN> float<SPAN class="punctuation token">,</SPAN> float<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">:</SPAN>

        <SPAN class="keyword token">try</SPAN><SPAN class="punctuation token">:</SPAN>
            tile <SPAN class="operator token">=</SPAN> self<SPAN class="punctuation token">.</SPAN>tiles<SPAN class="punctuation token">[</SPAN>self<SPAN class="punctuation token">.</SPAN>pos<SPAN class="punctuation token">]</SPAN>
        <SPAN class="keyword token">except</SPAN> IndexError<SPAN class="punctuation token">:</SPAN>
            <SPAN class="keyword token">raise</SPAN> StopIteration

        lower_l <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>tile<SPAN class="punctuation token">.</SPAN>xstart<SPAN class="punctuation token">,</SPAN> tile<SPAN class="punctuation token">.</SPAN>ystart<SPAN class="punctuation token">)</SPAN>
        ncols <SPAN class="operator token">=</SPAN> tile<SPAN class="punctuation token">.</SPAN>colend <SPAN class="operator token">-</SPAN> tile<SPAN class="punctuation token">.</SPAN>colstart
        rnrows <SPAN class="operator token">=</SPAN> tile<SPAN class="punctuation token">.</SPAN>rowend <SPAN class="operator token">-</SPAN> tile<SPAN class="punctuation token">.</SPAN>rowstart
        arr <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>RasterToNumPyArray<SPAN class="punctuation token">(</SPAN>self<SPAN class="punctuation token">.</SPAN>raster<SPAN class="punctuation token">,</SPAN> nodata_to_value<SPAN class="operator token">=</SPAN>np<SPAN class="punctuation token">.</SPAN>nan<SPAN class="punctuation token">,</SPAN>
                                       lower_left_corner<SPAN class="operator token">=</SPAN>lower_l<SPAN class="punctuation token">,</SPAN>
                                       ncols<SPAN class="operator token">=</SPAN>ncols<SPAN class="punctuation token">,</SPAN> nrows<SPAN class="operator token">=</SPAN>rnrows<SPAN class="punctuation token">)</SPAN>
        <SPAN class="comment token"># Loop over rows and extract X,Y,Z</SPAN>
        output_data <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN>
            <SPAN class="punctuation token">(</SPAN>j <SPAN class="operator token">*</SPAN> self<SPAN class="punctuation token">.</SPAN>cell_size <SPAN class="operator token">+</SPAN> tile<SPAN class="punctuation token">.</SPAN>xstart<SPAN class="punctuation token">,</SPAN> i <SPAN class="operator token">*</SPAN> self<SPAN class="punctuation token">.</SPAN>cell_size <SPAN class="operator token">+</SPAN> tile<SPAN class="punctuation token">.</SPAN>ystart<SPAN class="punctuation token">,</SPAN> value<SPAN class="punctuation token">)</SPAN>
            <SPAN class="keyword token">for</SPAN> i<SPAN class="punctuation token">,</SPAN> row <SPAN class="keyword token">in</SPAN> enumerate<SPAN class="punctuation token">(</SPAN>arr<SPAN class="punctuation token">)</SPAN>
            <SPAN class="keyword token">for</SPAN> j<SPAN class="punctuation token">,</SPAN> value <SPAN class="keyword token">in</SPAN> enumerate<SPAN class="punctuation token">(</SPAN>row<SPAN class="punctuation token">)</SPAN>
            <SPAN class="keyword token">if</SPAN> <SPAN class="operator token">not</SPAN> self<SPAN class="punctuation token">.</SPAN>no_data_val <SPAN class="operator token">==</SPAN> value <SPAN class="operator token">and</SPAN> <SPAN class="operator token">not</SPAN> np<SPAN class="punctuation token">.</SPAN>isnan<SPAN class="punctuation token">(</SPAN>value<SPAN class="punctuation token">)</SPAN>
        <SPAN class="punctuation token">]</SPAN>

        self<SPAN class="punctuation token">.</SPAN>pos <SPAN class="operator token">+=</SPAN> <SPAN class="number token">1</SPAN>
        <SPAN class="keyword token">return</SPAN> output_data‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<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>‍</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>

```

OlgaChesnokova

Hello,

I am trying to replace the holes on the Aster DEM with the values from the Ferranti DEM, and I get the same error as described in the first post. Therefore, I've decided to proceed with tiles. Now I receive the following error:

File "C:\Python27\ArcGIS10.1\Lib\site-packages\Pythonwin\pywin\framework\scriptutils.py", line 320, in RunScript
    debugger.run(codeObject, __main__.__dict__, start_stepping=1)
File "C:\Python27\ArcGIS10.1\Lib\site-packages\Pythonwin\pywin\debugger\__init__.py", line 60, in run
    _GetCurrentDebugger().run(cmd, globals,locals, start_stepping)
File "C:\Python27\ArcGIS10.1\Lib\site-packages\Pythonwin\pywin\debugger\debugger.py", line 655, in run
    exec cmd in globals, locals
File "D:\tests\aster_dem\aster_dem.py", line 53, in <module>
    new_raster = arcpy.CreateRasterDataset_management(arcpy.env.workspace, "aster_fix", cellsize, pixel_type, spatialReference)
File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy\management.py", line 11265, in CreateRasterDataset
    raise e
ExecuteError: ERROR 999999: Error executing function.
Failed to execute (CreateRasterDataset).

I am quite new to the python and arcpy, so the probability that I have some stupid mistakes inside is quite high. Any kind of help would be really appreciated!

code:
import os, sys
import string, math
import numpy
import arcpy
from arcpy.sa import *
# os.chdir(r"D:\tests\aster_dem")
arcpy.env.workspace = r"D:\tests\aster_dem"
arcpy.env.overwriteOutput = True

# input data
ferranti = arcpy.Raster(arcpy.env.workspace + "\\dhm25_ferranti_wgs84.tif")
aster = arcpy.Raster(arcpy.env.workspace + "\\e84_aster_switzerland_clip.tif")

# spatial analyst extension check
arcpy.CheckOutExtension("Spatial")

# difference
difference = aster - ferranti
difference.save(arcpy.env.workspace + "\\difference.tif")
threshold = -300

# data description
spatialReference=arcpy.Describe(ferranti).spatialReference
inRasterDesc = arcpy.Describe(ferranti)

# coordinates of the lower left corner
xmin = inRasterDesc.Extent.XMin
xmax = inRasterDesc.Extent.XMax
ymin = inRasterDesc.Extent.YMin
ymax = inRasterDesc.Extent.YMax

# cell size, raster size
rasMeanCellHeight = inRasterDesc.MeanCellHeight
rasMeanCellWidth = inRasterDesc.MeanCellWidth
nrows = inRasterDesc.Height
ncols = inRasterDesc.Width

cellsize = rasMeanCellHeight

# tiles
tile_size=(128,128)
ntiles = 2
pixel_type = "16_BIT_UNSIGNED"
RasterType = 'TIF'

new_raster = arcpy.CreateRasterDataset_management(arcpy.env.workspace, "aster_fix", cellsize, pixel_type, spatialReference)

for row in range(0,nrows,tile_size[1]):
    for col in range(0,ncols,tile_size[0]*ntiles):
        colend = min([ncols,col+tile_size[0]*ntiles-1])
        rowend = min([nrows,row+tile_size[1]-1])
        xstart = xmin + (cellsize*col)
        ystart = ymax - (cellsize*(rowend+1))
        print 'Processing cols',col,':',colend
        print 'Processing rows',row,':',rowend
        print 'Extracting lower left corner:',xstart,ystart
        ferranti_arr = arcpy.RasterToNumPyArray(ferranti, arcpy.Point(xmin,ystart), ncols, nrows)
        aster_arr = arcpy.RasterToNumPyArray(aster, arcpy.Point(xmin,ystart), ncols, nrows)
        difference_arr = arcpy.RasterToNumPyArray(difference, arcpy.Point(xmin,ystart), ncols, nrows)

        if difference_arr[row,col] < threshold:
            aster_arr[row,col] = ferranti_arr[row,col]

        out_raster = NumPyArrayToRaster(aster_arr, arcpy.Point(xmin,ystart), cellsize, cellsize)
        arcpy.Mosaic_management(out_raster,new_raster)

Thank you in advance!

Regards,
Olga

GraziaZulian

Thanks, I'm working on it now. I'll let you know

Luke_Pinner

Do I loop for rows and cols?

Depends on the type of raster you are using. Some rasters are tiled internally and it is much more efficient to access by tiles or blocks. Geodatabase rasters have 128x128 pixel tiles (by default), GRIDS are 256x4 (col then row), tiffs can be tiled or untiled (safest to process one row at a time).

Here is some code to process by tile (the commented out stuff is arcpy related that I can't test as I'm not at back work until January).

ncols=67000
nrows=58000
cellsize=2
xmin=0
ymax=58000
## ^^ Get all this info from your original raster

tile_size=(128,128) #GDB rasters=(128,128), tif=(ncols,1), GRID=(256,4)
ntiles = 4 #no. of tiles to process at a time

#arcpy.env.Extent=in_raster
#arcpy.env.outputCoordinateSystem=in_raster
#arcpy.env.cellSize=in_raster
#new_raster=arcpy.CreateRasterDataset(out_path, out_name, cellsize, pixel_type, 
#                                     spatial_reference)
for row in range(0,nrows,tile_size[1]):
    for col in range(0,ncols,tile_size[0]*ntiles):
        colend=min([ncols,col+tile_size[0]*ntiles-1])
        rowend=min([nrows,row+tile_size[1]-1])
        xstart=xmin+(cellsize*col)
        ystart=ymax-(cellsize*(rowend+1))
        print 'Processing cols',col,':',colend
        print 'Processing rows',row,':',rowend
        print 'Extracting lower left corner:',xstart,ystart
        #array=RasterToNumPyArray (in_raster, arcpy.Point(xmin,ystart), ncols, step)
        #Do your processing
        #out_raster=NumPyArrayToRaster(array, arcpy.Point(xmin,ystart), cellsize, cellsize)
        #arcpy.Mosaic_management(out_raster,new_raster)

GraziaZulian

Hi,
I need to go back to office as well to try.
I'll let you know tomorrow,
the question is:
Do I loop for rows and cols?
inside the loop I calculate

and then? do I create n rasters and than make the mosaic?
then I could delete the temporary rasters.... this could let me save space.

Could you give a look to the code that I'll post tomorrow?

Luke_Pinner

Sorry, I'm not at work so don't have access to ArcGIS to test an example. I think the mosaic tool will do what you need though:
http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Mosaic_To_New_Raster/001700000098000000/ or http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//001700000097000000

GraziaZulian

Thanks, but I can't figure out how to rebuild the output dataset.
Can you give me an example?

Luke_Pinner

You can specify what part of the raster to extract to a numpy array - RasterToNumPyArray(in_raster, lower_left_corner, ncols, nrows).
http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/NumPyArrayToRaster/000v00000130000000/

So you could loop through and process a few rows at a time.

GraziaZulian

Thanks,
I know that with the map algebra a can complete the same tasks, but I'm looking for a way to create less intermediate datasets and make my process less time consuming (I need 1 week to complete it). I'm working at continental scale with 25 m resolution data (The code I put was just a test).
Have you some suggestions?
Is there a way, for example to work with array by row?

Luke_Pinner

That's because your raster is too big to fit in memory. A 67000*58000 (1 band) 8 bit raster is 4GB, 16 bit is 7GB and 32 bit is 14GB...

If you must use numpy, you will have to clip out smaller sections/tiles of your raster, convert to numpy array, process each tile separately and then convert back to raster and mosaic.

Looking at the commented out section of your code, you don't need to use numpy at all, all those operations can be done in spatial analyst map algebra.

Nimm an der Unterhaltung teil

Angemeldete Mitglieder können Beiträge verfassen, Updates folgen und mehr. Neu hier? Registriere ein kostenloses Konto.

Schnellzugriffe

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