Available with Geostatistical Analyst license.

## Summary

Exports a 3D geostatistical layer created using the Empirical Bayesian Kriging 3D tool to a multidimensional Cloud Raster Format (*.crf file) raster dataset. Tools in the Multidimensional Analysis toolset of the Image Analyst toolbox are designed to work directly on multidimensional rasters and can identify the 3D nature of the data.

3D geostatistical layers store continuous 3D interpolation results and appear as 2D horizontal slices at a given elevation. The current elevation can be changed with the range slider, and the layer will calculate and render the predicted values of the new elevation. Rasters of interpolated predictions can be extracted at any elevation using the GA Layer To Rasters tool. The GA Layer 3D To Multidimensional Raster tool automates the process of extracting rasters at multiple elevations and stores them as a multidimensional raster dataset.

## Usage

The Input 3D geostatistical layer must be created using Empirical Bayesian Kriging 3D.

You can provide the desired elevations in either of the following ways:

- Provide elevation values using an iterator—Leave the Enter explicit elevation values parameter unchecked, and provide values to the Minimum elevation, Maximum elevation, and Elevation interval parameters. The iteration will start at the minimum and increase by the interval up to the maximum. For example, if you provide a minimum elevation of 0, a maximum of 100, and an interval of 20, the output will have six elevations: 0, 20, 40, 60, 80, and 100.
- Provide explicit elevation values—Check the Enter explicit elevation values parameter, and provide the desired elevations in the Elevation values parameter.

If the Elevation interval provided does not evenly divide the elevation range (Maximum elevation minus Minimum elevation), one elevation will be used that exceeds the maximum elevation. For example, if you provide a minimum elevation of 10, a maximum elevation of 80, and an elevation interval of 20, the output will have five elevations: 10, 30, 50, 70, and 90.

The default raster cell size of the Output multidimensional raster dataset parameter is the smaller of the output extent of the x and y dimensions, divided by 100. The default number of elevation slices is 10. This will create a multidimensional raster with an x and y resolution that is at least 10 times larger than the z-resolution. Consideration should be given to these values, particularly if you want approximately the same resolution in all three dimensions. Increasing the cell size or increasing the number of elevation slices will result in outputs that are closer to the same resolution in all dimensions.

The execution time of the tool is approximately proportional to the number of overall predictions that are required, and each raster cell of each elevation requires the calculation of a predicted value. By default, the execution time of the tool is approximately 60 percent longer than exporting a single elevation slice using the default options of the GA Layer To Rasters tool. Including Additional output types will not significantly increase the execution time because multiple output types can be predicted simultaneously.

If any Additional output types are provided, the output raster will be a multivariate multidimensional raster dataset with a different variable for each output type.

All tools in the Multidimensional Analysis toolset of the Image Analyst toolbox will have increased performance if transposes are built on the multidimensional raster. You can use the Build multidimensional transpose parameter to build them automatically, but this will increase the execution time of the tool. You can also use the Build Multidimensional Transpose tool to build transposes later.

## Syntax

GALayer3DToMultidimensionalRaster(in_3d_geostat_layer, out_multidimensional_raster, {cell_size}, {explicit_only}, {min_elev}, {max_elev}, {elev_interval}, {elev_values}, {elev_units}, {output_type}, {quantile_probability_value}, {additional_outputs}, {build_transpose})

Parameter | Explanation | Data Type |

in_3d_geostat_layer | The 3D geostatistical layer representing the model to be exported to a multivariate raster dataset. | Geostatistical Layer |

out_multidimensional_raster | The output raster dataset containing the results of exporting the geostatistical model. The output must be saved as a Cloud Raster Format file (*.crf). | Raster Dataset |

cell_size (Optional) | The cell size of the output multidimensional raster. | Analysis Cell Size |

explicit_only (Optional) | Specifies whether elevations will be provided as an explicit list or an iterator will be used. - EXPLICIT_VALUES —Elevation values will be provided as a list.
- NO_EXPLICIT_VALUES —Elevation values will be provided using an iterator. This is the default.
| Boolean |

min_elev (Optional) | The minimum elevation that will be used to start the iteration. | Double |

max_elev (Optional) | The maximum elevation that will be used to stop the iteration. | Double |

elev_interval (Optional) | The increment that the elevation will increase with each iteration. | Double |

elev_values [elev_values,...] (Optional) | The elevation values to export. | Double |

