nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Extraction de caractéristiques d'Arthur à partir de LiDAR, DEMs et imagerie

ArthurCrawford01

<\/HEAD>

26\/07\/2019 : Ajout de la méthode pour réaliser un inventaire d'arbres en 3D sans lidar ni hauteur (indice DBH).<\/P>

24\/10\/2018 : Ajout d'un lien et des instructions pour créer des bâtiments 3D à partir du lidar de Microsoft Building Footprints US.<\/P>

04\/10\/2018 : Ajout d'une mise à jour des outils dans le 3D mapping with Lidar Point Clouds<\/A> ci-dessous appelée ArthursLidarBuildingExtractionTools_10_04_2018.zip. Cela corrige des problèmes avec l'approche Unclassified pour la collecte de bâtiments à partir du lidar, un outil pour découper les bâtiments extraits selon les parcelles, et l'outil Jagged qui aide à identifier les bâtiments irréguliers résultant de la régularisation. Ces outils sont des améliorations que j'ai testées avec de nombreuses villes\/comtés pour l'extraction de bâtiments. Ces mises à jour seront finalement appliquées au 3D mapping with Lidar Point Clouds<\/A>. Pour créer des bâtiments 3D à partir des empreintes 2D extraites, veuillez utiliser la dernière version de la Local Government 3D Basemaps solution<\/A>.<\/P>

28\/09\/2017 : Veuillez consulter

3D mapping with Lidar Point Clouds<\/A>, ce sont des outils ArcGIS Pro, des données exemples et des exercices d'atelier issus de l'atelier IMF que nous avons développé pour le Forum Imagery and Mapping de cette année avant l'ESRI UC. Cela couvre l'extraction des empreintes de bâtiments à partir du lidar classifié pour les bâtiments, l'extraction des empreintes de bâtiments à partir du lidar non classifié pour les bâtiments, et les cartes de base 3D pour passer des empreintes au 3D et ajouter des arbres. Cela remplace la méthodologie présentée par Geoff Taylor dans Extracting Building Footprints from Classified LiDAR<\/SPAN><\/A> en utilisant certaines de mes techniques plus anciennes. <\/P>05\/12\/2016 : Ajout important aux suggestions de l'outil Regularize Building Footprint.<\/P>08\/12\/2016 : Changement de 1200 en 12000 dans les suggestions de l'outil Regularize Building Footprint.<\/P>21\/12\/2016 : Ajout de Clean Extracted Buildings Using Regularize Building Footprint Tool.<\/P>22\/12\/2016 : Correction d'une erreur dans les modèles pour Clean Extracted Buildings Using Regularize Building Footprint Tool due à un bug dans l'outil Regularize Building Footprint (entrée d'une classe d'entités vide).<\/P>03\/01\/2017 : Problème détecté avec Clean Extracted Buildings Using Regularize Building Footprint Tool ne fonctionnant pas s'il n'y a pas de grands bâtiments (plus de 25 000 pieds carrés), en cours d'investigation.<\/P><\/P><\/P>Ce blog sera dédié à l'extraction de caractéristiques depuis raster ou nuages de points vers vecteur. J'ai généré du contenu pour les World Topo Map Basemaps incluant les courbes de niveau, points cotés, végétation et même les bâtiments issus du Lidar. Dernièrement, je me suis concentré sur l'extraction des bâtiments et arbres à partir du Lidar et\/_ou imagerie. J'ai écrit un outil pour extraire les arbres du Lidar il y a 5 ans et il est publié ici<\/A>.<\/P>

<\/P>

<\/STRONG><\/SPAN><\/P>Instructions pour créer des bâtiments 3D à partir des

Microsoft US Building Footprints<\/A> :<\/STRONG><\/SPAN><\/P>

Microsoft a récemment publié un ensemble gratuit d'empreintes de bâtiments générées par deep learning couvrant les États-Unis d'Amérique. En soutien à ce travail remarquable et afin de rendre ces empreintes disponibles à la communauté ArcGIS, Esri a consolidé les bâtiments en une seule couche et les a partagés dans ArcGIS . <\/SPAN><\/SPAN>

Microsoft US Building Footprints<\/A> n'est pas parfait, mais offre une solution à une communauté sans empreintes de bâtiments et disposant du lidar pour créer une ville 3D. L'ancienneté des empreintes varie. <\/DIV>

1. Utilisez

