コメントは受け付けていません。Dataの測地系と投影法の変更の基本を学んでおきます。Developmentを参考にして、実践的に使えるPythonでGeoSpatialの使い方を学んでいくものです。code等に関しては、あくまでも参考としてください。なにがあっても責任とれませんので。PythonとGeoSpatialに関するlibraryが設定できているもとして進めていきます。UTM座標の投影法の変換GISのデータは多くの測地系や投影法があり、いつも自分の使いたい座標系ではないことが多いですが、この測地系と投影法の座標変換を学んでおくことは非常に有用です。ここでは、UTM座標系をWGS84の経度緯度座標系に変換する例を学びます。downloadします。folderに解凍したfileをcopyします。copyするfileはmiami.shpmiami.shxmiami.dbfmiami_datum_reference.txtの4つのfileです。codeとりあえず、Ediorで下記のcodeを入力して、保存します。folderにchangeProjection.pyで保存します。done.srcFile.Destroy()dstFile.Destroy()[/python]注意)上記codeで¥は全角ですので、半角にするかもしくは、バックスラッシュに置き換えてください。PyScriptでrunすると、何も起こりませんが、f:¥shapefile_data¥miami-reprojectedのfolderが作成されて、投影法がWGS84に変換されたfileができあがります。(projection定義)の4つのfileです。miami.prjをeditorで開いてみると、な感じでWGS84のprjが作成されていることが確認できます。codeの説明changeProjection.pyが何を行っているのかを順にみていきましょう。dstProjection)入力ShapefileのprojectionをsrcProjectionで定義します。SetUTM(17)でUTM17系とします。出力ShapefileのdstProjectionをsrcProjectionで定義します。SetWellKnownGeogCS('WGS84')で世界測地系とします。transformでCoordinateTransformationを定義します。feature.Destroy()srcLayerのfeature数でloopします。featureとgeometryに入力Shapefileの情報をcopyして、newGeometryのcloneを作成して、Transform(transform)で投影法の変換を行います。属性値については変更はありませんので、そのまま出力Shapefileに追加していきます。feature.SetGeometryで変換後のgeomertyを追加して、dstLayerに新しくfeatureを作成して追加していきます。srcFile.Destroy()dstFile.Destroy()実行後、QGISで確認してみるとの感じで投影変換されていることが確認できます。順番にみていくと理解できると思います。objectの引数が不明な場合は、PyScriptのインタプリタでhelp(osr)などを入力すると、使用方法がインタプリタ画面に出力されます。dataのShapefileの投影法変換を行いました。ExercisesでJGD2000の投影法を使って、Shapefileの変換をtyrしてください。次回は、測地系の変換をやってみます。Scriptの基本的なことはある程度理解している前提で今後も話を進めていきます。誤字、脱字、spell間違いや勘違いや翻訳間違いなど多々出てくると考えられます。それは違うじゃん!!とかいろんな意見をいただければと思います。Tutorialの目次に戻る。水文水理解析ソフトウェアの販売セール後1日です。水文水理解析ソフトウェアの販売PythonでMySQLを使ってみる。4]PythonでGUIを試してみる。第1弾ExcelでPythonを使ってみる。ExcelでPythonを動かす。3]PythonでExcelを読み書きする。PythonでGeocodingする。2]Pythonで水文水理解析をやってみる。AB-WebLog.com.。