地理情報システム

r.neighborsラスタ用のグラスアルゴリズムを使用して処理スクリプトを記述しようとしています。このアルゴリズムにGuiを使用する場合、GRASSレイヤーの範囲フィールドを空白のままにして、最小範囲を使用できます。逆に、Pythonコンソールで使用すると、このような「デフォルト」値を入力する方法がわかりません。 私は上を見https://docs.qgis.org/2.2/en/docs/user_manual/processing/console.htmlある程度パラメータのための4つの座標入力に必須であると考えられること。使用してみましたiface.mapCanvas().extent()が、返されたオブジェクトはパラメーターとして受け入れられず、座標を取得する方法が見つかりません。 xmin, ymin, xmax, ymaxPyQGISのレイヤーを取得してアルゴリズムに入力する方法はありますか？私はそれを完全に間違っていますか？

10 qgis  pyqgis  qgis-processing  extents  qgis-python-console 

私はgdalwarpが-multi複数のスレッドを使用するオプションを持っていることを見てきました。しかし、gdal_translateにはそのオプションがありません。gdal_translateはマルチスレッドをサポートしていますか？

10 gdal-translate 

QGISを使用していて、ポイントレイヤーから別のラインレイヤーのラインに多数の垂直線を作成できるスクリプトまたはプラグインを探しています。 これまで、MMQGISのハブ距離関数（ラインをポイントに変換し、ポイントを最も近いハブに接続する）とQGIS地理アルゴリズムのそれぞれのツールを使用してみました。どちらもうまくいきませんでした。どちらも2時間以上かかり、すべてのレイヤーの線、またはポイントに垂直または接続されていない線を作成します。 写真では、プロジェクトの現在のステータスを確認できます。垂直線は、ポイントから最も近い線まで延びる必要があります。最後に、ポイントと国境を結ぶ線との交点を使用して、2ポリゴンの深さである4辺ポリゴンのバッファを作成します。これを行う簡単な方法がある場合に備えて、これについて触れます。 垂線を作成する方法についていくつかの記事があることは知っていますが、どちらも私の問題を解決しませんでした。

10 qgis  pyqgis  line  perpendicular 

Openlayers（v4.0.1）の最新バージョンがGoogleマップをタイルレイヤーとしてサポートしているかどうかを知りたいです。これに関するドキュメントイオンは見つかりません。OpenlayersがGoogleマップをサポートしていない場合、これを行う方法があるかどうか誰かに教えてもらえますか？

10 javascript  openlayers  google-maps  google-maps-api 

Pythonスクリプトを使用して仮想レイヤーを作成することは可能ですか？たとえば、「道路」というレイヤーがあり、sql "SELECT * FROM road WHERE type = 'Expressway'"を実行したいとします。 これは可能でしょうか？私が参照できる例はありますか？

10 pyqgis  virtual-layer 

Leaflet 1.0.0ツールチップのデフォルトのスタイル、特にバブル/フレームを上書きしたいと思います。このためのオプションや方法はありません。CSSにフックするにはどうすればよいですか？レイヤーを個別に変更したいので、CSSで各レイヤーのツールチップを個別に選択する必要があります。

10 leaflet  css 

外部座標を取得するために、形の良いMultiPolygonをPolygonに変換する必要があります。私はそれをこのようにします： if poly.geometry.type == 'Polygon': x, y = poly.geometry.exterior.xy elif poly.geometry.type == 'MultiPolygon': allparts = [p.buffer(0) for p in poly.geometry] poly.geometry = shapely.ops.cascaded_union(allparts) x, y = poly.geometry.exterior.xy # here happens the error これは頻繁に成功しますが、次のエラーが発生するため、Polygonが明らかにMultiPolygonのままである場合もあります。 AttributeError: 'MultiPolygon' object has no attribute 'exterior' ただし、MultiPolygonのすべての部分がポリゴンであり、それ自体がMultiPolygonではないことを確認しました。 >>>>[p.type for p in poly.geometry] ['Polygon', 'Polygon'] これがなぜ起こるのか、そしてそれを修正する方法はありますか？ ポリゴンの穴にできますか？私はこのように見えます：

10 polygon  shapely  multipart  singlepart 

QGISには、ArcGISにあるものと同様のSWIPEツールはありますか。測位の問題について、サイドスキャンソナーモザイクと水深測定を比較しています。

10 qgis  tools  swipe 

私はマップをレンダリングするためにリーフレットを使用しています。マーカー付きのマップを作成しましたが、各マーカーにイベントリスナー「onClick」を実装する方法がわかりません。 私のコード var stops = JSON.parse(json); var map = new L.Map('map', { zoom: 12, minZoom: 12, center: L.latLng(41.11714, 16.87187) }); map.addLayer(L.tileLayer('https://api.tiles.mapbox.com/v4/{id}/{z}/{x}/{y}.png?access_token={accessToken}', { attribution: 'Map data', maxZoom: 18, id: 'mapbox.streets', accessToken: '-----' })); var markersLayer = new L.LayerGroup(); map.addLayer(markersLayer); //populate map from stops for (var i in stops) { L.marker(L.latLng(stops[i].Position.Lat, stops[i].Position.Lon), …