JSON to Features <\/>\ou Spatial ETL tool<\/> <\/>Import Building GeoJSON<\/><\/> ( ArcGIS Data Interoperability) pour importer les Microsoft US Building Footprints pour votre État. Une autre option serait d'utiliser la classe d'entités Microsoft US Building Footprints sur ArcGIS.com. Vous pouvez l'ajouter à ArcGIS Pro en utilisant l'URL : https://services.arcgis.com/P3ePLMYs2RVChkJx/arcgis/rest/services/MSBFP2/FeatureServer/0, sélectionnez les bâtiments que vous souhaitez utiliser et exportez-les vers une nouvelle classe d'entités. Si vous sélectionnez une grande zone, changez le paramètre In beyond en <None> pour voir les bâtiments à une échelle supérieure à 1:24 000.https://services.arcgis.com/P3ePLMYs2RVChkJx/arcgis/rest/services/MSBFP2/FeatureServer<\/><\/>2. Rendez-vous sur le The National Map TNM Download où, si cela échoue, sur le site United States Interagency Elevation Inventory pour vérifier si vous disposez du lidar sur la zone souhaitée et savoir où le télécharger.<\/><\/><\/> a. The National Map TNM Download cochez Elevation Source Data (3DEP) - Lidar, IfSAR puis Lidar Point Cloud.<\/> b. Naviguez sur la carte jusqu'à la zone que vous souhaitez télécharger depuis le TNM. puis cliquez sur Find Products, ajoutez au panier et téléchargez votre lidar. <\/> c. Si seul le format laz est disponible, utilisez alors l'outil Laszip.exe et ses instructions fournies par l'USGS pour convertir en format LAS.<\/> <\/><\/>3. Ouvrez le fichier las dans ArcGIS et identifiez la projection. Reprojetez une portion découpée des Microsoft Building Footprints dans la zone où se trouvent les fichiers las vers la même projection que ces derniers. Les unités de projection doivent être en pieds ou mètres.<\/><\/><\/>4. Si vous disposez de données Parcelles, vous pouvez utiliser Union avec les Microsoft Buildings puis appliquer une requête définie sur les bâtiments divisés pour sélectionner uniquement là où les bâtiments intersectent (FID Parcel and FID Buildings <> -1) pour exporter vers une nouvelle featureclass.<\/DIV><\/DIV>5. Installez les Local Government 3D Basemaps <\/A>et suivez les tâches pour utiliser les fichiers Microsoft Buildings\/las comme entrées. Assurez-vous d'utiliser le nouvel outil pour la segmentation des toits qui améliore grandement les bâtiments.<\/DIV><\/DIV>

<\/P>

Un exemple de ceci se trouve à :

http:\/\/www.arcgis.com\/home\/item.html?id=3f3de2298e27408689cb71b92340db7f<\/A> <\/P>

<\/P>

Voici un lien pour le voir en 3D :

Scene Viewer<\/A> <\/P>

<\/P>

Extraction semi-automatique des empreintes de bâtiments à partir d'images<\/SPAN><\/P> <\/SPAN>
<\/SPAN><\/P>Ci-joint l'outil pour ArcGIS Pro 1.3 pour Extraire les Bâtiments à partir d'Images. Vous aurez besoin de 3D Analyst et Spatial Analyst pour cet outil. Il prend des images segmentées avec l'outil Segment Mean Shift et converties en vecteurs. Vous sélectionnez ensuite les vecteurs, lancez l'outil. Il utilise un modèle pour convertir en bâtiments corrigés à partir des vecteurs raster bruts et les ajoute à une classe d'entités. Ci-dessous un exemple de vecteurs de segmentation sélectionnés puis traités par l'outil pour obtenir les empreintes des bâtiments. Je suggère de sélectionner plusieurs bâtiments puis d'exécuter l'outil. L'outil est joint en bas sous le nom ExtractBuildingsFromImageryPro1_3v4.zip. Après avoir décompressé les données, assurez-vous d'utiliser définir la projection sur la featureclass des bâtiments dans output.gdb pour correspondre à la projection de vos images.<\/P><\/P><\/P>Il y avait un problème avec les chemins qui ne passaient pas correctement. Je crois que c'est maintenant corrigé car l'espace de travail est défini comme in_memory. Si vous obtenez une erreur indiquant que le chemin est introuvable, ouvrez le modèle et changez les chemins dans le modèle pour C:\Users\student\Documents\ArcGIS\Projects\ImageryExtraction\ImageryExtraction.gdb\seg<\/SPAN> en in_memory\seg. Vous devrez faire cela pour tous les fichiers temporaires dans le modèle. Merci à ceux qui ont signalé cette erreur.<\/P><\/P><\/P>Téléchargez l'outil depuis mon blog.<\/LI>Décompressez l'outil.<\/LI>Obtenez des images haute résolution de la zone dont vous souhaitez extraire les bâtiments (6 pouces est idéal).<\/LI>Ouvrez ArcGIS Pro (1.3 ou supérieur) Carte<\/LI>Ajoutez vos images à la Carte<\/LI>Sous l'onglet Insertion, utilisez Ajouter Boîte à outils pour ajouter la boîte à outils ExtractBuildingsFromImageryPro1_3.tbx.<\/LI>Exécutez l'outil Imagery Segmentation to Vector sur vos images. Cela prend du temps (5 minutes pour un fichier de 349mg). Si quatre bandes, entrez 4 2 1 dans Band Indexes.<\/LI>Changez la symbologie de la sortie en rouge et placez-la au-dessus de l'image segment mean shift. Cela devrait ressembler à ceci :<\/LI>Désactivez l'image de segmentation et avec l'imagerie originale en dessous, sélectionnez les polygones des maisons que vous souhaitez convertir en polygones de bâtiment.
Pour la première fois, exécutez Définir Projection sur la featureclass des bâtiments pour correspondre à votre imagerie.
Ouvrez l'outil Segmentation Vectors to feature et exécutez-le après avoir défini vos paramètres.
Inspectez vos résultats :
< /OL >Outil Regularize Building Footprint - Nettoyez les entités brutes extraites du lidar ou des images.

<\/P>

<\/P>1. Pour utiliser cet outil, commencez par exécuter un Dissolve sur vos entités avec multipart désactivé.<\/SPAN><\/P>

