地理情報システム

データフレームのリストがあります。すべて同じ構造で、すべてポリゴンです。 ジオパンダを使用してそれらをすべて一緒に追加するにはどうすればよいですか？

15 python  geopandas 

QGISでは、道路名をできる限り行または列に揃えようとしています。Adobe Illustratorでマッピングプラグインを使用して実行できることは知っていますが、そのための$$はありません。 私はこのようなことを試みましたが、それは助けにならないようでした。 これは例です： 必ずしもそのようにする必要はありませんが、似たようなことをしたいと思います。

14 qgis  labeling 

QGIS 3をインストールして、Windows上でQGIS 2.xインストールと並行して実行できるかどうか知りたいですか？

14 qgis  installation  windows  qgis-3  compatibility 

Qgis 3.0の更新後、Qgis 3.0での処理スクリプトの記述に関する情報を見つけることは非常に困難になりました。 @Underdark（こちらを参照）は、スケルトンの基礎を提供しています。このコードは、テンプレート（Qgis 3.0.2）から新しいスクリプトを作成するときに、Qgisにも追加されたようです。 しかし、私のようなPython初心者がそのコードを変更する方法を理解するのを助ける方法を見つけることができませんでした。特に入力層と出力層については。 私の目標は、2つのラスターレイヤーとダブルを入力として、2つのレイヤーを出力するスクリプトを作成することです。 それを可能にするためにサンプルコードに必要な変更は何ですか？ Qgis 2.xの場合、次の構文を使用します。 ##Layer1=raster ##Layer2=raster ##myDouble=Double ##OutLayer1=output raster ##OutLayer2=output raster 私が理解していることから、次の手順で変更を行う必要がありますが、何を配置するのかわかりません。 def initAlgorithm(self, config=None): self.addParameter(QgsProcessingParameterFeatureSource( self.INPUT, self.tr("Input layer"), [QgsProcessing.TypeVectorAnyGeometry])) self.addParameter(QgsProcessingParameterFeatureSink( self.OUTPUT, self.tr("Output layer"), QgsProcessing.TypeVectorAnyGeometry)) 5月16日、Qgis python APIドキュメントがリリースされました。ただし、ここでの使用方法はまだ不明です。（Pythonの知識が不足している可能性があります）

14 qgis  pyqgis  qgis-3 

QGIS 2および3でベクトルポリゴンのレンガ壁スタイルを作成しようとしています。SVG塗りでレンガ壁のsvg描画をレンダリングできます（下の画像）。しかし、スタイルは無料の地質図の一部であるため、SVGを使用したくありません。SVGファイルは一般向けのデータの配布とレンダリングを複雑にします。デフォルトのQMLを介してSVGを保存すると、QMLファイル内のSVGファイルへの相対パスが強制され、SVGファイルはQGISシステムのSVGフォルダー内になければなりません。 誰でも次のようなポリゴンベクトルスタイルを作成しました。 レンダリングに外部SVGファイルは必要ありません。

14 qgis  style  qml 

上の米軍地図サービス（AMS C762） 1：アイスランド1948年から1951年製50,000マップ、ショーの度に、私は角度を疑い上部にマーキング規模があります。下の画像では、25°〜29°を示しています。それは左に増加します。これは、西経経度を指している可能性を示唆していますが、数値は正しくありません。マップは18°W–18°30'Wをカバーしています。それは、通常の方法でカラーの他の場所に示されている磁気偏角でもありません。シリーズ内の他のマップは同じ場所に同様の縮尺がありますが、数値が異なる場合があります。この角度スケールはどういう意味ですか？

14 maps  topography  angles  iceland 

sfギャップやスライバーを導入せずにポリゴンを単純化するにはどうすればよいですか？ たとえば、シェープファイルでは、次を使用しますrmapshaper::ms_simplify()。 library("pryr") library("rgdal") library("rmapshaper") download.file("https://borders.ukdataservice.ac.uk/ukborders/easy_download/prebuilt/shape/England_gor_2011.zip", destfile = "regions.zip") unzip("regions.zip") regions <- readOGR(".", "england_gor_2011") object_size(regions) # ~13MB regions <- ms_simplify(regions) object_size(regions) # < 1MB 私はsf::st_cast()、manページから次のように述べました： ジオメトリを別のタイプにキャスト：単純化するか、明示的にキャストします そして： 引数へ：文字; 対象タイプがない場合、簡素化が試行されます。xがsfgタイプ（つまり、単一のジオメトリ）の場合、指定する必要があります。 私がto行方不明になったとき、これは期待通りに機能しませんでした（それが本当であるには余りにも良いことを知っていました！）： library("sf") regions <- sf::read_sf("england_gor_2011.shp") object_size(regions) # ~13MB regions <- sf::st_cast(regions) object_size(regions) # Still 13MB 現在、でファイルを開きrgdal::readOGR()、簡略化して保存し、で再度読み込みsfます。 もっと良い方法はありますか？ rgeos::gSimplify() @skの提案はrgeos::gSimplify()、次の引数で指定された場合、トポロジを考慮した単純化（つまり、スライバを作成せずに単純化）を実行できます。 library("rgeos") regions_gSimplify …

