You have the right idea, you basically just need to transform the integration domain into a right triangle and use the change of variables theorem. If the triangle is composed of the segments ax+by=1, from (x
1, y
1) to (x
2, y
2) and cx+dy=1 from (x
2, y
2) to (x
3, y
3) (the third segment doesn't matter), then you get something like
integral[f(x,y) on triangle] = integral[|ad-bc| f(ax+by, cx+dy), {x, h(x
1, y
1), h(x
2, y
2)}, {y, h(x
2, y
2), h(x
3, y
3)} ]
where h(x,y)=(dx-by,-cx+ay)/(ad-bc) (just the inverse).
well, I might have missed some constant or something (can't be bothered to check everything

), but you get the idea.
[edit] Actually, on second thought a more straightforward approach might be to take the integral over the smallest square containing the triangle, and then just subtract out the three extra right triangle integrals to get the original triangle.