elev_units (Optional) | Specifies the measurement unit of the elevation values. - INCH —Elevation values are in inches.
- FOOT —Elevation values are in feet.
- YARD —Elevation values are in yards.
- MILE_US —Elevation values are in U.S. miles.
- NAUTICAL_MILE —Elevation values are in nautical miles.
- MILLIMETER —Elevation values are in millimeters.
- CENTIMETER —Elevation values are in centimeters.
- DECIMETER —Elevation values are in decimeters.
- METER —Elevation values are in meters.
- KILOMETER —Elevation values are in kilometers.
| String |

output_type (Optional) | Specifies the primary output type of the output multidimensional raster. The Additional output types parameter can be used to specify additional variables in the output multidimensional raster. For more information, see What output surface types can the interpolation models generate? - PREDICTION —A multidimensional raster of predicted values. This is the default.
- PREDICTION_STANDARD_ERROR —A multidimensional raster of standard errors of prediction.
- PROBABILITY —A multidimensional raster predicting the probability that a threshold is exceeded.
- QUANTILE —A multidimensional raster predicting the quantile of the predicted value.
| String |

quantile_probability_value (Optional) | If Output type is set to Quantile, use this parameter to enter the requested quantile. If Output type is set to Probability, use this parameter to enter the requested threshold value, and the probability that the threshold is exceeded will be calculated. Subtract this value from one to get the probability that the threshold is not exceeded. | Double |

additional_outputs [[output_type, quantile_probability_value],...] (Optional) | Specifies the output type and quantile or probability value for each additional output type. If multiple output types are provided, the output raster will be a multivariate raster dataset with a different variable for each output type. For more information, see What output surface types can the interpolation models generate? | Value Table |

build_transpose (Optional) | Specifies whether multidimensional transposes will be built on the output multidimensional raster. - BUILD_TRANSPOSE —Multidimensional transposes will be built on the output multidimensional raster.
- DO_NOT_BUILD_TRANSPOSE —Multidimensional transposes will not be built on the output multidimensional raster. This is the default.
| Boolean |

## Code sample

Interpolate 3D points and convert the output to a multivariate raster.

```
import arcpy
arcpy.ga.EmpiricalBayesianKriging3D("my3DLayer", "Shape.Z", "myValueField", "my3DGALayer")
arcpy.ga.GALayer3DToMultivariateRaster("my3DGALayer", r"C:\\gapydata\\outputMDRD.crf", 1000,
"NO_EXPLICIT_VALUES", 0, 500, 50, "","METERS",
"PREDICTION", "", [["PREDICTION_STANDARD_ERROR",""]],
"BUILD_TRANSPOSE")
```

Interpolate 3D points and convert the output to a multivariate raster with multiple output types and transposes.

```
# Name: GALayer3DToMDR_Example_02.py
# Description: Interpolates 3D points and exports to a multidimensional raster dataset.
# Requirements: Geostatistical Analyst Extension
# Author: Esri
# Import system modules
import arcpy
# Set local variables
in3DPoints = "C:/gapyexamples/input/my3DPoints.shp"
elevationField = "Shape.Z"
valueField = "myValueField"
outGALayer = "myGALayer"
# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")
# Execute Empirical Bayesian Kriging 3D
arcpy.ga.EmpiricalBayesianKriging3D(in3DPoints, elevationField, valueField, outGALayer)
# Export predictions and standard errors to multidimensional raster with tranposes.
# Set up variables
outMDRaster = r"C:\\gapydata\\outputMDRD.crf"
explicitVals = "NO_EXPLICIT_VALUES"
cell_size = 1000
min_elev = 0
max_elev = 500
elev_interval = 50
elev_list = ""
elev_units = "METERS"
out_type = "PREDICTION"
quan_value = ""
add_outputs = [["PREDICTION_STANDARD_ERROR",""]]
transpose = "BUILD_TRANSPOSE"
# Additionally output prediction standard errors.
arcpy.ga.GALayer3DToMultivariateRaster(outGALayer, outMDRaster, cell_size, explicitVals,
min_elev, max_elev, elev_interval, elev_list,
elev_units, out_type, quan_value, add_outputs,
transpose)
```

## Environments

## Licensing information

- Basic: Requires Geostatistical Analyst
- Standard: Requires Geostatistical Analyst
- Advanced: Requires Geostatistical Analyst