🎄teratailクリスマスプレゼントキャンペーン2024🎄』開催中!

\teratail特別グッズやAmazonギフトカード最大2,000円分が当たる!/

詳細はこちら
Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Q&A

解決済

2回答

1917閲覧

geopandasでshpファイルから読み込んだポイント座標(ジオメトリ)回転させたい。

melo_yuya

総合スコア16

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

0グッド

0クリップ

投稿2021/03/08 10:39

前提・実現したいこと

農地データと人工衛星のマイクロ波の照射方向の関係性について研究しています。
geopandasによってshpファイル(ポリゴン)を抽出しました。
そこでジオメトリのポイント座標を演算のために抽出し、
その座標を原点回りに回転させて表示しようと考えています。

実行環境

実行環境はWindows10、Annaconda(Spyder)を使用しています。

該当のソースコード

Python

1import geopandas as gpd 2import re 3 4#','のところで区切る 5def convert_coordinates(x): 6 s = re.sub('(^[^(]*(*)|()*)', '', x.wkt) 7 return [[float(x) for x in y.strip().split(' ')] for y in s.split(',')] 8 9#shpファイルの読み込み 10yoshi = gpd.read_file("yoshikawa_2019.shp") 11print(yoshi) 12 13# change the geometry 14yoshi_points = yoshi['geometry'] 15print(yoshi_points) 16 17#a度回転 18a = 10 19yoshi_rotate = yoshi_points.rotate(math.radians(a)) 20print(yoshi_rotate)

出力結果

以下が上記のコードの出力結果です。

Python

10 POLYGON Z ((273474.470 4123032.881 0.000, 2734... 21 POLYGON Z ((266669.419 4122033.251 0.000, 2666... 32 POLYGON Z ((266698.965 4122029.388 0.000, 2667... 43 POLYGON Z ((266659.482 4121927.719 0.000, 2666... 54 POLYGON Z ((266626.999 4121931.518 0.000, 2666... 6 73693 POLYGON Z ((266500.679 4121835.558 0.000, 2665... 83694 POLYGON Z ((266524.600 4121831.832 0.000, 2665... 93695 POLYGON Z ((266686.077 4121924.596 0.000, 2667... 103696 None 113697 None 12Name: geometry, Length: 3698, dtype: geometry 130 POLYGON Z ((273474.328 4123032.893 0.000, 2734... 141 POLYGON Z ((266669.263 4122033.226 0.000, 2666... 152 POLYGON Z ((266698.811 4122029.362 0.000, 2667... 163 POLYGON Z ((266659.326 4121927.703 0.000, 2666... 174 POLYGON Z ((266626.843 4121931.495 0.000, 2666... 18 193693 POLYGON Z ((266500.523 4121835.542 0.000, 2665... 203694 POLYGON Z ((266524.444 4121831.809 0.000, 2665... 213695 POLYGON Z ((266685.919 4121924.570 0.000, 2667... 223696 None 233697 None

質問

確かに座標は少し変わっているのですが、回転しているのか分からない状態です。
もし分かる方がいましたら、教えていただきたいです。

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答2

0

自己解決

地理参照系について調べたところ、赤道と子午線の交点であり、このコードで回転していることが分かりました。

投稿2021/06/29 02:40

melo_yuya

総合スコア16

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

0

plotしてみればいいんじゃないでしょうか。
plots

具体的なコードだと

python

1base=yoshi.plot(color='red') 2yoshi_rotate.plot(ax=base,color='blue')

(補足1)
rotateの角度の指定は、デフォルトではdegreesとなってるようです。
radiansでするには、use_radians=Trueのオプションをrotateメソッド呼び出し時に追加する必要があるみたいです。
geopandas.GeoSeries.rotate

python

1# radiansなら指定が必要 2yoshi_points.rotate(math.radians(a),use_radians=True) 3# degreesなら指定なし 4yoshi_points.rotate(a)

(補足2)
原点回りに回転ということであれば、originパラメータも必要では。

python

1yoshi_points.rotate(a,origin=(0, 0))

geopandas勉強中の為おかしなところがあったらすみません…

投稿2021/03/09 20:49

編集2021/03/09 21:42
xail2222

総合スコア1508

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.36%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問