14 r  simplify  sf 

私がやろうとしているのは、ポリゴンに100 mのバッファーを置くことですが、「カラスが飛ぶように」バッファーを固定幅にするのではなく、Lidar DEMの傾斜に追従させたいのです。 ArcMap（3D Analystを使用）があります。悲しいことに、ArcMapにはSpatialまたはNetwork Analyst Extensionはありません。 近接分析、経路距離分析、および異方性累積コストを検討しているインターネットによって、私はすべて好転しました。どういうわけか、時間とコストのグリッドを作成し、ポリゴンから100 mに到達するのにかかる時間を調べる必要があるように見えますが、それはむしろ逆に思えます。これらのツールの大部分は、2つのポイント間の最も簡単なルートを見つけるために作成されたようです。 誰もこれをやったことがありますか？

14 arcgis-desktop  distance  elevation  3d-analyst  cost-path 

リモートセンシングドメインにおけるピクセルベースとオブジェクトベースの分類の違いを明確に理解するのに苦労しており、このコミュニティの誰かが洞察を提供できることを期待しています。 私がこれまでに得た情報に基づいて、私の現在の理解はこれらの線に沿っています： ピクセルベースの分類： 分類はピクセルごとに行われ、その個々のピクセルで利用可能なスペクトル情報のみを使用します（つまり、局所性内のピクセルの値は無視されます）。この意味で、各ピクセルは分類アルゴリズムのトレーニング例であり、このトレーニング例はn次元ベクトルの形式になります。nは画像データのスペクトルバンドの数です。したがって、訓練された分類アルゴリズムは、画像内の個々のピクセルごとにクラス予測を出力します。 オブジェクトベースの分類： 分類は、各ピクセルが互いに関連しているときの各ピクセルの空間特性を考慮して、ローカライズされたピクセルのグループで行われます。この意味で、分類アルゴリズムのトレーニング例はピクセルのグループで構成され、トレーニングされた分類アルゴリズムはそれに応じてグループベースでピクセルのクラス予測を出力します。大まかな例では、画像を同じサイズのn個のセグメントに分割し、各セグメントにクラスを指定します（つまり、オブジェクトを含む/オブジェクトを含まない）。 この考えはこれらの用語の意味に関して正確ですか、それとも私が見逃したものがありますか？

14 remote-sensing  classification  machine-learning 

次の図のように、最長軸に垂直な（つまり、中間点の幅を横切る）中間点で約4000個のポリゴンを分割しようとしています。 理想的には、これを自動的に行い、各ポリゴンを手動で分割しないようにします。それぞれに描画できる最長の線を変換してポリゴンの中間点を抽出しましたが、この点を横切る線を自動的に描画する方法を決定するだけです。 ポリゴンの幅はさまざまであるため、特定の長さの幅の線を定義してポリゴンを分割するツールは、私が探しているものではありません。 何か案は？

14 arcpy  polygon  splitting 

USDA（National Map経由）は、米国のほとんどの州に4バンドNAIP ArcGIS "ImageServer"（サービス）を提供します。ここで説明するように、QGISのWMSサービスに正常に接続できます。これにより、デフォルトのRGB（1,2,3）バンドを使用してQGISに画像が読み込まれます。 ただし、CIR（バンド4、1、2）を表示するように画像を変更したいと思います。[スタイル]タブの[レンダリングタイプ]には、[シングルバンドカラーデータ]オプションのみがあり、ここの典型的なラスターの手順3.5で説明されている[マルチバンドカラー]オプションはありません。これはWMSサービスの制限のようです。 RESTエンドポイントを使用して、ArcGIS Desktopの同じサービスからこれを行うことができます。サービスのバンドを変更するArcGIS Online Webマップの例は、RESTエンドポイントを使用して利用できます。レイヤーの「...」に移動し、「画像表示」、「ユーザー定義レンダラー」の順に選択します。 更新： ここから、Chromeの開発者ツールとArcGIS.comマップサンプルをいじくり回してきましたが、小さな情報が見つかったと思います。画像の表示設定を変更した後、クエリパラメータが次のように設定されていることに気付きました。 https://gis.apfo.usda.gov/arcgis/rest/services/NAIP/Montana_2015_1m/ImageServer/exportImage?f=image&bandIds=3,0,1&bbox=... バンドIDを再度変更した後、このクエリパラメーターのバンドIDは1ベースではなく0ベースであることがわかりました。したがって、＆bandIds = 3,0,1は、CIR形式で画像を表示するために使用されるバンド4,1,2を実際に参照しています。そのため、QGISがこのパラメーターを受け入れるようにクエリ文字列を構成する方法を理解する必要があります。何かご意見は？

