地理情報システム

Pythonコンソールを使用して、ユーザー定義のプロジェクト変数の値を変更したい（プロジェクトの設定|変数で手動で編集できます）。QgsExpressionContextScopeクラスのsetVariable（）関数を追跡しましたが、プロジェクト設定の変数を実際に変更することに成功していません。これまでの私のコード： iface.mapCanvas().mapSettings().expressionContext().scope(0).setVariable('myvar',1) 私はさまざまな表現のコンテキストで迷子になっていると思います...

10 qgis  pyqgis 

私は、scipy-learn PythonパッケージのBirchアルゴリズムを使用して、10のセットで1つの小さな都市のポイントのセットをクラスター化しています。 私は次のコードを使用します： no = len(list_of_points)/10 brc = Birch(branching_factor=50, n_clusters=no, threshold=0.05,compute_labels=True) 私の考えでは、私は常に10ポイントのセットで終わります。今の場合、私はクラスターに650ポイントあり、n_clustersは65です。 しかし、私の問題は、しきい値が低すぎると、クラスターのアドレスが1つになり、非常に大きなしきい値（クラスターあたり40アドレス）になることです。 ここで何が悪いのですか？

10 python  clustering  scipy.spatial 

私の国の地図とさまざまなファイル（.tif、shp）をレイヤーとして表示するWebマップアプリケーションの作成に興味があり、Leafletの上に構築することにしました。 だから私が始めたとき、私はOpenStreetMapsからの標準マップを使用して、次のコードを使用しました： var map = L.map('map').setView([50, 9], 10); L.tileLayer('http://{s}.tile.osm.org/{z}/{x}/{y}.png', { attribution: '© <a href="http://osm.org/copyright">OpenStreetMap</a> contributors', }).addTo( map ); しかし、私の監督者は地図が「詳細すぎる」と言ったので、詳細が低く、さらに役立つかもしれない別のベースマップを探しています。 それで、私が使用できるOpenStreetMapsのような無料の地図の代替ソースはありますか？

10 openstreetmap  leaflet  basemap 

QGISを使用して保護地域の地図を作成しています。私はバージョン2.12を使用していますが、この質問はおそらくバージョン2.14にも適用されます。 このアトラスには、いくつかの興味のあるポイント（poi）オブジェクトが表示されます。駐車場、インフォメーションステーション、視点など、これらのいくつかはアイコンとして表示されます。これらが互いに近い場合は、アイコンをグループ化し、それらを（線上で）ポイオブジェクトの中心点の近くに配置することが適切です。現在、これはデスクトップパブリッシングツールを使用して行われ、次のようになります。 このプロセスをできる限りQGISアトラスシステムに取り入れたいと思います。私たちは多くの保護地域を持っているため、プロセスを自動化することで得られる利益は潜在的に大きくなります。 QGISでこれをどのように行うことができますか？ また、FMEを使用してPOIを1つのポイントに集約してコンテキストを確認する質問も参照してください。

10 qgis  cartography  point-of-interest  icon 

いくつかのフライトを示すデータを含むマップを作成しようとしていますが、グレートサークルアークを使用して、出発地と目的地を接続したいと考えています。 基本的に、私は有名なFacebookマップに似たようなことをしたいと思います。 私はこの投稿で提供された関数を使用しました：https : //gis.stackexchange.com/a/5205/442（つまり、このブログ記事：http : //anitagraser.com/2011/08/20/visualizing-global-connections /）そして線を引くことはできましたが、それらは極地で国際日付変更線を横切るだけでなく、束になっています。 @underdarkは、リンクされたブログポストでこれらの行を分割する必要があると述べていますが、PostGISでそれらを自動的に分割する方法がわかりません。 さらに、極付近の線の集まりも解決する必要があります。 これらのフライトの始点と終点のポイント位置がある場合、これらの両方をどのように行うのですか？

10 postgis  great-circle  geographiclib 

