I have a DEM of a study region; I want to modify the DEM in order to reflect the coastline outlook when the sea level was 5 meters higher than today. What I have done so far is to use 'Raster Calculator' with the quite simple following formula: DEM >= 5. Then, in the returned raster, I have set to 'hollow' the 1 values (representing cells with elevation >= 5m), and to 'blue' the 0 values. How and by which tools can I use the obtained raster as a mask to modify the original DEM. In other words, I want to mask out (i.e., remove) places with elevation below 5 m from the original DEM. I am using arcGIS 10.1 with advanced license.
A simpler workflow would be to use Con (SA):
outCon = Con(Raster("elevation") >= 5, 1, 0) arcpy.RasterToPolygon_conversion(outCon, "c:/output/NewLandSea.shp", "SIMPLIFY")
The Con tool is available interactively using Raster Calculator.
The reason for simplifying the raster to a binary (0/1) representation is for Raster to Polygon: a complex raster either will not convert (floating point is not suitable input) or possibly not finish due to high complexity, by making the raster binary the polygonization is simplified resulting in new land/sea polygons.