<\/P><\/P>2. Téléchargez et décompressez le fichier joint Clean Extracted Buildings Using Regularize Building Footprint Tool.zip<\/P>Exécutez d'abord l'outil "1. Regularize Building Footprints For Circle (Feet Projection)".<\/P><\/P>3. Exécutez ensuite l'outil "2. Regularize Building Footprints All Non-Circle (Feet Projection)".<\/P><\/P>4. Examinez d'abord votre sortie CircleBuildings, il ne devrait pas y en avoir beaucoup et vous devez vérifier s'ils sont réellement<\/P>des cercles. & nbsp ; Si ce n'est pas le cas, supprimez-les.& nbsp ; Cela devrait collecter des réservoirs, des bâtiments ronds, des kiosques, etc.& nbsp ; Si un bâtiment n'est pas complètement rond,<\/P>mais ressemble à un hexagone ou quelque chose de similaire, essayez d'utiliser le bouton généraliser dans l'édition avec un paramètre de 1,2 ou3 pour obtenir des côtés<\/P > droits comme un hexagone ou un octogone.& nbsp ; Ajoutez d'autres structures circulaires que vous voyez, souvent elles sont proches les unes des autres.& nbsp ; Enregistrez lorsque terminé.5. Ajoutez votre classe d'entités Buildings et sélectionnez par emplacement celles dont le centroïde se trouve dans les CircleBuildlings, puis supprimez-les.6. Examinez vos bâtiments.& nbsp ; Si vous avez des problèmes généraux d'extraction, veuillez examiner le type de bâtiment.& nbsp ; Allez dans l'outil "2.Regularize Building Footprints All Non-Circle (Feet Projection)" et ajustez ces paramètres en utilisant l'aide de l'outil Regularize Building Footprint.Des données d'exemple sont incluses ainsi que la sortie avant examen manuel.& nbsp & nbsp & nbsp Vous pouvez voir certains escaliers dans lesentités apparaître.& nbsp & nbsp Ceci est dû au fait que l'outil est uniquement configuré pour gérer les angles droits et diagonaux pour les types de grands et moyens bâtiments (Voir champ appelé Type).& nbsp & nbsp Le champ TotalProcess dans la classe d'entités Buildings vous montre les processus utilisés.& nbsp & nbsp Exemple:& nbsp & quot ;Regularized, Simplified, Simplified, Simplified, Simplified&quot & nbsp signifie qu'il a eu un outil Regularized Building Footprint avec une sortie Regularized, puis il a été simplifié en utilisant un outil Simplify Buildings avec des paramètres de simplification à2,4,6 et8.*& nbsp Problème détecté avec Clean Extracted Buildings Using Regularize Building Footprint Tool ne fonctionnant pas s'il n'y a pas de grands bâtiments (plus de25,000 pieds carrés), recherche d'une solution.Processus manuel décrit pour vous aider si vous souhaitez mieux comprendre ce que fait l'outil :Commencez par extraire les cercles de vos données.& nbsp & nbsp Je traite généralement un comté ou deux à la fois avec cela.Exécutez Circles d'abord sept fois – Regularize Building Footprint avec Tolérance2,4,5,6,8 et Tolérance12 et les requêtes de définition ci-dessous: minor-fareast; mso-font-kerning: 12.0pt; language: en-US; mso-style-textfill-type: solid; mso-style-textfill-fill-color: white; mso-style-textfill-fill-alpha: 100.0%;">Tolérance 1.6 avec les vecteurs bruts utilisant une requête de définition de Shape_Area <= 220 ET Shape_Area > 175

Tolérance 1.8 avec les vecteurs bruts utilisant une requête de définition de Shape_Area <= 240 ET Shape_Area > 220

Tolérance 2 avec les vecteurs bruts utilisant une requête de définition de Shape_Area <= 1400 ET Shape_Area > 240

Tolérance 4 avec les vecteurs bruts utilisant une requête de définition de Shape_Area <= 2500 ET Shape_Area > 1400

Tolérance&nbsp ;5 avec les vecteurs bruts utilisant une requête de définition de Shape_Area &lt ;=3000&nbsp ; ET Shape_Area &gt ;2500

Tolérance&nbsp ;6 avec les vecteurs bruts utilisant une requête de définition de Shape_Area &lt ;=8000&nbsp ; ET Shape_Area &gt ;3000

Tolérance&nbsp ;8 avec les vecteurs bruts utilisant une requête de définition de Shape_Area &lt ;=12000&nbsp ; ET Shape_Area &gt ;8000

Tolérance 12 avec les vecteurs bruts utilisant une requête de définition de Shape_Area &gt ;12000

Fusionnez la sortie avec BuildingsCircle puis appliquez une requête de définition de STATUS = 0. cela devrait fonctionner pour extraire automatiquement des cercles.

Je suggère d'examiner les cercles pour s'assurer qu'ils sont corrects, en supprimant ceux qui ne le sont pas.& nbsp;& nbsp;

Vous voudrez commencer à éditer sur les cercles fusionnés, sélectionner tout et utiliser la généralisation avec une entrée de 0,1 pied.

