Saya memiliki koordinat berikut
minx, maxx, miny ,maxy = 448262.080078, 450360.750122, 6262492.020081, 6262938.950073
Saya ingin membuat kotak persegi ukuran 1 m menggunakan python.
import math
minx,maxx,miny,maxy = 448262.080078, 450360.750122, 6262492.020081, 6262938.950073
size = 1
def set_bbox(minx, maxx, miny, maxy, distx, disty):
nx = int(math.ceil(abs(maxx - minx)/distx))
ny = int(math.ceil(abs(maxy - miny)/disty))
new_maxx = minx + (nx*distx)
new_miny = maxy - (ny*disty)
return ((minx, new_maxx, new_miny, maxy),ny,nx)
# shift the bottom (right - down)
coord, ny, nx = set_bbox(minx,maxx,miny,maxy,size,size)
# left-up origin
origin = coord[0],coord[3]
# number of tiles
ncell = ny*nx
python
polygon
vector-grid
Gianni
sumber
sumber
Jawaban:
Script berikut akan melakukan pekerjaan dengan GDAL dan Python:
sumber
Skrip Python ini menggunakan pyshp library, seperti yang disarankan oleh user16044:
Catatan: kisi persegi ukuran 1 m dengan tingkat yang sama dengan lapisan yang berisi sekitar 1 juta poligon dan karenanya kinerja skrip berkurang secara masuk akal.
sumber
Pertanyaan ini dijawab sejak waktu yang lalu, tetapi saya menambahkan solusi lain menggunakan pustaka rupawan dan fiona:
sumber
Jawaban untuk Membuat Shapnet jala jaring di QGIS? menunjukkan opsi buat kisi di kotak alat pemrosesan QGIS.
sumber