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; } } }

Extracción de características de Arthur a partir de LiDAR, DEMs e imágenes

ArthurCrawford01

<\/HEAD>

26\/7\/2019: Se añadió cómo hacer un inventario de árboles en 3D sin lidar ni altura (pista DBH).<\/P>

24\/10\/2018: Se añadió enlace e instrucciones para crear edificios 3D usando lidar de Microsoft Building Footprints US.<\/P>

4\/10\/2018: Se añadió una actualización a las herramientas en el 3D mapping with Lidar Point Clouds<\/A> abajo llamada ArthursLidarBuildingExtractionTools_10_04_2018.zip. Esto corrige problemas con el enfoque Unclassified para la colección de edificios a partir de lidar, una herramienta para recortar edificios extraídos a parcelas, y la herramienta Jagged que ayuda a identificar edificios dentados que resultan de la regularización. Estas herramientas son mejoras que he probado con muchas ciudades\/condados para la extracción de edificios. Estas actualizaciones eventualmente se aplicarán al 3D mapping with Lidar Point Clouds<\/A>. Para crear edificios 3D a partir de las huellas 2D extraídas, por favor use la última versión de la Local Government 3D Basemaps solution<\/A>.<\/P>

28\/9\/2017: Por favor vea

3D mapping with Lidar Point Clouds<\/A>, estas son ArcGIS Pro Tools, datos de ejemplo y ejercicios del taller del IMF que desarrollamos para el Imagery and Mapping Forum de este año antes del ESRI UC. Esto cubre Extracción de Huellas de Edificios desde Lidar Clasificado para edificios, Extracción de Huellas de Edificios desde Lidar No Clasificado para Edificios, y los 3D Basemaps para llevar las huellas a 3D y añadir árboles. Esto reemplaza la metodología que Geoff Taylor mostró en Extracting Building Footprints from Classified LiDAR<\/SPAN><\/A> usando algunas de mis técnicas más antiguas. <\/P>5\/12\/2016: Se añadieron extensamente sugerencias para la herramienta Regularize Building Footprint.<\/P>8\/12\/2016: Se cambió 1200 a 12000 en las sugerencias para la herramienta Regularize Building Footprint.<\/P>21\/12\/2016: Se añadió Clean Extracted Buildings Using Regularize Building Footprint Tool.<\/P>22\/12\/2016: Se corrigió error en modelos para Clean Extracted Buildings Using Regularize Building Footprint Tool debido a un bug en la herramienta Regularize Building Footprint (entrada vacía featureclass).<\/P>03\/01\/2017: Se encontró un problema con Clean Extracted Buildings Using Regularize Building Footprint Tool que no funciona si no hay edificios grandes (más de 25,000 pies cuadrados), investigando.<\/P><\/P><\/P>Este blog estará dedicado a la extracción de características desde raster o nubes de puntos a vector. He estado generando contenido para los World Topo Map Basemaps incluyendo Contornos, Elevaciones Puntuales, Vegetación e incluso edificios desde Lidar. Últimamente, me he enfocado en extraer edificios y árboles desde Lidar y/o imágenes. Escribí una herramienta para extraer Árboles desde Lidar hace 5 años y está publicada aquí<\/A>.<\/P>

<\/P>

<\/STRONG><\/SPAN><\/P>Instrucciones para crear Edificios 3D desde

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

Microsoft lanzó recientemente un conjunto gratuito de huellas de edificios generadas por aprendizaje profundo que cubren los Estados Unidos de América. En apoyo a este gran trabajo y para hacer estas huellas disponibles a la comunidad ArcGIS, Esri ha consolidado los edificios en una sola capa y los ha compartido en ArcGIS . <\/SPAN><\/SPAN>

Microsoft US Building Footprints<\/A> no son perfectos, pero ofrecen a una comunidad sin huellas de edificios y con lidar una forma de crear una ciudad 3D. La antigüedad de las huellas varía. <\/DIV>

1. Use

