[GIS] batch reclassify rasters with different elevation ranges

arcgis-10.0arcgis-desktopmodelbuilderreclassify

I have about 500 elevation rasters, and I would like to use modelbuilder with an iterator to reclassify each raster into distinct elevation ranges (0-1m, 1-2m, 2-3m, 3-4m, 4-5m). I believe my problem is that each raster has a different range of elevation, is there a work around this? I tried mosaic-ing all the rasters together but reclassify tool does not work on this mosaic.

Best Answer

Create a CSV file:

Low,High,Val
0,1,1
1,2,2
2,3,3
3,4,4
4,5,5

It is best to import it to a dbf, INFO or GDB table after entering the text values.

Then use Reclassify by Table, it doesn't care what the ranges of individual rasters are. According to the documentation, build the statistics first then run this tool (iteratively if you like). From_value_field is Low, To_value_field is High and Output_value_field is Val. Change the numbers in the Val column to suit your expected outputs, the values I have shown here are illustrative only.

Depending on what you want in your output you can specify Missing Values: DATA which will copy all values not listed in the table as they are or Missing Values: NODATA which will set all values not in the table to NODATA.

Related Solutions

[GIS] Reclassify many rasters in ModelBuilder

You can actually do this in model builder without python but it is finicky and a bit of a workaround.

create your model with the iterators for the first reclass and make the reclass table a parameter in the model.
right click on your model and set it to batch.
Now you have the reclass table as the parameter.
Use the + button to add the number of reclassses you have (20 different reclasses in this example - add twenty lines with the plus button).
Now in the twenty reclasses boxes on the batch screen load your remap tablesor set the remap in the box (it is kind of like CSV code in this box).

[GIS] Raster Calculator Iteration in ModelBuilder

The easiest solution is to replace the raster calculator with the "Plus" operator in the iterator.

EDIT: note that I wasn't able to reproduce the problem in ArcGIS 10.6:

("%Raster%") + 1

worked fine in my iterator ( with output =, for example, %Raster%_out )

Anyway, the best solution is to use a Python script

import arcpy 
from arcpy.sa import *

arcpy.env.workspace = r'your_path_to_directory'
for r in arcpy.ListRasters("*"): # list all rasters in workspace
    outRaster = Raster(r)+1 #perfoms addition
    outRaster.save(r[:-4] + "out.tif") #save output in tif format

Best Answer

Related Solutions

[GIS] Reclassify many rasters in ModelBuilder

[GIS] Raster Calculator Iteration in ModelBuilder

Related Question