ArcToolbox banner

Raster 4 - Dasymetric Calculations

A dasymetric map seeks to display statistical surface data by exhaustively partitioning space into zones where the zone boundaries reflect the underlying statistical surface variation. The process of dasymetric mapping is the transformation of data from a set of arbitrary source zones to a dasymetric map via the overlay of the source zones with an ancillary data set. In practice, dasymetric mapping is often considered a particular type of areal interpolation technique where source zone data are excluded from certain classes in a categorical ancillary data set. Dasymetric mapping is applicable to a wide variety of tasks where the user seeks to refine spatially aggregated data, for example in estimating local population characteristics in areas where only coarser, regional resolution census data are available. This tool addresses the design, implementation, validation, and application of a new ‘intelligent’ dasymetric mapping (IDM) technique that supports a variety of methods for characterizing the relationship between the ancillary data and underlying statistical surface. The technique is referred to as intelligent because an analyst may establish this relationship subjectively using their own domain knowledge, extract this relationship from the data using a novel empirical sampling technique, or combine the subjective and empirically-based methods.


Usage Tips

Here's a brief summary of the script process:

The script also produces two summary tables for your reference in the same output workspace as the dasymetric working table. The first, SamplingSummaryTable,contains statistics for the source population units regarded as representative of each ancillary class. The second, FinalSummaryTable, contains statistics for the final desnity values. A data dictionary text file that describes all of the attributes of all of the ouput tables should be found in the toolbox/script directory.


Command line syntax

DasymetricCalculations <Population_Working_Table> <Population_Count_Field> <Population_Area_Field> <Dasymetric_Working_Table> <Population_ID_Field> <Ancillary_Class_Field> <Combined_Area_Field> <Minimum_Sample> <Minimum_Sampling_Area> {Percent} {Preset_Table} <Preset_Field>

Parameters
Expression Explanation
<Population_Working_Table>

The standalone population table created in prior steps.

<Population_Count_Field>

The field in the population feature class that contains population count values.

<Population_Area_Field>

If this script is being used on raster data, this field should be set to "Count."

If this tool is being used with vector data, please make sure there is a field in the population feature class that contains area values that were calculated in this feature class' spatial reference. If the feature class is in a geodatabase, there should already be such a field. If the feature class is a shapefile, you must calculate this field yourself. Please see the ArcGIS help for more information on calculating the area of a shape.

<Dasymetric_Working_Table>

The standalone table, usually created in step two, that contains a row for each unique combination of population unit and ancillary class. This table will be used for the bulk of the dasymetric calculations.

If you re-use the same table more than once, the script will purge the old calculated values and start fresh.

<Population_ID_Field>

The field in the output working table that uniquely identifies the source units. This field name is usually derived from the name of the Population Raster.

<Ancillary_Class_Field>

Field containing ancillary classes to be used for redistributing the population within each source unit. The field may be of any type.

<Combined_Area_Field>

The name of the field in the output working table containing the raster cell count.

<Minimum_Sample>

This is the minimum threshold of representative source units required for an ancillary class to be considered sampled. We found that if only a small number (1-3) were sampled, those units were not often truly representative of others in that ancillary class, and the Intelligent Areal Weighting (IAW) method provided better results. Any class that is not preset or sufficiently sampled will be assigned a density using the IAW method.

<Minimum_Sampling_Area>

The minimum number of raster cells that may be considered "representative" of a source area. Increasing this number may eliminate some very small areas that might be considered anomalous outliers that could skew the class average.

{Percent}

This is the percent threshold value, in decimal notation, if Percent Area is specified above.

{Preset_Table}

If you wish, you may use this feature to manually "preset" the popluation density values of any of your ancillary classes with a priori values. For instance, if your ancillary dataset represents land cover, you may wish to preset the "Water" or "Glacier" classes to zero. To do so, create a standalone table with a field of the same name and type as the Ancillary Class Field above that contains only your preset classes, and add a field that contains their respective preset classes.

<Preset_Field>

The field in the preset table containing the preset population density values.

Scripting syntax

DasymetricCalculations (Population_Working_Table, Population_Count_Field, Population_Area_Field, Dasymetric_Working_Table, Population_ID_Field, Ancillary_Class_Field, Combined_Area_Field, Minimum_Sample, Minimum_Sampling_Area, Percent, Preset_Table, Preset_Field)

Parameters
Expression Explanation
Population Working Table (Required)

The standalone population table created in prior steps.

Population Count Field (Required)

The field in the population feature class that contains population count values.

Population Area Field (Required)

If this script is being used on raster data, this field should be set to "Count."

If this tool is being used with vector data, please make sure there is a field in the population feature class that contains area values that were calculated in this feature class' spatial reference. If the feature class is in a geodatabase, there should already be such a field. If the feature class is a shapefile, you must calculate this field yourself. Please see the ArcGIS help for more information on calculating the area of a shape.

Dasymetric Working Table (Required)

The standalone table, usually created in step two, that contains a row for each unique combination of population unit and ancillary class. This table will be used for the bulk of the dasymetric calculations.

If you re-use the same table more than once, the script will purge the old calculated values and start fresh.

Population ID Field (Required)

The field in the output working table that uniquely identifies the source units. This field name is usually derived from the name of the Population Raster.

Ancillary Class Field (Required)

Field containing ancillary classes to be used for redistributing the population within each source unit. The field may be of any type.

Combined Area Field (Required)

The name of the field in the output working table containing the raster cell count.

Minimum Sample (Required)

This is the minimum threshold of representative source units required for an ancillary class to be considered sampled. We found that if only a small number (1-3) were sampled, those units were not often truly representative of others in that ancillary class, and the Intelligent Areal Weighting (IAW) method provided better results. Any class that is not preset or sufficiently sampled will be assigned a density using the IAW method.

Minimum Sampling Area (Required)

The minimum number of raster cells that may be considered "representative" of a source area. Increasing this number may eliminate some very small areas that might be considered anomalous outliers that could skew the class average.

Percent (Optional)

This is the percent threshold value, in decimal notation, if Percent Area is specified above.

Preset Table (Optional)

If you wish, you may use this feature to manually "preset" the popluation density values of any of your ancillary classes with a priori values. For instance, if your ancillary dataset represents land cover, you may wish to preset the "Water" or "Glacier" classes to zero. To do so, create a standalone table with a field of the same name and type as the Ancillary Class Field above that contains only your preset classes, and add a field that contains their respective preset classes.

Preset Field (Required)

The field in the preset table containing the preset population density values.