JSON to Features <\/_A>o herramienta Spatial ETL<\/_SPAN> <\/_SPAN>Import Building GeoJSON<\/_A> ( ArcGIS Data Interoperability<\/_A>) para importar las Microsoft US Building Footprints<\/_A> de su Estado. Otra opción sería usar la clase de entidad Microsoft US Building Footprints<\/_A> en ArcGIS.com. Puede agregarla a ArcGIS Pro usando la url: https:\/services.arcgis.com/P3ePLMYs2RVChkJx/arcgis/rest/services/MSBFP2/FeatureServer/0<\/_A>, seleccione los edificios que desea usar y expórtelos a una nueva clase de entidad. Si está seleccionando un área grande, cambie el In beyond a <None> para ver los edificios a una escala mayor que 1:24,000.https:\/services.arcgis.com/P3ePLMYs2RVChkJx/arcgis/rest/services/MSBFP2/FeatureServer<\/_A><\/_DIV>

<\/_DIV>

2. Vaya al

The National Map TNM Download<\/_A> o, si eso falla, al sitio United States Interagency Elevation Inventory<\/_A> a verificar si tiene lidar sobre el área deseada y dónde descargarlo.<\/_DIV>

<\/_DIV>

The National Map TNM Download<\/_A> clic en <\/_SPAN>Elevation Source Data (3DEP) - Lidar, IfSAR y luego en Lidar Point Cloud.<\/_SPAN><\/_DIV> b. Vaya en el mapa al área que desea descargar del TNM. luego haga clic en Find Products, añada al Carrito y descargue su lidar. <\/_SPAN><\/_DIV> c. Si sólo está disponible el formato laz, entonces use <\/_SPAN>herramienta Laszip.exe e instrucciones del USGS<\/_A> para convertir a formato LAS.<\/_SPAN><\/_DIV> <\/_SPAN><\/_DIV>3. Abra el archivo las en ArcGIS e identifique la proyección. &nbsp.Reproyecte una porción recortada de Microsoft Building Footprints en el área donde están los archivos las a la misma proyección que los archivos las. &nbsp.Unidades de proyección deben estar en pies o metros.<\/_DIV><\/_DIV>4. Si tiene datos Parcel, puede usar Union con los Microsoft Buildings y luego poner una consulta de definición sobre los edificios divididos para seleccionar sólo donde los edificios intersectan (FID Parcel y FID Buildings <> -1) para exportar a una nueva featureclass.<\/DIV><\/DIV>5. Instale los Local Government 3D Basemaps <\/A>y siga las tareas para usar los archivos Microsoft Buildings\/las como entradas. Asegúrese de usar la nueva herramienta para Roof Segmentation que mejora mucho los edificios.<\/DIV><\/DIV>

<\/P>

Un ejemplo de esto está en:

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

<\/P>

Aquí hay un enlace para verlo en 3D:

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

<\/P>

