国土地理院の数値地図が試験期間を終え、いつの間にか、
基盤地図情報閲覧サービスとして発足している。 基盤情報として測量基準点、街区境界線、標高メッシュの情報が公開されている。
ただし、XML形式なので、javascriptやepsで利用するには解体してやる必要がある。
たとえば、基盤地図情報（縮尺レベル25000）で長野県を選択すると
FG-GML-20-AdmPt25000-20080331-0001.xml （市町村代表点）
FG-GML-20-AdmBdry25000-20080331-0001.xml （境界線）
FG-GML-20-AdmArea25000-20080331-0001.xml （市町村ポリゴン）
がダウンロードできる。中身はxmlEditor（windows）かXcode（mac）で閲覧できる。
市町村ポリゴンは、タグAdmAreaを単位として島や飛び地は別扱い。多角形はタグgml:exteriorが外周を、タグgml:interiorが空洞を表している。
PythonのTreeElementモジュールは、DOMを解析してくれる。
getで属性、findで子要素、getiteratorで子孫要素リストを返してくれる。名前空間を付ける必要がある。
from xml.etree.ElementTree import *
import codecs
import Polygon
csvOut = codecs.getwriter('utf-8')(file('areadata.csv', 'w'))
strLine = u"Pref\tCode\tName\tDirection\tArea\tPointArray\r"
csvOut.write(strLine)
for pref in xrange(1,48):
strpref = str(pref).zfill(2)
print strpref
unicode_file = codecs.open("FG-GML-" +strpref +"-05-Z001/FG-GML-" +strpref+ "-AdmArea25000-20080331-0001.xml", 'r', 'shift_jis').read().split('\n')
tmp_file = codecs.open('tmp.xml', 'w', 'UTF-8')
for line in unicode_file[2:]:
# １行目のDOC宣言を削除
tmp_file.write(line)
tmp_file.close()
elem = parse(open('tmp.xml', 'r'))
for e in elem.getiterator("{http://fgd.gsi.go.jp/spec/2008/FGD_GMLSchema}AdmArea"):
pre = unicode(strpref) + u"\t"+ e.find('{http://fgd.gsi.go.jp/spec/2008/FGD_GMLSchema}admCode').text +u"\t"+ e.find('{http://fgd.gsi.go.jp/spec/2008/FGD_GMLSchema}name').text
print pre
for ext in e.getiterator("{http://www.opengis.net/gml/3.2}exterior"):
for pos in ext.getiterator("{http://www.opengis.net/gml/3.2}posList"):
strLine = pre +u"\t1"
strPos = u""
pointList = list()
for coord in pos.text[:-1].split(u'\n')[1:]:
c = coord[1:].split(u' ')
strPos = strPos +u"\t" + c[0] +u"\t" + c[1]
pointList.append([float(c[0]), float(c[1])])
poly = Polygon.Polygon(pointList)
strLine = pre +u"\t1\t" + unicode(poly.area())
csvOut.write(strLine + strPos + u"\r")
for ext in e.getiterator("{http://www.opengis.net/gml/3.2}interior"):
for pos in ext.getiterator("{http://www.opengis.net/gml/3.2}posList"):
strLine = pre +u"\t1"
strPos = u""
pointList = list()
for coord in pos.text[:-1].split(u'\n')[1:]:
c = coord[1:].split(u' ')
strPos = strPos +u"\t" + c[0] +u"\t" + c[1]
pointList.append([float(c[0]), float(c[1])])
poly = Polygon.Polygon(pointList)
strLine = pre +u"\t1\t" + unicode(poly.area())
csvOut.write(strLine + strPos + u"\r")
」（邦訳は失われているので公正としての正義 再説
しかない）を批判した教授で、お陰で「正義論」改訂版が出た。「Justice as fairness」とか「無知のベール」という発想はアメリカでしか出てこなかっただろう。サンデル教授がどんな講義をするか楽しみだ。
