![\"Figure](\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/07\/Dialog.png\")
<\/a>The Con tool is available in the Conditional toolset in the Spatial Analyst toolbox. It is also available from the sa module in arcpy and can also be accessed using the Raster Calculator<\/a> tool.<\/p>\n Now that you have the basics of how Con works<\/a>, let\u2019s start the analysis.<\/p>\n Where are the best locations to grow my coffee?<\/span><\/strong><\/p>\n First, you are interested in locations where the elevations are greater than 1800 feet. If a location meets the criteria, you want to assign the location a \u201c1\u201d (True) and locations not meeting the criteria, you want to assign a \u201c0\u201d (False).<\/p>\n To do this in the Con tool dialog, specify your elevation raster (elevation.tif) as the Input condition raster and the value greater than 1800 feet as the condition in the Expression. To achieve the desired output, you will assign a \u201c1\u201d to the locations with elevations meeting this criterion and a \u201c0\u201d for elevations equal to or lower than 1800 (see Figure 1). To specify the input to the Con tool as a raster, create a raster object<\/a> from the input dataset. It is recommended to first create the raster object from the elevation dataset and then use the resulting object in the subsequent Con statement.<\/p>\n While the preceding statement will work fine, writing it in the following way is preferred, since it allows for some optimizations to be made on execution and thus should run faster:<\/p>\n Both these statements would give you the same output. Any location with an elevation greater than 1800 feet (the where_condition) meets your criteria and will be assigned \u201c1\u201d (the in_true_raster_or_constant), while locations equal to or below 1800 feet are assigned a \u201c0\u201d (the in_false_raster_or_constant).<\/p>\n Notice the second Con statement does not use the where_clause parameter explicitly, but specifies the condition directly to the in_conditonal_raster parameter instead.<\/p>\n Functionality: Defining ranges using Boolean evaluations (\u201c&\u201d) in a complex expression<\/em><\/strong><\/p>\n Once you have the minimum elevation set, you also need to identify the upper elevation limit, 3600 feet. An efficient way to identify the elevation range is to use a Boolean evaluation. You also realize that for subsequent analysis, you wish to know the elevation values for locations within the desired range, not just assign them a \u201c1\u201d, and you do not want to consider locations with elevations outside the range. Therefore, you must alter the above Con statement and replace ‘1’ with ‘elev’ to maintain the elevation values for those locations within the range, and eliminate the other locations from consideration by assigning them ‘NoData’ instead of ‘0’, as below:<\/p>\n In the statement, make note of the following: the condition is not a simple evaluation but a Boolean, the true condition is a raster and not a constant, and no false statement is specified. As a result, locations above 1800 feet and (\u201c&\u201d) below 3600 feet are assigned their elevation value using elev for the input_true_raster_or_constant parameter. By leaving the input_false_raster_or_constant parameter empty, you will assign locations equal to or below 1800 feet or equal to or above 3600 feet to NoData.<\/p>\n Functionality: When to use the where clause in Map Algebra<\/em><\/strong><\/p>\n Let\u2019s now incorporate another criterion in our analysis for identifying the best locations to grow coffee. Not only do you want elevations between 1800 and 3600 feet,<\/strong> you prefer the locations that are open and not forested within that range. You have a land use raster with a field called \u201cCategory\u201d which identifies the land use type for each location. To query a field other than \u201cValue\u201d in a raster, you must use the optional where clause.<\/p>\n In the above statement the condition is based on a field other than \u201cValue\u201d, thus, a where clause is necessary. As a result, all locations that are designated as \u201cOpenArea\u201d land use are assigned the elevation values from outElevRange and all other locations are assigned NoData. As you recall, in outElevRange, only locations with elevation values greater than 1800 and less than 3600 feet have values and all other locations are NoData.<\/p>\n Functionality: Creating a condition from multiple rasters<\/em><\/strong><\/p>\n The above types of conditions are the most common use for Con, but let’s explore a few additional conditions to address more complex problems. Continuing your analysis, <\/strong>since you prefer locations that you can access and that have water nearby, the best locations to grow coffee would be within a mile of a road and a half of mile of a stream.<\/p>\n You first need to create two rasters, one identifying how far each location is from the closest road and a second identifying how far each location is from the closest stream using the Spatial Analyst Euclidean Distance<\/a> tool.<\/p>\n In the Con statement above two different rasters are used in the conditional statement. As a result, only locations that are within a mile of road, within a half a mile (2640 feet) of a stream, have a land use of \u201cOpenArea\u201d, and are within elevations of 1800 and 3600 feet will be assigned their elevation values, while all other locations are assigned NoData.<\/p>\n Functionality: Using a nested Con<\/em><\/strong><\/p>\n Let’s complete the analysis.\u00a0 From the possible locations selected so far, you prefer flatter slopes with more southerly aspects for growing your coffee, but maintaining the southerly aspect is not mandatory.<\/p>\n First, two rasters are created identifying the slope and aspect for each location determined from the elevation surface using the Spatial Analyst Slope<\/a> and Aspect<\/a> tools.\u00a0Then, you will incorporate the slope and aspect preferences into the condition.<\/p>\n This Con statement is much more complex than the previous ones. In this statement, the in_true_raster_or_constant parameter for the first Con (the outer Con) is another Con (the inner Con, specifying southerly aspect) with its own condition and true and false statements. In the outer Con, the first condition determines if a location is less than 10 percent slope but also meets the other criteria identified above (being within the elevation range, open, within a mile from a road, and within a half a mile of a stream). If the location meets these criteria, it is then evaluated by a second condition, the inner Con, to see if the location faces south. If it does, then the location is assigned the elevation value (the in_true_raster_or_constant of the inner Con), if not, instead of assigning it NoData, it is assigned a \u201c0\u201d (the in_false_raster_or_constant of the inner Con) so that you can later identify these locations. Since the in_false_raster_or_constant for the outer Con is blank, any location that is steep and\/or does not meet the series of criteria defined above (the elevation, land use type, and distances criteria) will be assigned NoData. \u00a0Once you have mastered matching the conditions and the true and false statements within this complex Con statement, you will be able to unleash the power of Con to address many of your specific needs.<\/p>\n At this point, you know the locations that meet all your criteria. However, while selecting the area to expand your growing operation and to maintain a contiguous tract of land, you may need to relax one or more of the criteria (such as southerly aspect), and consider areas that are now assigned \u201c0\u201d, for your purchase.<\/p>\n Now get ready to plant!<\/p>\n"}],"authors":[{"ID":7161,"user_firstname":"Kevin","user_lastname":"Johnston","nickname":"Kevin Johnston","user_nicename":"kevinjohnston","display_name":"Kevin Johnston","user_email":"kjohnston@esri.com","user_url":"","user_registered":"2018-03-02 00:19:23","user_description":"Kevin has been a Product Engineer on the Spatial Analyst development team for over 29 years. He has a Master's in Landscape Architect from Harvard and a Ph.D. in environmental modeling from Yale. At Esri, Kevin\u2019s current focus is developing suitability and connectivity tools. He hopes the tools that he works on can help users make more informed decisions.","user_avatar":"import arcpy <\/code> from arcpy.sa import *<\/code> elev = Raster(\"elevation.tif\") # Creating a raster object from the elevation dataset<\/code> outElev1wc = Con(elev, 1, 0, \"Value > 1800\")<\/code><\/pre>\noutElev1 = Con(elev > 1800, 1, 0)<\/code><\/pre>\noutElevRange = Con(((elev > 1800) & (elev < 3600)), elev)<\/code><\/pre>\ninlanduse = Raster('landuse') # Creating a raster object from the landuse dataset<\/code> outElevLU = Con(inlanduse, outElevRange,'', \"Category = 'OpenArea'\")<\/code><\/pre>\ninDistRoads = EucDistance('roads')<\/code> inDistStreams = EucDistance('streams')<\/code> outElevDist = Con(((inDistRoads < 5280) & (inDistStreams < 2640)), outElevLU)<\/code><\/pre>\ninSlope = Slope(elev)<\/code> inAspect = Aspect(elev)<\/code> outRaster = Con((inSlope < 10) & (outElevDist > 0), Con((inAspect > 135) & (inAspect < 225), outElevDist, 0))<\/code><\/pre>\n
"},{"ID":7091,"user_firstname":"Sarmistha","user_lastname":"Chatterjee","nickname":"Sarmistha Chatterjee","user_nicename":"schatterjee17","display_name":"Sarmistha Chatterjee","user_email":"SChatterjee@esri.com","user_url":"","user_registered":"2018-03-02 00:19:17","user_description":"Sarmistha is a Product Engineer in Esri's raster analysis group, where she focuses on raster and scientific multidimensional data analysis. She is especially interested in combining GIS technology, spatial analysis, and software engineering to study changes for conservation and better decision-making. Before joining Esri in 2017, Sarmistha earned her PhD in Geography from the University of Delaware, where she researched fluvial geomorphology, water resources, and spatial modeling.","user_avatar":"
"}],"related_articles":"","card_image":false,"wide_image":false},"yoast_head":"\n