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

Point in Polygon Query

RichardAbram

Can the File Geodatabase API used to do a point in a polygon spatial query? If so does anyone have sample code?

Finde mehr Beiträge mit dem Tag

Accepted answers

VinceAngelo

The table::search method takes an envelope -- have you tried using the Querying
sample with an envelope where minx==maxx and miny==maxy? A pair of yin-
yang like shapes will tell the tale.

Worst case, you can make the envelope very small, and use a point-in-poly algorithm
to confirm that the center point is within the geometry (count the number of times a
line extended from the point to a point past the eastern [or western] edge of the shape
envelope crosses the segments [treating the second point as if it isn't part of the line,
so it isn't counted twice] -- odd crossing is within; even outside -- 20+ years ago I
solved the segments in the form Ax + By = C and used matrix algebra to locate the
point of intersection).

- V

All comments

EricONeal

Attached is a small file geodatabase with one featureclass. When I query with the point shown on the attached map(expanded slightly to make an envelope), I get three results. The correct one, David James, but also Tandy and Hamilton, which is obviously wrong. There should be only one row returned, but I am getting three.

Thanks

Code:

  static void Main(string[] args)
        {
            Geodatabase geodatabase = null;
            Table table = null;
            try
            {
                geodatabase = Geodatabase.Open(@C:\Temp\test.gdb);
                Esri.FileGDB.Point pt = new Esri.FileGDB.Point();

                double x = 1204123.4375;
                double y = 273148.87152777775;

                pt.x = x;
                pt.y = y;

                Envelope envelope = new Envelope();
                envelope.xMin = x - .00001;
                envelope.yMin = y - .00001;
                envelope.xMax = x + .00001;
                envelope.yMax = y + .00001;


                table = geodatabase.OpenTable("MetroCouncil");
                RowCollection spQueryRows = table.Search("*", "", envelope, RowInstance.Recycle);


                foreach (Row spQueryRow in spQueryRows)
                {
                    Console.WriteLine(spQueryRow.GetString("COUN_NAME").ToString());
                }

                Console.ReadLine();

            }

            catch (Exception ex)
            {
            }


            finally
            {
                table.Close();
                geodatabase.Close();
            }

        }

test.gdb.zip

CouncilPoint.jpg

EricONeal

THanks, I will provide a data and code sample in the morning when I am back at my PC.

VinceAngelo

It's only a bug if the three returned shapes do not overlap the requested envelope. Do they? Can you provide a FGDB data sample and the corners of the envelope with which you're querying?

ALL application code must handle the possibility that a point-in-poly query could return zero, one, or many features. With a simple envelope search, the probability of multiple return rows approaches unity.

Since there's no way to alter the design without including a full geometry processing library in the FGDBAPI DLL (and that doesn't seem likely), the best possible solution is going to be a simple helper function to test point-in-polygon (and even that will need to return a ternary response -- {"IN", "OUT", "ON"}).

- V

EricONeal

I am having the same issue. I have a point that I have converted to an extremely small envelope, but I get get 3 results back from my query and should only be getting one. This seems like a bug.

Anyone from ESRI reading these?

RichardAbram

Thanks for your response. I was hoping the search would return either 1 feature if it is inside it or 0 features if the point is not in the polygon.

When I do a table search (I am using the C# wrapper) setting the envelop min/max the same values as a point or even expanding the point to a small area I still get multiple results from the search.

I did find code for a winding number algorithm to detect whether the point is in or outside of a polygon. I reworked the c# code to deal with a MultipartShapeBuffer with parts and points. I now make the table search and then loop the results, usually 2 or 3 rows, to find which of the returned rows is the polygon I'm looking for.

VinceAngelo

Nimm an der Unterhaltung teil

Angemeldete Mitglieder können Beiträge verfassen, Updates folgen und mehr. Neu hier? Registriere ein kostenloses Konto.

Schnellzugriffe

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