Shape_Area<\/SPAN>"< 5000 sq. <\/SPAN>pieds <\/SPAN><\/SPAN><\/P> Right_ANGELS<\/SPAN> Tolérance2, <\/SPAN>Densification2, Précision 0,25<\/SPAN> <\/SPAN><\/SPAN><\/SPAN><\/P><\/P>•<\/SPAN>Moyen :"<\/SPAN>Shape_Area<\/SPAN>" >&#x3D ; 5000 ET <\/SPAN>"<\/SPAN>Shape_Area" <\/SPAN><&#x3D ; <\/SPAN>25000 <\/SPAN>Right_ANGLES_AND_DIAGONAL<\/SPAN> Tolérance3, <\/SPAN>Densification
Précision 0,25, Pénalité Diagonale 1,5<\/SPAN> <\/SPAN><\/SPAN><\/SPAN><\/P><\/P>Grand :<\/SPAN>"<\/SPAN><\/SPAN>Shape_Area<\/SPAN>">= <\/SPAN><\/SPAN>25000 <\/SPAN>Right_ANGLES_AND_DIAGONAL<\/SPAN><\/SPAN>
 Tole9rance 4, Densification 4, <\/SPAN> Pre9cision 0,25, Pe9nalite9 Diagonale 1,5<\/SPAN><\/SPAN><\/P><\/P><\/P><\/SPAN><\/P><\/P>Apr83s avoir ex85cut89 ces groupes, fusionnez le r89sultat avec BuildingsNonCircle et ex85cutez une dissolution sur les entit89s. . Ex85cutez une Eliminate Polygon Part sur les entit89s d'au moins 20 pieds carr89s et Eliminate contained part only. c'est pour se d89barrasser des trous int89rieurs dans les polygones nettoy89s. . Ex85cutez un Simplify Buildings avec une Tol89rance de simplification de 2 puis 85 nouveau avec cela comme entr89e avec une Tol89rance de simplification de 4.<\/SPAN><\/SPAN><\/P><\/P> <\/SPAN><\/P><\/P><\/P>Fusionnez vos b03timents circulaires et les BuildingsNonCirclesb4 pour obtenir vos b03timents finis. . Commencez 85 examiner les b03timents par rapport aux images lidar pour trouver des erreurs (il semble qu'il y aura toujours des erreurs).<\/P><\/P>Extraction des bâtiments et des arbres à partir du Lidar

Voici également un outil pour l'identification des bâtiments et des arbres à partir du Lidar développé par Joe McGlinchy et Roslyn Dunn.&amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;amp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a mp;a amp;s work in progress and only works with ArcGIS Desktop 10.3.1 currently.&;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;;amp;; amp collects building footprints and circle trees that can be run through the Regularize Building Footprint Tool above. Un autre outil que j'ai développé pour extraire les arbres utilise le Lidar et les images NAIP est ici.& nbsp & nbsp & nbsp Une version pour Pro est jointe (Trees From LIDAR and NAIP_Pro).Ajouter de la couleur aux toits des bâtiments pour l'affichage dans ArcGIS Pro (Scene Viewer）

Ajouter de la couleur aux toits dans ArcGIS Pro rend vraiment les scènes plus réalistes avec Local Government Scenes.& nbsp Il y a plus d'un an, j'ai développé un modèle pour ajouter des couleurs hexadécimales aux bâtiments comme champ à afficher dans ArcGIS Pro.& nbsp Avec le processus Local Government Scenes, vous pouvez utiliser un champ avec une couleur hexadécimale pour obtenir des couleurs ajoutées à l'affichage dans ArcGIS Pro en utilisant leurs couches.& nbsp & nbsp & nbsp & nbsp J'utilise souvent deux sources d'imagerie de différentes périodes et utilise l'une pour le toit, puis une autre pour la couleur latérale (ajoutez un autre champ appelé couleur latérale).& nbsp Bien que l'utilisation pour la couleur latérale soit idéale, c'est un moyen facile de colorer les bâtiments en masse. Instructions: 1． Découpez la feuille hors imagerie haute résolution à la zone où vous avez des bâtiments auxquels vous souhaitez ajouter de la couleur. 2． Rééchantillonnez à environ 3 pieds, cela permettra à l'outil Segment Mean Shift de fonctionner beaucoup plus rapidement. 3． Exécutez Segment Mean Shift sur l'imagerie sans feuilles avec les paramètres par défaut. 4． Exécutez Ajouter une couleur hexadécimale à partir de l'imagerie aux toits 10.4 ou ArcGIS Pro Ajouter une couleur hexadécimale aux toits (ci-joint ci-dessous). 5． Téléchargez les anciens fichiers de couche Local Government Scenes (ci-joint ci-dessous compressés) et décompressez-les. 6． Déplacez vos bâtiments vers la scène 3D dans Pro et cliquez sur propriétés pour définir l'altitude sur absolue. 7． Dans Pro, exécutez Appliquer la symbologie à partir de la couche en utilisant le fichier LOD2BuildingShellsFloors_feet.lyrx ou LOD2BuildingShells_meters.lyrx dans le répertoire OldLocalGovernmentScenesLayerFiles que vous avez décompressé.& nbsp La nouvelle version des Local Government Scenes n'a pas de couleur pour les toits ou les côtés des bâtiments, cela sera probablement modifié dans les futures versions.7． Définissez la couleur du toit sur votre champ de couleur hexadécimale dans la symbologie et appliquez. Ci-dessous se trouve St Louis en 3D avec la couleur ajoutée aux bâtiments.& nbsp L'Arche était la seule caractéristique dessinée.& nbsp Toutes les couleurs des bâtiments ont été appliquées via l'imagerie en utilisant le processus ci-dessus.& nbsp Les arbres ont été extraits en utilisant Trees From LIDAR (ici).& nbsp Les bâtiments ont été segmentés en utilisant l'outil Segment Mean Shift sur un DEM reclassifié découpé aux empreintes des bâtiments puis converti en polygones.& nbsp Cela donne aux bâtiments 3D plusieurs hauteurs pour différentes sections.& nbsp Ajouter une couleur hexadécimale à partir de l'imagerie aux toits a été exécuté sur des polygones individuels pour les différentes parties du toit.< IMG __jive_id = "211930" alt = "StLouis.jpg" class = "image-5 jive-image" src = "https://us.v-cdn.net/6038851/uploads/legacyfs/online/211930_StLouis.jpg" style = "width：808px;height：424px；" \/> Voici un autre exemple sur St Charles, County, MO.& nbsp Ils ont récemment rendu leurs données Open Data via MSDIS et le week-end suivant j'ai traité 160000 bâtiments en 3D en utilisant Local Government Scenes avec un modèle autour pour exécuter par lots(Scene Viewer).& nbsp Pour les couleurs latérales du bâtiment, j'ai utilisé cette formule pour calculer une couleur à partir des couleurs des toits. Les couleurs des côtés ne sont pas correctes, mais cela donne une vue plus réaliste.<\/P>

<\/P>

"#" + hex(int(((int(!HexColor![1:3],16))-random.randrange(30)+20)*-1))[-2:] + hex(int(((int(!HexColor![3:5],16))-random.randrange(30))*-1)+20)[-2:]+ hex(int(((int(!HexColor![5:7],16))-random.randrange(30))*-1)+20)[-2:]<\/P>

<\/P>

Une autre option est d'utiliser la vraie couleur des côtés des bâtiments. Pour ce faire, commencez par une recherche d'image sur un navigateur web pour des maisons dans la zone que vous traitez. Vous pouvez ensuite capturer les images, les géoréférencer côte à côte dans ArcMap et les mosaïquer ensemble. Exécutez Segment Mean Shift sur l'image mosaïquée puis créez une classe d'entités pour placer de faux bâtiments sur les zones dont vous souhaitez obtenir les couleurs. Exécutez l'outil Add_Hex_Color_From_RGB_Imagery_to_Roofs_10.4 sur ces bâtiments puis exportez-le vers un tableau. Regardez combien d'enregistrements se trouvent dans le tableau, ajoutez un champ aux bâtiments appelé join et calculez-le avec un nombre aléatoire égal au nombre d'enregistrements dans le tableau (import random dans le bloc de code, random.randint(1,100) , 100 étant le nombre d'enregistrements). Ensuite, ajoutez un champ à votre bâtiment appelé side color et calculez-le après avoir fait une jointure avec le tableau en utilisant objectid et le champ join. Utilisez ce champ comme entrée pour la symbologie des côtés des bâtiments. Un exemple est ci-dessous pour Fort Bend, TX, où nous avons extrait les empreintes des bâtiments à partir du lidar, les avons passés par les Local Government Scenes puis appliqué la couleur aux toits à partir des images et les couleurs des côtés de manière aléatoire.<\/P>

