OGR provider doesn't work with usgs waterdata endpoint

[C-Loftus]

Putting this here just for posterity

      - type: feature
        name: OGR
        data:
          source_type: OAPIF
          source: OAPIF:https://api.waterdata.usgs.gov/ogcapi/v0/
        id_field: id
        layer: monitoring-locations

[2025-12-15T14:39:47Z] {/Users/cloftus/github/WWDH/.venv/lib/python3.12/site-packages/pygeoapi/provider/ogr.py:379} ERROR - 'NoneType' object has no attribute 'AssignSpatialReference'
[2025-12-15T14:39:47Z] {/Users/cloftus/github/WWDH/.venv/lib/python3.12/site-packages/pygeoapi/api/__init__.py:590} ERROR - generic error (check logs)
Traceback (most recent call last):
  File "/Users/cloftus/github/WWDH/.venv/lib/python3.12/site-packages/pygeoapi/provider/ogr.py", line 363, in query
    result = self._response_feature_collection(
             ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/Users/cloftus/github/WWDH/.venv/lib/python3.12/site-packages/pygeoapi/provider/ogr.py", line 577, in _response_feature_collection
    json_feature = self._ogr_feature_to_json(
                   ^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/Users/cloftus/github/WWDH/.venv/lib/python3.12/site-packages/pygeoapi/provider/ogr.py", line 526, in _ogr_feature_to_json
    geom.AssignSpatialReference(None)
    ^^^^^^^^^^^^^^^^^^^^^^^^^^^
AttributeError: 'NoneType' object has no attribute 'AssignSpatialReference'

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/Users/cloftus/github/WWDH/.venv/lib/python3.12/site-packages/pygeoapi/api/itemtypes.py", line 546, in get_collection_items
    content = p.query(offset=offset, limit=limit,
              ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/Users/cloftus/github/WWDH/.venv/lib/python3.12/site-packages/pygeoapi/provider/ogr.py", line 380, in query
    raise ProviderGenericError(err)
pygeoapi.provider.base.ProviderGenericError: 'NoneType' object has no attribute 'AssignSpatialReference'
INFO:     127.0.0.1:57016 - "GET /collections/usgs-monitoring-locations/items?f=jsonld HTTP/1.1" 500 Internal Server Error

[webb-ben]

This is a known issue to do with how GDAL determines geom and the initial features from USGS containing null geom. Believe they have implemented a hack that intercepts this for QGIS but would not be surprised if its the source of this issue as well.

        QGIS_detected = query_helper.QGIS_detected
        LOGGER.debug(f"{QGIS_detected=}")

        # QGIS makes an initial test query for 10 records with no bbox to see if there's geometry
        # if this is that query, then we need to force at least one non-null geometry to be there
        hack_for_qgis_enabled_for_this_query = (
            self.hack_for_qgis_enabled and QGIS_detected and offset == 0 and limit == 10 and not bbox
        )
        LOGGER.debug(f"{hack_for_qgis_enabled_for_this_query=}")

[webb-ben]

I think you may be able to add query arg _QGIS=true to force the behavior