細密数値情報 （10m メッシュ土地利用）データの変換　QGIS読み込み

kanto.pl

#!/usr/local/bin/perl
#
#
#  細密数値情報 asc 変換  首都圏版
#
#
#   利用方法
#   Unix Linux 系     cat   <tdu file name> | perl kanto.pl > <output asc file name>
#   Win系            type    <tdu file name> | perl kanto.pl > <output asc file name>
#
#  利用例    type  s5_2010.tdu | perl kanto.pl  > s5_2010.asc
#
#  細密数値情報　首都圏のファイルを x,y,z 形式の ascii ファイル(asc ファイル)に変換する。
#　ただし　首都圏データのみ対応　　　（近畿圏、中部圏はたぶん原点座標が異なる）
#     http://www1.gsi.go.jp/geowww/LandUse/saimitsudr.html
#     http://www.gdal.org/frmt_xyz.html
#
#   2015.12.6    Yoichi.Kayama   yoichi.kayama@gmail.com
#

#             行単位ループ
while(<>) {

$filecode = substr( $_, 0, 4 );         #  ファイルコード
$linenum  = substr( $_, 4,  3 );        #  行番号


$xorg =-1 * ( substr($filecode, 0, 2 ) * 3 - 30 );    #  左上原点　X座標(測量の）　平面直角座標9系　Km単位
$yorg = substr( $filecode, 2, 2 ) * 4 -80;            #  左上原点　y座標(測量の)　平面直角座標9系　Km単位　　

    for ( $i=0; $i < 400; ++$i ) {                    #   デカルト座標系   Ｘ方向 カラム値の取得  400個
    
       $dcode = substr( $_, 7+($i*2), 2 );            #  値
       
       $xp = ($yorg * 1000) + ($i * 10 );             #   デカルト座標系   Ｘ座標　m単位　もしかすると　ここは -5　が必要かも
       $yp = ($xorg* 1000) - ($linenum * 10 );        #   デカルト座標系   Y座標　m単位　　　　　　　　こちらは　+5 ね
       
       print "${xp},${yp},${dcode}\n";                #   デカルト座標系  x,y,value    を出力

    }





}

tochiriyou.qml

<!DOCTYPE qgis PUBLIC 'http://mrcc.com/qgis.dtd' 'SYSTEM'>
<qgis version="2.12.0-Lyon" minimumScale="0" maximumScale="1e+08" hasScaleBasedVisibilityFlag="0">
  <pipe>
    <rasterrenderer opacity="1" alphaBand="-1" classificationMax="19" classificationMinMaxOrigin="User" band="1" classificationMin="0" type="singlebandpseudocolor">
      <rasterTransparency/>
      <rastershader>
        <colorrampshader colorRampType="DISCRETE" clip="0">
          <item alpha="255" value="1" label="山林・荒地等" color="#267500"/>
          <item alpha="255" value="2" label="田" color="#a3ff73"/>
          <item alpha="255" value="3" label="畑・その他の農地" color="#ffd37f"/>
          <item alpha="255" value="4" label="造成中地" color="#686868"/>
          <item alpha="255" value="5" label="空地" color="#ffff00"/>
          <item alpha="255" value="6" label="工業用地" color="#0070ff"/>
          <item alpha="255" value="7" label="一般低層住宅" color="#ffbee8"/>
          <item alpha="255" value="8" label="密集低層住宅地" color="#ff73df"/>
          <item alpha="255" value="9" label="中高層住宅地" color="#ff00c5"/>
          <item alpha="255" value="10" label="商業・業務地区" color="#ff0000"/>
          <item alpha="255" value="11" label="道路用地" color="#d7c29e"/>
          <item alpha="255" value="12" label="公園・緑地等" color="#70a800"/>
          <item alpha="255" value="13" label="その他の公共公益施設" color="#a900e6"/>
          <item alpha="255" value="14" label="河川・湖沼等" color="#00a5ff"/>
          <item alpha="255" value="15" label="その他" color="#828282"/>
          <item alpha="255" value="16" label="海" color="#002673"/>
          <item alpha="255" value="17" label="その他" color="#fffff0"/>
          <item alpha="255" value="18" label="ダミーコード" color="#000000"/>
          <item alpha="255" value="19" label="ダミーコード" color="#000000"/>
        </colorrampshader>
      </rastershader>
    </rasterrenderer>
    <brightnesscontrast brightness="0" contrast="3"/>
    <huesaturation colorizeGreen="128" colorizeOn="0" colorizeRed="255" colorizeBlue="128" grayscaleMode="0" saturation="0" colorizeStrength="100"/>
    <rasterresampler maxOversampling="2"/>
  </pipe>
  <blendMode>0</blendMode>
</qgis>

細密数値情報のファイルをＸＹＺ形式のテキストに変換するスクリプトです。
とりあえず首都圏分しか対応していません。
原点座標を変えれば近畿と中部も対応可能だと思います

http://www1.gsi.go.jp/geowww/LandUse/saimitsudr.html

ＸＹＺ形式にしたファイルはＱＧＩＳのラスタレイヤの追加で読み込み、ラスタレイヤとして利用できます。
この時首都圏のデータはデータの座標系を平面直角座標９系（EPSG;2451)で指定する必要があります。（近畿、中部は別の座標系です）

tochiriyou.qml　は読み込んだファイルに対して適用できるスタイル例です。

国土地理院ＷＥＢページでのデータ閲覧

変換したＸＹＺファイルを平面直角座標９系指定でＱＧＩＳに読み込む

各点の値別に色を指定する

http://www1.gsi.go.jp/geowww/LandUse/saimitsu_shuto.pdf　　首都圏データの解説

凡例のラベルを指定

このように表示できます