Calculate visibility of astronomical bodies using GDAL/PyEphem and GeoTIFF