Como criar uma grade dentro de um polígono orientado com o lado mais longo?

Eu tenho usado um processo chamado "Criar gratícula de polígono" para criar uma grade dentro de um polígono com atributos de linha e coluna atribuídos. É uma ferramenta SAGA que eu tenho executado no QGIS. Ele gera uma grade com base na extensão do polígono de entrada com a célula executando N a S e E a W. Seria melhor para minha análise se a direção da grade fosse paralela à borda mais longa do polígono. Para referência, anexei imagens da saída automatizada atual e uma saída que eu alterei manualmente para alinhar a grade com o polígono. As grades originalmente cobriam toda a extensão do polígono, mas foram cortadas para mostrar apenas as partes da célula que estão dentro.

A imagem acima é o que a ferramenta coloca por padrão.

Esta é uma imagem de uma grade que foi editada manualmente para ajustar os parâmetros necessários.

Minha pergunta é: existe uma maneira de produzir grades alinhadas, como a da segunda imagem, de maneira automatizada? Vou precisar gerar grades para alguns milhares de polígonos.

Se você criar uma grade de vetor Vector -> Research Tools -> Vector Grid, poderá especificar o ângulo de rotação.

QGIS 2.18

Como o @jbgramm comentou a resposta fornecida pelo @neogeomat, não vemos a Rotation angleopção Vector Grid na ferramenta QGIS 2.18 (2.18.15-17 verificada pelo @jbgramm e o mesmo com o 2.18.18).

Isso nos deixa com a função GRASS v.mkgrid .

O ângulo é medido a partir do norte, no sentido anti-horário. Portanto -30, a figura produz uma grade orientada para N30E.

Eu também encontrei esse problema no QGIS 3 e encontrei esta solução no estouro de pilha

Basicamente, a idéia seria aplicar o ângulo no polígono no qual a extensão é definida antes de criar a grade. Se o seu polígono não for um retângulo, você precisará criar uma camada a partir da extensão do seu polígono antes e girá-la. Em seguida, você pode criar a grade de acordo com essa nova extensão e, em seguida, gire o polígono e a grade de volta à extensão original do polígono. Tudo isso enquanto garante que a mesma coordenada x, y seja usada como ponto de ancoragem nas duas camadas.

#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'})