地理情報システム

コードをspパッケージから新しいsfパッケージに移行しています。以前のコードでは、SpatialDataFrame（censimentoMap）とSpatialPointDataFrame（indirizzi.sp）のポリゴンがあり、以下の手順で各ポイントのポリゴンセルID（Cell110）を取得しました。 points.data <- over(indirizzi.sp, censimentoMap[,"Cell110"]) 実際に、2つのsfオブジェクトを作成しました。 shape_sf <- st_read(dsn = shape_dsn) shape_sf <- st_transform(x=shape_sf, crs=crs_string) そして indirizzi_sf = st_as_sf(df, coords = c("lng", "lat"), crs = crs_string) そして、上記の命令に相当するsfを探しています... ids<-sapply(st_intersects(x=indirizzi_sf,y=shshape_sfpeCrif), function(z) if (length(z)==0) NA_integer_ else z[1]) cell_ids <- shape_sf[ids,"Cell110"]

16 r  overlay  sp  sf 

ArcGISのフィーチャレイヤーをGeoPackageに直接変換することはできますか？ ArcGIS Desktop 10.5のArcToolboxの[変換ツール/ジオパッケージへ]に表示されるオプションは、[ジオパッケージへのラスターの追加]のみです。「機能の追加」などはありません。 で、このブログの記事、ほぼ3年前から、ESRIは、「10.2.1で、または10.2.2のArcGISデスクトップと、あなたは空のGeoPackageを作成することができますし、それに特徴データをコピーすることによってGeoPackageを移入」と主張しました。しかし、私はどのように見ることができません。ArcCatalogには、GeoPackageファイルを作成するオプションはありません。

16 arcgis-desktop  convert  arccatalog  geopackage  arcgis-10.5 

いくつかのポリゴン（境界）のクリーンアップを行っていますが、選択したポリゴンのセットまたは単一のポリゴンにノード/頂点を追加できるかどうか疑問に思っています。 ところで、QGISトポロジ編集ツールは非常にシンプルです。非常にシンプルですが、非常に正確であり、そのような手作業を行うことから多くの作業を奪います。 QGISには、選択したノードを削除する機能があります-驚くべきことには、余分なノードをクリーンアップするときにトポロジが保持され、さらに改善されます。素晴らしいもの！

16 qgis  editing  vertices 

postgisにラインストリングテーブルとポイントテーブルがあります。 任意のポイントに最も近い線を知っています。私が知っておく必要があるのは、その行のどの「サイド」がポイントであるかです。与えられたポイントからライン（ライン上の最も近いポイント）に垂直なラインを作成し、座標を比較することでそれを行う必要があると思いますが、それを行う方法は正確にはわかりません。線が変化するため、方向が変わります。 私は自分の仕事を説明するために写真を作りました。 線自体は黒で、その方向は緑の矢印で示されています。ポイントテーブルに「サイド」列を追加する必要があります。そのため、赤いポイントの値は「right」、青いポイントの値は「left」になります。 誰かがポイントの「サイド」値を計算するSQLコードの例を提供できますか？

16 postgis  point  linestring 

QGISのラスター内のセル値のカウントの概要を取得するにはどうすればよいですか？レイヤーのプロパティに要約統計とヒストグラムを提供するために、QGISがこれを行っていることを知っています。要約統計もヒストグラム画像も必要ではなく、個々の値のカウントが必要です。 ある時点で、このための特別な機能を追加する試みがあったようです（機能＃4321、http：//hub.qgis.org/issues/4321）。 私は本当に明白な何かを見逃していますか？

16 qgis  raster 

Pythonのgeopyパッケージには、Great CircleとVincentyの公式という 2つの距離測定技術が備わっています。 >>> from geopy.distance import great_circle >>> from geopy.distance import vincenty >>> p1 = (31.8300167,35.0662833) # (lat, lon) - https://goo.gl/maps/TQwDd >>> p2 = (31.8300000,35.0708167) # (lat, lon) - https://goo.gl/maps/lHrrg >>> vincenty(p1, p2).meters 429.16765838976664 >>> great_circle(p3, p4).meters 428.4088367903001 違いはなんですか？どの距離測定が望ましいですか？

16 python  distance  great-circle  geopy  vincenty-formulae 

Pythonで作成されたQGISプラグインの自動テストを作成するためのアドバイスを探しています。 過去にPyUnit（unittestモジュール）を使用してPythonスクリプトのテストを作成しましたが、GUIを使用するアプリケーションではテストを実行したことがありません。PyQt4.QTestを使用してQtウィジェットでユニットテストを行う方法を説明するページを見つけました（http://www.voom.net/pyqt-qtest-example） QGIS内から実行するように設計されたウィジェットを使用します。 PyQGISドキュメントの「テスト」のセクションは特にありません。 私がこれまでに持っているものは： 分離されたモジュールまたは関数でデータの実際の処理を維持し、それらのユニットテストを記述します。 QTestを使用してUIの基本的なテストを実行します。 QGIS内からプラグインを使用するときにすべてが一緒に保持されるように。 もっと良い方法はありますか？

16 pyqgis  qgis-plugins  unit-testing 

