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

How to get the MIN and Max from a while loop - python

DanielAbera

Hi Eveyone,

I got the following python code to get the geometry(X, Y, Z) of a polyline feature. It can print all the values. But what I am looking is that, how do we get the MIN and MAX values for Z -value which is an elevation value in my case; and pass it as a variable. I am new to python and having trouble to figure it out. The code is as follows. Thanks in advance for your help.

Daniel A.

# Import native arcgisscripting module # import arcgisscripting  # Create the geoprocessor object # gp = arcgisscripting.create(9.3)  infc = gp.GetParameterAsText(0)    # Identify the geometry field # desc = gp.Describe(infc) shapefieldname = desc.ShapeFieldName  # Create search cursor # rows = gp.SearchCursor(infc) row = rows.Next()  # Enter while loop for each feature/row # while row:     # Create the geometry object     #     feat = row.GetValue(shapefieldname)      # Print the current multipoint's ID     #     print "Feature " + str(row.getvalue(desc.OIDFieldName)) + ":"     partnum = 0      # Count the number of points in the current multipart feature     #     partcount = feat.PartCount      # Enter while loop for each part in the feature (if a singlepart feature     # this will occur only once)     #     while partnum < partcount:         # Print the part number         #         print "Part " + str(partnum) + ":"         part = feat.GetPart(partnum)         pnt = part.Next()         pntcount = 0          # Enter while loop for each vertex         #         while pnt:             # Print x,y coordinates of current point             #             print pnt.x, pnt.y, pnt.z             pnt = part.Next()             pntcount += 1              # If pnt is null, either the part is finished or there is an              #   interior ring             #             if not pnt:                  pnt = part.Next()                 if pnt:                     print "Interior Ring:"         partnum += 1      row = rows.Next()

Find more posts tagged with

Accepted answers

jskinner

Did you create the empty list? It's hard to see the bold in the pasted code before. Here is the subset of it:

while row:     list = []

All comments

FerialAssmani

Lot of thanks JSkinn3 😉

jskinner

Yes, that is correct. The MIN value will be the first value in the list, i.e. list[0].

FerialAssmani

Hi JSkinn3, Thank you it work 10/10
Please for the lower case (min value), is it correct to write ?
*************
dict = {}

for v in centerDisplacement.values:
dict[v.nodeLabel] = v.data[0]

#sort dictionary
list = sorted(dict, key=dict.get)
max = list[-1]
min = list[0]

print 'Max displacement = ', dict[max]
print 'For the Node label = ', max

print 'Min displacement = ', dict[min]
print 'For the Node label = ', min
*************
Thanks a lot

jskinner

Hi Ferial,

You could accomplish this using a dictionary:

dict = {}

for v in centerDisplacement.values:
    dict[v.nodeLabel] = v.data[0]    

#sort dictionary
list = sorted(dict, key=dict.get)
max =  list[-1]

print 'Max displacement = ', dict[max]
print 'For the Node label = ', max

FerialAssmani

Dear all, For my case I am doing a simple scripting with Python for a Abaqus odb reading, so i want to extract the max value and its Node label of my list (reading via a loop)
#
for v in centerDisplacement.values:
    print 'Node label          = ', v.nodeLabel
    print 'X displacement = ', v.data[0]
#
-----------------------------
Node label          = 1
X displacement =0.000
Node label          = 2
X displacement =1.050
Node label          = 3
X displacement =2.306
--------------------------------
So, I want to print at next

Max displacement = 2.306
For the Node label= 3

Thanks

DanielAbera

Thanks alot it works !!!
I just made a little change to get just the min and max, so i put the code in the main loop as follows.

# Identify the geometry field
#
desc = gp.Describe(infc)
shapefieldname = desc.ShapeFieldName

# Create search cursor
#
rows = gp.SearchCursor(infc)
row = rows.Next()

