Using Stations and Provincies data sets from PDOK, creating layers for them in QGIS, creating a virtual layer from Provincies (where naam ="Utrecht"), clipping the Stations Layer using the virtual layer for Utrecht and exporting the clipped layer - with CRS set to Mercator - EPSG:3857 - to GeoJSON