J-SHIS - 4-5 条件付超過確率 地図データ

4-5 条件付超過確率　地図データ
条件付超過確率地図は、想定した地震が発生した場合に、個々の地点での揺れがある震度を上回る確率を示した地図です。（「用語集[http://www.j-shis.bosai.go.jp/glossary]」より）

「項目」の「条件付超過確率」の「地図データ」をクリックします。

条件を次のように設定します。

圧縮形式：zip
ファイルフォーマット：シェープファイル
地震種別：主要活断層

断層を選択します。

震源断層：布田川断層帯宇土区間

「ダウンロード開始」ボタンをクリックします。

ダウンロードしたファイルを解凍します。

user@deb8-vmw:~$ cd ダウンロード
user@deb8-vmw:~/ダウンロード$ ls
---
C-V2-F013002-SHAPE.zip
---

user@deb8-vmw:~/ダウンロード$ unzip C-V2-F013002-SHAPE.zip 
Archive:  C-V2-F013002-SHAPE.zip
   creating: F013002/
  inflating: F013002/C-V2-F013002-MAP-SHAPE-CASE1.dbf  
  inflating: F013002/C-V2-F013002-MAP-SHAPE-CASE1.shx  
  inflating: F013002/C-V2-F013002-MAP-SHAPE-CASE1.prj  
  inflating: F013002/C-V2-F013002-MAP-SHAPE-CASE1.shp

user@deb8-vmw:~/ダウンロード$ mv F013002/C-V2-F013002-MAP-SHAPE-CASE* ../mapdata/jshis/
user@deb8-vmw:~/ダウンロード$ cd ../mapdata/

user@deb8-vmw:~/mapdata/$ ogrinfo jshis
INFO: Open of `jshis'
      using driver `ESRI Shapefile' successful.
---
8: C-V2-F013002-MAP-SHAPE-CASE1 (Polygon)
---

user@deb8-vmw:~/mapdata/$ ogrinfo -summary jshis C-V2-F013002-MAP-SHAPE-CASE1
INFO: Open of `F013002'
      using driver `ESRI Shapefile' successful.

Layer name: C-V2-F013002-MAP-SHAPE-CASE1
Geometry: Polygon
Feature Count: 666351
Extent: (129.081250, 31.189583) - (132.331250, 34.345833)
Layer SRS WKT:
GEOGCS["JGD2000",
    DATUM["Japanese_Geodetic_Datum_2000",
        SPHEROID["GRS_1980",6378137,298.257222101]],
    PRIMEM["Greenwich",0],
    UNIT["Degree",0.017453292519943295]]

CODE: String (11.0)
AVE_SI: Real (9.5)
I45_PS: Real (9.5)
I50_PS: Real (9.5)
I55_PS: Real (9.5)
I60_PS: Real (9.5)

「データダウンロード」の「項目」の下にある「規約集」から、属性値は次の通りです。

属性一覧
属性名  型    書式 説明
CODE   String 11  250m メッシュコード(日本測地系)
AVE_SI Double 9.5 地表の計測震度(平均値)
I45_PS Double 9.5 震度 5 弱以上となる確率
I50_PS Double 9.5 震度 5 強以上となる確率
I55_PS Double 9.5 震度 6 弱以上となる確率
I60_PS Double 9.5 震度 6 強以上となる確率

user@deb8-vmw:~/mapdata/$ cd ../mapfile/
user@deb8-vmw:~/mapfile$ vim jshis-4930.map

MAP
 NAME jshis-4930_map #任意の名前
 STATUS ON
 SIZE 600 600
# EXTENT 130.000000 32.666667 131.000000 33.333333
 Extent 129.081250 31.189583 132.331250 34.345833
 #ogrinfoで得た情報から

---
 LAYER
  NAME i60ps #任意の名前
  TYPE POLYGON #ogrinfoで得た情報から
  STATUS ON
  DATA "F013002/C-V2-F013002-MAP-SHAPE-CASE1.shp"
  #シェイプファイルの相対位置
  OPACITY 80
  CLASSITEM "I60_PS"

  CLASS
   NAME "50 to 100" #任意の名前
   EXPRESSION ([I60_PS] > 0.5 AND ([I60_PS] <= 1)
   STYLE
    COLOR 102 0 0
   END
  END

  CLASS
   NAME "40 to 50"
   EXPRESSION ([I60_PS] > 0.4 AND [I60_PS] <= 0.5 )
   STYLE
    COLOR 255 51 0
   END
  END

  CLASS
   NAME "30 to 40"
   EXPRESSION ([I60_PS] > 0.3 AND [I60_PS] <= 0.4 )
   STYLE
    COLOR 255 102 51
   END
  END

  CLASS
   NAME "20 to 30"
   EXPRESSION ([I60_PS] > 0.2 AND [I60_PS] <= 0.3 )
   STYLE
    COLOR 255 153 102
   END
  END

  CLASS
   NAME "10 to 20"
   EXPRESSION ([I60_PS] > 0.1 AND [I60_PS] <= 0.2 )
   STYLE
    COLOR 255 204 153
   END
  END

  CLASS
   NAME "1 to 10"
   EXPRESSION ([I60_PS] > 0.01 AND [I60_PS] <= 0.1 )
   STYLE
    COLOR 255 255 153
   END
  END

  CLASS
   NAME 0 to 1"
   EXPRESSION ([I60_PS] > 0 AND [I60_PS] <= 0.01 )
   STYLE
    COLOR 153 255 153
   END
  END
 END
---

http://localhost/cgi-bin/mapserv?mode=browse&template=openlayers&layer=dansoparam&layer=i60ps&layer=kumamoto&layer=saga&layer=fukuoka&layer=nagasaki&map=/home/user/mapfile/jshis-4930.map

「震度6強以上となる確率」