Debianでka-MapのメモプラスOpenLayers とGeoExt: 12 基盤地図情報の数値標高モデル

基盤地図情報の数値標高モデルから標高画像データを作成します。

12-1 ダウンロード
「基盤地図情報（数値標高モデル）」の「10mメッシュ（標高）」の左側の＋をクリックし、「東京都（神奈川県、埼玉県、山梨県）」をチェックします。
「選択して次へ」ボタンをクリックします。
「ダウンロードファイルリスト」の「ダウンロード」ボタンをクリックしてダウンロードしました。
（tokyo_tokosen、 kanagawa_tokosen、 saitama_tokosen、 yamanashi_tokosen フォルダを作成しそれぞれにダウンロードしました。）
FG-GML-*****-DEM10B.zip

12-2 標高画像データの作成
「基盤地図情報閲覧コンバートソフト」
1 ツールバーの「新規作成」ボタンをクリックします。
2 「新規プロジェクトの作成」ダイアログで、
2-1「このプロジェクトのタイトル」：tokyo_tokosen
2-2「読み込むファイル」の「追加」ボタンをクリックします。
2-3 基盤地図情報データのあるフォルダを「ファイルの場所」で指定し、「ファイルの種類」を「基盤地図情報（JPGIS）圧縮ファイル（*.zip）」を選び、表示されたファイルを全部選んで「開く」ボタンをクリックします。
2-4「大量のデータに対して専用ファイルを構成して利用する」をチェックします。
2-5 「OK」ボタンをクリックします。
3 諸島部を含む東京都の地図が表示されます。（表示までに時間がかかります。）
表示部分と配色を変えました。
3-1 メニューバーの「表示」->「縮尺設定」をクリックします。
3-2 「表示縮尺設定」の「画面上でのおおよその縮尺」で「500000」、「あとで表示位置を設定する」にチェックをして「OK」ボタンをクリックします。
3-3 ポインタの先の表示範囲を東京に合わせてクリックします。
3-4 メニューバーの「設定」->「表示設定」をクリックします。
3-5 「表示設定」の「DEM表示設定」ボタンをクリックします。
3-6 現在の設定を「ファイルへ保存」ボタンをクリックして保存します。
3-7 「ペグの追加」で
2000m 茶色
100m 黄色
0m 黄緑
-4m 水色
と設定し「データの標高地で最適化」ボタンをクリックします。
3-8 現在の設定を「ファイルへ保存」ボタンをクリックして保存します。
tokosen.hcrm
3-9 「OK」ボタンをクリックします。
3-10 「表示設定」の「閉じる」ボタンをクリックします。
3-11 「ファイル」->「上書き保存」で保存します。
tokyo_tokosen フォルダに保存しました。

「表示縮尺設定」の「画面上でのおおよその縮尺」で「50000」まではシャギーが目立ちませんでした。

「ファイル」->「表示画面を画像データとして保存」で必要な部分をドラッグして選択し、「ワールドファイルも生成する」と「高解像度」にチェックをつけて保存します。
tokyo_tokosen.tif tokyo_tokosen.tfw

12-3 標高画像データの投影法の変換

基盤地図情報閲覧コンバートソフト FGDV で作成した標高画像データの情報を見てみます。


