To promote urban sustainability, many cities are adopting bicycle-friendly policies, leveraging GPS trajectories as a vital data source. However, the inherent errors in GPS data necessitate a critical preprocessing step known as map-matching. Due to GPS device malfunction, road network ambiguity for cyclists, and inaccuracies in publicly accessible streetmaps, existing map-matching methods face challenges in accurately selecting the best-mapped route. In urban settings, these challenges are exacerbated by high buildings, which tend to attenuate GPS accuracy, and by the increased complexity of the road network. To resolve this issue, this work introduces a map-matching method tailored for cycling travel data in urban areas. The approach introduces two main innovations: a reliable classification of road availability for cyclists, with a particular focus on the main road network, and an extended multi-objective map-matching scoring system. This system integrates penalty, geometric, topology, and temporal scores to optimize the selection of mapped road segments, collectively forming a complete route. Rotterdam, the second-largest city in the Netherlands, is selected as the case study city, and real-world data is used for method implementation and evaluation. Hundred trajectories were manually labelled to assess the model performance and its sensitivity to parameter settings, GPS sampling interval, and travel time. The method is able to unveil variations in cyclist travel behavior, providing municipalities with insights to optimize cycling infrastructure and improve traffic management, such as by identifying high-traffic areas for targeted infrastructure upgrades and optimizing traffic light settings based on cyclist waiting times.