ラスターのフットプリントを示すシェープファイルを作成しますか？

TIFF形式で約1,000枚の衛星画像があり、ラスタのインデックスとして機能するシェープファイルを作成したいと思います。これはラスターカタログに似ていますが、ラスターカタログを作成したくありません。

私が予見できるいくつかの障害は、画像がジオリフレンシー化されているため、形状が長方形ではないことです（データ領域について話している）。

明確にするために、ポリゴンは、長方形のラスタ全体ではなく、ラスタの非ゼロ（または非データ）ピクセルのみをカバーする必要があります。これまでの回答のほとんどは、データと非データピクセルをカバーする長方形のポリゴンを提供します。

私のイメージ

私が調べたツールによって得られた結果（ラスターcatlog、さまざまなArcscript、回答の1つで与えられたカスタムPythonスクリプトなど）：

私が望む結果：

QGISにはImage Boundaryと呼ばれるプラグインがあります。それは素晴らしいツールです。このツールには、たとえば、衛星画像の黒い縁を省く「有効なピクセル」のオプションがあります。

次のコードは、入力ラスターを取得し、その範囲を取得して、その範囲をポリゴンフィーチャクラスに挿入します。

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ウィンドウでそれを実行することができますin_rasterし、fcそうのように：

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

どこr_extent既存のポリゴンフィーチャクラスです。次に、コードをコピーして実行します。私はこれを得る：

これにはgdaltindexを使用できます：http ://www.gdal.org/gdaltindex.html

ただし、画像と同じ参照システムに四角形（4 + 1ポイントなど）を作成します。しかし、それは本当に問題なのだろうか：あなたの画像はどれくらい大きいのか？

QGISで「画像境界」ツールを試しましたが、一部のラスタで失敗します。

私のために働いた別のオプションは次のとおりです。

gdal_translate -co compress = lzw -b 1 -ot byte -scale 1 1 inp.tif out.tif

これにより、inp.tifの値がゼロ以外のすべての場所で値1のラスタが作成されます。

2. 次に、ArcMapで「ラスタからポリゴン」ツールを使用して、ポリゴン（シェープファイル）を取得できます。

出力シェープファイルには、画像の実際の境界に沿って、複数の多角形の頂点を持つ非四辺形の境界があります。

スクリプト化されたソリューションに興味がある場合は、gdalコマンドnearblackを試してください。フラグ-whiteを使用して、空白で区切られたデータも取得できます。出力をポリゴン化して、ラスタデータの非長方形のフットプリントを生成できます。

##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の世界では、3D Analyst \ Conversion \ From Raster \ Raster Domainツールで実現できると思います。

ArcGIS 3D Analystには、ラスタデータセットのデータ部分のポリゴンまたはポリラインフットプリントを作成し、NoDataセルを無視するラスタドメインと呼ばれるツールがあります。私にとっては素晴らしい仕事でした。

空間参照が割り当てられたラスタカタログを作成し、ラスタカタログをArcMapに追加してから、「フットプリントのエクスポート」を実行できます。各ラスタデータセットの複製を作成しないように、「管理されていない」ラスタカタログを作成するオプションがあることを忘れないでください。