Select to view content in your preferred language

3D Analyst performance - Beta 1

2165
3
02-01-2010 11:25 AM
BrianQuinn
Occasional Contributor
I haven't posted this until now - perhaps because it's not a problem.

Over the past year we've compiled a gob of data into our countywide Terrain Dataset:
- photogrammetric 1/4800 (10-ft) and 1/1200 (2-ft) contours and breaklines
- FEMA map modernization LiDAR bare-earth filtered points and breaklines
- NCALM high-res LiDAR along select seismic hazards

For best use with Terrain Dataset, the contours were vertex densified to 5-foot (horizontal) intervals along each contour line; several iterations of QA have repaired problem contours and stream-following water line features.  The FEMA LiDAR points were retained in original cloud positions, while we had to consume the NCALM data in their gridded product form, to reset Z reference from an up-to-the minute, research-grade ellipsoid into WGS84-NAVD88 (Geoid 2003).  Once into the geoid, the grid was converted to point-Z features and poured into the Terrain Dataset.

By the time we were ready to try out Beta 1, we had 960 million point-Z features in about 30 feature classes, water polyline-Z and hard break polyline-Z features, and a clipping polygon all set to go.  Under 9.3.1 we'd actually ground out some countywide 30cm grids, but it took more than 20 tiles to get the job done without crashing ArcMap.

Almost on a whim, we re-imaged an existing GIS workstation with Windows Server 2003 OS and chose to install 9.4 Beta 1 for its version of ArcGIS.  The OS change was to accomodate larger-sized seamless mosaics with ERDAS Imagine, where mosaics larger than about 70 GB tend to fail in non-server flavors of Windows.

About three weeks ago, we set up the Beta 1 system with a torture test: build the terrain, and then grind a single seamless countywide 40cm gridded surface, at pyramid level 0, from the Terrain Dataset that was approaching a giga-point of content.  The first reflexive attempt to run a 30-cm gridding produced a scratch output file of 133 GB, but we were not confident that ArcMap was running properly, so we revised the grid to 40cm, and 77 GB of output, and watched just over one Xeon core of processor loading running, and running, and running...

=> after 135 CPU hours and 7.5 calendar days, it finished.  The Beta 1 did hang when trying to create pyramids, but the output ".img" grid was serviceable in ERDAS, where we calculated statistics and pyramids.

As far as I can tell, the 3D analyst "Terrain to Raster" function has really grown some teeth since 9.3.1;  thank you, Beta crew!  The 40cm DSM product has been used already for some hydrology applications.
Tags (1)
0 Kudos
3 Replies
EricRice
Esri Regular Contributor
Great Test!  Thank you for your feedback.  The pyramid issue is known to occur when the pyramid file is going to exceed 4gb by itself.  We are working on a solution to that problem.  As you noted, the base pixels should be good to go.
0 Kudos
ClaytonCrawford
Esri Contributor
Thanks for the report. We're glad you've seen improvement in robustness. Are you using linear or natural neighbor interpolation during rasterization?

Regards,
Clayton Crawford
ESRI Software Products
0 Kudos
ClaytonCrawford
Esri Contributor
Was background processing enabled for GP when executing TerrainToRaster? There is a serious performance hit if that's the case - at least with beta 2 and probably beta 1. You may want to try another (smaller) example and see if you notice a difference when turning it off.

In general, terrain related GP tools are not ready for prime time when it comes to background processing. We're hoping to have this addressed by pre-release.

I'm curious about some of the data you're using and the terrain schema definition. Send me an e-mail if you have the time to describe this in more detail -> ccrawford @ esri dot com

Regards, Clayton
0 Kudos