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?

Encontre mais publicações marcadas com

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

Participe da Conversa

Membros conectados podem postar, seguir atualizações e mais. Novo aqui? Registre uma conta gratuita.

Links Rápidos

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