Strahlerのストリームの順序付け方法に従っています。 この方法によると、緑色のストリームは4次ですか？

10 stream-order  strahler 

QGISでPythonスクリプトを作成し、いくつかのデータ（この例ではshp + tif）を取得し、（ファイルから）tamplate Map Composerを使用して、作成したコンポジションをpng画像にエクスポートします。 プログラミングの経験はほとんどありません（基本的なPythonの知識はありません）。いくつかのコードスニペットをググって、それらを一緒に動作させるようにしています。以前に回答したいくつかの質問からMap Composerコードを取得しました 。Pythonを使用して（表示レイアウトの何も変更せずに）Print / Map QGIS composerビューをPNG / PDFとして保存しますか？ データとMap Composerテンプレート（単一のマップと凡例項目が定義されている）をなんとかロードしましたが、エクスポートpng画像には空のマップフレームがあります（フレーム内にベクター/ラスターデータはありません）。ただし、凡例のアイテムは問題ありません。 このコードを機能させるための助けはありますか？ from qgis.core import * import qgis.utils from PyQt4 import QtCore, QtGui from qgis import core, gui # ADD VECTOR LAYER data_folder = "D:/QGIS/dane/" granica = "granica_SZ VI_UTM34.shp" granica_name = granica[0:-4] granica = data_folder …

10 qgis  python  print-composer  png 

継続的に接続されるはずのラスターをいくつかインポートしましたが、その代わりに、それらのラスター間に境界付きグリッドが作成されます。これは個々のラスター間のキャリブレーションの問題だと思うかもしれませんが、データのように実行する方法はわかりません。 すべてのラスターの全体がよりスムーズになるように、誰かが以前にこのようなことを経験していて、おそらく可能な解決策を知っているでしょう。 これらはすべて地理的に接続されているはずです。したがって、正方形の「境界線」の問題。 WindowsでQGISを使用しています。

10 qgis  raster  windows  smoothing 

一見シンプルなツールをWebAppBuilderマップに組み込むことを検討しています。 いくつかの基本： コンテンツはArcGIS 10.3.1で作成されます コンテンツはキャッシュサービスとしてArcGIS 10.3.1 Serverに公開されます マップはArcGIS Onlineで開発されています マップはWebAppBuilderにプッシュされます 現在、標準のウィジェットが使用されています（開発者バージョンのWebAppBuilderがあるため、必要に応じてコードを操作してカスタム作成できます） タスク：WebAppBuilderでインタラクティブに機能の選択機能を複製します。 これまで、ユーザーは「空間フィルターを使用して機能を制限する」を使用して機能をインタラクティブに選択でき、次に「ユーザー定義領域に触れる機能のみ」を使用して、さまざまなクエリウィジェットを使用してみました（下の画像を参照）。 ただし、これは私が求めているもののバリエーションであるため、成功したとは言えません。ユーザーがジオメトリ（ポリゴン、ライン、ポイント）を選択して、マップ上でハイライト表示できるようにしたいだけです。 おまけとして、私はユーザーに選択を強調表示する色を制御できるようにしたいと思っています-選択オプションと同様です。（下の画像を参照）

10 arcgis-desktop  arcgis-server  arcgis-10.3  arcgis-online  arcgis-web-appbuilder 

を使用してポリゴン間の小さな「断片」を排除する方法はありRますか？望ましいソリューションはSpatialPolygonsDataFrame、ポリゴン間の共有境界が一致する新しいものを作成します。RArcMapやQGISではなく、を使用するソリューションに特に興味があります。 そもそもなぜこれらのギャップが存在するのかについての説明も聞きたいです。 これは、私が使用している空間データの再現可能な例です。 library(rgdal) library(sp) library(tigris) library(magrittr) library(leaflet) library(gplots) # This project will use WGS 84 projected coordinate system crs_proj <- CRS("+init=epsg:4326") # These are the FIPS codes of the specific block groups in my study area sel <- c("530330079005", "530330079001", "530330079004", "530330085002", "530330085003", "530330086003", "530330087003", "530330085001", "530330090001", "530330091001", "530330091002", …