<\/P>

Fort Bend County, TX, en 3D:<\/P>

http:\/\/arcg.is\/2efCU6r<\/A><\/P>

<\/P>

St. Charles en 3D

Scene Viewer<\/A><\/P>

<\/P>

Prendre votre inventaire d'arbres et l'afficher en 3D dans ArcGIS Pro:<\/STRONG><\/SPAN><\/P><\/P>Il y a quelque temps, un parc local m'a demandé quoi faire avec leur nouvel inventaire d'arbres destiné à leur conseil qui a été livré sous forme de classe d'entités ponctuelles avec genre. J'ai suggéré de le mettre en valeur en utilisant ArcGIS Pro en 3D. Ci-dessous se trouvent les instructions pour ajouter une hauteur aux points d'arbres à partir du Lidar puis calculer un champ pour le genre à partir du nom scientifique. En calculant la hauteur et un champ avec juste le genre, vous pouvez afficher facilement les arbres en 3D dans Pro. Voici un lien vers tous les arbres disponibles:<\/P>

Types de genres d'arbres pris en charge014ArcGIS Pro | ArcGIS for Desktop<\/A><\/P>

Avec les arbres, avoir cela basé sur le genre signifie qu'un érable argenté (Acer saccharinum) sera affiché sur le terrain simplement comme "Acer" pour piloter la symbologie. Avec le lien ci-dessus, vous pouvez voir que Acer sera affiché comme un érable à sucre. C'est généralement assez proche pour donner une bonne représentation du type d'arbre.<\/P>

<\/P>

En regardant la recherche, il semble que l'obtention de la hauteur d'un arbre à partir du DBH ne fonctionne pas très bien. Donc, je suggérerais d'utiliser le Lidar ou la hauteur si elle est dans l'inventaire des arbres. Le seul problème que j'ai rencontré était les arbres nouvellement plantés que le lidar n'avait pas ou qui étaient plus petits lorsque le lidar a été réalisé.<\/P>

<\/P>

