def intersect(self, ra, dec, radius, inclusive=True, convertCoordinates=True):
<EOL>if convertCoordinates == True:<EOL><INDENT>converter = unit_conversion(<EOL>log=self.log<EOL>)<EOL>ra = converter.ra_sexegesimal_to_decimal(<EOL>ra=ra<EOL>)<EOL>dec = converter.dec_sexegesimal_to_decimal(<EOL>dec=dec<EOL>)<EOL><DEDENT>if inclusive:<EOL><INDENT>inc = <NUM_LIT:1><EOL><DEDENT>else:<EOL><INDENT>inc = <NUM_LIT:0><EOL><DEDENT>return super(HTM, self).intersect(ra, dec, radius, inc)<EOL>
*return IDs of all triangles contained within and/or intersecting a circle centered on a given ra and dec*

        **Key Arguments:**
            - ``ra`` -- RA of central point in decimal degrees or sexagesimal
            - ``dec`` -- DEC of central point in decimal degrees or sexagesimal
            - ``radius`` -- radius of circle in degrees
            - ``inclusive`` -- include IDs of triangles that intersect the circle as well as those completely inclosed by the circle. Default *True*
            - ``convertCoordinates`` -- convert the corrdinates passed to intersect. Default *True*
            -

        **Return:**
            - ``trixelArray`` -- a numpy array of the match trixel IDs

        **Usage:**

            To return the trixels overlapping a circle with a 10 arcsec radius centred at 23:25:53.56, +26:54:23.9

            .. code-block:: python

                overlappingTrixels = mesh16.intersect(
                    ra="23:25:53.56",
                    dec="+26:54:23.9",
                    radius=10 / (60 * 60),
                    inclusive=True
                )

            Or to return the trixels completing enclosed by a circle with a 1 degree radius centred at 23:25:53.56, +26:54:23.9

            .. code-block:: python

                overlappingTrixels = mesh16.intersect(
                    ra="23:25:53.56",
                    dec="+26:54:23.9",
                    radius=1,
                    inclusive=False
                )