Tento blog se vyvíjí, právě začíná a bude růst s časem:<\/STRONG><\/P><\/P>V posledních několika letech jsem publikoval lidarové barevné bodové mraky jako doplňky k 3D městům, jako například oblast St. Louis<\/A>, které jsem vytvořil pomocí 3D Basemap řešení a extrahováním budov z lidaru. Sean Morrish napsal skvělý blog<\/A> o mechanice publikování lidaru. V poslední době jsem se zabýval použitím barevných lidarových bodových mraků jako efektivního a relativně levného způsobu tvorby 3D Basemapů. Má to několik výhod, protože většina států, okresů a měst již má k dispozici lidar, většina má vysoce kvalitní snímky a všichni mají NAIP snímky, které jsou potřeba k vytvoření těchto scén. <\/P><\/P>Lidarové scény okresu Manitowoc:<\/STRONG><\/P>Některé okresy, města a státy to už dělají. Okres Manitowoc, WI., má 3D LiDAR Point Cloud<\/A> a také vytváří scény<\/A> s těmito daty. Okres Manitowoc také vytvořil skvělou story mapu ukazující, jak jejich lidar je použit jako barevný lidarový bodový mrak <\/SPAN>zde<\/A> a velmi doporučuji se na ni podívat. <\/SPAN><\/P><\/A><\/P>StoryMap ukazuje ve videu, jak zachytit vertikální, horizontální a přímé vzdálenosti LiDAR povrchů země na jakémkoli místě v celém okrese. Jak získat detailní měření malých LiDAR povrchových prvků, jako je hloubka příkopů k silnici. <\/P><\/P>Měřte vzdálenosti v lidarových bodových mracích vzhledem k LiDAR povrchům země pomocí střech domů jako příkladu pro zjištění výšky budovy.<\/P><\/P>Zde je jedna ze scén vrstev budov Manitowoc, kde měli kreativní nápad použít stejný barevný lidar k vytvoření falešných stran budov tím, že ukázali lidar klasifikovaný jako budovy několikrát za sebou, každou s mírně nižší nadmořskou výškou změnou offsetu v symbologii. Dále v blogu ukazuji, jak to udělat.<\/P>
<\/P>
Lidarové scény okresu Manitowoc:<\/STRONG><\/P>Některé okresy, města a státy to už dělají. Okres Manitowoc, WI., má
StoryMap ukazuje ve videu, jak zachytit vertikální, horizontální a přímé vzdálenosti LiDAR povrchů země na jakémkoli místě v celém okrese. Jak získat detailní měření malých LiDAR povrchových prvků, jako je hloubka příkopů k silnici. <\/P>
Měřte vzdálenosti v lidarových bodových mracích vzhledem k LiDAR povrchům země pomocí střech domů jako příkladu pro zjištění výšky budovy.<\/P>
Zde je jedna ze scén vrstev budov Manitowoc, kde měli kreativní nápad použít stejný barevný lidar k vytvoření falešných stran budov tím, že ukázali lidar klasifikovaný jako budovy několikrát za sebou, každou s mírně nižší nadmořskou výškou změnou offsetu v symbologii. Dále v blogu ukazuji, jak to udělat.<\/P>
Klobouk dolů Bruce Riestererovi z okresu Manitowoc, který to všechno sestavil před odchodem do důchodu, včetně nápadu použít body budov vícekrát s různými barvami pro zobrazení stran budov a nyní pracuje pro soukromý sektor, jeho novou práci najdete na
Leaf On High Resolution:<\/SPAN><\/SPAN><\/P><\/SPAN><\/P>Leaf On NAIP:<\/SPAN><\/P><\/SPAN><\/P>Leaf Off High Resolution:<\/SPAN><\/P><\/SPAN><\/P>Někdy mohou být barevné podzimní snímky před opadem listí skvělé pro scény, ale záleží na tom, co chcete ukázat.<\/SPAN><\/P><\/P>Způsob přidání stromů s větvemi a body jako listy:<\/STRONG><\/P>Zde je test, který jsem provedl pomocí symboliky Dead Tree k reprezentaci kmene a větví, zatímco listy pocházejí z barevného lidarového zobrazení. Můj nástroj Trees from LIdar nebo 3D Basemap Solution tree tools může vytvořit stromy z lidarových dat a pak mohou být body symbolizovány symbolem Dead Tree. Tento obrázek byl vytvořen v ArcGIS Pro, ne ve scéně.<\/P><\/P><\/P>Město Klaipėda 3D viewer (Litva):<\/STRONG><\/P>Město
Zde je další pohled s menším přiblížením:<\/P>
A ještě jeden s dalším oddálením:<\/P>
Některé výše uvedené scény používají NAIP snímky místo snímků s vyšším rozlišením. Proč, když je dostupné vyšší rozlišení? Často jsou snímky s vysokým rozlišením poněkud šikmé, zatímco NAIP je obvykle sbírán z větší výšky. S vyšší výškou získáte méně naklonění budov na snímcích a často rozestupy lidarových bodů nepodporují vyšší rozlišení snímků. V těchto případech je často preferován NAIP, ale vytvořte bodová mračna s několika las soubory, abyste viděli rozdíly.<\/P>
Zde je jeden nad Fairfax County, VA, používající jejich publikovaný lidar colorized (z jejich
Níže je práce, se kterou jsem pomohl státu Kentucky. Výsledek je v některých ohledech podobný práci neo-impresionistů Georgese Seurata a Paula Signaca, kteří průkopnili malířskou techniku nazvanou Pointillism.<\/P>
Kentucky lidar kolorovaný:<\/P>
S iluzí přidaných zdí a ztluštěnou střechou (Cumberland, KY):<\/P>
Zde je další pohled ukazující strany budovy:<\/P>
Je to iluze se stranami budov, barvy neodpovídají skutečné straně budov a rozdíly v intenzitě použitého lidarového snímání neukazují skutečné rozdíly na stranách jako vnímaná okna výše. Jako většina iluzí nebo reprezentací, klame oči, aby si myslely, že jsou to strany budov. Někdy, jako většina iluzí nebo reprezentací, to nefunguje ze všech úhlů, jak tato budova s nízkým počtem bodů:<\/P>
Celkově si myslím, že reprezentativní strany vám pomáhají identifikovat, co je budova a co není. Barvy stran budov nebudou odpovídat realitě, ani nezobrazují skutečná okna nebo skutečné dveře. Dělat to obvykle znamená velmi nákladný a časově náročný proces po vytvoření 3D modelů budov. Stále můžete generovat 3D budovy pomocí Local Government 3D Basemap Solution pro ArcGIS Pro a mít je poloprůhledné a používat pro výběr a analýzu bez jakéhokoli zbarvení. Také můžete použít rpks aplikované na budovy k přesnějšímu znázornění 3D modelů, ale to by opět byla reprezentace bez velmi podrobných dat. <\/P>
Jak snadno získat Lidar pro velkou oblast:<\/STRONG> Zde je
Základy Colorized Lidar Point Clouds, vytváření balíčků scénových vrstev a vizualizace:<\/STRONG><\/P>Zde je
@ArthurCrawford I wanted to share the St. Louis imagery web app linked in the blog (https://www.arcgis.com/home/webscene/viewer.html?webscene=6ece8b4050294139969fc4c10ad6b9b1) with my students in a Geovisualization and Cartography class. When we log in to our org we can see it, but we get a few layers that are "service not found" or "You do not have permission to access this resource". They are:
St Louis Iconic Waterfront Features 3D
City of St. Louis Lidar Dec 2016 Colorized
CityofStLouis_6_27_2018
Was trying to show them how the Colorize LAS tool worked and our computers were not happy (end of semester and not great graphics cards) but the colorized lidar was not showing up in this. 🙂 Thanks in advance.
Best, @shannonhwhite (not sure why its not defaulting to my wm-gis org account but either way I will get it).
Hi @ArthurCrawford,
Very cool looking scenes!
Can you colorize a point cloud scene layer with imagery?
Thanks,
Simon
At this point whatever I can display nicely online is what I will go with.
Hi Rickey,
Is your intent to use the multipatches for analysis and use a colorized point cloud for rendering? I did something like that for Chicago using the LOD1/LOD2 buildings with a transparency. With the 3D trees in as multipatches as per the Local Government 3D Basemap, you could do some nice 3D analysis like line of sight or flood estimations with it. With both set to a 100% transparency, then use the point cloud for rendering quickly.
Arthur
Arthur Crawford
This was made from lidar points using LAS building Multipatch. I have tried different Simplification Tolerance but I have found one that works.
Hi Rickey, good to hear from you. How's Ashland theses days? There is not a real good way to do that right now. It's something we have asked for, several universities and companies are working on. The goal would be to take a tin building and then flatten the roofs, superstructure and the sides of the buildings. I assume this is a building generated from a tin approach. There is some ideas of taking the planes generated to extend to the outside of the building polygon. Have you tried to segment roof tool with the Local Government 3D Basemap Solution?
Arthur CrawfordIs there a way to smooth multipatch?
Thanks, Arthur. I did download the NAIP locally and it matched the area of the point cloud, so I'm not sure what the issue was.
However, I performed the same colorization using ESRI's NAIP service and it worked perfectly!
Thanks again.
Hello John,
It looks like you are using NAIP as an imagery source. Did you download the NAIP locally and is it matching the area of the point cloud?
Arthur Crawford - Esri
Thanks Arthur!
I was able to resolve my vertical problem. I will seek out the Imagery Showcase booth (virtually) with other questions I have.
Speaking of which, I get to the Colorize Las stage of your workflow and no matter what adjustments I make to the tool settings, I still just end up with a grey colorized lasd. I'm thinking I must be missing something really basic but the solution hasn't been forthcoming! See attached.
What I did was saw the vertical projection in meters, but the one list was in feet. Here's the text in that projection file: USGS_LPC_CA_Wildfires_B1_2018_w1927n1436_LAS_2019.prj
PROJCS["NAD_1983_2011_Contiguous_USA_Albers",GEOGCS["GCS_NAD_1983_2011",DATUM["D_NAD_1983_2011",SPHEROID["GRS_1980",6378137.0,298.257222101,AUTHORITY["EPSG",7019]],AUTHORITY["EPSG",1116]],PRIMEM["Greenwich",0.0,AUTHORITY["EPSG",8901]],UNIT["Degree",0.0174532925199433,AUTHORITY["EPSG",9102]],AUTHORITY["EPSG",6318]],PROJECTION["Albers",AUTHORITY["Esri",43007]],PARAMETER["False_Easting",0.0,AUTHORITY["Esri",100001]],PARAMETER["False_Northing",0.0,AUTHORITY["Esri",100002]],PARAMETER["Central_Meridian",-96.0,AUTHORITY["Esri",100010]],PARAMETER["Standard_Parallel_1",29.5,AUTHORITY["Esri",100025]],PARAMETER["Standard_Parallel_2",45.5,AUTHORITY["Esri",100026]],PARAMETER["Latitude_Of_Origin",23.0,AUTHORITY["Esri",100021]],UNIT["Meter",1.0,AUTHORITY["EPSG",9001]],AUTHORITY["EPSG",6350]],VERTCS["NAVD88 (height) (ftUS)_1",VDATUM["North_American_Vertical_Datum_1988",AUTHORITY["EPSG",5103]],PARAMETER["Vertical_Shift",0.0,AUTHORITY["Esri",100006]],PARAMETER["Direction",1.0,AUTHORITY["Esri",100007]],UNIT["Meter",1.0,AUTHORITY["EPSG",9001]]]
Where the "Meter" is was originally "Feet", so I changed it. You can do this manually or in ArcGIS Pro by modifying the projection and saving it. I created a simple model with an iterator that would run through the las files and copy the projection file so it would match each las file.
It's looks to me yours might be in meters, but need to be set to feet, from the image you shared. Basically, the opposite of what I did with the Redlands lidar. You can create a projection file for one tile, change it from meters to feet and see if it works.
Here's one that might work, copy and paste below into a txt file. Then rename with las name and change from .txt file to .prj.
PROJCS["NAD_1983_2011_UTM_Zone_18N",GEOGCS["GCS_NAD_1983_2011",DATUM["D_NAD_1983_2011",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-75.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0],AUTHORITY["EPSG",6347]],VERTCS["NAVD88_height_(ftUS)",VDATUM["North_American_Vertical_Datum_1988"],PARAMETER["Vertical_Shift",0.0],PARAMETER["Direction",1.0],UNIT["Foot_US",0.3048006096012192],AUTHORITY["EPSG",6360]]
I'll be on the question and answer team for the Imagery Showcase booth this year for Elevation and Lidar if you want to ask some questions in person.
Arthur Crawford
Living Atlas Team/Community Maps
Hi Arthur - really interesting blog.
I have a question that is semi-related to your blog entry -- more specifically related to your #3 step towards the end of the blog: 3. Looked at the projection in the reports in the meta data downloaded and found it was not supported because it used meters instead of feet, modified a projection and then copied it using a simple model builder tool to give a projection file to each LAS file.
When you say "modified a projection", what exactly are you doing in ArcGIS Pro? I have a las file with a known projection and vertical coordinate system (NAD 1983 (2011) UTM Zone 18N, NAVD88 height - Geoid12A (metres)) but when I generate the las dataset it floats above the surface. This suggests to me that I have in an incorrect vertical reference so I'm trying to figure out how to modify the projection or reproject the lasd to get the stars to align.
See below.
Any thoughts?
Hi Trish, thanks for the comment. I added the location to it.
Hi! Interesting blog. I noticed you have Barnegat Bay Scene Viewer listed under the examples of some colorized lidar to look at. Barnegat Bay is in New Jersey, for those who are curious!
Přihlášení členové mohou přispívat, sledovat aktualizace a další. Nový zde? Zaregistrujte si bezplatný účet.
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.