Examinez les points, vous voudrez probablement reprojeter les points dans la même projection que votre lidar. <\/SPAN><\/LI>La première étape était de télécharger le fichier LAS de la zone
 (format lidar des données). <\/SPAN><\/LI>La deuxième étape était de créer un jeu de données las avec les fichiers las
 .<\/SPAN><\/LI>Ensuite, nous avons utilisé LAS Point as Statistics as Raster en utilisant la méthode Z_Range
 pour obtenir un MNT des hauteurs des arbres.<\/SPAN><\/LI>J'ai exécuté un outil Focal Statistics avec Max de cercle 3 pour
 obtenir la propagation du point le plus haut afin qu'il intersecte le point.<\/SPAN><\/LI>Exécuté Add Surface information pour obtenir le Z du
 raster vers les points de l'inventaire des arbres, c'est la hauteur de l'arbre.<\/SPAN><\/LI>Ajoutez un champ appelé Genus et calculez le champ en utilisant
 le nom scientifique. ArcGIS nécessite
 l'entrée du genre seulement, donc j'ai utilisé ce script vb dans 10.4 pour faire cela à partir du nom scientifique en python (ArcGIS Pro) !BOTANICAL!.split(' '<\/> ,1)[0] ou <\/>vb script Split([BOTANICAL] , " ",2)(0) <\/>
Ouvrez ArcGIS Pro Scene et ajoutez les arbres comme

Ajouter des couches prédéfinies<\/>
Arbres réalistes utilisant ce champ Genus, vous obtenez une symbologie plus correcte
pour les arbres et le champ hauteur. Cela ne fonctionne pas avec tous vos
arbres, mais beaucoup sont pris en charge.& nbsp;
& nbsp ;& nbsp ;& nbsp ;& nbsp ;Encore une fois, vous pouvez voir ceux pris en charge à ce lien : Types de genres d'arbres pris en charge014ArcGIS Pro | ArcGIS for Desktop<\/>

Ajoutez des images et vérifiez vos arbres pour vous assurer qu'ils ont l'air corrects.& nbsp ;& nbsp ;Une erreur courante est que la hauteur peut être en mètres ou en pieds, mais j'avais choisi la mauvaise unité.< /LI>Voici quelques données avec lesquelles je l'ai essayé depuis Amherst, MA avec les données sources et un sous-ensemble traité attaché appelé SubsetOfAmherstMATreeInventory2010.zip :& nbsp ;& nbsp ;& nbsp ;& nbsp ;Inventaire des arbres : < A href = "http://www.arcgis.com/home/item.html?id=4eaa4f8ceecf47cc89ea5726ddbcaad7" title = "http://www.arcgis.com/home/item.html?id=4eaa4f8ceecf47cc89ea5726ddbcaad7" rel = "nofollow noopener noreferrer">http://www.arcgis.com/home/item.html?id=4eaa4f8ceecf47cc89ea5726ddbcaad7< / A >& nbsp ;& nbsp ;& nbsp ;& nbsp ;Lidar : < A href = "http://www.arcgis.com/apps/Legend/main/index.html?appid=d6226fcfd1844e9fb45303e5aa41b730" title = "http://www.arcgis.com/apps/Legend/main/index.html?appid=d6226fcfd1844e9fb45303e5aa41b730" rel = "nofollow noopener noreferrer">http://www.arcgis.com/apps/Legend/main/index.html?appid=d6226fcfd1844e9fb45303e5aa41b730< / A >< IMG __jive_id = "212238" class = "jive-image image-8" src = "https://us.v-cdn.net/6038851/uploads/legacyfs/online/212238_pastedImage_0.png" style = "max-height :900px ; max-width :1200px ;" \/ >Un autre à regarder est < A href = "http://www.arcgis.com/home/item.html?id=aa306931082342048396d5501aca1f31" rel = "nofollow noopener noreferrer">Tower Grove Park < / A >à St. Louis, MO.& nbsp ;& nbsp ;& nbsp ;La scène comprend plus de 6 000 arbres rendus.& nbsp ;Certains arbres n'apparaissent pas avec le genre correct et s'affichent en blanc.< IMG __jive_id = "332516" alt = "Tower Grove Park" class = "image-22 jive-image" src = "https://us.v-cdn.net/6038851/uploads/legacyfs/online/332516_TowerGrovePark.JPG" style = "width :845px ; height :544px ;" \/ >< A href = "https://arcg.is/1qOqbf" rel = "nofollow noopener noreferrer">Faire passer votre inventaire d'arbres en 3D sans lidar ?< / A >& nbsp;(indice DBH)J'ai reçu certaines données d'inventaire d'arbres du cimetière Bellefontaine ici à St. Louis, MO.& nbsp ; & nbsp ;J'ai utilisé une jointure avec des arbres similaires au nom scientifique du genre Le Missouri Botanical Garden a réalisé pour leurs arbres un rendu en 3D. Esri prend en charge environ 60 genres d'arbres pour le rendu 3D. Ensuite, j'ai utilisé le lidar pour obtenir une formule simple afin de convertir le DBH (Diamètre à hauteur de poitrine) en hauteur d'arbre. Ce n'est pas précis, mais cela a rendu les arbres assez bien. Il existe des formules très complexes que je pourrais étudier davantage pour lier le DBH à la hauteur, ce qui pourrait permettre à de nombreux inventaires d'arbres de devenir facilement 3D. Le tableau pour calculer un genre similaire est ci-dessous avec les autres téléchargements appelés MissouriBotanicalGardenSubstitutesForTreeGenus.txt.<\/P>

<\/P>