user@deb-vmw:~/mapdata$ gdalinfo tokyo_tokosen.tif 
Driver: GTiff/GeoTIFF
Files: tokyo_tokosen.tif
Size is 3255, 1632
Coordinate System is:
GEOGCS["JGD2000",
    DATUM["Japanese_Geodetic_Datum_2000",
        SPHEROID["GRS 1980",6378137,298.2572221010002,
            AUTHORITY["EPSG","7019"]],
        AUTHORITY["EPSG","6612"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4612"]]
Origin = (138.922483001119701,35.925224852022623)
Pixel Size = (0.000315613541764,-0.000275756703337)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_DOCUMENTNAME=PC-Mapping
  TIFFTAG_SOFTWARE=PCM(25000)
  TIFFTAG_XRESOLUTION=300
  TIFFTAG_YRESOLUTION=300
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( 138.9224830,  35.9252249) (138d55'20.94"E, 35d55'30.81"N)
Lower Left  ( 138.9224830,  35.4751899) (138d55'20.94"E, 35d28'30.68"N)
Upper Right ( 139.9498051,  35.9252249) (139d56'59.30"E, 35d55'30.81"N)
Lower Right ( 139.9498051,  35.4751899) (139d56'59.30"E, 35d28'30.68"N)
Center      ( 139.4361440,  35.7002074) (139d26'10.12"E, 35d42'0.75"N)
Band 1 Block=3255x1 Type=Byte, ColorInterp=Red
Band 2 Block=3255x1 Type=Byte, ColorInterp=Green
Band 3 Block=3255x1 Type=Byte, ColorInterp=Blue

作成された標高画像データは、EPSG:4326 になっています。
国土数値情報のマップファイル nippon_dnli_tokyo_pgis.map に次のようにレイヤを追加します。


---
 LAYER
   NAME tokyo_tokosen
   GROUP tokyo_tokosen
   TYPE RASTER
   STATUS ON
   DATA "../mapdata/tokyo_tokosen.tif"
 END

 LAYER
   NAME tokyo_kuiki
---

これを gdalwarp コマンドで投影法を EPSG:2456 に変換します。


user@deb-vmw:~/mapdata$ gdalwarp -s_srs EPSG:4326 -t_srs EPSG:2456 tokyo_tokosen.tif tokyo_tokosen2.tif
Creating output file that is 3252P x 1806L.
Processing input file tokyo_tokosen.tif.
0...10...20...30...40...50...60...70...80...90...100 - done.


user@deb-vm:~/mapdata$ gdalinfo tokyo_tokosen2.tif 
Driver: GTiff/GeoTIFF
Files: tokyo_tokosen2.tif
Size is 3252, 1806
Coordinate System is:
PROJCS["JGD2000 / Japan Plane Rectangular CS XIV",
    GEOGCS["JGD2000",
        DATUM["Japanese_Geodetic_Datum_2000",
            SPHEROID["GRS 1980",6378137,298.2572221010002,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6612"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4612"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",26],
    PARAMETER["central_meridian",142],
    PARAMETER["scale_factor",0.9999],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","2456"]]
Origin = (-279323.341314172896091,1104680.310845412779599)
Pixel Size = (29.000127704889998,-29.000127704889998)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( -279323.341, 1104680.311) (138d54'18.51"E, 35d55'29.20"N)
Lower Left  ( -279323.341, 1052306.080) (138d55'23.93"E, 35d27'12.17"N)
Upper Right ( -185014.926, 1104680.311) (139d56'57.51"E, 35d56'49.82"N)
Lower Right ( -185014.926, 1052306.080) (139d57'40.91"E, 35d28'31.40"N)
Center      ( -232169.134, 1078493.196) (139d26'4.86"E, 35d42'4.71"N)
Band 1 Block=3252x1 Type=Byte, ColorInterp=Red
Band 2 Block=3252x1 Type=Byte, ColorInterp=Green
Band 3 Block=3252x1 Type=Byte, ColorInterp=Blue

ワールドファイルも変更します。
元のワールドファイルは次のようになっています。

user@deb-vmw:~$ vim mapdata/tokyo_tokosen.tfw
0.00031561354176
0
0
-0.0002757567033
138.922640807891
35.925086973671

１行目と４行目は、gdalinfo で出力したデータの pixel サイズ、５行目と６行目は、値は少し違いますが、Origin または Upper Left の値みたいです。

変換した標高画像データのワールドファイルは次のようにました。

user@deb-vmw:~$ vim mapdata/tokyo_tokosen2.tfw
29.0001277048
0
0
-29.0001277048
-279323.3413141728
1104680.3108454127

基盤地図情報のマップファイル nippon_bmi_tokyo_pgis.map に次のようにレイヤを追加します。


---
 LAYER
   NAME tokyo_tokosen
   GROUP tokyo_tokosen
   TYPE RASTER
   STATUS ON
   DATA "../mapdata/tokyo_tokosen2.tif"
 END

 LAYER
   NAME tokyo_kukaku
---

12-4 同じ方法で神奈川県、埼玉県、山梨県の標高画像データを作成しました。

12-4-1 神奈川県


user@deb-vmw:~/mapdata$ gdalinfo kanagawa_tokosen.tif 
Driver: GTiff/GeoTIFF
Files: kanagawa_tokosen.tif
Size is 2436, 1815
Coordinate System is:
GEOGCS["JGD2000",
    DATUM["Japanese_Geodetic_Datum_2000",
        SPHEROID["GRS 1980",6378137,298.2572221010002,
            AUTHORITY["EPSG","7019"]],
        AUTHORITY["EPSG","6612"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4612"]]
Origin = (138.883886386362491,35.686520927277449)
Pixel Size = (0.000389353672356,-0.000318376607750)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_DOCUMENTNAME=PC-Mapping
  TIFFTAG_SOFTWARE=PCM(25000)
  TIFFTAG_XRESOLUTION=300
  TIFFTAG_YRESOLUTION=300
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( 138.8838864,  35.6865209) (138d53'1.99"E, 35d41'11.48"N)
Lower Left  ( 138.8838864,  35.1086674) (138d53'1.99"E, 35d 6'31.20"N)
Upper Right ( 139.8323519,  35.6865209) (139d49'56.47"E, 35d41'11.48"N)
Lower Right ( 139.8323519,  35.1086674) (139d49'56.47"E, 35d 6'31.20"N)
Center      ( 139.3581192,  35.3975942) (139d21'29.23"E, 35d23'51.34"N)
Band 1 Block=2436x2 Type=Byte, ColorInterp=Red
Band 2 Block=2436x2 Type=Byte, ColorInterp=Green
Band 3 Block=2436x2 Type=Byte, ColorInterp=Blue


user@deb-vmw:~/mapdata$ gdalwarp -s_srs EPSG:4326 -t_srs EPSG:2451 kanagawa_tokosen.tif kanagawa_tokosen2.tif
Creating output file that is 2446P x 1825L.
Processing input file kanagawa_tokosen.tif.
0...10...20...30...40...50...60...70...80...90...100 - done.


user@deb-vm:~/mapdata$ gdalinfo kanagawa_tokosen2.tif 
Driver: GTiff/GeoTIFF
Files: kanagawa_tokosen2.tif
Size is 2446, 1825
Coordinate System is:
PROJCS["JGD2000 / Japan Plane Rectangular CS IX",
    GEOGCS["JGD2000",
        DATUM["Japanese_Geodetic_Datum_2000",
            SPHEROID["GRS 1980",6378137,298.2572221010002,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6612"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4612"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",36],
    PARAMETER["central_meridian",139.8333333333333],
    PARAMETER["scale_factor",0.9999],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","2451"]]
Origin = (-86551.229141840289230,-34363.557355713906873)
Pixel Size = (35.350898701874428,-35.350898701874428)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  -86551.229,  -34363.557) (138d52'37.52"E, 35d41'11.28"N)
Lower Left  (  -86551.229,  -98878.947) (138d53'2.14"E, 35d 6'17.99"N)
Upper Right (     -82.931,  -34363.557) (139d49'56.70"E, 35d41'24.95"N)
Lower Right (     -82.931,  -98878.947) (139d49'56.72"E, 35d 6'31.37"N)
Center      (  -43317.080,  -66621.252) (139d21'23.24"E, 35d23'54.80"N)
Band 1 Block=2446x1 Type=Byte, ColorInterp=Red
Band 2 Block=2446x1 Type=Byte, ColorInterp=Green
Band 3 Block=2446x1 Type=Byte, ColorInterp=Blue

user@deb-vmw:~/mapdata$ vim kanagawa_tokosen.tfw
0.00038935367236
0
0
-0.0003183766078
138.884081063199
35.6863617389736

user@deb-vmw:~/mapdata$ vim kanagawa_tokosen2.tfw
35.3508987018744
0
0
-35.350898701874
-86551.229141840
-34363.557355713

12-4-2 埼玉県


user@deb-vmw:~/mapdata$ gdalinfo saitama_tokosen.tif 
Driver: GTiff/GeoTIFF
Files: saitama_tokosen.tif
Size is 2682, 1512
Coordinate System is:
GEOGCS["JGD2000",
    DATUM["Japanese_Geodetic_Datum_2000",
        SPHEROID["GRS 1980",6378137,298.2572221010002,
            AUTHORITY["EPSG","7019"]],
        AUTHORITY["EPSG","6612"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4612"]]
Origin = (138.673237951228458,36.304917546304509)
Pixel Size = (0.000478243821535,-0.000388520895571)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_DOCUMENTNAME=PC-Mapping
  TIFFTAG_SOFTWARE=PCM(25000)
  TIFFTAG_XRESOLUTION=300
  TIFFTAG_YRESOLUTION=300
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( 138.6732380,  36.3049175) (138d40'23.66"E, 36d18'17.70"N)
Lower Left  ( 138.6732380,  35.7174740) (138d40'23.66"E, 35d43'2.91"N)
Upper Right ( 139.9558879,  36.3049175) (139d57'21.20"E, 36d18'17.70"N)
Lower Right ( 139.9558879,  35.7174740) (139d57'21.20"E, 35d43'2.91"N)
Center      ( 139.3145629,  36.0111957) (139d18'52.43"E, 36d 0'40.30"N)
Band 1 Block=2682x2 Type=Byte, ColorInterp=Red
Band 2 Block=2682x2 Type=Byte, ColorInterp=Green
Band 3 Block=2682x2 Type=Byte, ColorInterp=Blue


user@deb-vmw:~/mapdata$ gdalwarp -s_srs EPSG:4326 -t_srs EPSG:2451 saitama_tokosen.tif saitama_tokosen2.tif
Creating output file that is 2692P x 1526L.
Processing input file saitama_tokosen.tif.
0...10...20...30...40...50...60...70...80...90...100 - done.


user@deb-vmw:~/mapdata$ gdalinfo saitama_tokosen2.tif 
Driver: GTiff/GeoTIFF
Files: saitama_tokosen2.tif
Size is 2692, 1526
Coordinate System is:
PROJCS["JGD2000 / Japan Plane Rectangular CS IX",
    GEOGCS["JGD2000",
        DATUM["Japanese_Geodetic_Datum_2000",
            SPHEROID["GRS 1980",6378137,298.2572221010002,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6612"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4612"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",36],
    PARAMETER["central_meridian",139.8333333333333],
    PARAMETER["scale_factor",0.9999],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","2451"]]
Origin = (-104962.096120072616031,34455.375603438791586)
Pixel Size = (43.109948452022763,-43.109948452022763)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( -104962.096,   34455.376) (138d39'52.55"E, 36d18'17.40"N)
Lower Left  ( -104962.096,  -31330.406) (138d40'23.94"E, 35d42'43.25"N)
Upper Right (   11089.885,   34455.376) (139d57'24.58"E, 36d18'37.74"N)
Lower Right (   11089.885,  -31330.406) (139d57'21.27"E, 35d43'3.15"N)
Center      (  -46936.106,    1562.485) (139d18'45.48"E, 36d 0'46.63"N)
Band 1 Block=2692x1 Type=Byte, ColorInterp=Red
Band 2 Block=2692x1 Type=Byte, ColorInterp=Green
Band 3 Block=2692x1 Type=Byte, ColorInterp=Blue

user@deb-vmw:~/mapdata$ vim saitama_tokosen.tfw
0.00047824382154
0
0
-0.0003885208956
138.673477073139
36.3047232858567

user@deb-vmw:~/mapdata$ vim saitama_tokosen2.tfw
43.1099484520227
0
0
-43.109948452022
-104962.09612007
34455.3756034387

12-4-3 山梨県


user@deb-vmw:~/mapdata$ gdalinfo yamanashi_tokosen.tif 
Driver: GTiff/GeoTIFF
Files: yamanashi_tokosen.tif
Size is 1677, 1602
Coordinate System is:
GEOGCS["JGD2000",
    DATUM["Japanese_Geodetic_Datum_2000",
        SPHEROID["GRS 1980",6378137,298.2572221010002,
            AUTHORITY["EPSG","7019"]],
        AUTHORITY["EPSG","6612"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4612"]]
Origin = (138.109818475870782,35.999445358876471)
Pixel Size = (0.000666855886739,-0.000545231105732)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_DOCUMENTNAME=PC-Mapping
  TIFFTAG_SOFTWARE=PCM(25000)
  TIFFTAG_XRESOLUTION=300
  TIFFTAG_YRESOLUTION=300
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( 138.1098185,  35.9994454) (138d 6'35.35"E, 35d59'58.00"N)
Lower Left  ( 138.1098185,  35.1259851) (138d 6'35.35"E, 35d 7'33.55"N)
Upper Right ( 139.2281358,  35.9994454) (139d13'41.29"E, 35d59'58.00"N)
Lower Right ( 139.2281358,  35.1259851) (139d13'41.29"E, 35d 7'33.55"N)
Center      ( 138.6689771,  35.5627152) (138d40'8.32"E, 35d33'45.77"N)
Band 1 Block=1677x3 Type=Byte, ColorInterp=Red
Band 2 Block=1677x3 Type=Byte, ColorInterp=Green
Band 3 Block=1677x3 Type=Byte, ColorInterp=Blue


user@deb-vmw:~/mapdata$ gdalwarp -s_srs EPSG:4326 -t_srs EPSG:2450 yamanashi_tokosen.tif yamanashi_tokosen2.tif
Creating output file that is 1686P x 1607L.
Processing input file yamanashi_tokosen.tif.
0...10...20...30...40...50...60...70...80...90...100 - done.


user@deb-vmw:~/mapdata$ gdalinfo yamanashi_tokosen2.tif 
Driver: GTiff/GeoTIFF
Files: yamanashi_tokosen2.tif
Size is 1686, 1607
Coordinate System is:
PROJCS["JGD2000 / Japan Plane Rectangular CS VIII",
    GEOGCS["JGD2000",
        DATUM["Japanese_Geodetic_Datum_2000",
            SPHEROID["GRS 1980",6378137,298.2572221010002,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6612"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4612"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",36],
    PARAMETER["central_meridian",138.5],
    PARAMETER["scale_factor",0.9999],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","2450"]]
Origin = (-35560.814174073195318,183.648549885500898)
Pixel Size = (60.468102128726208,-60.468102128726208)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  -35560.814,     183.649) (138d 6'19.99"E, 36d 0'3.62"N)
Lower Left  (  -35560.814,  -96988.592) (138d 6'35.36"E, 35d 7'30.46"N)
Upper Right (   66388.406,     183.649) (139d14'10.94"E, 35d59'57.82"N)
Lower Right (   66388.406,  -96988.592) (139d13'42.25"E, 35d 7'24.84"N)
Center      (   15413.796,  -48402.472) (138d40'12.14"E, 35d33'48.96"N)
Band 1 Block=1686x1 Type=Byte, ColorInterp=Red
Band 2 Block=1686x1 Type=Byte, ColorInterp=Green
Band 3 Block=1686x1 Type=Byte, ColorInterp=Blue

user@deb-vmw:~/mapdata$ vim yamanashi_tokosen.tfw
0.00066685588674
0
0
-0.0005452311057
138.110151903814
35.9991727433236

user@deb-vmw:~/mapdata$ vim yamanashi_tokosen2.tfw
60.4681021287262
0
0
-60.468102128726
-35560.814174073
183.648549885500

Debianでka-MapのメモプラスOpenLayers とGeoExt

ブログアーカイブ

自己紹介

2011年4月20日水曜日

12 基盤地図情報の数値標高モデル

0 件のコメント:

Debianでka-Mapのメモ プラスOpenLayers とGeoExt

ブログ アーカイブ