# Enter while loop for each feature/row
list = []
while row:
   
    # Create the geometry object
    #
    feat = row.GetValue(shapefieldname)

    # Print the current multipoint's ID
    #
    
    partnum = 0

    # Count the number of points in the current multipart feature
    #
    partcount = feat.PartCount

    # Enter while loop for each part in the feature (if a singlepart feature
    # this will occur only once)
    #
    while partnum < partcount:
        # Print the part number
        #
        
        part = feat.GetPart(partnum)
        pnt = part.Next()
        pntcount = 0

        # Enter while loop for each vertex
       
        while pnt:
            # Print x,y coordinates of current point
            ##print pnt.x, pnt.y, pnt.z
            list.append(pnt.z)
            pnt = part.Next()
            pntcount += 1

            # If pnt is null, either the part is finished or there is an 
            #   interior ring
            #
            if not pnt: 
                pnt = part.Next()
                if pnt:
                    print "Interior Ring:"
        partnum += 1
    row = rows.Next()
list.sort()
min = list[0]
max = list[-1]
print min, max
   

del row, rows

DanielAbera

Looks like I forgot to include the append portion:

while pnt:
    # Print x,y coordinates of current point
    ##print pnt.x, pnt.y, pnt.z
    list.append(pnt.z)
    pnt = part.Next()
    pntcount += 1

So, the entire code I was able to get working is:

# Import native arcgisscripting module
#
import arcgisscripting

# Create the geoprocessor object
#
gp = arcgisscripting.create(9.3)

infc = gp.GetParameterAsText(0)
 

# Identify the geometry field
#
desc = gp.Describe(infc)
shapefieldname = desc.ShapeFieldName

# Create search cursor
#
rows = gp.SearchCursor(infc)
row = rows.Next()

# Enter while loop for each feature/row
#
while row:
    list = []
    # Create the geometry object
    #
    feat = row.GetValue(shapefieldname)

    # Print the current multipoint's ID
    #
    print "Feature " + str(row.getvalue(desc.OIDFieldName)) + ":"
    partnum = 0

    # Count the number of points in the current multipart feature
    #
    partcount = feat.PartCount

    # Enter while loop for each part in the feature (if a singlepart feature
    # this will occur only once)
    #
    while partnum < partcount:
        # Print the part number
        #
        print "Part " + str(partnum) + ":"
        part = feat.GetPart(partnum)
        pnt = part.Next()
        pntcount = 0

        # Enter while loop for each vertex
        #
        while pnt:
            # Print x,y coordinates of current point
            ##print pnt.x, pnt.y, pnt.z
            list.append(pnt.z)
            pnt = part.Next()
            pntcount += 1

            # If pnt is null, either the part is finished or there is an 
            #   interior ring
            #
            if not pnt: 
                pnt = part.Next()
                if pnt:
                    print "Interior Ring:"
        partnum += 1
    list.sort()
    min = list[0]
    max = list[-1]
    print min, max

    row = rows.Next()

del row, rows

Thanks alot it works !!!
I just made a little change to get just the min and max, so i put the code in the main loop as follows.

jskinner

Yes, i used it. May be can be more specific where i should put the append, because there are 3 while loops and which one to append ( pnt or pnt.z ).

thanks again

Looks like I forgot to include the append portion:

while pnt:
    # Print x,y coordinates of current point
    ##print pnt.x, pnt.y, pnt.z
    list.append(pnt.z)
    pnt = part.Next()
    pntcount += 1

So, the entire code I was able to get working is:

# Import native arcgisscripting module
#
import arcgisscripting

# Create the geoprocessor object
#
gp = arcgisscripting.create(9.3)

infc = gp.GetParameterAsText(0)
 

# Identify the geometry field
#
desc = gp.Describe(infc)
shapefieldname = desc.ShapeFieldName

# Create search cursor
#
rows = gp.SearchCursor(infc)
row = rows.Next()

