37 - 陸前高田市の地図の表示 6 - 数値標高モデル 2

37-6 陸前高田市の基盤地図情報の数値標高モデルをデータベースに追加
37-6-1 陸前高田市の基盤地図情報シェイプファイルを保存
陸前高田市の基盤地図情報のデータを追加します。
a 「設定」 -> 「プロジェクト設定」をクリックします。

b 「プロジェクト設定」ダイアログボックスの「追加読み込み」ボタンをクリックします。

c 「ファイルを開く」ダイアログボックスで陸前高田市の基盤地図情報のデータを選択して「開く」ボタンをクリックします。

d 「プロジェクト設定」ダイアログボックスの「表示設定」ボタンをクリックします。

e 「表示設定」ダイアログボックスで
「形状の表示・非表示」：行政区画（輪郭線）　のみ
「名称等の表示・非表示」：（すべてOFF）
「閉じる」ボタンをクリック







f 「プロジェクト設定」ダイアログボックスの「OK」ボタンをクリックします。
g ウィンドウを小さくして陸前高田市だけ表示されるようにします。
（画像を拡大すると時間がかかりました。）
 「エクスポート」 -> 「標高メッシュをシェープファイルへ出力」をクリックします。

h 「標高メッシュデータのシェープファイルデータへの変換」ダイアログボックスで
「（おおむね）画像が表示されている要素のみを出力」をチェック
「出力先ファイル」欄の右側のフォルダアイコンをクリック

i 「ファイルの保存先」ダイアログボックスで
「保存する場所」を選択
「ファイル名」: rikuzentakata5m を入力
「保存」ボタンをクリック



j 「標高メッシュデータのシェープファイルデータへの変換」ダイアログボックスの「OK」ボタンをクリックします。

37-6-2 データベースに追加
rikuzentakata5m.dbfのタイトルが日本語なのでアルファベットにします。
LibreOffice のCalc で開きます。
（私の場合は、Debian では修正後の保存ができなかったので、Windows で修正しました。もしくは、PostgresQL に登録してから変更します。）
標高 -> elevation
user@debian7-vmw:~/mapdata/rikuzentakata$ cd ../
user@debian7-vmw:~/mapdata$ ogrinfo -summary rikuzentakata rikuzentaka5m.shp

INFO: Open of `rikuzentakata'
      using driver `ESRI Shapefile' successful.
FAILURE: Couldn't fetch requested layer rikuzentaka5m.shp!
user@debian7-vmw:~/mapdata$ ogrinfo -summary rikuzentakata rikuzentaka5m
INFO: Open of `rikuzentakata'
      using driver `ESRI Shapefile' successful.

Layer name: rikuzentaka5m
Geometry: 3D Point
Feature Count: 14733499
Extent: (55997.639472, -118897.330308) - (79068.143237, -96565.386996)
Layer SRS WKT:
PROJCS["JGD2000_Japan_Zone_10",
    GEOGCS["GCS_JGD_2000",
        DATUM["Japanese_Geodetic_Datum_2000",
            SPHEROID["GRS_1980",6378137.0,298.257222101]],
        PRIMEM["Greenwich",0.0],
        UNIT["Degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["False_Easting",0.0],
    PARAMETER["False_Northing",0.0],
    PARAMETER["Central_Meridian",140.833333333333],
    PARAMETER["Scale_Factor",0.9999],
    PARAMETER["Latitude_Of_Origin",40],
    UNIT["Meter",1.0]]
�W��: Real (6.2)

user@debian7-vmw:~/mapdata$ cd rikuzentakata/
user@debian7-vmw:~/mapdata/rikuzentakata$ shp2pgsql -W sjis rikuzentaka5m.shp rikuzentaka5m > rikuzentaka5m.sql
Shapefile type: PointZ
Postgis type: POINT[4]
user@debian7-vmw:~/mapdata/rikuzentakata$ psql -d nippon_bmi -f rikuzentaka5m.sql
user@debian7-vmw:~/mapdata/rikuzentakata$ psql nippon_bmi
psql (9.1.9)
"help" でヘルプを表示します.

nippon_bmi=> \dt

                   リレーションの一覧
 スキーマ |          名前            |    型    | 所有者 
----------+-------------------------+----------+--------
---
 public   | rikuzentaka5m            | テーブル | user
 public   | rikuzentakata_choaza     | テーブル | user
 public   | rikuzentakata_doro       | テーブル | user
 public   | rikuzentakata_dorokozo   | テーブル | user
 public   | rikuzentakata_kaigan     | テーブル | user
 public   | rikuzentakata_kenchiku   | テーブル | user
 public   | rikuzentakata_kido       | テーブル | user
 public   | rikuzentakata_kukakusen  | テーブル | user
 public   | rikuzentakata_suigai     | テーブル | user
 public   | rikuzentakata_suiiki     | テーブル | user
 public   | rikuzentakata_suikozo    | テーブル | user
 public   | rikuzentakata_suikozosen | テーブル | user
---
(34 行)

(END):q
nippon_bmi=> \d rikuzentaka5m

                              テーブル "public.rikuzentaka5m"
  カラム  |        型        |                           修飾語                            
----------+------------------+-------------------------------------------------------------
 gid      | integer          | not null default nextval('rikuzentaka5m_gid_seq'::regclass)
 標高     | double precision | 
 the_geom | geometry         | 
インデックス:
    "rikuzentaka5m_pkey" PRIMARY KEY, btree (gid)
CHECK 制約:
    "enforce_dims_the_geom" CHECK (st_ndims(the_geom) = 4)
    "enforce_geotype_the_geom" CHECK (geometrytype(the_geom) = 'POINT'::text OR the_geom IS NULL)
    "enforce_srid_the_geom" CHECK (st_srid(the_geom) = (-1))

~
(stdin):q

nippon_bmi=> ALTER TABLE rikuzentaka5m RENAME COLUMN 標高 to elevation;
ALTER TABLE
nippon_bmi=> \d rikuzentaka5m

                              テーブル "public.rikuzentaka5m"
  カラム   |        型        |                           修飾語                            
-----------+------------------+------------------------------------------------------------
 gid       | integer          | not null default nextval('rikuzentaka5m_gid_seq'::regclass)
 elevation | double precision | 
 the_geom  | geometry         | 
インデックス:
    "rikuzentaka5m_pkey" PRIMARY KEY, btree (gid)
CHECK 制約:
    "enforce_dims_the_geom" CHECK (st_ndims(the_geom) = 4)
    "enforce_geotype_the_geom" CHECK (geometrytype(the_geom) = 'POINT'::text OR the_geom IS NULL)
    "enforce_srid_the_geom" CHECK (st_srid(the_geom) = (-1))

~
(stdin):q

nippon_bmi=> \q