地理情報システム

PostgreSQLデータベースに通り（道路の中心線）とポリゴンテーブルを作成しています。サンプルシナリオは次のとおりです。 問題： 通りの周りの50メートルのバッファーと両側の最も近い建物のポリゴンの交点で、通りに沿った平行線を計算する必要があります。望ましい出力シナリオは次のとおりです。 私が試したこと： 私のアプローチは： 1) Generate 50m buffer around street layer 2) get the intersection of buffer and polygons 3) Compute the distance 4) Draw offset curves (parallel lines) at both sides of street layer 5) Merge both curves to get parallel lines at the intersection これが私の試みです： WITH street_buffer …

10 postgis  postgresql 

座標について 31.96212, -103.004715 UTMコンバーターは、UTM座標がであることを示し13/R/FRます。 コンバーターの例はこちら：http : //www.rcn.montana.edu/resources/converter.aspx しかし、それらの多くはあり、それらはこれらの座標に対して同様の答えを出します。 同時に、ここのSentinel-2データセットhttp://sentinel-s2-l1c.s3-website.eu-central-1.amazonaws.com/#tiles/13/R/ FRサブディレクトリが見つかりません。 グーグルでこの場所はここにあります： そして、私が見るSentinel画像ブラウザで同じ場所を見つけると、そのタイルは異なります これを意味し13/S/FR、すなわち同じUTMと四角が、異なるバンド。 これはどのようにして可能ですか？ 更新 Sentinel-2タイルを使用したKMLはS、指定された場所のタイルも報告します アップデート2 この写真によると ここから取得すると、FR正方形はSUTMゾーンの半分とRゾーンの半分に配置されます。明らかに、ほとんどの自動コンバーターはこの四角をRゾーンに割り当てますが、Sentinel-2はそれをSゾーンに割り当てます。 ここに真実はありますか？ アップデート3 ここから取得した単純なPythonコード/gis//a/224994/32207 bandVals = "CDEFGHJKLMNPQRSTUVWXX" lon = 31.96212 lat = -103.004715 zone = int(lat + 186.0) / 6 if (lon >= 84.0): band = 'Y' if (lat < 0.0) else …

10 coordinate-system  utm  sentinel-2  mgrs 

シェープファイルをディゾルブしているとき： ogr2ogr "b:\out.shp" "b:\Tele.shp" -dialect sqlite -sql "SELECT ST_Union(geometry), nev, jar FROM Tele GROUP BY jar" 次のようなエラーメッセージが表示されました。 警告1：1つまたは複数の文字をUTF-8からISO-8859-1に正しく変換できませんでした。 Tele.cpgに "UTF-8"が含まれているため、理解できません これを試すとき： ogr2ogr "b:\out.shp" "b:\Tele.shp" -dialect sqlite -sql "SELECT ST_Union(geometry), nev, jar FROM Tele GROUP BY jar" ENCODING=UTF-8 エラーメッセージは次のとおりです。 警告1：-sqlと組み合わせてレイヤー名は無視されます。エラー1：フィールド名 'nev'の作成に失敗しました：「UTF-8」に変換できません それで、私の奇妙なハンガリーの性格を保つために私は何をすべきですか？

10 ogr2ogr  utf-8  encode 

（ウェブ）マップを作成します。これは、いくつかの空港周辺の飛行機の範囲を示しています。 飛行機の射程を考慮してバッファを計算してみました。ここで結果を見ることができます。 しかし、飛行機はまっすぐなルートをとらず、カーブが短いために飛行するため、結果が間違っていることに気づきました。 短い曲線で範囲を計算する方法はありますか？

10 qgis  buffer 

ポイントが重複している場合、「ポイントの変位」と呼ばれる、ポイントの周囲にそれらの多くを自動的に個別に表示できるようにするこのプロパティがあります。しかし、それはラインでは機能しません、それでも私はそのようなことを達成するためにかなり概念的に実現可能であるように見えます： 実際にはすべて同じ場所にあるさまざまな回線を確認する必要があります（私は通信ネットワークで作業しています）。現時点で私が見る唯一の方法は、上の図のように実際に異なる線を作成することであり、それにより空間的な間違いを引き起こします。 QGIS 2.14を使用しています。

10 qgis  symbology  style  visualisation  qgis-custom-function 

QGISではright-click、ツールバー領域のどこかを使用して、下のスクリーンショットのようにパネルをオン/オフに切り替えることができます（たとえば、GPS情報パネル）。 これをプログラムで行う必要があるので、PyQGISからパネルのオン/オフを切り替える方法は？QgsInterfaceクラスにaddDockWidgetand removeDockWidgetメソッドがあることがわかりましたが、これがこれを達成する正しい方法であるかどうか本当にわかりませんでした...

10 pyqgis  panel 

私はこの状況でここで少し問題を抱えています。 「NumeroCamp」に「Text」に数字が表示されるようにしたいだけです。 それを助けることができるフィールド電卓の機能はありますか？

10 qgis  field-calculator  string 

マルチポイント入力機能から溶解バッファーを作成する必要があります。以下の例では、入力テーブルに4つのフィーチャが含まれています。フィーチャ#2は2つのポイントジオメトリで構成されています。バッファを作成した後、4つのポリゴンジオメトリを取得します。 結果をグループ化する方法はありますか？ポイント#1とのバッファ#2はディゾルブされ、単一のマルチポリゴンフィーチャである必要があります（a）。 これまでに行ったこと： -- collect all buffers to a single multi-polygon feature -- dissolve overlapping polygon geometries CREATE TABLE public.pg_multibuffer AS SELECT row_number() over() AS gid, sub_qry.* FROM (SELECT ST_Union(ST_Buffer(geom, 1000, 8))::geometry(MultiPolygon, /*SRID*/) AS geom FROM public.multipoints) AS sub_qry; 編集： -- create sample geometries CREATE TABLE public.multipoints ( gid serial …

