Raster ayak izlerini gösteren shapefile oluşturmak?

Tiff formatta yaklaşık 1.000 uydu görüntüsüm var ve rasterlere bir indeks görevi yapacak bir shapefile oluşturmak istiyorum. Bu raster kataloguna benzer bir şey, ancak raster katalog oluşturmak istemiyorum.

Görebildiğim bazı engeller, görüntünün coğrafi olarak güçlendirilmiş olmasıdır, bu yüzden onlar şekil olarak dikdörtgen değillerdir (veri alanı hakkında konuşuyorum).

Açıklığa kavuşturmak için, poligonun yalnızca rasterin sıfır olmayan (veya nodata olmayan) piksellerini kaplamasını ve tüm dikdörtgen rasterleri içermesini istemiyorum. Şimdiye kadar verilen cevapların çoğu, veriyi ve veri olmayan pikselleri kapsayan dikdörtgen bir çokgen vermektedir.

Benim görselim

İnceledim araçları tarafından verilen sonuç (raster catlog, çeşitli Arcscripts, cevaplardan birinde verilen özel Python betiği gibi):

İstediğim sonuç:

QGIS'de Image Boundary ( Görüntü Sınırı) adlı bir eklenti var . Bu harika bir araçtır. Bu araç içerisinde, örneğin bir uydu görüntüsünün siyah kenarlarını atlayacak olan "Geçerli Piksel" seçeneği vardır.

Aşağıdaki kod bir giriş rasterini alır, boyutunu alır ve bu dereceyi çokgen bir özellik sınıfına ekler:

import arcpy

r = arcpy.Raster(in_raster)
point = arcpy.Point()
array = arcpy.Array()

corners = ["lowerLeft", "lowerRight", "upperRight", "upperLeft"]

cursor = arcpy.InsertCursor(fc)
feat = cursor.newRow()
for corner in corners:    
    point.X = getattr(r.extent, "%s" % corner).X
    point.Y = getattr(r.extent, "%s" % corner).Y
    array.add(point)
array.add(array.getObject(0))
print len(array)
polygon = arcpy.Polygon(array)
feat.shape = polygon
cursor.insertRow(feat)
array.removeAll()
del feat
del cursor

ArcMap Python penceresinde ayarlayarak in_rasterve fcbenzeri şekilde çalıştırabilirsiniz:

>>> fc = 'r_extent'
>>> in_raster = 'CaliCoast10mNED_HavCurvedPCS'

nerede r_extentbir varolan poligon featureclass olduğunu. Sonra kodu kopyala ve çalıştır. Bunu anladım:

Bunun için gdaltindex kullanabilirsiniz: http://www.gdal.org/gdaltindex.html

Bununla birlikte, görüntülerle aynı referans sisteminde hala dikdörtgenler (örneğin 4 + 1 puan) yaratacaktır. Ama bunun gerçekten bir sorun olup olmadığını merak ediyorum: resimleriniz ne kadar büyük?

QGIS’de “Image Boundary” aracını denedim, ancak bazı rasterlerle başarısız oluyor.

İşe yarayan bir diğer seçenek ise: - 1. gdal_translate komut satırı yardımcı programını kullanarak giriş rasterini tek bir dijital numaraya dönüştürün.

gdal_translate -co sıkıştır = lzw -b 1 -ot byte -scale 1 1 inp.tif

Bu, inp.tif in sıfır olmayan değerlere sahip her yerde 1 değerine sahip raster yaratacaktır.

2. Ardından ArcMap'te çokgen elde etmek için "çokgen raster" aracını kullanabilirsiniz (şekildosyası).

Çıktı şekil dosyası, görüntünün gerçek sınırını izleyen çoklu poligon köşeleriyle birlikte dörtgen olmayan sınırlara sahip olacaktır.

Komut dosyası ile yazılan bir çözümle ilgileniyorsanız, yakındaki gdal komutunu deneyin . Bayrak -beyazı da beyaz boşlukla sınırlanmış verileri almak için kullanabilirsiniz. Sen edebilirsiniz Çokgenselleştirmek kafes veri dikdörtgen olmayan ayak izi oluşturmak için çıkış.

##run nearblack
shellcmds = "nearblack -white -o " + outpath + " " + inpath
os.system(shellcmds)

##run polygonize
shellcmds = "gdal_polygonize.py " + outpath + " -f" + """ "ESRI Shapefile"  """ + outpath2 + " nearblack"
os.system(shellcmds)

ArcGIS dünyasında bunu 3D Analyst \ Conversion \ Raster \ Raster Domain aracını kullanarak elde edebileceğinizi düşünüyorum .

ArcGIS 3D Analyst'te Raster Domain adı verilen ve raster veri kümesinin veri bölümlerinin çokgen veya çok satırlı bir alan oluşturduğu ve NoData hücrelerini yok sayan bir araç var. Benim için harika çalıştı.

Uzamsal bir referans atanmış bir tarama kataloğu oluşturabilir, tarama kataloğunu ArcMap'e ekleyebilir ve ardından "Ayakizi Dışa Aktar" komutunu kullanabilirsiniz. Unutmayın "UNMANAGED" raster kataloğu oluşturma seçeneği, böylece her raster veri kümesinin kopyasını oluşturmazsınız.