Rasterisation de nuages de points

Bonjour,
j'ai un problème que je n'arrive pas à résoudre.
Je possède deux nuages de points dont j'aimerais comparer leurs altitudes. Des outils existent mais ils interpoles et ça ne me convient pas.
Je souhaiterais donc les rasteriser, c'est à dire créer une grille dont chaque pixel représente un point. Il faut donc créer une sorte d'orthophoto dont chaque pixel possède les coordonnées des points correspondant, et la valeur du pixel correspond à la valeur d'altitude. Si aucun point ne correspond à un pixel, celui ci prend comme valeur 0.
J'arrive à le faire en passant par Qgis, mais j'aimerai automatiser le traitement soit uniquement avec Python, soit en insérant un scripts Python dans Qgis. Les deux rasters doivent avoir la même emprise. La suite du programme est déjà rédigée (comparaison des deux rasters).
Merci par avance pour votre aide.