A road network consists of multiple objects that deteriorate over time with different speeds of deterioration. In order to provide an adequate level of service over time, these objects will eventually require interventions. As road managers are trying, in general, to maximize the benefit obtained from the road network, it is in their interest to determine intervention programs, which consist of the grouping of interventions in work zones. The determination of optimal intervention programs is relatively complicated when considering single objects, but it becomes even more so when considering multiple objects embedded within a network. The objects to be included in the work zones at each time interval depend on many factors, such as the interventions to be executed on the objects, the maximum allowable length of the work zones, the traffic configurations to be used in the work zones and the available financial resources. Although some initial research in this area has been conducted, none has determined the optimal set of work zones on large infrastructure networks in a geographical information system (GIS) framework, which is necessary in the world of modern infrastructure management. In the work presented in this paper, a GIS-based program was developed to determine optimal intervention programs for large infrastructure networks using a linear optimization model, which can be linked directly to a GIS. The model includes constraints on the amount of available resources, on the length of the work zone, and on the distance between two work zones. A constraint-constructing algorithm is used in order to set up the latter two constraints. The program is illustrated by determining the optimal set of work zones for an example road network similar to the one in the canton of Wallis, Switzerland, including more than 2,000 bridges, tunnels, and road sections.