10 polygon  leaflet  r  tiger  slivers 

PythonとArcpyを使用してシェープファイルのサイズを取得することは可能ですか？もしそうなら、どうですか？

10 arcpy  shapefile  file-size 

QGISバージョン2.8.2でpostGISポイントデータセットを分析しています。合計約200kのフィーチャです。加重平均値を計算するには、逆距離加重補間手法を使用しています。ただし、出力に表示される「高い」値は、意味のある空間パターンの基礎ではなく、少数の外れ値の結果であることに懸念があります。 SAGA IDWプラグイン-または別のIDWツール-に、ラスターサーフェスが加重平均を計算するために必要な最小ポイント数を設定する方法はありますか？

10 qgis  saga  inverse-distance-weighted 

ArcGIS Proタスクは、マップフレームで動作するPython Toolboxツールをサポートしていますか？ 私が尋ねる理由は、私が次のことをしようとするときに立ち往生していることです： ArcGIS Pro 1.1.1を起動します 新しいプロジェクトを作成します-私は自分のTestProjectを呼び出し、C：\ Tempに配置しました プロジェクトペインを使用して、Natural Earthから世界の国々のシェープファイルがある場所にフォルダー接続を追加します（） ne_10m_admin_0_countries.shpをマップにドラッグアンドドロップして、ne_10m_admin_0_countriesというレイヤーを作成します 新しいレイアウトを挿入する-A3ランドスケープを使用しました レイアウトに新しいマップフレームを挿入する プロジェクトペインで、TestProjectフォルダーに新しいPythonツールボックスを作成します-私はTestPYTと呼んでいます プロジェクトペインでTestPYTを右クリックして編集します 以下のコードでコードを置き換えて、Python Toolboxにチリとスイスと呼ばれる2つのツールを提供します。 スクリプトを保存してPythonツールボックスを更新し、2つの新しいツールを確認します チリツールを実行して、レイアウトにマップフレームを表示してチリにズーム スイスツールを実行して、レイアウト上のマップフレームを表示しますスイスにズームします 新しいタスクアイテムを挿入する タスクアイテム1に新しいタスクを挿入し、チリと呼びます チリタスクで新しいステップを挿入し、それをZoom to Chileと呼びます ステップ動作の場合、それを自動で非表示にします [アクション]タブで、コマンド/ジオプロセシングをジオプロセシングツールとして設定して、チリツールを選択しようとします OKを選ぶとくっついてしまうようです 完了をクリックするとツールが「失われる」ようです 具体的には、作成しようとしているのは、チリにズームまたはスイスにズームするためにクリックできる2つのタスクを含むワークフローですが、上記のステップ19で立ち往生しています。 私が全体的にやろうとしているのは、ArcPy（ArcGIS 10.xアーキテクチャ用）のPythonアドインツールバーに相当するArcPy（ArcGIS Pro用）を見つけて、それらの国にズームするための2つのボタン（チリとスイス）を備えていることです。 私はこの手順を数回実行しましたが、ある場合にはチリとスイスのツールをタスクとして固定することができましたが、それでもそれらはマップフレームと相互作用していなかったようです（エラーなし、変更なし） Pythonツールボックスから実行したときにツールが問題なく機能し続けたとしても、実行時にマップフレームに表示されていたものと同じです）。 これは、Pythonツールボックス（TestPYT）にコピーして貼り付けるコードです。 import arcpy class Toolbox(object): def __init__(self): """Define the toolbox (the name of the …

10 arcpy  arcgis-pro  python-toolbox  tasks 

パッケージgstat、variogram()関数を使用して、Rでセミバリオグラムを作成しました。モデルの残差に空間的自己相関があるかどうかを確認したい（glmmを使用して、数kmから900 km間隔で離れたサイト間で、生息地の関数としての種の存在量）。 私の単位はkmです。したがって、私の解釈では、空間的自己相関が「問題」でなくなるまで、範囲は100 kmを少し超えていると解釈されます。ナゲットがこんなに高いように見える理由を誰かが説明できるかどうか疑問に思っていますか？これは、同じような場所でも比較的大きな差があることを意味しますか？または、この波状のバリオグラムは、より一般的な形状になるまでラグ数とラグ距離を調整する必要があることを意味しますか？ もう少し詳しく調べるために、私variog()はパッケージgeoR の関数も使用し、を使用breaks=seq(0,100,10)して、より近い距離を（同じ点と同じモデルの残差を使用して）調べようとしました。これは、最も近い点がより異なっていることを示していますが、これも意味がありません。多分これは空間的自己相関がないことを示しており、私のモデルはすでにこれを説明しています。 私はこの優れた情報源「涙のないジオスタット」を見つけました。51ページで、バリオグラムのフィッティングについていくつかの良いアドバイスがあります。このアドバイスにより、私の最初のものは正しい範囲を持っているようです。これは最初の質問に戻ります-これをどのように解釈しますか？

10 r  interpolation  spatial-statistics  autocorrelation  variogram 

私は、地域に多数の現役軍人がいることと、暴力犯罪のレベルの高低が空間的に相関しているかどうかを調べています。つまり、大規模な軍事基地を取り巻く地域は、軍事基地の近くにない地域よりも、平均して、暴力が激しい/少ないのでしょうか？ 次の2つのデータセットを使用しています。 （1）米国本土の軍事基地とそれに対応する部隊レベルの一連のポイントデータ： （2）町/都市ごとの暴力犯罪率に関する全国的な一連のデータ： 「ベース」機能が各ベースで部隊レベルを与える、ある種の重力ベースのモデルを探しているような気がします。したがって、大規模な兵士の存在は、より広い領域に影響を及ぼし、重心（つまり、GISレイヤー内のポイントの位置）近くでより強い影響を及ぼします。 概念的には、次のようになると思います。 この図では、X、Y、Zは軍事基地を表しています。a、b、c、dはそれぞれ都市を表しています（各都市の属性テーブルには暴力率フィールドがあります）。 基底の周りの勾配は影響の場を表し、基底の重心から離れるにつれて指数関数的に減少します。部隊の存在が大きくなると、影響範囲が大きくなり（しきい値の最大距離がいくらか大きくなります）、また、基地の近くのエリアに比べて、センターの近くの影響が強くなります。 各都市には、影響半径が存在する周囲のすべてのベースからのすべての「力」ベクトルの大きさの合計に基づいてスコアが割り当てられます。たとえば、私の図では、City aのスコアが0になります。基地の半径の外側。市bは、ベースXによってのみ影響を受けます。シティcはベースZの影響のみを受け、XはZよりもはるかに大きいベースであるため、そのスコアはbよりも低くなります。最後に、City dはXとYの両方のベースの半径内にあります、両方のベースからの影響の大きさの合計に基づいてスコアを受け取ります。次に、都市のより高いスコアとより高い暴力率との間に相関関係があるかどうかを確認します。 私はさまざまな重力ベースのモデル（ハフモデルなど）を調査してきましたが、QGIS / Pythonまでは見つけることができず、上で説明したものを実装する方法がよくわかりません...誰かが提案をしていますか？このため？他の分野でこの種の分析を行ったことがありますか？ したがって、TLDRは次のとおりです。 この種の質問にはどのような統計手法を使用できますか？ これを行うことができるQGISに組み込まれた（またはプラグインとして利用可能な）ツールはありますか？ QGISにこのようなものがない場合、この種の分析を実行できるPythonライブラリはありますか？

10 qgis  spatial-statistics  statistics  correlation