Extracción semi-automática de huellas de edificios a partir de imágenes<\/SPAN><\/P> <\/SPAN>
<\/SPAN><\/P>Adjunto está la herramienta para ArcGIS Pro 1.3 para Extraer Edificios De Imágenes. Necesitará 3D Analyst y Spatial Analyst para esta herramienta. Toma imágenes segmentadas con la herramienta Segment Mean Shift y convertidas a vectores. Luego selecciona los vectores, ejecuta la herramienta. Usa un modelo para convertir en edificios corregidos a partir de los vectores ráster sin procesar y anexarlos a una clase de entidad. A continuación hay un ejemplo de vectores de segmentación seleccionados y luego procesados con la herramienta para obtener huellas de edificios. Sugiero seleccionar varios edificios y luego ejecutar la herramienta. La herramienta está adjunta al final como ExtractBuildingsFromImageryPro1_3v4.zip. Después de descomprimir los datos, asegúrese de usar definir proyección en la featureclass de edificios en el output.gdb para que coincida con la proyección de sus imágenes.<\/P><\/P><\/P>Hubo un problema con las rutas que no se procesaban correctamente. Creo que ahora está solucionado ya que el espacio de trabajo se define como in_memory. Si obtiene un error por ruta no encontrada, abra el modelo y cambie las rutas en el modelo para C:\Users\student\Documents\ArcGIS\Projects\ImageryExtraction\ImageryExtraction.gdb\seg<\/SPAN> a in_memory\seg. Tendrá que hacer esto para todos los archivos temporales en el modelo. Gracias a quienes reportaron este error.<\/P><\/P><\/P>Descargue la herramienta desde mi blog.<\/LI>Descomprima la herramienta.<\/LI>Obtenga imágenes de alta resolución del área de la cual desea extraer edificios (6 pulgadas es lo mejor).<\/LI>Abra ArcGIS Pro (1.3 o superior) Mapa<\/LI>Agregue sus imágenes al Mapa<\/LI>En la pestaña Insertar, use Agregar Caja de Herramientas para añadir la caja ExtractBuildingsFromImageryPro1_3.tbx.<\/LI>Ejecute la herramienta Imagery Segmentation to Vector contra sus imágenes. Esto toma tiempo (5 minutos para un archivo de 349mg). Si es cuatro bandas, ingrese 4 2 1 en los Índices de Banda.<\/LI>Cambie la simbología del resultado a rojo y colóquelo encima de la imagen segment mean shift. Debería verse así:<\/LI>Apague la imagen de Segmentación y con las imágenes originales debajo, seleccione polígonos para casas que desea convertir en polígonos de edificio. <\/SPAN><\/LI>Por primera vez, ejecute Definir Proyección en la featureclass de edificios para que coincida con sus imágenes.
Abra la herramienta Segmentation Vectors to feature y ejecútela después de configurar sus parámetros.
Inspeccione sus resultados :
< /OL >Herramienta Regularize Building Footprint - Limpie las características extraídas sin procesar extraídas del lidar o imágenes.< SPAN style = "color :#303030 ;letter-spacing :-0.05pt ;font-family :Arial ;font-size :12pt ;mso-ascii-font-family :Arial ;mso-fareast-font-family : '\uff2d\uff33 \uff30\u30b4\u30b7\u30c3\u30af' ;mso-bidi-font-family :Arial ;mso-ascii-theme-font :major-latin ;mso-fareast-theme-font :minor-fareast ;mso-color-index :1 ;mso-font-kerning :12.0pt ;language :en-US ;mso-style-textfill-type :solid ;mso-style-textfill-fill-themecolor :text1 ;mso-style-textfill-fill-color :white ;mso-style-textfill-fill-alpha :100.0 % ;">Los edificios extraídos sin procesar del Lidar y Raster necesitan ser limpiados y la herramienta Regularize Building Footprint hace un buen trabajo, pero es difícil saber qué ingresar en los parámetros para la herramienta.& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;& nbsp;. He estado trabajando en modelos para hacerlo y ahora tengo< SPAN style = "color:#303030 ;letter-spacing :-0.05pt;font-family :Arial;font-size :12pt;mso-ascii-font-family :Arial;mso-fareast-font-family :'\uff2d\uff33 \uff30\u30b4\u30b7\u30c3\u30af';mso-bidi-font-family :Arial;mso-ascii-theme-font :major-latin;mso-fareast-theme-font:minor-fareast;mso-color-index :1;mso-font-kerning :12.0pt;languag e=en-US;mso-style-textfill-type=solid;mso-style-textfill-fill-themecolor=text1;mso-style-textfill-fill-color=white;m so-style-textfill-fill-alpha=100.0%;"> añadido un modelo para ejecutar la herramienta Regularize Building Footprint contra edificios extraídos por raster llamado Clean Extracted Buildings Using Regularize Building Footprint Tool. <\/P><\/P>Tiempo para ejecutar ambas operaciones (Advertencia: La herramienta usará la mayor parte de la potencia de cálculo de su sistema):<\/SPAN><\/P><\/P>Intel Core i7 - 4600 CPU @ 2.10 GHz 2.70 GHz U con 8GB de memoria (64bit) - 25 minutos para 10,000 features<\/SPAN><\/P>Intel Xeon CPU E5 -1620 0 @ 3.60GHz 3.60GHz con 32GB de memoria (64bit) - 8 minutos para 10,000 features<\/SPAN><\/P><\/P>Abajo en blanco están los edificios extraídos originales sin procesar del lidar y Aqua Green es el CircleBuildings extraído. Los rojos son la salida de Buildings. Los edificios bajo el circlebuildings deberán borrarse<\/SPAN><\/P><\/P><\/P>1. Para usar esta herramienta, primero ejecute un Dissolve en sus features con multipart desactivado.<\/SPAN><\/P><\/P><\/P>2. Descargue y descomprima el archivo adjunto Clean Extracted Buildings Using Regularize Building Footprint Tool.zip<\/P>Ejecute primero la herramienta "1. Regularize Building Footprints For Circle (Feet Projection)" dentro. <\/P><\/P>3. Ejecute a continuación la herramienta "2. Regularize Building Footprints All Non-Circle (Feet Projection)".<\/P><\/P>4. Revise primero su salida CircleBuildings, no deberían ser tantos y debe verificar si realmente son<\/P>círculos. Si no, elimine. Esto debería recoger tanques, edificios redondos, gazebos, etc. Si un edificio no es completamente redondo,<\/P>pero es como un hexágono o algo similar, intente usar el botón generalizar en edición con un valor de 1,2 o3 para que los lados<\/P>queden rectos como un hexágono u octágono. . Añada estructuras circulares adicionales que vea, a menudo están cerca unas de otras. . Guarde cuando<\/P>termine.<\/P><\/P>5. Añada su featureclass Buildings y seleccione por ubicación aquellos cuyo centroide esté en CircleBuildlings, luego elimine.<\/P><\/P>6. Revise sus buildings. . Si tiene problemas generales de extracción, por favor revise el tipo de edificio. . Vaya a la herramienta "2.<\/P> < P style = " margin p x *& nbsp;. Issue found with Clean Extracted Buildings Using Regularize Building Footprint Tool & nbsp;. not working if there are no large buildings (over25 ,000 sq. feet), looking into how to fix this.< \/ P > :< \/ P > Ejecute Circles primero siete veces – Regularize Building Footprint con Tolerancia2 ,4 ,5 ,6 ,8 , y Tolerancia12 y< \/ span >< \/ P > las consultas de definición abajo:< \/ span >< \/ P > 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%;">Tolerancia 1.6 con los vectores sin procesar usando una consulta de definición de Shape_Area <= 220 Y Shape_Area > 175<\/SPAN><\/P>Tolerancia 1.8 con los vectores sin procesar usando una consulta de definición de Shape_Area <= 240 Y Shape_Area > 220<\/SPAN><\/P>Tolerancia 2 con los vectores sin procesar usando una consulta de definición de Shape_Area <= 1400 Y Shape_Area > 240<\/SPAN><\/P>Tolerancia 4 con los vectores sin procesar usando una consulta de definición de Shape_Area <= 2500 Y Shape_Area > 1400<\/SPAN><\/P>Tolerancia 5 con los vectores sin procesar usando una consulta de definición de Shape_Area <= 3000 Y Shape_Area > 2500<\/SPAN><\/P>Tolerancia&nbsp ;6 con los vectores sin procesar usando una consulta de definición de Shape_Area &lt ;=8000&nbsp ;Y Shape_Area &gt ;3000<\/SPAN><\/P>Tolerancia&nbsp ;8 con los vectores sin procesar usando una consulta de definición de Shape_Area &lt ;=12000&nbsp ;Y Shape_Area &gt ;8000<\/SPAN><\/P>Tolerancia 12 con los vectores sin procesar usando una consulta de definición de Shape_Area &gt ;12000<\/SPAN><\/SPAN><\/P> <\/SPAN><\/P>Fusiona la salida con BuildingsCircle y luego aplica una consulta de definición de STATUS = 0. . Esto debería funcionar para extraer círculos automáticamente.<\/SPAN><\/P>Sugiero revisar los círculos para asegurarse de que son correctos, eliminando cualquiera que no lo sea.& nbsp;. <\/SPAN><\/P>Querrás comenzar a editar en los círculos fusionados, selecciona todo y usa la generalización con un valor de entrada de o .1 pie.<\/SPAN><\/P> <\/SPAN><\/P>Esto convierte los círculos perfectos en otros con vértices que se usan como entrada para Local Government Scenes.<\/SPAN><\/P><\/P>Selecciona los edificios extraídos sin procesar que no son círculos usando Seleccionar por Ubicación.& nbsp;. Ejecuta un Disolver con Crear características multipartes desactivado.& nbsp;. <\/SPAN><\/P> <\/P><\/P>Divide en tres grupos por tamaño a continuación y ejecuta Regularizar huella del edificio con los valores a continuación:<\/SPAN><\/P><\/P>•<\/SPAN> Pequeño:"<\/SPAN> Shape_Area"< 5000 sq. pies