# Enter while loop for each feature/row
#
while row:
    list = []
    # Create the geometry object
    #
    feat = row.GetValue(shapefieldname)

    # Print the current multipoint's ID
    #
    print "Feature " + str(row.getvalue(desc.OIDFieldName)) + ":"
    partnum = 0

    # Count the number of points in the current multipart feature
    #
    partcount = feat.PartCount

    # Enter while loop for each part in the feature (if a singlepart feature
    # this will occur only once)
    #
    while partnum < partcount:
        # Print the part number
        #
        print "Part " + str(partnum) + ":"
        part = feat.GetPart(partnum)
        pnt = part.Next()
        pntcount = 0

        # Enter while loop for each vertex
        #
        while pnt:
            # Print x,y coordinates of current point
            ##print pnt.x, pnt.y, pnt.z
            list.append(pnt.z)
            pnt = part.Next()
            pntcount += 1

            # If pnt is null, either the part is finished or there is an 
            #   interior ring
            #
            if not pnt: 
                pnt = part.Next()
                if pnt:
                    print "Interior Ring:"
        partnum += 1
    list.sort()
    min = list[0]
    max = list[-1]
    print min, max

    row = rows.Next()

del row, rows

DanielAbera

minval, maxval = None, None
cur = gp.SearchCursor(infc)
for row in iter(cur.next, None):
   colval = row.getValue(field_name)
   if minval is None or colval < minval:
      minval = colval
   if maxval is None or colval > maxval:
      maxval = colval

del row
del cur

Hi Jason,

Thanks for the reply

can you rewrite this so that it can be used to extract shape field values ( X, Y, & Z)

Thanks,

DanielAbera

Did you create the empty list? It's hard to see the bold in the pasted code before. Here is the subset of it:
while row:
    list = []

Yes, i used it. May be can be more specific where i should put the append, because there are 3 while loops and which one to append ( pnt or pnt.z ).

thanks again

JasonScheirer

minval, maxval = None, None
cur = gp.SearchCursor(infc)
for row in iter(cur.next, None):
   colval = row.getValue(field_name)
   if minval is None or colval < minval:
      minval = colval
   if maxval is None or colval > maxval:
      maxval = colval

del row
del cur

jskinner

Did you create the empty list? It's hard to see the bold in the pasted code before. Here is the subset of it:

while row:     list = []

DanielAbera

Hi Jake,

Thanks for the reply. I did try that. I get an error message that float object has no attribute append. How did u do the append. What i did is ptn.z.append(list).

jskinner

Hi Daniel,

Are you looking to return only the MIN and MAX Z values? If so, you can create a list, append the Z values to the list, sort the list and then print out the min and max values. Ex: (lines in bold)

# Import native arcgisscripting module
#
import arcgisscripting

# Create the geoprocessor object
#
gp = arcgisscripting.create(9.3)

infc = gp.GetParameterAsText(0)
 

# Identify the geometry field
#
desc = gp.Describe(infc)
shapefieldname = desc.ShapeFieldName

# Create search cursor
#
rows = gp.SearchCursor(infc)
row = rows.Next()

# Enter while loop for each feature/row
#
while row:
    list = []
    # Create the geometry object
    #
    feat = row.GetValue(shapefieldname)

    # Print the current multipoint's ID
    #
    print "Feature " + str(row.getvalue(desc.OIDFieldName)) + ":"
    partnum = 0

    # Count the number of points in the current multipart feature
    #
    partcount = feat.PartCount

    # Enter while loop for each part in the feature (if a singlepart feature
    # this will occur only once)
    #
    while partnum < partcount:
        # Print the part number
        #
        print "Part " + str(partnum) + ":"
        part = feat.GetPart(partnum)
        pnt = part.Next()
        pntcount = 0

        # Enter while loop for each vertex
        #
        while pnt:
            # Print x,y coordinates of current point
            ##print pnt.x, pnt.y, pnt.z
            pnt = part.Next()
            pntcount += 1

            # If pnt is null, either the part is finished or there is an 
            #   interior ring
            #
            if not pnt: 
                pnt = part.Next()
                if pnt:
                    print "Interior Ring:"
        partnum += 1
    list.sort()
    min = list[0]
    max = list[-1]
    print min, max

    row = rows.Next()

Join the Conversation

Signed in members can post, follow updates, and more. New here? Register a free account.

Quick Links

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