Le portail francophone de la géomatique

Bonjour. Je cherche à satisfaire une demande : représenter les zones Lambert 1 à 4 dans un logiciel SIG bureautique (par ex QGIS), avec en fond de carte les données IGN GEOFLA, par exemple, afin de produire une belle carte au format A0 avec les principales villes, les départements, etc.. que je puisse placarder dans mon bureau.

Après pas mal de recherches, je n'ai trouvé que des cartes peu précises, à petite échelle, représentant sommairement les 4 zones, essentiellement pour un besoin visuel sommaire.
Donc si je ne trouve pas de carte "toute faite", je prend le problème "mathématiquement", et je cherche donc un moyen de représenter dans un SIG l'ensemble des points contenus dans une emprise délimitée par les latitudes N et S de chaque zone Lambert (lat = cte), et le contour de la France aux extrémités E et O, et d'exporter le polygone résultant en shape.
Mathématiquement, les délimitations N et S correspondent à un courbe paramétrée très simple : Lat = constante, qui correspond aux parallèles extrêmes de la zone) mais je ne vois pas comment exploiter "Lat= cte" dans un SIG (QGIS) (possibilité de représenter des points définis par une courbe paramétrée ?).

Ou sinon, je ne suis pas récalcitrant au Python ni à GDAL, donc si qqn a déjà planché là-dessus, je suis à l'écoute.

Ou alors, existe-t-il un tableau des villes Française, avec des données attributaires (Lat, Long) et l'information d'appartenance à la zone Lambert (1, 2, 3, 4) ?
Bref, la question me semblait banale au départ, mais je coince. 
Merci d'avance pour vos contributions.
Alexis GAILLOT, ERT TECHNOLOGIES, Géomaticien - projeteur FTTH