Spreadsheet layer

[YI]

[YI]

お世話になっております。QGIS3.18.

ShapefileをCSVでexportしてエクセルで編集して、またspreadsheet layerでimportすると、上記のようにシンボロジーの位置がおかしくなってしまいます。ミラーイメージのような感じになってしまっているようです。

Coordinates XとYも指定して、CRSも合わせているのですがどうしても適切な位置に反映されません。

どうしたらよいでしょうか？

2021年4月19日月曜日 16:37:53 UTC+2 YI:

[sakuya_izayoi]

YIさま

えっと、

①下の図のポリゴン重心計算した後に、CSVで出力して

②エクセルで編集して

③編集したCSVをQGISに取り込んだ

という流れでしょうか？

自分がやった感じですとミラーリングはしませんでした。

サンプルを添付します。（JGD2000の7系で作りました。）

ちょっと判断が付きにくいため、私と同じようなサンプルファイルを添付していただけると回答しやすいと思います。

申し訳ありませんが添付することはできますでしょうか？

よろしくお願いします。

2021年4月19日月曜日 23:41:44 UTC+9 YI:

[福岡]

YI 様

こんにちは

直感での投稿ですので、間違っていたらご容赦ください。m(__)m

地名をみると、南半球の地名のように思いました。

左右逆ということは、北半球のものとして地物が配置されたということのように思います。

2021年4月20日火曜日 9:42:11 UTC+9 sakuya_izayoi:

[ありた]

こんにちは

 

Coordinates XとYも指定して、CRSも合わせているのですがどうしても適切な位置に反映されません。

とありますが、具体的な状況や操作内容を記述してください。

0. XYZ Tiles の OpenStreetMap などを表示し、元データが正しい位置にあるか確認

1. 読み込んだ Shapefile の状況

　・レイヤプロパティの一般情報タブを開き「座標参照系 (CRS)」

　・フィールド計算機を開き

　　　　新しいフィールドを作るにチェック

　　　　仮想フィールド作成にチェック

　　　　出力する属性の名前：wkt

　　　　フィールド型：テキスト (string)

　　　　式： geom_to_wkt($geometry)

　・属性テーブルを開き wkt 列を確認（画像キャプチャ等）

2. エクスポートした際のパラメータ

　・座標参照系 (CRS)

　・レイヤオプションの GEOMETRY や SEPARATOR

3. エクスポートした CSV ファイル

　・メモ帳などで開く（画像キャプチャ等）

4. エクスポートした CSV ファイルを編集したあとの状態

　・メモ帳などで開く（画像キャプチャ等）

5. QGIS で CSV ファイルを読み込んだときのパラメータ

　・ファイル形式の区切り文字

　・ジオメトリ定義： X 属性、 Y 属性

　・ジオメトリの CRS

6. 読み込んだ CSV の状況

　・レイヤプロパティの一般情報タブを開き「座標参照系 (CRS)」

　・フィールド計算機を開き

　　　　新しいフィールドを作る

　　　　仮想フィールド作成

　　　　出力する属性の名前：wkt

　　　　フィールド型：テキスト (string)

　　　　式： geom_to_wkt($geometry)

　・属性テーブルを開き wkt 列を確認（画像キャプチャ等）

正しく設定できていれば正しく表示されるはずですので、どこかで不適切な設定が

あったか、もしくは CSV を編集した際にデータがおかしくなったかです。

CSV の編集を行わずに、そのまま QGIS に読み込んだ場合はいかがでしょうか。

[福岡]

YI 様

こんにちは

メールの方にご連絡頂きましたが、他の方へもわかるように掲示板の方へ投稿いたします。

なお、ご連絡頂いた内容ですが、地物はやはり南半球のものだったということ、その関連での対処方法はあるのでしょうかということでした。

情報が不足していて確実にこれだというようにお答えできませんが、私が疑っている内容を記載いたします。

さて私が想像しているのは、あらかじめ属性テーブルに記載されていた座標値をそのまま使っているのではないかということです。

属性テーブルのY座標が南緯の数値を正しいものであったとしても、プラス値で記載されていた場合、投稿なさったような現象が発生するように思いました。

この場合、単純に解決するのにはマイナス値にして取り込めば良いのですが、ちょっとズボラすぎるような感じです。

できれば

１　元になるShapeファイルが正しい位置に表示されているかについて、背面にWeb地図を表示して正しい位置に投影されているかを確認する。

２　フィールド計算機を使ってXY座標を自分で取得する。

３　csvデータは、取得した座標値を使ってQGISに取り込む。

４　取り込んだcsvデータも背面にWeb地図を表示して正しい位置に投影されているかを確認する。

ということはしておいた方が良いと思います。

