Làm cách nào để tạo lưới bên trong đa giác được định hướng với cạnh dài nhất?

Tôi đã và đang sử dụng một quá trình gọi là "Tạo đa giác Graticule" để tạo lưới trong một đa giác với các thuộc tính hàng và cột được gán. Nó là một công cụ SAGA mà tôi đã chạy trong QGIS. Nó tạo ra một lưới dựa trên phạm vi của đa giác đầu vào với ô chạy N đến S và E đến W. Sẽ tốt hơn cho phân tích của tôi nếu hướng của lưới song song với cạnh dài nhất của đa giác. Để tham khảo tôi đã đính kèm hình ảnh của đầu ra tự động hiện tại và đầu ra mà tôi đã thay đổi bằng tay để căn chỉnh lưới với đa giác. Các lưới ban đầu bao phủ toàn bộ phạm vi của đa giác nhưng đã được cắt bớt để chỉ hiển thị các phần tế bào bên trong.

Hình trên là những gì công cụ đặt ra theo mặc định.

Đây là hình ảnh của lưới được chỉnh sửa thủ công để phù hợp với các tham số cần thiết.

Câu hỏi của tôi là, có cách nào để tạo ra các lưới được căn chỉnh như cái trong hình ảnh thứ hai theo kiểu tự động không? Tôi sẽ cần tạo lưới cho vài nghìn đa giác.

Nếu bạn tạo lưới vector từ Vector -> Research Tools -> Vector Grid, bạn có thể chỉ định góc xoay.

Quốc tế 2.18

Như @jbgramm đã nhận xét câu trả lời được cung cấp bởi @neogeomat, chúng tôi không thấy Rotation angletùy chọn trong công cụ Vector Grid trong QGIS 2.18 (2.18.15-17 được kiểm tra bởi @jbgramm, và tương tự với 2.18.18).

Điều đó cho chúng ta chức năng GRASS v.mkgrid .

Góc được đo từ phía bắc, ngược chiều kim đồng hồ. Vì vậy, -30trong hình tạo ra một lưới định hướng đến N30E.

Tôi cũng gặp vấn đề này trong QGIS 3 và tìm thấy giải pháp này trong stack stack

Về cơ bản, ý tưởng sẽ là áp dụng góc trên đa giác mà phạm vi được xác định trước khi tạo lưới. Nếu đa giác của bạn không phải là hình chữ nhật, bạn sẽ cần tạo một lớp từ phạm vi đa giác của bạn trước đó, sau đó xoay nó. Sau đó, bạn có thể tạo lưới theo phạm vi mới này và sau đó xoay đa giác của bạn và lưới trở lại phạm vi Đa giác ban đầu. Tất cả điều này trong khi đảm bảo cùng tọa độ x, y được sử dụng làm điểm neo trong cả hai lớp.

#Define extent of Polygon
ext = QgsVectorLayer('path_to_polygon.shp', '', 'ogr' ).extent()
xmin = ext.xMinimum()
xmax = ext.xMaximum()
ymin = ext.yMinimum()
ymax = ext.yMaximum()
coords = "%f,%f,%f,%f" %(xmin, xmax, ymin, ymax)

#Define The angle of rotation. Change value to yours
azimut = 70.043

#define anchor point for rotation
anchor = "%f, %f" % (xmin, ymax)

#define x and y spacing of grid. Update to your desired spacing.
x = 3
y = 6

#create new polygon from extent
processing.run("native:extenttolayer", {'INPUT':coords,'OUTPUT':'Path_to_Output.shp'})

#Rotate Extent
processing.run("native:rotatefeatures", {'INPUT': 'Path_to_extent_Polygon.shp','ANGLE': azimut,'ANCHOR':anchor + '[EPSG:4326]','OUTPUT': 'Path_to_rotated_extent.shp'})

#Define extent of Rotated Polygon 
ext1 = QgsVectorLayer('Path_to_Rotated_Extent.shp', '', 'ogr' ).extent()
xmin1 = ext1.xMinimum()
xmax1 = ext1.xMaximum()
ymin1 = ext1.yMinimum()
ymax1 = ext1.yMaximum()
coords1 = "%f,%f,%f,%f" %(xmin1, xmax1, ymin1, ymax1)

#Create grid 
processing.run("qgis:creategrid", {'TYPE':0,'EXTENT': coords1 +'[EPSG:4326]','HSPACING':x,'VSPACING':y,'HOVERLAY':0,'VOVERLAY':0,'CRS':'EPSG:4326','OUTPUT': 'Path_to_grid.shp'})

#Rotate Grid to original extent
processing.run("native:rotatefeatures", {'INPUT': 'path_to_grid.shp','ANGLE': - 
azimut,'ANCHOR':rotate + '[EPSG:4326]','OUTPUT': 'path_to_rotated_grid.shp'})

# Clip Grid to Original Polygon
processing.run("native:clip", {'INPUT':'path_to_rotated_grid.shp','OVERLAY': 
'path_to_original_Polygon.shp','OUTPUT':'path_to_final_grid.shp'})