Right_ANGELS Tolerancia2, Densification2, Precisión 0.25

Medio:"Shape_Area">= <\/SPAN><\/SPAN>25000 <\/SPAN>Right_ANGLES_AND_DIAGONAL<\/SPAN><\/SPAN>
Tolerancia 4, Densificación 4, <\/SPAN> Precisión 0.25, Penalización Diagonal 1.5<\/SPAN><\/SPAN><\/P><\/P><\/P><\/SPAN><\/P><\/P>Después de ejecutar estos grupos, combine el resultado con BuildingsNonCircle y ejecute disolver sobre las entidades. eje una Eliminación de Parte de Polígono en las entidades con al menos 20 pies cuadrados y elimine solo la parte contenida. &nbspEsto es para eliminar los agujeros interiores en los polígonos limpiados. &nbspEjecute una Simplificación de Edificios con tolerancia de simplificación 2 y luego nuevamente con eso como entrada con una tolerancia de simplificación 4.<\/SPAN><\/SPAN><\/P><\/P> <\/SPAN><\/P><\/P><\/P>Combine sus edificios circulares y los BuildingsNonCirclesb4 para obtener sus edificios terminados. &nbspComience a revisar los edificios contra imágenes lidar para encontrar errores (siempre parece haber errores).<\/P><\/P>Extracción de Edificios y Árboles de Lidar<\/SPAN><\/P>< \/ P > Aquí; también hay una herramienta para la Identificación de Edificios y Árboles a partir de Lidar desarrollada por Joe McGlinchy y Roslyn Dunn. És un trabajo en progreso y actualmente solo funciona con ArcGIS Desktop 10.3.1. Él recopila huellas de edificios y árboles circulares que pueden ser procesados mediante la herramienta Regularize Building Footprint mencionada arriba.< \/ P > Otra herramienta para extraer árboles que desarrollé, que utiliza Lidar e imágenes NAIP, está

