Points d'observation (ArcGIS Spatial Analyst)
Récapitulatif
Identifie les points d'observation visibles depuis chaque emplacement de la surface raster.
Utilisation
-
La détermination de points d'observation est un processus nécessitant d'importantes ressources de traitement par ordinateur. Le temps de traitement dépend de la résolution. Pour les études préliminaires, vous pouvez utiliser une taille de cellule plus grossière afin de réduire le nombre de cellules dans l'entrée. Utilisez le raster de résolution maximale lorsque les résultats finaux sont prêts à être générés.
-
Si le raster en entrée comporte du bruit indésirable causé par des erreurs d'échantillonnage, vous pouvez lisser le raster avec un filtre passe-bas, comme l'option Moyenne de Statistiques focales, avant d'exécuter cet outil.
-
La visibilité de chaque centre de cellule est déterminée par la comparaison de l'angle d'altitude vis-à-vis du centre de cellule avec l'angle d'altitude vis-à-vis de l'horizon local. L'horizon local est calculé en fonction du MNT intermédiaire entre le point d'observation et le centre de cellule actuel. Si le point se trouve au-dessus de l'horizon local, il est considéré comme visible.
Un raster en sortie facultatif AGL (Above Ground Level, au-dessus du niveau du sol) est fourni par l'outil. Chaque cellule du raster en sortie AGL enregistre la hauteur minimale devant être ajoutée à cette cellule pour la rendre visible par au moins un observateur.
Lorsque les entités d'observation en entrée contiennent plusieurs observateurs, la valeur en sortie est le minimum des valeurs AGL provenant de chaque observateur.
Si le raster en entrée doit être ré-échantillonné, la technique bilinéaire sera utilisée. Un raster en entrée peut être ré-échantillonné par exemple lorsque le système de coordonnées en sortie, l'étendue ou la taille de cellule est différent(e) de celui ou celle en entrée.
Syntaxe
| Paramètre | Explication | Type de données |
in_raster |
Raster de surface en entrée. | Raster Layer |
in_observer_point_features |
Classe d'entités points qui identifie les emplacements des points d'observation. Le nombre maximal de points autorisés est 16. | Feature Layer |
z_factor (Facultatif) | Nombre d'unités x, y au sol dans une unité z de surface. Le facteur z ajuste les unités de mesure des unités z lorsqu'elles sont différentes des unités x,y de la surface en entrée. Les valeurs z de la surface en entrée sont multipliées par le facteur z lors du calcul de la surface finale en sortie. Si les unités x,y et les unités z utilisent les mêmes unités de mesure, le facteur z est égal à 1. Si les valeurs x,y et les valeurs z sont exprimées dans des unités de mesure différentes, le facteur z doit être défini comme approprié, sinon les résultats sont incorrects. Par exemple, si les unités z sont des pieds et les unités x,y sont des mètres, vous devez utiliser un facteur z égal à 0,3048 pour convertir les unités z de pieds en mètres (1 pied = 0,3048 mètre). | Double |
curvature_correction (Facultatif) |
Permet de corriger en intégrant la courbure de la Terre.
| Boolean |
refractivity_coefficient (Facultatif) |
Coefficient de réfraction de la lumière visible dans l'air. La valeur par défaut est 0,13. | Double |
out_agl_raster (Facultatif) |
Raster AGL en sortie. Le raster AGL obtenu est un raster dans lequel la valeur de chaque cellule correspond à la hauteur minimale devant être ajoutée à une cellule (qui, sinon, ne serait pas visible) pour la rendre visible par au moins un observateur. Les cellules qui étaient déjà visibles auront la valeur 0 dans ce raster en sortie. | Raster |
Valeur renvoyée
| Nom | Explication | Type de données |
| out_raster |
Raster en sortie. La sortie identifie avec précision les points d'observation visibles depuis chaque emplacement de la surface raster. | Raster |
Exemple de code
Cet exemple identifie avec précision les points d'observation visibles depuis chaque emplacement de la surface raster.
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outObsPoints = ObserverPoints("elevation","observers.shp", 1, "CURVED_EARTH", 0.13)
outObsPoints.save("C:/sapyexamples/output/outobspnt01")
Cet exemple identifie avec précision les points d'observation visibles depuis chaque emplacement de la surface raster.
# Name: ObserverPoints_Ex_02.py
# Description: Identifies exactly which observer points are visible
# from each raster surface location.
# Requirements: Spatial Analyst Extension
# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Set environment settings
env.workspace = "C:/sapyexamples/data"
# Set local variables
inRaster = "elevation"
inObsPoints = "observers.shp"
zFactor = 1
useEarthCurv = "CURVED_EARTH"
refractionVal = 0.13
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute ObserverPoints
outObsPoints = ObserverPoints(inRaster, inObsPoints, zFactor,
useEarthCurv, refractionVal)
# Save the output
outObsPoints.save("C:/sapyexamples/output/outobspnt02")