topic Re: Adding a feature class to terrain in python in Python Questions

Adding a feature class to terrain in python

BethanyMarshall — Tue, 01 Mar 2016 17:38:04 GMT

Hello all,

I am having a problem with the adding feature classes to the terrain in my python code. I am trying to go from las files to a terrain. The code I am using is listed below: The LAS_FD is the Las feature dataset

arcpy.env.workspace = LAS_FD

fcList = arcpy.ListFeatureClasses()

params = "fcList shape mass points <None> 0 48 <none> false <none> <none>"

print "Adding feature classes to the terrain"

arcpy.AddFeatureClassToTerrain_3d (in_terrain, params)

Thanks in advance for any help!

Bethany

Re: Adding a feature class to terrain in python

AdrianWelsh — Tue, 01 Mar 2016 17:40:21 GMT

What is the error message that you are getting?

Re: Adding a feature class to terrain in python

BethanyMarshall — Tue, 01 Mar 2016 17:46:27 GMT

File "C:\Python_ArcGIS\DTM.py", line 95, in <module>

arcpy.AddFeatureClassToTerrain_3d (in_terrain, params)

Error Info:

Object: Error in executing tool

Re: Adding a feature class to terrain in python

DanPatterson_Retired — Tue, 01 Mar 2016 17:46:52 GMT

where is this defined? LAS_FD

Re: Adding a feature class to terrain in python

DanPatterson_Retired — Tue, 01 Mar 2016 17:47:49 GMT

go to the Results window, and get the whole message, there should be more reported, including paths inputs options etc.

Re: Adding a feature class to terrain in python

BethanyMarshall — Sat, 11 Dec 2021 23:44:04 GMT

print "Importing modules..."
import os
import arcpy, sys, traceback, exceptions
try:
    #STEP 1: SET OVERWRITE OUTPUT TO TRUE
    arcpy.env.overwriteOutput = True
    print "Overwrite Output = True"


    #STEP 2: CHECK OUT EXTENSION
    #STEP 3: ACTIVATE THE SPATIAL ANALYST EXTENSION
    print "The code is importing modules"
    if arcpy.CheckExtension("3D"):
        if arcpy.CheckOutExtension("3D"):


            #STEP 4: SET THE WORKSPACE
            arcpy.env.workspace = r"C:\Python_ArcGIS\LAS"
            thePath = r"C:\Python_ArcGIS\LAS"
            LAS = r"C:\Python_ArcGIS\LAS\LAS.gdb"


            # Create the file gdb that will store the feature dataset
            if arcpy.Exists("LAS.gdb"):
                arcpy.Delete_management("LAS.gdb")
            arcpy.management.CreateFileGDB(thePath, "LAS.gdb")
            gdb = '{0}/{1}'.format(thePath, LAS)


            #Create Feature Dataset
            # CreateFeatureDataset_management (out_dataset_path, out_name, {spatial_reference})
            if arcpy.Exists("LAS_FD"):
                arcpy.Delete_management("LAS_FD")
            arcpy.management.CreateFeatureDataset(LAS, "LAS_FD", "C:\Python_ArcGIS\SpatialReference.gdb\NAD_1983_UTM_Zone_15N")
            fd = '{0}/{1}'.format(LAS, "LAS_FD")
            LAS_FD = r"C:\Python_ArcGIS\LAS\LAS.gdb\LAS_FD"


            # Create a LOOP for to reach all the files in the folder
            for root, dirs, files in os.walk(thePath):
                print "**********************"
                print "ROOT: " + root
                for aFile in files:
                    print root + "\\" + aFile


                        #Creating multipoint features from LAS
                        # LASToMultipoint_3d (input, out_feature_class, average_point_spacing, {class_code}, {return}, {attribute}, {input_coordinate_system}, {file_suffix}, {z_factor}, {folder_recursion})
                    extension1 = os.path.splitext(aFile)[0]
                    print extension1
                    LAS_File = extension1 + ".las"
                    if aFile == LAS_File:
                        MultiName = "DTM_MTP_" + aFile
                        arcpy.LASToMultipoint_3d (aFile, MultiName, 0.7, 2)
                        print MultiName


                        extension = os.path.splitext(MultiName)[0]
                        print extension
                        Multi_FC = extension + ".shp"


                        #Move Shapefiles to Geodatabase
                        #FeatureClassToGeodatabase_conversion (Input_Features, Output_Geodatabase)
                        arcpy.FeatureClassToGeodatabase_conversion (Multi_FC, LAS_FD)


                        #Create Terrain Dataset for each file
                    else:
                        print "No More Files to Process: Move to Terrain"
                        print "Creating Terrain"
                        #Create Terrain
                        # CreateTerrain_3d (in_feature_dataset, out_terrain_name, average_point_spacing, {max_overview_size}, {config_keyword}, {pyramid_type}, {windowsize_method}, {secondary_thinning_method}, {secondary_thinning_threshold})
                        arcpy.CreateTerrain_3d (LAS_FD, "DTM_Terrain", 0.7)


                        print "Adding terrain pyramid levels"
                        # Add Terrain Pyramid Level
                        # AddTerrainPyramidLevel_3d (in_terrain, pyramid_level_definition, {pyramid_type})
                        pyramid_level = ["1.5 1000", "3 2500", "6 5000", "12 10000", "24 20000", "48 50000"]
                        in_terrain = LAS_FD + "\\" + "DTM_Terrain"
                        print in_terrain
                        arcpy.AddTerrainPyramidLevel_3d (in_terrain, 'WINDOWSIZE', pyramid_level)


                        print "Creating a list of feature classes"
                        arcpy.env.workspace = LAS_FD
                        fcList = arcpy.ListFeatureClasses()
                        params = "fcList shape mass points <None> 0 48 <none> false <none> <none>"


                        print "Adding feature classes to the terrain"
                        #Add feature class to terrain
                        # AddFeatureClassToTerrain_3d (in_terrain, in_features)
                        # [[in_features, height_field, SF_type, group, min_resolution, max_resolution, overview, embed, embed_name, embed_fields, anchor],...]
                        arcpy.AddFeatureClassToTerrain_3d (in_terrain, params)