10 javascript  leaflet 

パッケージを使用して、Rを使用したインタラクティブマップに関する投稿を読みましたleaflet。 この記事では、著者は次のようなヒートマップを作成しました。 X=cbind(lng,lat) kde2d <- bkde2D(X, bandwidth=c(bw.ucv(X[,1]),bw.ucv(X[,2]))) x=kde2d$x1 y=kde2d$x2 z=kde2d$fhat CL=contourLines(x , y , z) m = leaflet() %>% addTiles() m %>% addPolygons(CL[[5]]$x,CL[[5]]$y,fillColor = "red", stroke = FALSE) 私はbkde2D関数に精通していないので、このコードをシェープファイルに一般化できるかどうか疑問に思っていますか？ 各ノードにヒートマップで表現したい特定の重みがある場合はどうなりますか？ leafletRのマップでヒートマップを作成する他の方法はありますか？

10 r  leaflet  heat-map 

大量の画像コレクションをオンザフライでモザイクするためにvrtを使用する場合、個々のラスターまたは最終的な.vrt、あるいはその両方に概要（ピラミッド）を追加する方が良いですか？（「より良い」とは、最終的な画像の使用と閲覧にスムーズなエクスペリエンスを提供しながら、必要以上にストレージを使用しないことの間のバランスポイントを見つけることを意味します。）

10 vrt  image-pyramids  raster-catalog 

このhttp://tinyurl.com/odfbanuのような本当に単純なシェープファイルを少しラスタライズする必要があります。これは、米国のシェープファイルクーティン郡です。GDALRasterizeLayerはすべてのポリゴンをラスターに書き込みませんか？しかし、Giopandasやfionaを使用して、tiffの記述部分にrastraioを使用する方法があるかどうか疑問に思っていました。 したがって、私の目標は、ラスタライズして、共通の値（例ではLSAD）を共有するすべてのポリゴンに値を割り当てることです。 だから私はスレッドでshongololoに触発されたコードの最初を書きました： Pythonで属性に基づいてポリゴンをディゾルブします（シェイプリー、フィオナ）？。 from geopandas import GeoDataFrame name_in = 'cb_2013_us_county_20m.shp' #Open the file with geopandas counties = GeoDataFrame.from_file(name_in) #Add a column to the Geodataframe containing the new value for i in range (len(counties)): LSAD = counties.at[i,'LSAD'] if LSAD == 00 : counties['LSAD_NUM'] == 'A' elif LSAD == 03 : …

10 python-2.7  fiona  geopandas 

QGIS APIを使用してPythonスクリプトを作成しています。postgresデータベーステーブルから結果を取得する必要があります。テーブルには、ジオメトリフィールドとその他の非ジオメトリフィールドがあります。 非ジオメトリフィールドでSQLクエリを実行し、コードで結果を取得したいと考えています。これを行う方法？そうするクラスはありますか？ QgsDataSourceURIクラスについては知っていますが、私が知る限り、ジオメトリフィールドからの結果しか得られません。

10 qgis  postgresql  pyqgis 

私は通常、静的な地形図を作成します。これはほとんどが印刷されます。OSMとラスターファイル（主に異なるDEM）をよく使用しています。私がマッピングしているエリアが20.000km²を超えることはめったにありません。QGIS、SAGA、GRASSなどの排他的なオープンソースソフトウェアを使用しています。 シェープファイルの問題は知っていますが、処理は簡単です。そして、私はpostgisやSpatiaLiteにはあまり詳しくありません。私はpostgisを試しましたが、それを処理するのに多くの問題がありましたが、SpatiaLiteの方が簡単に思えます。 どう思いますか：シェープファイルを使い続けるか、SpatiaLiteやpostgisなどのデータベースに切り替えますか？私の場合の長所と短所は何ですか？

10 postgis  shapefile  spatialite 

これが繰り返し質問である場合は申し訳ありませんが、ここで検索したところ、回答が見つかりませんでした。Openpaths.ccからエクスポートしたJSONデータがあり、緯度と経度の値が含まれています。ただし、データはGeoJSON形式ではないため、OGR2OGRで読み取ることができません。私のデータをGeoJSONに変換する方法について誰かが正しい方向に向けてくれますか？最終的にはリーフレットで表示したいです。データは次のようになります。 [ { "lon": -73.97, "device": "iPhone3,3", "version": "1.1", "t": 1381167616, "lat": 40.66, "alt": 67, "os": "6.1.3" }, { "lon": -73.96, "device": "iPhone3,3", "version": "1.1", "t": 1381171200, "lat": 40.66, "alt": 45, "os": "6.1.3" } ]

10 leaflet  geojson  ogr2ogr  json