コメントは受け付けていません。dataの表示(抽出)の基本を学んでおきます。Developmentを参考にして、実践的に使えるPythonでGeoSpatialの使い方を学んでいくものです。code等に関しては、あくまでも参考としてください。なにがあっても責任とれませんので。PythonとGeoSpatialに関するlibraryが設定できているもとして進めていきます。地理空間情報を扱う場合、多くはShapefileで提供されることが多いですが、Shapefilesは不適当であるか、能率が悪いことがあります。libraryで地理空間情報を扱ったり、保存したりできる独立したFormatを持つ必要があります。TextあるいはWKTと呼ばれます。WKTは、point、lineあるいはpolygonのような地理空間情報のobjectのコンパクトなテキストベースの記述です。バチカン市国のWKTはで表されます。WKTは、geometryがどんなに複雑になってもテキストベースです。WKTは、投影法、データおよび(または)座標系を包含する空間的基準を表わすためにも使用することができます。ここに、osgeoがあります。WKTに変換されたWGS84データを使用して、測地系を表わすosr.SpatialReferenceオブジェクトは、として表現されます。referenceに変換するとができます。libraryをつかってタイとミャンマーの国境のintersectionをやってみます。fileをclickしてdownload、その後解凍しておきます。私の開発環境では、f:¥shapefile_dataのfolderを作成して、解凍したfileをcopyしておきます。codeとりあえず、Ediorで下記のcodeを入力して、保存します。folderにcalcCommonBorders.pyで保存します。donefeature.Destroy()dstFile.Destroy()[/python]注意)上記codeで¥は全角ですので、半角にするかもしくは、バックスラッシュに置き換えてください。PyScriptでrunすると、何も起こりませんが、f:¥shapefile_data¥common-borderのfolderが作成されて、タイとミャンマーの接する国境のShapefileができあがります。(projection定義)の4つのfileです。border.prjをeditorで開いてみると、な感じでWGS84のprjが作成されていることが確認できます。QGISで確認すると、でタイとミャンマーの接する国境が作成されています。codeの説明changeDatum.pyが何を行っているのかを順にみていきましょう。libraryを使います。manualを参考にしてください。shapefile.GetLayer(0)Shapefileのlayerはただ一つのlayerですから、GetLayer(0)となります。shapely.wkt.loads(geometry.ExportToWkt())タイとミャンマーのfeatureをそれぞれthailandとmyanmarにgeometry.ExportToWkt()で追加します。thailand.intersection(myanmar)intersectionでcommonBorderを抽出します。spatialReference)spatialReferenceを”WGS84”として出力Shapefileを定義します。osgeo.ogr.Feature(dstLayer.GetLayerDefn())feature.SetGeometry(osgeo.ogr.CreateGeometryFromWkt(wkt))dstLayer.CreateFeature(feature)commonBorderをshapely.wkt.dumpsでwtkにdumpします。featureを定義して、wktからGeometryを作成してfeatureにsetします。dstLayerの出力Shqpefileにfeatureを作成します。feature.Destroy()dstFile.Destroy()順番にみていくと理解できると思います。dataの表示について、WKTの基本を学びました。manualのshapelyのexampleを実行して、shapelyでできることを確認すると良いと思います。Dataの保存をやってみます。Scriptの基本的なことはある程度理解している前提で今後も話を進めていきます。誤字、脱字、spell間違いや勘違いや翻訳間違いなど多々出てくると考えられます。それは違うじゃん!!とかいろんな意見をいただければと思います。Tutorialの目次に戻る。水文水理解析ソフトウェアの販売セール後1日です。水文水理解析ソフトウェアの販売PythonでMySQLを使ってみる。4]PythonでGUIを試してみる。第1弾ExcelでPythonを使ってみる。ExcelでPythonを動かす。3]PythonでExcelを読み書きする。PythonでGeocodingする。2]Pythonで水文水理解析をやってみる。AB-WebLog.com.。