I created contours by using SRTM30 data. Now I need to find area between two contours like in between 0 m to 50 m depth and 50 m to 100 m.
You might also want to refer to the create Field Contours page in which you will find a comment about script modifications to allowed for customized contour intervals.
One of these builds a "filled polygon contour" from a raster dem. Not exactly what you are looking for but may be useful.
Here is a vector-based suggested workflow to compute the area for each contour interval in your data set. You might want to look for some "easier" to implement solutions in other sources - here is one.
As to the suggested work flow:
- Start with converting your contours to polygons (Feature to Polygon tool) and add a field named: "Within" of an integer type. This next phase should be based on a model or a python script. It is rather simple, thus model would do the job. Build a model that will
- Iterate rows in your new polygon layer
- for each row - create a temporary layer using Make Feature Layer; use the Get Field Value tool to extract the contour value of the temporary layer into a model variable.
- Use select layer by location tool to select from your polygons layer the features that are WITHIN_CLEMENTINI your temporary layer.
- Use select by attributes tool and define REMOVE_FROM_SELECTION as your selection_type. Use an expression that will unselect (remove) all selected features that ther contour value distance from your row contour value is bigger than the contours interval. For example: Let's say that you have contours every 15 meters. Than you should write an expression that will select all [Temporary_Contour_Value]-[Contour]>15
- Now you have the contour the "biggest"/"highest" contour within your row. Use get Field value again to extract the selection ID, you might need to write it as temporary layer. Than write the extracted ID value to the within field in your iterated row. You might consider to use the Calculate field tool.
- Last phase of each loop is to CLEAR_SELECTION using Select layer by attributes.
Now after the model have run you have, for each row, a field which tells you which contour polygon is within it (by ID). All you have to do is compute shape are (you can use the python expression !shape.area! in the field calculator - it is a double/float type). Make a copy of this layer and join it to itself when the PK is your within field and the FK is your original ID. Use field calculator again to compute the area within two contours, bu expression Area_1 (of origianl layer) - Area_2 (of the joint layer).
You can also incorparate all steps into one model. It is just a workflow, since I havn't tried it since I don't have arcmap on my computer. Some difficulties might come up with the "get field value" variables - so you might want to use .py after all (with the same workflow).