##                        print "Building Terrain"
##                        #Build Terrain
##                        # BuildTerrain_3d (in_terrain, {update_extent})
##                        arcpy.BuildTerrain_3d (in_terrain, "NO_UPDATE_EXTENT")
##
##                        #terrain to raster
##                        # TerrainToRaster_3d (in_terrain, out_raster, {data_type}, {method}, {sample_distance}, {pyramid_level_resolution})
##                        arcpy.TerrainToRaster_3d ("Terrain", in_terrain, "INT", "NATURAL_NEIGHBORS")
##                        print "All files are done"
        else:
                print "The 3D is not check out."


        arcpy.CheckInExtension("3D")
    else:
            print "The 3D extension is not present."




except:


    tb = sys.exc_info()[2]
    tbinfo = traceback.format_tb(tb)[0]
    pymsg = "PYTHON ERRORS:\nTraceback Info:\n" + tbinfo + "\nError Info:\n     " + str(sys.exc_info()[1])


    msgs = "ARCPY ERRORS:\n" + arcpy.GetMessages(2) + "\n"


    arcpy.AddError(msgs)
    arcpy.AddError(pymsg)


    print msgs
    print pymsg


    arcpy.AddMessage(arcpy.GetMessages(1))
    print arcpy.GetMessages(1)

Re: Adding a feature class to terrain in python

DanPatterson_Retired — Tue, 01 Mar 2016 17:52:00 GMT

And what did the Results window return compared to what you expected in your script?

Using the Results window—Help | ArcGIS for Desktop

ADDENDUM

code formatting would help others

Posting Code blocks in the new GeoNet

Code Formatting... the basics++

Re: Adding a feature class to terrain in python

BethanyMarshall — Tue, 01 Mar 2016 17:55:10 GMT

I am not doing this in ArcGIS. I am working on the script in pyscripter. I think this is the full error message.

PYTHON ERRORS:

Traceback Info:

File "C:\Python_ArcGIS\DTM.py", line 95, in <module>