aquí< \/ A >. &nbspUna versión para Pro están adjunta (Trees From LIDAR and NAIP_Pro).< \/ P >A&ntildeadiendo color a los techos de los edificios para mostrar en ArcGIS Pro (Scene Viewer< \/ A >)< \/ STRONG >< \/ SPAN >< \/ P >< \/ P > A&ntildeadir color a los techos en ArcGIS Pro realmente hace que las escenas se vean mas realistas con Local Government Scenes.&nbspHace mas de un año desarroll&eacute un modelo para agregar colores hexadecimales a los edificios como un campo para mostrar en ArcGIS Pro. &nbspCon el proceso Local Government Scenes, puede usar un campo con un color hexadecimal para agregar colores para mostrar en ArcGIS Pro usando sus capas. a menudo uso dos fuentes de im&aacutegenes de diferentes tiempos y uso una para el techo y otra para el color lateral (agregue otro campo llamado color lateral).&nbspAunque usarlo para el color lateral es ideal, es una forma fác il de colorear edificios en masa.< \/ P >< \/ P > Instrucciones:< \/ P >< \/ P > 1 . Recorte la imagen sin hojas de alta resolución al &aacuterea donde tiene edificios a los que desea agregar color.< \/ P > 2 . Remuestree a aproximadamente 3 pies, esto permitir&aacute que la herramienta Segment Mean Shift funcione mucho mas r&aacutepido.< \/ P > 3 . Ejecute Segment Mean Shift en la imagen sin hojas con par&aacutemetros predeterminados.< \/ P > 4 . Ejecute Agregar Color Hexadecimal desde Im&aacutegenes a Techos 10 .4 o ArcGIS Pro Agregar Color Hexadecimal a Techos (adjunto abajo).< \/ P > 5 . Descargue los archivos de capa Old Local Government Scenes (adjuntos abajo comprimidos) y descompr&iacutemalos.< \/ P > 6 . Mueva sus edificios a la escena 3D en Pro y haga clic en propiedades para establecer la elevación como absoluta.< \/ P > 7 . En Pro, ejecute Aplicar Simbología desde Capa usando LOD2BuildingShellsFloors_feet.lyrx o LOD2BuildingShells_meters.lyrx en el directorio OldLocalGovernmentScenesLayerFiles que descomprimio.& nbspLa nueva versión de Local Government Scenes no tiene color para techos o lados del edificio, esto probablemente cambiar&aacute en versiones futuras.< \/ P >7 .& nbspEstablezca el color del techo a su campo de color hexadecimal en la simbología y aplique.< \/ P >< \/ P > A continuación se muestra St.Louis en 3D con color agregado a los edificios.& nbspEl Arco fue la &uacutenica caracter&iacutestica dibujada.& nbspTodos los colores del edificio se aplicaron mediante im&aacutegenes utilizando el proceso anterior.& nbspLos árboles se extrajeron usando Trees From LIDAR (aquí< \/ A >).& nbspLos edificios se segmentaron usando Segment Mean Shift sobre un DEM reclasificado recortado a las huellas del edificio y luego convertido a pol&iacutegonos.& nbspEsto da a los edificios 3D alturas multiples para diferentes secciones.& nbspAgregar Color Hexadecimal desde Im&aacutegenes a Techos se ejecut&oacute sobre pol&iacutegonos individuales para las diferentes partes del techo.< \/ P >< 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;"><\ /P>

