多边形 shapefile 中的垂直线

Question

我正在使用多边形 shapefile 的单个输入来解决 GIS 问题。

考虑一个不规则的多边形。我想以相等的间距在多边形的范围内绘制垂直线。

我打算如何继续：

1. 识别边界框（使用 PyShp 完成）

2. 以等间距绘制平行于边界框左边缘的垂直线（如何？）

3. 将线条剪裁到多边形的范围（如何，不使用 ArcPy？）

注意：它们只能是垂直的，并且不是刻度线。另外，我不打算使用 ArcPy，而是打算在 Python (2.7) 中完成编码，因为这部分代码需要进入从 PyQt 生成的工具中。

Answer 1

• 在左边缘找到描述直线的点/顶点 - (x1, y1), (x2, y2)
• 为 x 值添加一个常数 - (x1+k, y1), (x2+k, y2)
• 在新的 x 值 - (x1+k, y3), (x2+k, y4) 处找到多边形上的 y 值
• 在这两点之间画线。

Answer 2

终于找到了我的问题的代码.. !!因此回答它...谢谢您的投入..

Ipath = raw_input("Enter the input file :- ")
Opath = raw_input("Enter the output directory :- ")
Ipath = Ipath.replace("\\", "/") # Python requirement for paths
Opath = Opath.replace("\\", "/")
copyfile(str(Ipath) + ".prj", str(Opath) + "/" + "Out_Lines" + ".prj") # Copying projection file


sf = shapefile.Reader(str('Input Path'))
shapes = sf.shapes()
Box = shapes[0].bbox
Spc = input("Enter the grid spacing :- ") # Grid Spacing read

x_min = Box[0] # Save the coordinates of the right-bottom, left-top bounding box
y_min = Box[1]
x_max = Box[2]
y_max = Box[3]

A_bbox = [x_min, y_min] # First assignment of coordinates
B_bbox = [x_max, y_max]
C_bbox = [x_min, y_max]
D_bbox = [x_max, y_min]

w = shapefile.Writer(shapefile.POLYLINE) # Shapefile writer
w.line(parts = [[A_bbox, C_bbox]])
w.field('Path number', 'C', '50') 
w.record(str(1)) # Writes the first line, that is the left 'side' of the bounding box

# Increasing the X coordinate to generate a line at a specified spacing 
i = 2
while (A_bbox[0] <= x_max):
    A_bbox = [A_bbox[0] + Spc, A_bbox[1]]
    C_bbox = [C_bbox[0] + Spc, C_bbox[1]]
    w.line(parts = [[A_bbox, C_bbox]])
    w.record(str(i))
    i = i+1

w.save(str(Opath) + "/" + "Out_Lines")

这会将结果保存在 shapefile 中。

作为上述问题的延续，问题的解决方案可在Clipping Line shapefiles within extent of Polygon shape 获得。我认为这组问题现在可以被视为已回答和结束。

感谢大家的帮助。