Tendance (3D Analyst)
Récapitulatif
Interpole une surface raster à partir de points à l’aide d’une méthode de tendance.
Utilisation
-
La surface interpolée se complexifie au fur et à mesure de l'accroissement de l'ordre polynomial. Un ordre polynomial supérieur ne génère pas nécessairement la surface la plus précise ; cela dépend des données.
-
Le fichier EQM facultatif en sortie contient les informations sur l'erreur EQM (racine carrée moyenne) de l'interpolation. Ces informations permettent de déterminer la meilleure valeur à utiliser pour l'ordre polynomial, en modifiant la valeur d'ordre jusqu'à obtenir l'erreur QM la plus faible. Reportez-vous à la rubrique Fonctionnement de l'outil Tendance pour plus d'informations sur le fichier EQM.
-
Pour l'option LOGISTIC du Type de régression, le champ de valeur z d'entités ponctuelles en entrée doit avoir des codes de zéro (0) et un (1).
-
Certains jeux de données en entrée peuvent avoir plusieurs points avec les mêmes coordonnées x,y. Si les valeurs des points à l'emplacement commun sont les mêmes, elles sont considérées comme étant en double et n'ont aucune incidence sur la sortie. Si les valeurs sont différentes, elles sont considérées comme étant des points coïncidents.
Les différents outils d'interpolation peuvent gérer cette condition de données différemment. Par exemple dans certains cas, le premier point coïncident détecté est utilisé pour le calcul, tandis que dans d'autres cas c'est le dernier point détecté qui est utilisé. Cela peut entraîner des valeurs inattendues pour les emplacements du raster en sortie. La solution consiste à préparer vos données en supprimant ces points coïncidents. L'outil Collect Events de la boîte à outils Outils de statistiques spatiales sert à identifier tous points coïncidents de vos données.
Syntaxe
| Paramètre | Explication | Type de données |
in_point_features |
Les entités points en entrée contenant les valeurs z à interpoler dans un raster de surface. | Feature Layer |
z_field |
Champ contenant une valeur de hauteur ou de magnitude pour chaque point. Il peut s'agir d'un champ numérique ou du champ Forme, si les entités ponctuelles en entrée contiennent des valeurs Z. Si la régression est du type Logistique, les valeurs du champ peuvent uniquement être 0 ou 1. | Field |
out_raster |
Raster de surface interpolé en sortie. | Raster Layer |
cell_size (Facultatif) |
Taille des cellules qui sera utilisée pour la création du raster en sortie. Il s'agit de la valeur dans l'environnement s'il est explicitement défini, ou de la plus petite largeur ou hauteur de l'étendue des entités ponctuelles en entrée, dans la référence spatiale en entrée, divisée par 250. | Analysis Cell Size |
order (Facultatif) |
Ordre de la tendance polynomiale. Cette valeur doit être un entier compris entre 1 et 12. Une valeur de 1 ajuste une surface plane au niveau des points, alors qu'une valeur supérieure ajuste une surface plus complexe. La valeur par défaut est 1. | Long |
regression_type (Facultatif) |
Type de régression à effectuer.
| String |
out_rms_file (Facultatif) |
Nom du fichier texte en sortie dans lequel enregistrer les informations sur l'erreur QM et le khi carré de l'interpolation. Ce fichier doit posséder une extension ".txt". | File |
Exemple de code
Cet exemple entre un fichier de formes ponctuel et interpole la surface en sortie en tant que raster TIFF.
import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.Trend_3d("ca_ozone_pts.shp", "ozone",
"C:/output/trendout.tif", 2000, 2, "LINEAR")
Cet exemple entre un fichier de formes ponctuelles et interpole la surface en sortie en tant que raster Grid.
# Name: Trend_3d_Ex_02.py
# Description: Interpolate a series of point features onto a
# rectangular raster using a trend technique.
# Requirements: 3D Analyst Extension
# Import system modules
import arcpy
from arcpy import env
# Set environment settings
env.workspace = "C:/data"
# Set local variables
inPointFeatures = "ca_ozone_pts.shp"
zField = "ozone"
outRaster = "C:/sapyexamples/output/trendout02"
cellSize = 2000.0
PolynomialOrder = 2
regressionType = "LINEAR"
# Check out the ArcGIS 3D Analyst extension license
arcpy.CheckOutExtension("3D")
# Execute Trend
arcpy.Trend_3d(inPointFeatures, zField, outRaster, cellSize,
PolynomialOrder, regressionType)