Aqu&iacute hay otro ejemplo sobre St.Charles, County, MO.& nbspRecientemente hicieron sus datos Open Data a través de MSDIS y el siguiente fin de semana proces&eacute160,000 edificios a 3D usando Local Government Scenes con un modelo alrededor para ejecutar por lotes(

Scene Viewer<\ /A>).& nbspPara los colores laterales del edificio us&eacute esta f&oacutermula para calcular un color a partir de los colores del techo. Los colores laterales no son correctos, pero le da una vista más realista.<\/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>

Otra opción es usar el color real de los lados de los edificios. Para hacer esto, comienza con una búsqueda de imágenes en un navegador web para casas en el área que estás trabajando. Luego puedes capturar las imágenes, georreferenciarlas lado a lado en ArcMap y unirlas en un mosaico. Ejecuta Segment Mean Shift en la imagen mosaica y luego crea una clase de entidad para poner edificios falsos sobre las áreas de las que quieres obtener colores. Ejecuta la herramienta Add_Hex_Color_From_RGB_Imagery_to_Roofs_10.4 en esos edificios y luego expórtala a una tabla. Mira cuántos registros hay en la tabla, agrega un campo a los edificios llamado join y calcúlalo con un número aleatorio igual al número de registros en la tabla (importa random en el bloque de código, random.randint(1,100) , siendo 100 el número de registros). Luego agrega un campo a tu edificio llamado side color y calcúlalo después de hacer una unión con la tabla usando el objectid y el campo join. Usa este campo como entrada para la simbología de los lados de los edificios. Un ejemplo está abajo para Fort Bend, TX, donde extraímos las huellas de los edificios del lidar, las procesamos con Local Government Scenes y luego aplicamos color a los techos desde imágenes y los colores laterales aleatoriamente.<\/P>

<\/P>

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

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

<\/P>

St. Charles en 3D

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

<\/P>

Llevando tu Inventario de Árboles y Mostrándolo en 3d en ArcGIS Pro:<\/STRONG><\/SPAN><\/P><\/P>Hace un tiempo un parque local me preguntó qué hacer con su nuevo Inventario de Árboles para su junta que fue entregado como una clase de entidad puntual con género. Sugerí mostrarlo usando ArcGIS Pro en 3D. A continuación están las instrucciones para agregar una altura a los puntos de árboles desde Lidar y luego calcular un campo para el Género a partir del nombre científico. Calculando la altura y un campo solo con el Género, puedes mostrar los árboles en 3D en Pro fácilmente. Aquí hay un enlace a todos los árboles disponibles:<\/P>

Tipos de géneros de árboles soportados013ArcGIS Pro | ArcGIS for Desktop<\/A><\/P>

Con los árboles, tenerlo basado en Género significa que un Arce Plateado (Acer saccharinum) se mostrará en el campo solo como "Acer" para manejar la simbología. Con el enlace anterior, puedes ver que Acer se mostrará como un árbol Arce Azucarero. Generalmente es lo suficientemente cercano para dar una buena representación del tipo de árbol.<\/P>

<\/P>

Revisando investigaciones, parece que obtener la altura de un árbol desde el DBH no funciona muy bien. Por eso, sugeriría usar el Lidar o la altura si está en el Inventario de Árboles. El único problema que encontré fue con árboles recién plantados que el lidar no tenía o eran más pequeños cuando se hizo el lidar.<\/P>