背景地図はXYZTilesを使うのが一般的かと思うのですが、入門当初ならQuickMapServicesというプラグインを使うのが簡単なのでお勧めしています。

次の情報の中の「5.Webマップに重ねてみよう」でプラグインについて説明していますので、よろしければご覧ください。

参照座標系のざっくりとしたお話し

https://qiita.com/fukuoka1960/items/2318898e47ec1e749e1c

フィールド計算機では「$x」と「$y」という関数で、座標値を得ることができます。

次の情報の中の「フィールド計算機で座標を取得する」という項での説明が参考になると思います。

QGISでラインをポイントに変換してｘｙ座標を取り出す

https://qiita.com/ishiijunpei/items/a90385cca43ddcb91428

あくまでも原因を想像した上での投稿ですので、見当違いの内容でしたらお読み捨てくださいますようお願いいたします。

2021年4月20日火曜日 21:16:05 UTC+9 ありた:

[ありた]

メールにて添付資料込みで返信をいただきました。

いただいた GIS ファイルや操作画面画像を確認させていただいたところ、

行った操作はおそらく下記のとおりだと推察されました。

1. 街（？）ごとのポリゴンデータがある

　　・ EPSG:4326

　　・ coordx という属性があり 28° 52' 1,518" E のような形式の文字列

　　・ coordy という属性があり 2° 35' 31,395" S のような形式の文字列

2. ポリゴンデータを CSV 形式でエクスポートした

　　・レイヤオプションの GEOMETRY はデフォルト（指定なし）とし、

　　　位置情報の付加は行わない

3. CSV データを編集

　　・今回の内容には無関係

4. Spreadsheet Layers プラグインを用いて CSV を読み込み

　　・ジオメトリフィールドとして coordx / coordy を指定

ということだと思います。

（プラグインを使用したときは、その旨を明記してください。）

この上記4の工程で、指定したジオメトリフィールドの書式が NESW の

サフィックス付きの度分秒形式であったため、正しく認識できず、認識できた

箇所までの (28, 2) すなわち「東経28度　北緯2度」のように判定された

ものだと思います。

解消法としては、座標値の属性情報は、人が見て判別しやすい列であっても

構いませんが、それとは別にに、 28.86709 や -2.59205 のような度単位の

機械読み取り用の列を別途用意し、読み込む際にそれを座標値として

指定させます。

また、他の方法として、点の位置が既存の属性データ（都市の中心部等の

任意の地点）でなく、ポリゴンから機械的に算出したものでよければ

1. ポリゴンに対し、プロセシングツールの「重心」、「内部保証点」、

　　「到達不能極」などを実行し、代表点を算出する

2. 算出した点に対し、 CSV のエクスポートを行う。その際にオプションの

　　レイヤプロパティのジオメトリとして AS_XY を指定します

3. 必要に応じ CSV データを編集します

4. CSV を読み込む際に、ジオメトリフィールドとして X / Y を指定

という方法もあります。

[YI]

皆様

ご丁寧なご説明誠にありがとうございます。

おっしゃる通りCoordinates表記がQGISとCSVで微妙に変わっていました。

表記を度単位に変えようと思いますが、

Coordinates converterのようなサイトもしくは一括変換するためのExcelのformulaのようなものはご存じではないでしょうか？

もしくはほかのやり方がありますでしょうか？

下記添付のように１対1で変えるものは見つかったのですが個体数が増えるとそれはかなり時間的にも難しくなってしまいそうです。

ご教示いただけますと幸いです。

YI

2021年4月23日金曜日 14:29:21 UTC+2 ありた:

[福岡]

YI 様

こんにちは

「度＋分/60＋秒/360」で計算できると思いますので、Excelで式を組み立てることで対応できると思います。

ただ、ありたさんの情報をみると文字列でデータが整備されているので、文字関数を使って必要な値を抜き出す必要があります。

また例示されたものでは、秒の小数点をカンマにしているようでしたので、これを置き換える必要もありそうです。

文字列の関数としては「LEFT」、「RIGHT」、「MID」、「FIND」、「SUBSTITUTE」を使っています。

使い方については、検索すると出てくると思います。

ちょっと注意が必要なのは文字列は「"」で囲みますが、「"」を囲んだときの表現は「"""」ではなく「""""」になります。

サンプルがあった方が良いと思いましたので、Excelファイルを添付します。

状況に応じて、適宜修正してお使い頂ければと思います。

2021年4月26日月曜日 21:39:42 UTC+9 YI:

[YI]

福岡様、皆様

解決できました。

本当にありがとうございました。

今後ともよろしくお願いします。

YI

2021年4月27日火曜日 2:50:41 UTC+2 福岡: