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

Calculating rectangular extends for polygon grids

WadeGivens

Looking for a little from the guru's out there on a task I'm trying to accomplish. I have a vector file I created using the fishnet tool that consists of hundreds of rectagular "grids". I'd like to loop through the records and find the uuper and lower Y, and left and right X values for each grid and write them to the appropriate new columns in the dataset (UpperY, LowerY, LeftX, RightX)

Any ideas on how to go about doing this?

Thanks again,

Najděte více příspěvků označených jako

Accepted answers

jskinner

Here is an example that may help:

fc = "fishnet"  rows = arcpy.UpdateCursor(fc) for row in rows:     geom = row.shape     row.LowerX = geom.extent.XMin     row.UpperX = geom.extent.XMax     row.LowerY = geom.extent.YMin     row.UpperY = geom.extent.YMax     rows.updateRow(row)  del row, rows

All comments

VAMSHIKRISHNA

Hi,

when i tried to give the field calculator expression in Quantum GIS it is throwing error like "An error occured while evaluating the calculation string"..Please provide me solution.

Thanks in advance,
Vamshi

Hi Wade,

one way you could go is with field calculator:
1. create fields for centroidX/centroidY of the grids, calculate geometry for them.
2. create fields for UpperY, LowerY, LeftX, RightX

in field calculator you can then calculate each field:
UpperY = centroidY + (gridheight*0.5)
LowerY = centroidY - (gridheight*0.5)
LeftX = centroidX - (gridwidth*0.5)
RightX = centroidX + (gridwidth*0.5)

WadeGivens

Here is an example that may help:

fc = "fishnet"

rows = arcpy.UpdateCursor(fc)
for row in rows:
    geom = row.shape
    row.LowerX = geom.extent.XMin
    row.UpperX = geom.extent.XMax
    row.LowerY = geom.extent.YMin
    row.UpperY = geom.extent.YMax
    rows.updateRow(row)

del row, rows

Thanks a bunch! Worked like a charm!

RaphaelR

Going to try the python code above

that´s what i´d do too (if i had known it was that straight forward, thanks Jake!) 🙂

WadeGivens

Hi Wade,

one way you could go is with field calculator:
1. create fields for centroidX/centroidY of the grids, calculate geometry for them.
2. create fields for UpperY, LowerY, LeftX, RightX

in field calculator you can then calculate each field:
UpperY = centroidY + (gridheight*0.5)
LowerY = centroidY - (gridheight*0.5)
LeftX = centroidX - (gridwidth*0.5)
RightX = centroidX + (gridwidth*0.5)

Thought about doing something similar to this after the first reply, but a little worried as to what the results would be if I want them in decimal degrees. Going to try the python code above first after my computer finishes up its processing task.

Thanks!

RaphaelR

Hi Wade,

one way you could go is with field calculator:
1. create fields for centroidX/centroidY of the grids, calculate geometry for them.
2. create fields for UpperY, LowerY, LeftX, RightX

in field calculator you can then calculate each field:
UpperY = centroidY + (gridheight*0.5)
LowerY = centroidY - (gridheight*0.5)
LeftX = centroidX - (gridwidth*0.5)
RightX = centroidX + (gridwidth*0.5)

jskinner

Here is an example that may help:

fc = "fishnet"  rows = arcpy.UpdateCursor(fc) for row in rows:     geom = row.shape     row.LowerX = geom.extent.XMin     row.UpperX = geom.extent.XMax     row.LowerY = geom.extent.YMin     row.UpperY = geom.extent.YMax     rows.updateRow(row)  del row, rows

DarrenWiens2

Hmmm... You can get X and Y coordinates of the rectangle centroid from Calculate Geometry (right-click column, choose Calculate Geometry). You can get width and height value from the Minimum Bounding Geometry tool (group by Object ID, envelope shape). Join output back to fishnet. Then, you should be able to subtract or add half the width and height to get your desired values. Warning: totally untested.

edit: I just realized this is in the Python forum. I guarantee there is a Python solution, not that I know it off hand. Check out Reading Geometry. Also, looks relatively simple to get to an Extent object which has the values you're looking for (see the example).

Připojte se k diskusi

Přihlášení členové mohou přispívat, sledovat aktualizace a další. Jste tu noví? Zaregistrujte si bezplatný účet.

Rychlé odkazy

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