10 postgis  buffer  multipart 

グラフィカルモデラーとメインのQGISアプリケーションを頻繁に切り替え、モデラーウィンドウが開いているときにQGISを操作できないのはなぜですか。モデルを保存してモデラーを閉じたり、再度開いたりすると困るので... モデラーがモーダルになるより深い理由はありますか？ この動作をオフにする簡単な方法はありますか（Pythonコンソール）？

10 qgis  qgis-modeler  modal 

SQLiteデータベースにMultiLineStringZMレイヤーがあり、頂点のメジャーまたはM値を視覚化しようとしています。QGISでこれを行う方法に関する情報を調べてみましたが、収集できたすべてのことは、これをラインストリングレイヤーから直接行うことはできず、ポイントを別のレイヤーに抽出する必要があることです。 マルチラインストリングレイヤーの頂点を表すマルチポイントレイヤーを作成するためにVector-> Geometry Tools-> Extract nodesを使用しましたが、プロセスは頂点のm値を失います。m値をポイントの属性として保存するか、何か他の方法でm値を保持する必要がありますか？ 内部的には、ラインストリングをポイントシェープファイルに変換するコマンドラインツールがあり、m値が各ポイントの属性として保存されています。これを使用して、頂点にm値が割り当てられていることを確認しました。それは私がしなければならない場合ですが、可能であれば、これをQGISの内部で直接実行できればよいでしょう。 編集 -私が上で言ったことを繰り返しますが、GDALライブラリを使用する私が探している結果を達成できるコマンドラインツールがあるという事実を強調します。 PyQGISは私が探している答えではありません。組み込みツール、QGIS用にすぐに作成できるプラグイン、またはMultiLineStringZMまたはLineStringZMジオメトリからm値を抽出（作成/生成ではなく）して視覚化できる完全なスクリプトを探しています。

10 qgis  linestring  linear-referencing  extract  m-values 

印刷するマップを作成しようとしていますが、ポリゴンのサイズが大きく異なります。これは参照マップになるので、すべてのポリゴンにラベルを付ける必要がありますが、一部の名前はポリゴンに収まりません。そのような場合は、代わりに機能IDを使用し、IDを名前にリンクする属性テーブルを用意します。 私がやりたいのは、マップで使用するラベルと属性テーブルに含める行を選択するために使用できるある種の式です。このようなisLabelledものは、マップの縮尺、ラベルの長さ、フォントサイズなどを考慮に入れることができます。

10 qgis  labeling  print-composer  qgis-custom-function 

地質図を一般化しています。最初のステップでは、削除、拡大、集約などの小さなポリゴンを扱います。 上記の操作を適用する前に、細長いポリゴンやより多くの丸いポリゴンなどの形状分析から始めたいと思います。これは対処する重要な要素になる可能性があるためです。 地質図は非常に複雑な形状をしており、可能な限り考慮および維持する必要があります。 ジオメトリに基づいてポリゴンの形状を識別する方法はありますか？ 基本的にはArcGIS、Pythonを使用しています。

10 arcgis-desktop  arcpy  shapefile  geology 

初心者はここでGISと格闘しています。私は、郡のウェブサイトの郡のウェブサイトにあるシェープファイルを使用して、ミルウォーキー市の病棟を作成しようとしています。私はいくつかの成功を収めてここのスレッドをフォローしています。私のコードは与えます： from pyproj import Proj, transform # wisconsing EPSG:32054 # epsg:4326 is for the entire world, wgs 84...not obvious inProj = Proj(init='epsg:32054') outProj = Proj(init='epsg:4326') x1,y1 = 2560131.496875003, 406816.434375003 x2,y2 = transform(inProj,outProj,x1,y1) print(x2,y2) 出力付き、 -65.70220967836329 43.08590211722421 問題はこれが間違っていることです。ミルウォーキーの経度/緯度は-87.863984および42.920816です。 次に、シェープファイル全体に対してプログラムでこれを行うにはどうすればよいですか。これをベースマップにプロットしたいと思います。このスレッドをフォローしようとすると、エラーコードが表示されます： with fiona.open("ward2012/ward.shp") as shp: ori = Proj(init='epsg:32054' ), dest= Proj(init='EPSG:4326',preserve_units=True) with fiona.open('ward2012/MKE_wards_lat_lon.shp', …

10 python  shapefile  fiona  pyproj 

QGISラスタースタイルダイアログを使用して、変更されたラスター値を表示します。具体的には、ラスターの値の平方根を使用したグラデーションを見たいのですが。 ラスター計算機を使用してラスター値を変換することは可能であることを理解していますが、可能であれば、その手順をスキップしてスタイリング中に数式を適用し、元のラスターの変更されたコピーを作成しないようにします。 このタイプの操作は、色/幅などのデータ定義値を使用して、ベクターファイルで可能です。ラスタでこれを行うことは可能ですか？

10 qgis  raster  style 

PostGISにポリゴンレイヤーがあります。QGISを使用してフィーチャを追加しようとすると、編集を保存するとエラーが発生します。 これは、「結果バッファがありません」というエラー通知です。 私はここで同様のケースを見つけました：機能の追加中にPostGISエラー：結果バッファーがありません。しかし、それは私にそれを解決する答えを与えません。何か案は ？

10 qgis  postgis