<\/P>

Revisa los puntos, probablemente querrás reproyectar los puntos al mismo sistema que tu lidar. <\/SPAN><\/LI>El primer paso fue descargar el archivo LAS del área
 (formato lidar de datos). <\/SPAN><\/LI>El segundo paso fue crear un dataset las con los archivos las
 <\/SPAN><\/LI>Luego usamos LAS Point as Statistics as Raster usando el método Z_Range
 para obtener un DEM con las alturas de los árboles.<\/SPAN><\/LI>Ejecuté una herramienta Focal Statistics con Max círculo 3 para
 obtener la extensión del punto más alto para que intersecte el punto.<\/SPAN><\/LI>Ejecuté Add Surface information para obtener la Z del
 raster a los puntos del Inventario de Árboles, esta es la altura del árbol.<\/SPAN><\/LI>Agrega un campo llamado Genus y calcula el campo usando
 el nombre científico. ArcGIS necesita
 solo el género como entrada, así que usé este script vb en 10.4 para hacerlo desde el nombre científico en python (ArcGIS Pro) !BOTANICAL!.split(' '<\/>),1)[0] o <\/>vb script Split([BOTANICAL] , " ",2)(0) <\/>
Abre ArcGIS Pro Scene y agrega los árboles como

Agregar capas predefinidas<\/>
Árboles realistas usando ese campo Genus, obtienes una simbología más correcta
para los árboles y el campo altura. Esto no funciona con todos tus
árboles, pero muchos están soportados.& nbsp;
& nbsp;. Nuevamente, puedes ver cuáles están soportados en este enlace: Tipos de géneros de árboles soportados013ArcGIS Pro | ArcGIS for Desktop<\/>

Agrega imágenes y revisa tus árboles para asegurarte que se vean correctos.& nbsp;. Un error común es que la altura puede estar en metros o pies, pero elegí la incorrecta. Aquí hay algunos datos con los que lo probé desde Amherst, MA con los datos fuente y un subconjunto procesado adjunto llamado SubsetOfAmherstMATreeInventory2010.zip : & nbsp;. Inventario de Árboles : < A href = El Missouri Botanical Garden ha hecho para sus árboles para renderizar en 3D. Esri soporta alrededor de 60 géneros de árboles para renderizado 3D. Luego utilicé el lidar para obtener una fórmula simple para convertir el DBH (Diámetro a la Altura del Pecho) a la altura del árbol. No es exacto, pero renderizó los árboles bastante bien. Existen fórmulas muy complejas que podría investigar más a fondo para vincular el DBH con la altura, esto podría permitir que muchos inventarios de árboles se conviertan fácilmente en 3D. La tabla para calcular un género similar está abajo con las otras descargas llamadas MissouriBotanicalGardenSubstitutesForTreeGenus.txt.<\/P>

<\/P>

Escena mostrando árboles hechos en 3D usando una fórmula aplicada al DBH ((DBH x 2)+20). 5,150 de los 9,153 géneros de árboles fueron soportados usando los 63 géneros <\/SPAN>tipos soportados por Esri ArcGIS Pro<\/A>. 2,120 árboles no coincidieron con el género soportado para renderizado, pero se usó una tabla suministrada por el Missouri Botanical Garden para géneros de árboles similares. Los árboles en blanco no fueron soportados ni mapeados para la simbología de género, esto incluyó 1,883 árboles y arbustos. Para calcular el género, usé el campo SPP (Nombre Científico) con <\/SPAN> !SPP!.split(<\/SPAN>' '<\/SPAN>,1)[0]<\/SPAN><\/P>Scene Viewer<\/A> - Asegúrate de activar el lidar para ver qué tan cerca o no están los árboles en altura.<\/P>

<\/A><\/P>

Arthur Crawford - Esri<\/P>

<\/P>

Agradecimientos especiales al Condado de St. Charles, Monica Balestreri, MSDIS, Michael Owen, Town of 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 y los demás que me han ayudado.<\/P><\/BODY><\/HTML>

Encuentra más publicaciones etiquetadas con

Comentarios

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

Únete a la Conversación

Los miembros registrados pueden publicar, seguir actualizaciones y más. ¿Nuevo aquí? Regístrate gratis.

Enlaces Rápidos

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