Rastere büyük şekil dosyası

Büyük bir şekil dosyam var (1 gb) ve rasterleştirmem gerekiyor. Şunu zaten denedim.

1.) Bunu kullanarak v.in.ogrhata mesajıyla GRASS'a aktarın: ERROR: G_realloc: break_polygons.c adresinde 498240036 bayt ayrılamadı.: 188

2.) İkinci fikrim PostGIS kullanmaktı. Şekil dosyasını içe aktarın, ızgaranın x, y konumunda yeniden örnekleyin ve ardından bu noktaları dışa aktarın ve xyz'den bir ızgara oluşturun. Şekil dosyalarını (çokgenler ve noktalar) başarıyla içe aktardım, ancak 300k noktalı kesişen 1 milyon çokgen çok yavaş görünüyor. Aşağıdaki PostGIS doygunluğunu kullandım, belki iyileştirme için yer var.

select polygons.land_id,grid.geom from grid,polygons where grid.geom && polygons.geom and within(grid.geom,polygons.geom)

3.) simplify()PostGIS'te de kullanmaya çalıştım. Ama birçok küçük çokgene kaybettim (yani sadece küçük çokgenlerle kaplı bazı alanlar geçersiz hale geldi).

Herhangi bir fikir büyük mutluluk duyacağız.

Bu kadar büyük bir şekil dosyasıyla kullanmadığım halde gdal_rasterize'yi deneyebilirsiniz , bu yüzden GRASS ile aynı sorunlara sahip olabilirsiniz. Aşağıdaki gibi bir şey (GDAL> = 1.8.0 ile) çalışması gerektiğini düşünüyorum:

gdal_rasterize -a AN_ATTRIB -l THE_LAYER -a_nodata -9999 -a_srs EPSG:27700 -co TILED=YES -tr 10 10 -ot Float32 src.shp dest.tif

Tabii ki, kaynak şekil dosyanıza bağlı olarak bazı seçeneklerle oynamak zorunda kalacaksınız. En önemli parametre bir pikselin çözünürlüğünü belirten -tr'dir; onsuz, kendinizi çok büyük bir tarama ile bulabilirsiniz ...

GRASS ile yapışmak istiyorsanız, rasterleştirme için daha küçük bir boyut ayarlamayı deneyin ve işlemi yönetilebilir parçalara ayırın, ardından rasterleri bir tanesine mozaikleyin.

İlk niyetinizle ilgili olarak, daha fazla RAM veya takas alanı olan bir makinede v.in.ogr komutunu çalıştırmayı deneyebilir misiniz ?

Değilse, bir çok dosyaya bölebilirsiniz. Yeniden birleştirmeden önce her birini rasterleştirin.

Bir rasterleştirme yaparken daha fazla kontrole ihtiyacınız varsa, perrygeo'nun başlık altında GDAL kullanan, ancak örtüşen özellikler veya tek başına mümkün olanların ötesinde koşullu değerlendirmeler eklemek için kullanılabilen poly_density.py betiğine göz atın gdal_rasterize.

Karındeşen Jack gibi parçalara ayır. Postgis veri dilimleri ihracat çim gelen ve dönüştürmek.

İşte PostGIS'te bir vektör tablosunun rasterleştirilmesi hakkında bir rehber ...