Scène montrant des arbres rendus en 3D en utilisant une formule appliquée au DBH ((DBH x 2)+20). 5 150 des 9 153 genres d'arbres ont été pris en charge en utilisant les 63 genres <\/SPAN>types pris en charge par Esri ArcGIS Pro<\/A>. 2 120 arbres ne correspondaient pas au genre pris en charge pour le rendu, mais un tableau fourni par le Missouri Botanical Garden a été utilisé pour des genres d'arbres similaires. Les arbres en blanc n'étaient pas pris en charge ni cartographiés pour la symbologie du genre, cela incluait 1 883 arbres et arbustes. Pour calculer le genre, j'ai utilisé le champ SPP (Nom scientifique) avec <\/SPAN> !SPP!.split(<\/SPAN>' '<\/SPAN>,1)[0]<\/SPAN><\/P>Scene Viewer<\/A> - Assurez-vous d'activer le lidar pour voir à quel point les arbres sont proches ou non en hauteur.<\/P>

<\/A><\/P>

Arthur Crawford - Esri<\/P>

<\/P>

Remerciements spéciaux au comté de St. Charles, Monica Balestreri, MSDIS, Michael Owen, ville d'Amherst (MA), Dan Hedges, Joe McGlinchy, Roslyn Dunn, Thomas Maurer, Hua Wei, Jeff Liedtke, J. D. Overton, Khalid Duri, David Altergott, Tower Grove Park, Andy Watson et aux autres qui m'ont aidé.<\/P><\/BODY><\/HTML>

Trouver plus de publications taguées avec

Réponses

RobertFord1

Hi Arthur!

Do you have a source for McGlinchy and Dunn?

Thanks!

ArthurCrawford01

Hi Robert,

By source, do you mean if they both work here at Esri? If you have any questions, I contact them to get a response.

Thanks,

Arthur

RobertFord1

Thanks Arthur!

ArthurCrawford01

Questions are welcome on these tools or processes, we are trying to improve these. If you run into issues, we would like to know.

PhilLarkin1

Arthur-

Thanks for the informative post. What zipfile contains the Regularize Building Footprint Tool? If you have not published it yet, when do you expect to do so?

ArthurCrawford01

Hi Phil, the tool is a model that goes around the Regularize Building Footprint tool that is currently in ArcGIS 10.4 or Pro 1.3 or 1.4. I offered my best advice above on how to use it by dividing buildings into different sizes. I plan on adding a the model built around it, but am still doing testing with versions of it to get the best outcome.

PhilLarkin1

Ah, I see. I thought you might have published an improved version. Thanks again for this post.

I ran the Version 3 of McGlinchy/Dunn and then ran your suggestions through the Regularize Building Footprint tool without first running the Circle method. Here are my results: http://arcg.is/2c4B75R

ArthurCrawford01

Hi Phil,

The building footprints and trees would not show up in the link provided. Could you send me the results (zipped) or share your layers/build a new map?

Arthur Crawford

acrawford@esri.com

enricobonansea

Really great Work Arthur!

running the "segmentation vector to feature" to extract buildings the tool give us always the same error regarding a schema problem

ERROR 000466: C:\Users\enribona\Documents\ArcGIS\Projects\segmentation\segmentation.gdb\Buildingstemp_EliminatePolyg0 does not match the schema of target c:\users\enribona\documents\arcgis\projects\segmentation\Segmentation.gdb\Buildings
Failed to execute (Append).

To run the tool I first created a Building feature class in the segmentatin.gdb with tid, shape, gridcode, shapearea, shapelenght

Thank you very much

Enrico

ArthurCrawford01

Hello Enrico,

I replaced the ExtractBuildingsFromImageryPro1_3.zip with a new one. It has sample data, a geodatabase with sample output and an empty building featureclass to use for your output. I also am adding instruction in a work document. Thanks to Enrico Bonsansea and Leilei Duan for finding this issue. Please let me know if you run into any other problems.

Arthur Crawford

acrawford@esri.com

enricobonansea

Hi Arthur and thank you for the update. I followed the tuorial using my image but I received a different error while i run the "segmentation o vector " tool.

The error is:

Start Time: martedì 30 agosto 2016 17:14:15
Executing (Copy Features): CopyFeatures seg041053sub C:\Users\student\Documents\ArcGIS\Projects\ImageryExtraction\ImageryExtraction.gdb\seg # # # #
Start Time: martedì 30 agosto 2016 17:14:16
ERROR 000210: Cannot create output C:\Users\student\Documents\ArcGIS\Projects\ImageryExtraction\ImageryExtraction.gdb\seg
Failed to execute (Copy Features).

but i havent a "student" user on my PC

Thank you for your Help!!

enrico

ArthurCrawford01

Hello Enrico,

I think when I fixed the tool it change the paths from in_memory to the path for the project. I think it is fixed. Let me know if this works for you now. Being in_memory as storage speeds up the tool.

Thanks Enrico,

Arthur

enricobonansea

Hello Arthur, I got this new error maybe due to the same path problem

thank you

Enrico

Executing (Copy Features): CopyFeatures Seg037_062sub C:\Users\student\Documents\ArcGIS\Projects\ImageryExtraction\ImageryExtraction.gdb\seg # # # #
Start Time: mercoledì 31 agosto 2016 11:48:33
ERROR 000210: Cannot create output C:\Users\student\Documents\ArcGIS\Projects\ImageryExtraction\ImageryExtraction.gdb\seg
Failed to execute (Copy Features).
Failed at mercoledì 31 agosto 2016 11:48:33 (Elapsed Time: 0,41 seconds)
Failed to execute (Model1).
Failed at mercoledì 31 agosto 2016 11:48:34 (Elapsed Time: 1,76 seconds)

