NDVI（MOD13Q1）MODISシリーズの5：NDVI処理フロー

序文：（個人的なアドバイス）

1.通常はデータ量が比較的多いため、データ処理作業を実行します。そのため、パソコンでは、基本データ、成長過程データ、結果データを分類し、見つけやすくする必要があります。

2.重要なデータをバックアップします。データ処理の過程で、間違いによって重要なデータが失われるため、重要なデータが絶版にならないようにしてください。したがって、圧縮ファイルを生成または保存することをお勧めします。

3.データ量が多く、バッチ処理が必要な場合。テンプレートドキュメントまたはコードは、詳細エラーの発生を防止する必要があります。完了するのは簡単で、時間がかかり、非効率的です。

処理ツール：ArcMap

処理データ：MOD13Q1のNDVI完成品（.tif）ファイル

調査地域のベクターファイル（河南省を例にとります）

処理フロー：

1.スタディエリアベクトルファイルをインポートします（中空処理を行います）。

2. .tifファイルをArcMapにインポートします（主にコンピューターの耐性に基づいて、マルチリードにすることができます）。

ピラミッドを構築するかどうか：はい

3.バッチ定義プロジェクション（MRTはバッチ処理でそれを実行しましたが、二重の保険であり、このステップは簡単です）（ArcMapには多数または単一または複数のバッチプロジェクションメソッドがあり、この記事では1つだけ紹介します）

レイヤーをダブルクリック

座標系をクリックします

投影座標系をクリック

UTMをクリックします

WGS1984をクリックします。

北半球

この記事では、WGS 1984 UTM ZONE 49Nを使用しています。

[適用]をクリックしてから[OK]をクリックします。（よく使用される座標がある場合は、右クリックして座標をクリックし、フォルダーのお気に入りを選択します）

4.マスクごとに抽出（クロップも機能しますが、マスクごとに抽出することをお勧めします）（このステップは1つずつしか処理できません。Pythonバッチ処理はまだ更新されていません）

次の図に示すように、ArcToolboxをクリックしてArcToolboxを呼び出し、マスクで抽出します。ダブルクリックしてマスクで抽出します。

OKをクリック

便宜上、元の.tifファイルは削除できます。結果は次のとおりです。

このステップでコードが必要な場合、コードメソッドは次のとおりです。

同じマスクファイルをバッチカットします。

インポートArcPy
 から ArcPy をインポートENV
 から arcpy.sa インポート * ＃設定環境 
arcpy.env.workspace = R " ：FAN InterpolationlKrigO E " 
ラスタの。= ArcPy ListRasters（" * "、" GRID "）＃GETファイルリスト＃の出力パス 
outRasterPath = " E：EANInterpolationKrigO_ extractlKrigO_。" ras in rasの場合：
     print outRasterPathtras ＃必要なライセンスを確認します







    
    Arcpy。CheckOutExtension（" Spatial " ）   
    
    ＃ローカル変数
    inRaster ras 
    inMaskData = " E：FUJIANModelBoundary.shp "  ＃マスク 
    膜文件outRaster-outRasterPath-ras + " .tif " 
    
    ＃プロセス：マスクによる抽出
    arcpy.gp。ExtractByMask_ sa（inRaster、inMaskData、outRaster）

print（" OK "）