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

Hello everybody I need some help about python scripting to create a square buffer around points feature class ?

MahmoudMdookh

If there is a problem by giving me the script , I just need some explanation about how can I start scripting the square buffer?

I'm a beginner python programmer, I'm looking forward to be professional in python , and I hope that any one will help me !

Thanks in advance...

Find more posts tagged with

Accepted answers

XanderBakker

... and if you want to create an output that contains the attributes of the input, you can try this:

import arcpy, os

def main():
    arcpy.env.overwriteOutput = True

    # inputs
    fc_in = r'D:\Xander\EPM\Datos pruebas\gdb\datos erase.gdb\points03'
    width = 0.005 # with of the output rectangle
    height = 0.005 # height of the output rectangle

    # output
    fc_out = r'D:\Xander\EPM\Datos pruebas\gdb\datos erase.gdb\rectangles3'

    # create empty output featureclass
    sr = arcpy.Describe(fc_in).spatialReference
    out_ws, out_name = os.path.split(fc_out)
    arcpy.CreateFeatureclass_management(out_ws, out_name, "POLYGON", fc_in, "DISABLED", "DISABLED", sr)

    # make fields list
    flds_use = correctFieldList(arcpy.ListFields(fc_in))

    # create insert cursor
    with arcpy.da.InsertCursor(fc_out, flds_use) as curs_out:

        # loop through input features
        with arcpy.da.SearchCursor(fc_in, flds_use) as curs_in:
            for row_in in curs_in:
                polygon = createRectangleFromPoint(row_in[0].firstPoint, width, height, sr)
                row_out = buildRow(row_in, flds_use, polygon)
                curs_out.insertRow(row_out)

def createRectangleFromPoint(pnt, width, height, sr):
    arrPnts = arcpy.Array()
    pnt2 = arcpy.Point(pnt.X - width, pnt.Y - height)
    arrPnts.add(pnt2)
    pnt2 = arcpy.Point(pnt.X - width, pnt.Y + height)
    arrPnts.add(pnt2)
    pnt2 = arcpy.Point(pnt.X + width, pnt.Y + height)
    arrPnts.add(pnt2)
    pnt2 = arcpy.Point(pnt.X + width, pnt.Y - height)
    arrPnts.add(pnt2)
    pnt2 = arcpy.Point(pnt.X - width, pnt.Y - height)
    arrPnts.add(pnt2)
    return arcpy.Polygon(arrPnts, sr)

def correctFieldList(flds):
    flds_use = ['Shape@']
    fldtypes_not = ['Geometry', 'Guid', 'OID']
    for fld in flds:
        if not fld.type in fldtypes_not:
            flds_use.append(fld.name)
    return flds_use

def buildRow(row_in, flds_use, polygon):
    try:
        lst_in = list(row_in)
        lst_in[0] = polygon
        return tuple(lst_in)
    except Exception, e:
        print("buildRow exception: {0}".format(e))
        return None

if __name__ == '__main__':
    main()

All comments

AndrewL

Thanks for your help! I used "width" twice because I am just creating squares.

XanderBakker

Hi Andrew Louchios ,

