1

My data consists of several points each with xy coordinates and height.

I used spatial analyst tools -> interpolation -> natural neighbor to get an interpolation map for the height distribution, separated in categories like 10m-20m, 20m-30m and so on. This works fine.

Now, I would like to get the area of those categories. For example, I want to see that the interpolated area in the 10m-20m range is 200m².

3
  • 2
    Use reclass in spatial analyst. You'll end up with integer raster containing values 1,2,..n zones. Convert it to polygon using raster to polygon. If output is feature class it will have shape_area field. Use summary statistics to find totals of area per zone
    – FelixIP
    Commented Jul 13, 2015 at 4:58
  • I converted it to a raster using 'reclassify' and then to a poly successfully - it now looks like contour lines (height) on a map. But the next step does not work. I go to 'Analysis Tools -> Statistics -> Summary Statistics' and choose the new polygon. But in 'statistic field' I only got 'FID, ID and GRIDCODE' but no 'area' field. Why is there no area field?
    – user136036
    Commented Jul 13, 2015 at 8:22
  • Just found out: You have to use 'Spatial Statistic Tools -> Utilities -> Calculate Areas' first on this poly. Then you can use 'Summary Statistics' on the new poly and select the area field.
    – user136036
    Commented Jul 13, 2015 at 8:30

1 Answer 1

1

When calculating areas, be aware in which Coordinate Reference System (CRS) data is defined and pay attention if used tools are honoring the CRS (as opposed to project data on the fly and then, calculate areas).

Also, verify if the calculation method is aligned with your goal. For example, a projected coordinate system will yield planar areas (the area in a cartesian system), and if it is an equal area projection the result will be exact, otherwise, approximate. A geographic coordinate system will yield geodetic areas (the area based on a geoid/ellipsoid).

First, get the height classes from the raster resulted from the interpolation analysis:

  • Use the Reclassify tool to convert continuous values of height into desired classes (bins): 10m-20m, 20m-30m, etc. Set the remap parameter equal to RemapRange, i.e., each value within a given range (class) will be remapped to a common value.

From here, there are two main strategies:


  1. Calculate areas directly from the raster;

If CRS is in a projected coordinate system, one can count how many pixels exist in each class and multiply them by the pixel area. See Measuring area of raster classes?.


  1. Convert the Raster to Polygon and then, calculate areas (for each polygon):

There are several tools (some of them use different methods) one can use to calculate area from polygons:

See the following example from Mahr (2015) - Calculating Geodesic Areas in ArcMap with Field Calculator:

!Shape!.getArea("GEODESIC","SQUAREKILOMETERS")

After calculating areas from polygons, add up (sum) their areas within each height class using the Summary Statistics tool.

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.