arcpy.AddFeatureClassToTerrain_3d (in_terrain, params)

Error Info:

Object: Error in executing tool

Re: Adding a feature class to terrain in python

DanPatterson_Retired — Tue, 01 Mar 2016 18:02:44 GMT

In that case, I would get rid of the try except blocks and interpret the error message directly, it is returning pretty useless information as you can tell. Some people prefer them, but I generally find them useless with working with Arc* stuff

I will defer to someone else ... who can get your try-except block to reveal useful information if you prefer to keep them.

Re: Adding a feature class to terrain in python

BethanyMarshall — Tue, 01 Mar 2016 18:05:54 GMT

I ran it in Arc and this is the message I received. The add terrain pyramid level succeeded in the geoprocessing results. It does not even have the add feature class to terrain in the results.

PYTHON ERRORS:

Traceback Info:

File "<string>", line 84, in <module>

Error Info:

Object: Error in executing tool

Re: Adding a feature class to terrain in python

AdrianWelsh — Tue, 01 Mar 2016 18:36:43 GMT

Is line 84 this one?

arcpy.AddTerrainPyramidLevel_3d (in_terrain, 'WINDOWSIZE', pyramid_level)

Maybe try using double quotes instead of single quotes?

Re: Adding a feature class to terrain in python

BethanyMarshall — Tue, 01 Mar 2016 18:44:24 GMT

In ArcGIS, line 84 is this one:

arcpy.AddFeatureClassToTerrain_3d (in_terrain, params)

I am thinking it must be a problem with my params, but I don't know what is wrong.

params = "fcList shape mass points <None> 0 48 <none> false <none> <none>"

Re: Adding a feature class to terrain in python

AdrianWelsh — Tue, 01 Mar 2016 18:58:50 GMT

Hmmm, maybe it has to do with the params. Maybe try changing shape to SHAPE and the "<none>" values to "<None>" to stay with the naming conventions. Sometimes Python is picky about upper and lower case.

Though it is possible that the ListFeatureClasses list may be causing the problem too.

I was looking at the example here:

Add Feature Class To Terrain—Help | ArcGIS for Desktop

and it seems that the feature class used for the AddFeatureClassToTerrain_3d is always one feature class versus a list.

Re: Adding a feature class to terrain in python

DanPatterson_Retired — Tue, 01 Mar 2016 20:27:54 GMT

The documentation has a lot of preconditions which need to be met and confirmed

as indicated, I would check the documentation that Adrian listed. Add Feature Class To Terrain—Help | ArcGIS for Desktop Then I would do a manual run through the process in ArcMap (ie not script) using the same inputs as you want to do with your script. Once that works, you should compare the Results window to your code...in fact you can steal the results window outputs and use it to formulate the base of your code.

Re: Adding a feature class to terrain in python

AdrianWelsh — Wed, 02 Mar 2016 18:20:09 GMT

Bethany, did you find where the error was happening?

Re: Adding a feature class to terrain in python

BethanyMarshall — Thu, 03 Mar 2016 19:47:44 GMT

I actually changed the shape to SHAPE and the none to None values and I still got the error. I have tried to only input one feature class for the in_feature, but it doesn't seem to make a difference if it is one feature class or a list of feature classes. I have used the arc help, but their code is using a TIN and I am not using a TIN so their inputs are a little different.

I can't seem to figure out what the problem is with my code.

Re: Adding a feature class to terrain in python

DanPatterson_Retired — Thu, 03 Mar 2016 21:57:24 GMT

if shape is supposed to refer to the field in a featureclass/shapefile, it would be proper case as in Shape as would None

Re: Adding a feature class to terrain in python

AdrianWelsh — Thu, 03 Mar 2016 22:08:52 GMT

Dan, that is interesting because the documentation in that link shows it as caps:

However, further down below, it has Shape like you listed:

That would be nice if it were as easy as a case-sensitive issue... Let's hope!

Re: Adding a feature class to terrain in python

DanPatterson_Retired — Thu, 03 Mar 2016 22:11:58 GMT

tHat is sO conFusIng.... mayBe iT WaS caseFrEe DAy