You cannot define the width like you did and as you can notice in the error you cannot mix numeric and text values. You will probably have to do something like this:

    fld_width <SPAN class="operator token">=</SPAN> <SPAN class="string token">"name of your width field"</SPAN>
    fld_height <SPAN class="operator token">=</SPAN> <SPAN class="string token">"name of your height field"</SPAN>

    <SPAN class="comment token"># ...</SPAN>

    <SPAN class="keyword token">with</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>fc_in<SPAN class="punctuation token">,</SPAN> <SPAN class="punctuation token">(</SPAN>flds_use<SPAN class="punctuation token">,</SPAN> fld_width<SPAN class="punctuation token">,</SPAN> fld_height<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">as</SPAN> curs_in<SPAN class="punctuation token">:</SPAN>
        <SPAN class="keyword token">for</SPAN> row_in <SPAN class="keyword token">in</SPAN> curs_in<SPAN class="punctuation token">:</SPAN>
            width <SPAN class="operator token">=</SPAN> row_in<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN>
            height <SPAN class="operator token">=</SPAN> row_in<SPAN class="punctuation token">[</SPAN><SPAN class="number token">2</SPAN><SPAN class="punctuation token">]</SPAN>
            polygon <SPAN class="operator token">=</SPAN> createRectangleFromPoint<SPAN class="punctuation token">(</SPAN>row_in<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">.</SPAN>firstPoint<SPAN class="punctuation token">,</SPAN> width<SPAN class="punctuation token">,</SPAN> height<SPAN class="punctuation token">,</SPAN> sr<SPAN class="punctuation token">)</SPAN>
<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

In the beginning you will have to define the names of your width and height field and in the search cursor you will have to include the field names in the list of fields. Inside the loop you will have to read the values and use them in the CreateRectangleFromPoint function.

AndrewL

Is it possible to use this code using a width and height attribute field in the input file for different widths and heights? Thank you.

I tried the following but get an error: pnt2 = arcpy.Point(pnt.X - width, pnt.Y - height)
TypeError: unsupported operand type(s) for -: 'float' and 'str'

width = "float(!newWidth!)"

height = "float(!newHeight!)"

XanderBakker

... and if you want to create an output that contains the attributes of the input, you can try this:

import arcpy, os

def main():
    arcpy.env.overwriteOutput = True

    # inputs
    fc_in = r'D:\Xander\EPM\Datos pruebas\gdb\datos erase.gdb\points03'
    width = 0.005 # with of the output rectangle
    height = 0.005 # height of the output rectangle

    # output
    fc_out = r'D:\Xander\EPM\Datos pruebas\gdb\datos erase.gdb\rectangles3'

    # create empty output featureclass
    sr = arcpy.Describe(fc_in).spatialReference
    out_ws, out_name = os.path.split(fc_out)
    arcpy.CreateFeatureclass_management(out_ws, out_name, "POLYGON", fc_in, "DISABLED", "DISABLED", sr)

    # make fields list
    flds_use = correctFieldList(arcpy.ListFields(fc_in))

    # create insert cursor
    with arcpy.da.InsertCursor(fc_out, flds_use) as curs_out:

        # loop through input features
        with arcpy.da.SearchCursor(fc_in, flds_use) as curs_in:
            for row_in in curs_in:
                polygon = createRectangleFromPoint(row_in[0].firstPoint, width, height, sr)
                row_out = buildRow(row_in, flds_use, polygon)
                curs_out.insertRow(row_out)

def createRectangleFromPoint(pnt, width, height, sr):
    arrPnts = arcpy.Array()
    pnt2 = arcpy.Point(pnt.X - width, pnt.Y - height)
    arrPnts.add(pnt2)
    pnt2 = arcpy.Point(pnt.X - width, pnt.Y + height)
    arrPnts.add(pnt2)
    pnt2 = arcpy.Point(pnt.X + width, pnt.Y + height)
    arrPnts.add(pnt2)
    pnt2 = arcpy.Point(pnt.X + width, pnt.Y - height)
    arrPnts.add(pnt2)
    pnt2 = arcpy.Point(pnt.X - width, pnt.Y - height)
    arrPnts.add(pnt2)
    return arcpy.Polygon(arrPnts, sr)

def correctFieldList(flds):
    flds_use = ['Shape@']
    fldtypes_not = ['Geometry', 'Guid', 'OID']
    for fld in flds:
        if not fld.type in fldtypes_not:
            flds_use.append(fld.name)
    return flds_use

def buildRow(row_in, flds_use, polygon):
    try:
        lst_in = list(row_in)
        lst_in[0] = polygon
        return tuple(lst_in)
    except Exception, e:
        print("buildRow exception: {0}".format(e))
        return None

if __name__ == '__main__':
    main()

XanderBakker

Below the code based on the thread mentioned by Dan Patterson‌ adapted to create squares or rectangles. Change lines 6 and 11 to point to your input and output (will be created) and change the sizes on lines 7 and 8. For some more snippets with Python have a look at: https://community.esri.com/docs/DOC-1927

import arcpy
def main():
    arcpy.env.overwriteOutput = True

    # inputs
    fc = r'D:\Xander\Rectangles\data.gdb\points'
    width = 0.005 # with of the output rectangle
    height = 0.005 # height of the output rectangle

    # output
    fcout = r'D:\Xander\Rectangles\data.gdb\rectangles'

    sr = arcpy.Describe(fc).spatialReference
    features = []
    with arcpy.da.SearchCursor(fc, ("SHAPE@")) as cursor:
        for row in cursor:
            rectangle = createRectangleFromPoint(row[0].firstPoint, width, height, sr)
            features.append(rectangle)

    # write to output
    arcpy.CopyFeatures_management(features, fcout)

def createRectangleFromPoint(pnt, width, height, sr):
    arrPnts = arcpy.Array()
    pnt2 = arcpy.Point(pnt.X - width, pnt.Y - height)
    arrPnts.add(pnt2)
    pnt2 = arcpy.Point(pnt.X - width, pnt.Y + height)
    arrPnts.add(pnt2)
    pnt2 = arcpy.Point(pnt.X + width, pnt.Y + height)
    arrPnts.add(pnt2)
    pnt2 = arcpy.Point(pnt.X + width, pnt.Y - height)
    arrPnts.add(pnt2)
    pnt2 = arcpy.Point(pnt.X - width, pnt.Y - height)
    arrPnts.add(pnt2)
    return arcpy.Polygon(arrPnts, sr)

if __name__ == '__main__':
    main()

Venkata_RaoTammineni

Hi,

try this Buffers With Points - A Tool for Creating Rectangular Buffers | GISNuts.com

Thanks and Regards,

Venkat

DanPatterson_Retired

See Xander Bakker's diamond buffer script...a little bit of fiddling and you can turn a diamond into a square

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.