意匠の時空間的分析

基礎な手法を用いて日本における登録意匠の時空間的分析などを探検する目的で、このブログによって、手法などについてシェアしたり、意見交換したりします。

佐賀県と福岡県の市区町村ペンタゴンをダウンロードし、Rを用いて表示する手順

目的:2010年3月29日現在の市区町村シェープファイルをダウンロードし、Rにインポートと表示をすること

  1. ダウロード

    URL http://nlftp.mlit.go.jp/ksj/index.html

    f:id:ttcrossroads09:20180927140151p:plain

    ここにクリック

    URL http://nlftp.mlit.go.jp/ksj/gml/datalist/KsjTmplt-N03-v2_3.html

    f:id:ttcrossroads09:20180927140755p:plain

    3.ファイルの選択 http://nlftp.mlit.go.jp/ksj/gml/cgi-bin/download.php (1)

    f:id:ttcrossroads09:20180927143829p:plain

    次へ

    アンケートのご協力のお願い

    回答する

    4.国土数値情報利用約款

    上記のことを理解しました。-> はい

    5.ダウンロード

  2. ZIPファイルを展開する

    佐賀県の例)
    KS-META-N03-100329_41-g.xml
    N03-100329_41-g.xml
    N03-100329_41-g_AdministrativeBoundary.dbf
    N03-100329_41-g_AdministrativeBoundary.shp
    N03-100329_41-g_AdministrativeBoundary.shx
  3. Rへのインポート

    コード 備考
    library(rgdal) rgdal
    setwd(".../R Projects/Saga-Fukuoka") シェープファイルのフォルダー
    dsnFukuoka <- c(".../R Projects/Saga-Fukuoka/N03-100329_40-g_AdministrativeBoundary.dbf")
    ogrInfo(dsnFukuoka) チェック
    Source: "...\R Projects\Saga-Fukuoka\N03-100329_40-g_AdministrativeBoundary.dbf", layer: "N03-100329_40-g_AdministrativeBoundary"
    Driver: ESRI Shapefile; number of rows: 666
    Feature type: wkbPolygon with 2 dimensions
    Extent: (130.0322 32.99997) - (131.1906 34.24996)
    LDID: 19
    Number of fields: 7
    name type length typeName
    1 N03_001 4 8 String
    2 N03_002 4 10 String
    3 N03_003 4 10 String
    4 N03_004 4 14 String
    5 N03_005 4 1 String
    6 N03_006 4 1 String
    7 N03_007 4 5 String
    fukuoka <- readOGR(dsn=dsnFukuoka)
    plot(fukuoka)

    f:id:ttcrossroads09:20180927154733p:plain

    dsnSaga <- c(".../R Projects/Saga-Fukuoka/N03-100329_41-g_AdministrativeBoundary.dbf")
    saga <- readOGR(dsn=dsnSaga)
    plot(saga)

    f:id:ttcrossroads09:20180927155403p:plain