XMLを配列にする部分を関数化するに当たって、関数内でエラーが発生した場合、どうすりゃいいの？ 「例外とか投げるのかな？」とググったら見つかった。 8. エラーと例外 — Python 3.6.1 ドキュメント XML内に<gml:tupleList>が無い場合にExceptionを投げて、例外があっても</gml:tuplelist>…

配列をMatplotlibで図にできるのが便利なので、2次メッシュ全体を1つの図に表示できたらいいなぁと配列の結合でググった。 numpy.concatenate — NumPy v1.13 Manual やはりサイズというか結合する方向の項目数を合わせないといけないので、まずは横に結合し…

# coding: utf-8 import xml.etree.ElementTree as ET import numpy as np import matplotlib.pyplot as plot import sys GEO_DIR = "FG-GML-{0:04d}-{1:02d}-DEM5A" GEO_XML = "FG-GML-{0:04d}-{1:02d}-{2:02d}-DEM5A-20161001.xml" if len(sys.argv) < 4: …

基盤地図情報～数値標高モデルを利用する - Qiita このページを参考にNumPyを使ってみた。 # coding: utf-8 import xml.etree.ElementTree as ET import numpy as np import sys GEO_DIR = "FG-GML-{0:04d}-{1:02d}-DEM5A" GEO_XML = "FG-GML-{0:04d}-{1:02d…

# coding: utf-8 import xml.etree.ElementTree as ET import sys GEO_DIR = "FG-GML-{0:04d}-{1:02d}-DEM5A" GEO_XML = "FG-GML-{0:04d}-{1:02d}-{2:02d}-DEM5A-20161001.xml" if len(sys.argv) < 4: print("xmlparse.py mesh1 mesh2 mesh3") exit(-1) mesh…

ググっていたらこんな情報があった。 qiita.com ふむふむ、5340は次のような意味になるのか。 意味 値 緯度 53 / 1.5 = 35.3333 経度 40 + 100 = 140

100個のXMLファイルを個別に処理するにしても、その度にソースコードを書き換えるのは面倒なので、引数で指定したい。 引数 引数にアクセスするにはsysをインポートして、配列sys.argvで参照する。 # coding: utf-8 import sys print(sys.argv) argv.pyとし…

PythonのXMLパーサーをググると、これが見つかりました。 20.5. xml.etree.ElementTree — ElementTree XML API — Python 3.6.1 ドキュメント まずはこのXMLをパースして<gml:tupleList>を見つけて、中のテキストを取り出す（だけの）プログラム(xmlparse.py)を書いてみまし</gml:tuplelist>…

ダウンロードしたファイルは、FG-GML-5340-22-DEM5A.zipというファイルでした。 解凍すると、FG-GML-5340-22-DEM5Aというディレクトリの中に次のような100個のXMLファイルができました。 FG-GML-5340-22-nn-DEM5A-20161001.xml (nn=00〜99) 中身は、次のよう…

まずは国土地理院の標高データを手に入れないと… 基盤地図情報ダウンロードサービス ここの「数値標高モデル」というのがそれらしい。 ダウンロードするには利用者登録しないといけないので、メールアドレスを登録して「ファイル選択」から地図上で欲しい地…

実は結構前に本を買っていて、ピラミッドを作るところまで読んで放っておいた。 やっぱり、自分でやりたい事がないとやる気にならない。 で、ようやくやりたい事ができた。 「国土地理院の標高データから実際の地形をマイクラ上に再現する」という事。 今後…

マイクラでPythonを使えるようにするのは私は本を見てやったので、それをそのまま紹介するわけにもいかないのでググって見つかったリンクを紹介しておきます。 www.timeless-education.com P.S. 2020/05/15 上記サイトでもリンク先が無効になっていたりした…