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

Raster tif count pixels of certain values

AnneTrampe

Hi,

I have a raster tif with different values and want to count the pixels in the images, but for certain value classes

for example count the pixels of value: 1-3, 3-15, 15-50, 50-200, 200-400, 400-1000

Any idea, how this could been done with arcgis tools. I am not really good in scripting anything, so would prefer to us only a few tools. Also I have to proceed for 25 images, so I am looking for a simple and fast solution

Find more posts tagged with

Accepted answers

XanderBakker

You will probably have to reclassifiy the rasters See Reclassify—Help | ArcGIS Desktop

This will create a new raster of type integer that will have an attribute table which will provide the number of pixels for each class.

To do this on 25 rasters you can use Model Builder see What is ModelBuilder?—Help | ArcGIS Desktop or a few lines of Python code; see What is Python?—Help | ArcGIS Desktop

All comments

XanderBakker

You will probably have to reclassifiy the rasters See Reclassify—Help | ArcGIS Desktop

This will create a new raster of type integer that will have an attribute table which will provide the number of pixels for each class.

To do this on 25 rasters you can use Model Builder see What is ModelBuilder?—Help | ArcGIS Desktop or a few lines of Python code; see What is Python?—Help | ArcGIS Desktop

TimothyHales

Check out the Lookup tool - Lookup—Help | ArcGIS Desktop Keep in mind that it requires an integer raster.

-Credit to Raster Master Jeffrey Swain‌

Being more of a feature guy myself, I would have said to convert the raster to points and then pull the values from the output feature class.

XanderBakker

Actually, the LookUp tool will not give you the result you are looking for. It will create a raster based on an attribute value of the integer raster. This will not provide you the pixel count fro each class that you specified. Converting the raster to point, could potentially a huge set of points and will still require you to classify the values or define where clauses and do counts. To do this on 25 point featureclasses (derived from the rasters) would still require some programming or the use of Model Builder.

Although you said that you are not very experienced in programming, below an example Python script that could do the job. Just adjust the folder on line 6 where you have stored you TIFFs:

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">main</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">import</SPAN> arcpy
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">import</SPAN> os

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># folder with input rasters (TIF)</SPAN>
&nbsp;&nbsp;&nbsp; ws <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\GeoNet\CompareDEMs\DEM'</SPAN>&nbsp; <SPAN class="comment token">### change this folder</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># definition of the classes</SPAN>
&nbsp;&nbsp;&nbsp; classes <SPAN class="operator token">=</SPAN> <SPAN class="string token">"1 3 1;3 15 2;15 50 3;50 200 4;200 400 5;400 1000 6"</SPAN>

&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> ws
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>
&nbsp;&nbsp;&nbsp; ws_tmp <SPAN class="operator token">=</SPAN> <SPAN class="string token">'IN_MEMORY'</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Spatial"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># create a list of rasters</SPAN>
&nbsp;&nbsp;&nbsp; lst_rasters <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># loop through the list of rasters</SPAN>
&nbsp;&nbsp;&nbsp; cnt <SPAN class="operator token">=</SPAN> <SPAN class="number token">0</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> raster_name <SPAN class="keyword token">in</SPAN> lst_rasters<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; cnt <SPAN class="operator token">+=</SPAN> <SPAN class="number token">1</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># classify the raster</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; raster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws<SPAN class="punctuation token">,</SPAN> raster_name<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; out_ras <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_tmp<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"reclass{0}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>cnt<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>Reclassify_sa<SPAN class="punctuation token">(</SPAN>raster<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"VALUE"</SPAN><SPAN class="punctuation token">,</SPAN> classes<SPAN class="punctuation token">,</SPAN> out_ras<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"NODATA"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># print the resulting statistics</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"\n"</SPAN><SPAN class="punctuation token">,</SPAN> raster_name
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; PrintStats<SPAN class="punctuation token">(</SPAN>out_ras<SPAN class="punctuation token">)</SPAN>


<SPAN class="keyword token">def</SPAN> <SPAN class="token function">PrintStats</SPAN><SPAN class="punctuation token">(</SPAN>ras<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># get the statistics and print them</SPAN>
&nbsp;&nbsp;&nbsp; fld_val <SPAN class="operator token">=</SPAN> <SPAN class="string token">'Value'</SPAN>
&nbsp;&nbsp;&nbsp; fld_cnt <SPAN class="operator token">=</SPAN> <SPAN class="string token">'Count'</SPAN>
&nbsp;&nbsp;&nbsp; flds <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">(</SPAN>fld_val<SPAN class="punctuation token">,</SPAN> fld_cnt<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; curs <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>ras<SPAN class="punctuation token">,</SPAN> flds<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> row <SPAN class="keyword token">in</SPAN> curs<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; value <SPAN class="operator token">=</SPAN> row<SPAN class="punctuation token">.</SPAN>getValue<SPAN class="punctuation token">(</SPAN>fld_val<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; count <SPAN class="operator token">=</SPAN> row<SPAN class="punctuation token">.</SPAN>getValue<SPAN class="punctuation token">(</SPAN>fld_cnt<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> value<SPAN class="punctuation token">,</SPAN> count


<SPAN class="keyword token">if</SPAN> __name__ <SPAN class="operator token">==</SPAN> <SPAN class="string token">'__main__'</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; main<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>‍</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>

GünterDörffel

Hi,

making use of the Raster Attribute Table ( a rather old concept) would be nice here ... but fails when trying to export it. Based on Xanders code here is a modified version that exports the actual values to an excel file (after writing them to an in_memory table). In addition, the ListRasters has been modified to only read tifs - as in the original question ... AND the resulting table in excel will have two additional fields that reflect the re-classification done in the "classes" variable by Xander.

Hope this helps. I added it as a python file ...

reclass_and_count.py.zip

Join the Conversation

Signed in members can post, follow updates, and more. New here? Register a free account.

Quick Links

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