ArthurCrawford01

Hello Enrico,

I added a new version of the tool with the workspace as a variable. I believe this should fix the problem of in_memory being substituted with C:\Users\student\Documents\ArcGIS\Projects\ImageryExtraction\ImageryExtraction.gdb. I tested in on another system and it appears to work. Let me know if this works for you now.

Thanks again for you testing of the tool, it's still in prototype.

Arthur

enricobonansea

Hello Arthur , thank you.
I just tried the new version and the path error disappear, the "segmentation to feature" process finish ok.

However at the end no building layer is added in the contents.
I added the building feature class manually but it dosn't show me any polygon.
Opening the attribute table I see that new records are added but they seems empty (lenght , area are equal to zero).

I think that the problem coudl be due to difference in spatial reference system between image and feature class in output.gdb

My images infact are in WGS84(utm) while the geodatabase is in NAD.

I tried to projected the Seg037_062sub and Building feature classes in geographic (4326), then i deleted oll the records in the new feature classes

With these feature classes I re-executed the tools and now i have buildings!

I'm not sure about the process but it works...

Thank You
Enrico

ArthurCrawford01

Hello Enrico,

I'm going to add that you need to use define projection on the output building featureclass in the output.gdb.

Thanks,

Arthur

JasonGillow1

What extansions does this tool require? I go to add the tool and I get a message "Tool is not Licensed".

What to do?

ArthurCrawford01

Hello Jason,

I believe you need Spatial Analyst (for Segment Mean Shift) and 3D Analyst (for Regularize Building Footprint).

Thanks,

Arthur Crawford - Esri

JasonGillow1

So to at least create building footprints you need both extensions then? I just have spatial analyst - would there be a workaround to get the footprints at least?

Thanks

ArthurCrawford01

Hello Jason,

For using this process that is semi-automatic feature extraction you do need both. If you have just spatial analyst, you can get rough polygons by bypassing the Regularize Building Footprints in the model. The polygons will not have the right angles and will have rough edges. You might want to ask your account manager at Esri for a temporary license to test it to see if it is worth it in saving you time on extraction.

I know St. Charles County (Missouri) that has been using this tool increased the number of licenses so more could use this process on extracting building footprints.

To make sure we are talking about the same tool, is ExtractBuildingsFromImageryPro1_3v4.zip the one you are trying to use?

Thanks,

Arthur Crawford - Esri

JasonGillow1

Yes that is the tool.

JoanBiediger

This looks like a great tool. Would it be possible to create a similar process to extract clouds on high resolution 4 band imagery?

MarieDucharmeFEC_0-1761075037627.png

ArthurCrawford01

Hi Joan,

I would likely use the method in Understanding segmentation and classification. Then you could convert to polygon. Using this tool is good for building extraction as is.

GigiB910

Hi Arthur-

I have both licenses installed (3D Analyst is a trial), but I'm still getting the "Tool is not licensed" message when trying to run "Segmentation Vectors to feature". Any idea what the issue may be? The license is active under the Project tab. Thanks!

MeherChelbi1

Hi Arthur,

I would like to test your method to extract building from imagery but I'm operating on ArcMap 10.3 so I can't run the tool as it is.Can you please tell us what are the parameters that you use in the image segmentation: Do you use the Segment mean shift tool and if so what are the parameters you are using (Spectral detail, Spatial detail and Minimum segment size in pixel)

Thank you,

ArthurCrawford01

Hello Meher,

If you are current on your license, you should be able to use ArcGIS Pro to do this, it installs separately. If you use ArcMap 10.3, the settings I used are the defaults which are Spectral detail 15.5, Spatial detail 15 and Minimum segment size in pixel 20. Originally, I did this tool in 10.3. There were issues with the way it interfaced with the client that caused me to use Pro.

Thanks,

Arthur Crawford - Esri

ArthurCrawford01

Hello Gigi,

I would call the Esri help desk, it might be a issue with not having Advanced ArcGIS Pro.

Thanks,

Arthur Crawford - Esri

StephenMorgan1

Hi Arthur, this looks fantastic!

Are the Lidar tools available for Pro 1.3? I downloaded the "3DMappingWorkshop_IMF2017.zip" and tried using them in Pro 1.3, but the tools fail (cannot open/read script/script not found) and I assume this is because of versioning issues (Pro and Python). Are the extract buildings from Lidar tools in lower versions of python/pro? I would love it use these! Unfortunately I cannot upgrade to 2.0 because of work.

ArthurCrawford01

Hello Stephen,

Sorry, but the tools were developed for 2.0. I'm not sure how much work it would take to get it at 1.3. You might try the ArcGIS for Personal Use Program to try it at home and then see if the results can convince your management to upgrade. You will require Spatial Analyst and 3D Analyst extensions.

Arthur Crawford - Esri Content

AkhilChhibber2

I used your model i.e. the first tool "Segmentation to vector" and I am facing the error as shown in the screenshot. This is the error which I am facing while using the tool "Segmentation to vector"

How should I fix this error and how should I proceed

My tif image is in Projected coordinate system i.e. "WGS_1984_UTM_Zone_44N"

Please let me know regarding any possible solution for this error

Rejoindre la conversation

Les membres connectés peuvent publier, suivre les mises à jour, et plus encore. Nouveau ici ? Inscrivez-vous gratuitement.

Liens rapides

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.