GSI_DEM
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
単語検索
|
最終更新
|
ヘルプ
]
開始行:
[[GIS/RS一般]]
*国土地理院のXML標高データをGeoTiffに変換する [#rf3221b1]
+ダウンロードした圧縮データを解凍し、変換するxmlファイル...
+株式会社エコリスの[["基盤地図情報標高DEM変換ツール">http...
# "変換結合.vbs"を開く~
# 投影法を "緯度経度"~
# xmlが格納されているフォルダを指定~
# Nodata(海の部分)は"-9999"または"0"と指定~
# ファイル名がメッシュコードのtif画像が出力される
+JGD2000緯度経度のままですべてのDEMをモザイク
--コマンドプロンプトを起動しカレントディレクトリをtif画像...
cd /d D:\DEM\GSI5m
--.tif拡張子のついたディレクトリ内のすべてのファイルのフ...
dir /b /s *.tif>list.txt~
--OSGeo4wシェルからgdal_merge.pyを使用してモザイク処理を...
# --optfile でlist.txtを読み込むことで、処理対象ファイル...
# -init でNodata(海の部分)は"-9999"または"0"と指定~
python C:\OSGeo4W\bin\gdal_merge.py -o D:/DEM/GSI5m/merge...
+投影変換~
-- gdalwarpを使用してJGD2000緯度経度からJGD2000UTM54帯の...
# -tr で解像度を指定（ここでは5m×5m）~
# -dstnodata でNodata値を"-9999"または"0"と指定 -s_srs...
gdalwarp -s_srs EPSG:4612 -t_srs EPSG:3100 -dstnodata -99...
+(ArcGIS10で利用する場合）ピラミッドレイヤと統計情報の作成
-- 作成したすべてのファイルに対してピラミッドレイヤと統計...
# 重たいのでArcGISのGUIは開かずにpythonでarcgisの機能を...
---すべてのプログラム→ArcGIS→python2.6→IDLE(Python GUI)を...
---*.py拡張子を付与したファイル名で保存する。
---以下の内容を打ち込む(workspaceの場所は適宜変えてくださ...
# arcpyをインポート
import arcpy
# ファイルが格納されているフォルダを指定する
workspace = "D://DEM//GSI5m//"
# ピラミッドレイヤと統計情報の構築
arcpy.BuildPyramidsandStatistics_management(workspace, "...
# 正常に終了したら"finish"と表示させる
print "finish"
---保存したのちrun→Run Moduleで実行する
終了行:
[[GIS/RS一般]]
*国土地理院のXML標高データをGeoTiffに変換する [#rf3221b1]
+ダウンロードした圧縮データを解凍し、変換するxmlファイル...
+株式会社エコリスの[["基盤地図情報標高DEM変換ツール">http...
# "変換結合.vbs"を開く~
# 投影法を "緯度経度"~
# xmlが格納されているフォルダを指定~
# Nodata(海の部分)は"-9999"または"0"と指定~
# ファイル名がメッシュコードのtif画像が出力される
+JGD2000緯度経度のままですべてのDEMをモザイク
--コマンドプロンプトを起動しカレントディレクトリをtif画像...
cd /d D:\DEM\GSI5m
--.tif拡張子のついたディレクトリ内のすべてのファイルのフ...
dir /b /s *.tif>list.txt~
--OSGeo4wシェルからgdal_merge.pyを使用してモザイク処理を...
# --optfile でlist.txtを読み込むことで、処理対象ファイル...
# -init でNodata(海の部分)は"-9999"または"0"と指定~
python C:\OSGeo4W\bin\gdal_merge.py -o D:/DEM/GSI5m/merge...
+投影変換~
-- gdalwarpを使用してJGD2000緯度経度からJGD2000UTM54帯の...
# -tr で解像度を指定（ここでは5m×5m）~
# -dstnodata でNodata値を"-9999"または"0"と指定 -s_srs...
gdalwarp -s_srs EPSG:4612 -t_srs EPSG:3100 -dstnodata -99...
+(ArcGIS10で利用する場合）ピラミッドレイヤと統計情報の作成
-- 作成したすべてのファイルに対してピラミッドレイヤと統計...
# 重たいのでArcGISのGUIは開かずにpythonでarcgisの機能を...
---すべてのプログラム→ArcGIS→python2.6→IDLE(Python GUI)を...
---*.py拡張子を付与したファイル名で保存する。
---以下の内容を打ち込む(workspaceの場所は適宜変えてくださ...
# arcpyをインポート
import arcpy
# ファイルが格納されているフォルダを指定する
workspace = "D://DEM//GSI5m//"
# ピラミッドレイヤと統計情報の構築
arcpy.BuildPyramidsandStatistics_management(workspace, "...
# 正常に終了したら"finish"と表示させる
print "finish"
---保存したのちrun→Run Moduleで実行する
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