14 qgis  configuration  image-service 

ラベルを付ける2つの機能が互いに近すぎる場合、基本的に3つのオプションがあります。デフォルトでは、QGISはラベルの1つを省略します。問題のレイヤーまたはグローバルラベリング設定のいずれかで「すべてのラベルを表示」プロパティを有効にすると、重複するラベルが印刷されます。 ラベルの位置を保存するデータ列を割り当てていない限り、ラベルの1つをより良い位置に移動することで競合を解決できます。 私の問題は、ラベルの省略もラベルの重複も適切な解決策とは考えていないことです。競合するすべてのラベルを体系的に反復して、ラベルの1つを移動するか、おそらくそれらの1つを非表示にすることで、競合を手動で解決する方法がありますか？ 配置は縮尺と地図の範囲に依存しているように見えるため、これは印刷コンポーザの特定の地図に固有のものであると思われます。

14 qgis  labeling 

ボロノイダイアログで「バッファ」が何をしているのか、つまり数字が何を表しているのか説明できますか？

14 qgis  buffer  voronoi-thiessen 

RのHistDataパッケージ（https://r-forge.r-project.org/R/?group_id=574）には、1854年のロンドンでのコレラの発生に関するJohn Snowのマップに関連するデータセットがあります。ウォルター・トブラーの監督の下で慎重にデジタル化された権威です。これらのデータセットの詳細については、John Mackenzieがhttp://www1.udel.edu/johnmack/frec480/cholera/cholera2.htmlで説明しています。 残念ながら、死亡、ポンプ、および道路の座標は、他のGISアプリケーションまたはRのマッピングソフトウェア（空間パッケージ、ggmapなど）に適したマップ座標ではなく、任意の座標系を使用します でhttp://freakonometrics.hypotheses.org/19213アーサー・シャルパンティエはからスノーデータのバージョンでggmapを使用しています http://www.rtwilson.com/downloads/SnowGIS_v2.zip。Cholera_Deaths.shpただし、このファイルには、私が記録した578人ではなく、489人の死亡のみがリストされていHistData::Snow.deathsます。 1つのアイデアは、平均と（x、y）座標の標準偏差との関係を見つけて線形に再スケーリングすることですが、おそらくもっと良い方法がありますか？ ここに私が今まで試したものがあります > data(Snow.deaths, package="HistData") > D <- Snow.deaths[,2:3] > colMeans(D) x y 13.03312 11.69721 > var(D) x y x 3.8150987 0.3802654 y 0.3802654 2.7213828 Cholera_deathsファイルを読む > folder <- "C:/Dropbox/R/data/Snow/SnowGIS_v2/SnowGIS" > library(maptools) > deaths <- readShapePoints(file.path(folder, "Cholera_Deaths")) > head(deaths@coords) coords.x1 coords.x2 0 529308.7 181031.4 …

14 r  coordinate-system  scale 

QGIS 2.14.5 LTR（および最近のQGIS 3.0）では、各値の範囲に目盛り付きレンダラーとベクトルフィールドマーカーを使用して美しいマップデザインを設定しました（矢印の長さは、0.001から〜 5m / s）。 これはうまくいきましたが、伝説の矢印を短くする設定を逃したのではないかと思います。「通常の」幅のレイヤーウィンドウでは、クラス値は表示されません。 単位：ミリメートル 私の最初の推測は、それがsthを持っているということでした。矢印の長さ、線幅の単位をミリメートルからマップ単位に切り替えることに関係していますが、これは凡例に影響しません。 単位：マップ単位 更新 問題はQGIS 3.0でも残っています。 これを視覚化するスクリーンキャストを作成しました：https : //youtu.be/oxJ3tGXFB18 プロパティ/シンボルで、シンボルは正しくレンダリングされます。 レイヤーを複製するか、レイヤーパネルの幅が変わらない限りスタイルに変更を適用すると、レイヤーリストにマーカーが正しく表示されることに注意してください。

14 qgis  legend  markers  qgis-3  vector-field