手順は東京都のデータと同じです。
データの内容を少し見てみます。
11-1 国土数値情報
座標系は都道府県に関係なく同じです。
user@deb-vmw:~/deb-win/dnli_shp/kanagawa$ ogrinfo -summary N03-100329_14_EC01.shp N03-100329_14_EC01
INFO: Open of `N03-100329_14_EC01.shp'
using driver `ESRI Shapefile' successful.
Layer name: N03-100329_14_EC01
Geometry: Polygon
Feature Count: 610
Extent: (118.915767, 35.128496) - (119.835839, 35.672897) <-1
Layer SRS WKT: <-2
(unknown)
PRN: String (8.0)
SUN: String (1.0)
CON: String (8.0)
CN2: String (10.0)
AAC: String (5.0)
1 緯度経度が度(degree)で表示
2 WKT(Well Known Text)で座標系が unknown (分からない)とき、特別な記載がないときは、EPSG:4326 で表示してみます。
11-2 基盤地図情報
11-2-1 神奈川県
AdmArea.shp のシェイプファイルの内容
user@deb-vmw:~/deb-win/bmi_shp/kanagawa$ ogrinfo -summary AdmArea.shp AdmArea
INFO: Open of `AdmArea.shp'
using driver `ESRI Shapefile' successful.
Layer name: AdmArea
Geometry: Polygon
Feature Count: 118
Extent: (-83341.671850, -96635.492871) - (-3373.197019, -36080.147641)
Layer SRS WKT:
PROJCS["JGD2000_Japan_Zone_9", <-1
GEOGCS["GCS_JGD_2000",
DATUM["Japanese_Geodetic_Datum_2000",
SPHEROID["GRS_1980",6378137.0,298.257222101]], <-2
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"], <-3
PARAMETER["False_Easting",0.0],
PARAMETER["False_Northing",0.0],
PARAMETER["Central_Meridian",139.833333333333], <-4
PARAMETER["Scale_Factor",0.9999],
PARAMETER["Latitude_Of_Origin",36], <-5
UNIT["Meter",1.0]]
ID: String (6.0)
UUID: String (27.0)
PRESENCE: Integer (8.0)
FINISHED: Integer (8.0)
ORGGILVL: String (5.0)
ORGMDID: String (1.0)
CATEGORY: String (4.0)
TYPE: String (6.0)
NAME: String (24.0)
CODE: Integer (5.0)
1 座標系 JGD2000(新日本測地系)ゾーン9
2 楕円体 GRS_1980
3 図法 横メルカトール図法 (tmerc)
4 中央子午線 139.83... (lon_0)
5 原点緯度 36 (lat_0)
を表しています。
/usr/share/proj/epsg から、JGD2000 に関係するコードをみてみます。
user@debian:/usr/share/proj$ grep -A1 JGD2000 epsg
# JGD2000
<4612> +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs <>
--
...
# JGD2000 / Japan Plane Rectangular CS IX
<2451> +proj=tmerc +lat_0=36 +lon_0=139.8333333333333 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
...
上は、経緯度表示、下は、平面直角座標系の定義のIXです。
これから、基盤地図情報の神奈川県のEPSGコードは 2451になります。
11-2-2 埼玉県
AdmArea.shp のシェイプファイルの内容
user@deb-vmw:~/deb-win/bmi_shp/saitama$ ogrinfo -summary AdmArea.shp AdmArea
INFO: Open of `AdmArea.shp'
using driver `ESRI Shapefile' successful.
Layer name: AdmArea
Geometry: Polygon
Feature Count: 82
Extent: (-101187.409166, -27323.089993) - (6047.540603, 31688.480017)
Layer SRS WKT:
PROJCS["JGD2000_Japan_Zone_9", <-1
GEOGCS["GCS_JGD_2000",
DATUM["Japanese_Geodetic_Datum_2000",
SPHEROID["GRS_1980",6378137.0,298.257222101]], <-2
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"], <-3
PARAMETER["False_Easting",0.0],
PARAMETER["False_Northing",0.0],
PARAMETER["Central_Meridian",139.833333333333], <-4
PARAMETER["Scale_Factor",0.9999],
PARAMETER["Latitude_Of_Origin",36], <-5
UNIT["Meter",1.0]]
ID: String (6.0)
UUID: String (27.0)
PRESENCE: Integer (8.0)
FINISHED: Integer (8.0)
ORGGILVL: String (5.0)
ORGMDID: String (1.0)
CATEGORY: String (4.0)
TYPE: String (6.0)
NAME: String (22.0)
CODE: Integer (5.0)
1 座標系 JGD2000(新日本測地系)ゾーン9
2 楕円体 GRS_1980
3 図法 横メルカトール図法 (tmerc)
4 中央子午線 119.83... (lon_0)
5 原点緯度 36 (lat_0)
を表しています。
これから、基盤地図情報の埼玉県のEPSGコードは 2451になります。
11-2-3 山梨県
AdmArea.shp のシェイプファイルの内容
user@deb-vmw:~/deb-win/bmi_shp/yamanashi$ ogrinfo -summary AdmArea.shp AdmArea
INFO: Open of `AdmArea.shp'
using driver `ESRI Shapefile' successful.
Layer name: AdmArea
Geometry: Polygon
Feature Count: 37
Extent: (-28933.909054, -92258.075411) - (57472.335133, -3130.562875)
Layer SRS WKT:
PROJCS["JGD2000_Japan_Zone_8", <-1
GEOGCS["GCS_JGD_2000",
DATUM["Japanese_Geodetic_Datum_2000",
SPHEROID["GRS_1980",6378137.0,298.257222101]], <-2
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"], <-3
PARAMETER["False_Easting",0.0],
PARAMETER["False_Northing",0.0],
PARAMETER["Central_Meridian",138.5], <-4
PARAMETER["Scale_Factor",0.9999],
PARAMETER["Latitude_Of_Origin",36], <-5
UNIT["Meter",1.0]]
ID: String (6.0)
UUID: String (27.0)
PRESENCE: Integer (8.0)
FINISHED: Integer (8.0)
ORGGILVL: String (5.0)
ORGMDID: String (1.0)
CATEGORY: String (4.0)
TYPE: String (6.0)
NAME: String (26.0)
CODE: Integer (5.0)
1 座標系 JGD2000(新日本測地系)ゾーン8
2 楕円体 GRS_1980
3 図法 横メルカトール図法 (tmerc)
4 中央子午線 138.5 (lon_0)
5 原点緯度 36 (lat_0)
# JGD2000 / Japan Plane Rectangular CS VIII
<2450> +proj=tmerc +lat_0=36 +lon_0=138.5 +k=0.9999 +x_0=0
平面直角座標系の定義のVIIIです。
これから、基盤地図情報の山梨県のEPSGコードは 2450になります。
■参考 grep -A1 JGD2000 epsg の内容
user@deb-vmw:/usr/share/proj$ grep -A1 JGD2000 epsg
# JGD2000
<4612> +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs <>
---
# JGD2000 / Japan Plane Rectangular CS I
<2443> +proj=tmerc +lat_0=33 +lon_0=129.5 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS II
<2444> +proj=tmerc +lat_0=33 +lon_0=131 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS III
<2445> +proj=tmerc +lat_0=36 +lon_0=132.1666666666667 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS IV
<2446> +proj=tmerc +lat_0=33 +lon_0=133.5 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS V
<2447> +proj=tmerc +lat_0=36 +lon_0=134.3333333333333 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS VI
<2448> +proj=tmerc +lat_0=36 +lon_0=136 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS VII
<2449> +proj=tmerc +lat_0=36 +lon_0=137.1666666666667 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS VIII
<2450> +proj=tmerc +lat_0=36 +lon_0=138.5 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS IX
<2451> +proj=tmerc +lat_0=36 +lon_0=139.8333333333333 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS X
<2452> +proj=tmerc +lat_0=40 +lon_0=140.8333333333333 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS XI
<2453> +proj=tmerc +lat_0=44 +lon_0=140.25 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS XII
<2454> +proj=tmerc +lat_0=44 +lon_0=142.25 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS XIII
<2455> +proj=tmerc +lat_0=44 +lon_0=144.25 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS XIV
<2456> +proj=tmerc +lat_0=26 +lon_0=142 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS XV
<2457> +proj=tmerc +lat_0=26 +lon_0=127.5 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS XVI
<2458> +proj=tmerc +lat_0=26 +lon_0=124 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS XVII
<2459> +proj=tmerc +lat_0=26 +lon_0=131 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS XVIII
<2460> +proj=tmerc +lat_0=20 +lon_0=136 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / Japan Plane Rectangular CS XIX
<2461> +proj=tmerc +lat_0=26 +lon_0=154 +k=0.9999 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
---
# JGD2000 / UTM zone 51N
<3097> +proj=utm +zone=51 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / UTM zone 52N
<3098> +proj=utm +zone=52 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / UTM zone 53N
<3099> +proj=utm +zone=53 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / UTM zone 54N
<3100> +proj=utm +zone=54 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
# JGD2000 / UTM zone 55N
<3101> +proj=utm +zone=55 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs <>
---
0 件のコメント:
コメントを投稿