矢量GIS数据视线？

Question

我有一种情况，我只关心几个矢量层和二维视线。我知道视线通常是在栅格数据上执行的，因为典型的用途是地形。因为这个措辞含糊其辞，而且几乎毫无用处，下面是我的情况：

我有一个多边形海岸线矢量shapefile，一个放置在水中某处的“源”点，以及一个表示围绕“源”点的大半径圆的“缓冲”多边形层。我只对缓冲区多边形中“看得见”源点的部分感兴趣。在下面的图像中，红色的点是源，橙色的多边形是用海岸线裁剪的缓冲区，黄色的多边形是我感兴趣的。即使这样也不像我希望的那样好。

图片来源：http://i.stack.imgur.com/IKBLv.png

我想自动化我现在使用的过程(相当耗时)，并且更喜欢使用python/numpy/scipy/OGR/GRASS而不是ESRI的东西。

你知道如何沿着这条线追踪并检查“可见性”吗？我可以光栅化所有的东西，并在GRASS中使用传统的径向视线脚本，但这似乎是一种处理内存中保存的太多数据并对像素运行检查的方法，我们知道这些像素不会在几个矢量的交点上产生冲突。我想尽可能的轻巧，同时保持尽可能高的精确度。

Answer 1

如何(迭代地)考虑点和海岸线shapefile中的每个点之间的线？如果它与“陆地”多边形相交(在陆地上交叉)，那么岸上的那个点是不可见的。获取所有可见的点，并使用它们形成可见区域的新多边形。