En Uzun Tarafa Yönlendirilmiş Bir Çokgenin İçinde Izgara Nasıl Oluşturulur?

Atanmış satır ve sütun öznitelikleri olan bir çokgen içinde bir ızgara oluşturmak için "Çokgen Graticule Oluştur" adlı bir işlem kullanıyorum. QGIS'te çalıştırdığım bir SAGA aracı. N ile S ve E ile W arasında çalışan hücre ile giriş çokgeni boyutuna göre bir ızgara oluşturur. Izgara yönünün çokgenin en uzun kenarına paralel olması analizim için daha iyi olur. Referans için geçerli otomatik çıktının görüntülerini ve ızgarayı çokgene hizalamak için manuel olarak değiştirdiğim bir çıktıyı ekledim. Izgaralar başlangıçta çokgenin tüm alanını kapladı, ancak yalnızca içerideki hücre parçalarını göstermek için kırpıldı.

Yukarıdaki resim, araç çıkışının varsayılan olarak koyduğu şeydir.

Bu, gerekli parametrelere uyacak şekilde manuel olarak düzenlenen bir ızgaranın görüntüsüdür.

Benim sorum, ikinci görüntüdeki gibi otomatikleştirilmiş ızgaralar üretmenin bir yolu var mı? Birkaç bin çokgen için ızgara oluşturmam gerekecek.

'Den bir vektör ızgarası oluşturursanız Vector -> Research Tools -> Vector Grid, dönüş açısını belirleyebilirsiniz.

QGIS 2.18

@Jbgramm @ neogeomat tarafından verilen cevaba yorum yaptığında, QGIS 2.18'de (2.j.1gram-17 @jbgramm tarafından kontrol edilen ve 2.18.18 ile aynı) Vector Grid aracında bir Rotation angleseçenek görmüyoruz .

Bu da bizi GRASS v.mkgrid işleviyle terk ediyor .

Açı kuzeyden saat yönünün tersine ölçülür. Yani -30resimde N30E'ye yönelik bir ızgara üretiliyor.

Ayrıca QGIS 3'te bu problemle karşılaştım ve bu çözümü yığın taşmasında buldum

Temel olarak fikir , ızgarayı oluşturmadan önce açının tanımlandığı çokgene uygulanmasıdır . Çokgeniniz bir dikdörtgen değilse, daha önce çokgeninizin boyutundan bir katman oluşturmanız ve ardından döndürmeniz gerekir. Ardından ızgarayı bu yeni dereceye göre oluşturabilir ve daha sonra çokgeninizi ve ızgarayı orijinal Çokgen kapsamına döndürebilirsiniz. Bütün bunlar , her iki katmanda da sabit x, y koordinatının bağlantı noktası olarak kullanıldığından emin olur .

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