featureCollection要素を親とするproperties要素を持つことは有効ですか？ geojson.orgによると、これは有効です。 { "type": "FeatureCollection", "features": [ { "type": "Feature", "geometry": {"type": "Point", "coordinates": [102.0, 0.5]}, "properties": {"prop0": "value0"} } ] } しかし、私はそれが有効であるか、これを持っていることが無効であるかどうかを見つけることができません： { "type": "FeatureCollection", "properties" : { "description" : "This is the geometry for..." } "features": [ { "type": "Feature", "geometry": {"type": "Point", "coordinates": [102.0, 0.5]}, "properties": {"prop0": …

16 geojson 

閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか？この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 2ヶ月前に閉店。 GeoTIFF形式の多数のファイルを含むフォルダーがあります。 これらの各イメージに対してPRJファイルとTFWファイルの両方を生成する必要があります。 これを行う方法はありますか？

16 geotiff-tiff 

与えられたgeopandas GeoDataFrameポリゴンのシリーズを含む、私は私のリストにある各機能のキロ平方の領域を取得したいと思います。 これはかなり一般的な問題であり、過去に通常の提案された解決策を使用するようになっているshapelyとpyproj（例えば、直接こことここ）。 これを純粋に行う方法はありgeopandasますか？

16 python  area  shapely  geopandas  pyproj 

QGISの起動中、スプラッシュスクリーンの下部に「ロードされたプラグインを復元しています」のようなステータスメッセージが表示されます。 私は、startup.py機能を使用しており、この機能を使用して、起動スクリプトのどの部分が現在実行されているかをユーザーに通知します。 スプラッシュスクリーンにこの情報を表示することは可能ですか？ 編集1： 回避策として、起動時に自分のスプラッシュスクリーンを使用することができました： from qgis.gui import * from qgis.utils import * from qgis.core import * from PyQt4.QtGui import * from qgis.PyQt.QtCore import QSettings, Qt import time template=QgsApplication.qgisSettingsDirPath() + "python/" app=QgsApplication.instance() splash_pix = QPixmap(template+'splashscreen.png') splash = QSplashScreen(splash_pix, Qt.WindowStaysOnTopHint) splash.setWindowFlags(Qt.WindowStaysOnTopHint | Qt.FramelessWindowHint) splash.setEnabled(False) splash.setMask(splash_pix.mask()) progressBar = QProgressBar(splash) progressBar.setMaximum(10) progressBar.setGeometry(0, splash_pix.height() …

15 qgis  pyqgis  pyqt4 

PyQGISでは、レイヤーパネル（または目次）のグループを太字にするなど、さまざまな設定を変更できます。 from PyQt4.QtCore import QSettings QSettings().setValue("/Qgis/legendGroupsBold", True) 設定は、オプションまたは入力から確認できるように書き込まれます。 QSettings().value("/Qgis/legendGroupsBold") >>> u'true' しかし、それは適用されません。私は使用すると思った QSettings().sync() 助けになるかもしれませんが、残念ながらそうではありません。適用する唯一の方法は、GUI（Settings > Options）に移動し、OKQGISをクリックするか再起動することです。 PyQGISを介してこれを行う方法はありますか？ コンテキスト： これをプロジェクトファイルにマクロとして追加して、グループ名のみが太字になり、レイヤー名が太字にならないように（美的目的で）します。したがって、プロジェクトがロードされたときにQSettingsが適用されると理想的です（プロジェクトが閉じられたときに初期設定が再適用されます）。

15 pyqgis  table-of-contents  macro 

MySQLのドキュメントでは ST_Distance_Sphere 計算では、球状の地球と構成可能な半径を使用します。オプションの半径引数はメートル単位で指定する必要があります。省略した場合、デフォルトの半径は6,370,986メートルです。radius引数が存在するが正ではない場合、ER_WRONG_ARGUMENTSエラーが発生します。 PostGISはのドキュメントで述べてST_Distance_Sphereいます（ドキュメントはもう正確ではありませんが） 球状の地球と6370986メートルの半径を使用します。 どこからデフォルトの6,370,986メートルを取得しましたか？ WGS84によると、長軸半径は6,378,137.0 mです。現在、平均半径を使用しているPostGISは、基本的に6371008を使用しています。 コードを見る #define WGS84_MAJOR_AXIS 6378137.0 #define WGS84_INVERSE_FLATTENING 298.257223563 #define WGS84_MINOR_AXIS (WGS84_MAJOR_AXIS - WGS84_MAJOR_AXIS / WGS84_INVERSE_FLATTENING) #define WGS84_RADIUS ((2.0 * WGS84_MAJOR_AXIS + WGS84_MINOR_AXIS ) / 3.0) つまり -- SELECT 6378137.0 - 6378137.0 / 298.257223563; WGS84_MINOR_AXIS = 6356752.314245179498 -- SELECT ( 2.0 * 6378137.0 + …

15 postgis  postgresql  mysql  mysql-spatial 

QGISでは、形状の塗りつぶしをラベル値のパターンの繰り返し（つまり、アイテムの属性）に設定したいと思います。 これは可能ID = 1,2,3ですか、下のスケッチに似ています（を持つ図形の場合）？

15 qgis  labeling  style  fields-attributes  fill 

QGIS 2.14を使用して、ルールベースのラベル付けがあることを確認してください。大文字の町のラベル（CRAWLEYなど）を表示し、小文字の町のラベル（Faygateなど）は表示しない 属性テーブルは、名前以外のさまざまな町を区別しません。これを行う方法はありますか？

15 qgis  labeling  qgis-2.14