地理情報システム

グループ内のレイヤーパネルには多数のレイヤーがありますが、それらを表示したくありません（かなりたくさんあるので、乱雑に見えます）。これらの層は、他のいくつかの層とも結合されています。 これらのレイヤーをユーザーから見えないようにすることは可能ですか（つまり、レイヤーパネルには表示されません。レイヤーの可視性を設定することではありません）または、おそらくグループを非表示にしますか？ 私が見つけた最も近いものはこの投稿でしたが、より直接的なアプローチがあることを期待しています： QGIS凡例にいくつかのレイヤーのみを表示しますか？

13 qgis  pyqgis  table-of-contents 

QGISを使用しているユーザーがアクセスできるようにPostGISデータベースを設定するためのマニュアルに従うことを試みましたが、現時点では十分に理解できずに進行しますか？ 私の状況; 私は、緑のオブジェクト（樹木など）のデジタル化のためのPostGISデータベースを作成し、クエリを実行する研修生です。たとえば、ストリートリストと画像ファイルを属性テーブルにリンクします。 「シェープファイルスケルトン」（現在空のオブジェクトなし、属性フレームのみ）があり、現在および将来のジオリファレンスオブジェクト（タイプごとに異なるシェープファイル）が貼り付け/作成されます。これらのスケルトンシェープファイルをインポートできると信じていますデータベース？ 状況; 私たちは、ファイルごとにバックアップサーバーを備えた独自のLANとインターネット接続を備えた別の建物で働いています。したがって、ローカルファイルアクセスとインターネット接続を介して共有するファイルがあります。 データベースのセットアップは、PostGISおよびPostgreSQLの最新の.exe実行ファイル（postgresql-9.6.2-3-windows-x64およびpostgis-bundle-pg96x64-setup-2.3.2を使用して、リモートアクセス経由で管理者によって行われました。 -1）64ビットマシンの場合。管理者からログイン情報が記載されたテキストファイルが提供されました。ここでのトレーニングマニュアルは、Ubuntuを扱っているため、かなり混乱しています。 書き留めたログイン情報をどこで使用するかわからないので、行き詰まっていますか？pgAdminを使用する必要があると思いますか？

13 qgis  postgis  postgresql  installation  windows 

layer-propertiesを介して手動で各レイヤーを有効にする代わりに、新しいレイヤーを追加した直後に、レイヤーウィンドウにデフォルトで透明度スライダーを表示することは可能ですか？

13 qgis  transparency  user-interface 

線形ネットワークの大規模なセットがあり、ネットワークに沿って互いに最も離れている各ネットワークの両端を見つけたいと思います（下の画像では、DからKです）。この問題に対するブルートフォースソリューションは、発信元の各ペアのネットワークに沿った最短パスを計算することですが、数千のエンドを持つ数百のネットワークがあるため、可能な各パスの計算は非常に重いです。 総当たりを使用せずにこれを計算する最適な方法はありますか？いくつかの巧妙なルールに基づいていくつかのポイントを除外できますか？ 編集：私の質問を明確にするために、@ Alex Tereshenkovが言及した最長パスの図を追加しました。黒いパスは、最長パスアルゴリズム（頂点を繰り返さない最長パス）の結果です。私の場合、あなたはいずれかの文字からネットワークに入り、できるだけ早く別の文字に運転する必要があると想像してください。参加するのに最も難しいのはどの2文字ですか？

13 algorithm  network 

私は、小さなセグメントと、それらが最も関連している可能性が高い大きなセグメントとを一致させようと試みています。 これが私が持っているデータの典型的な例です： ここでは、セグメント652を198969に一致させる必要がありますが、711および707は何にも一致しません。 さまざまな方法、特にハウスドルフ距離（ここの回答に基づいて）を探しました。PostGISを使用して計算しましたが、奇妙な結果が得られます：最短距離は707と198985の間であり、652はたとえば198985よりも198969までの距離が大きいです（必要に応じてクエリと結果を追加できます）。 Hausdorffは実際にこれを解決する正しい方法ですか？他のアプローチはありますか？言及したパラメーター（距離、方位など）の一連のチェックを作成することを考えましたが、エッジケースまたはそれらがどれだけしきい値に達しているかなどを処理するために条件の束をすべて追加する必要があります向かい合っています。 更新：許容できる妥協案と思われる方法を見つけました： 最初に<->、10メートル以内にある（PostGIS オペレーターを使用して）一致させようとしている青いセグメントから最も近い10個の黒いセグメントを見つけます。 次に、それぞれの黒いセグメントの青いセグメントの端に最も近い点を見つけて（を使用してST_ClosestPoint）新しいセグメントを作成し、長さが青いセグメントの90％未満の結果を除外します（セグメントはそうではありません）直面している、またはベアリングの差が〜20°以上であること） 次に、距離とハウスドルフ距離がある場合、ソートされた最初の結果を取得します。 微調整を行う必要があるかもしれませんが、今のところ許容できる仕事をしているようです。いくつかのエッジケースを見逃した場合に実行する他の方法や追加のチェックを探しています。

13 algorithm  gis-principle 

閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか？この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 7か月前に閉鎖されました。 サイズが1GB（.tif）で、幅と高さが94000x71680のイメージがあります。この画像を20000X20000タイルにチャンクして、処理できるようにします。 これどうやってするの？

13 geotiff-tiff  tiles  digital-image-processing 

8つの列と1670万件のレコードを持つテーブルがあります。列に対してif-else方程式のセットを実行する必要があります。UpdateCursorモジュールを使用してスクリプトを記述しましたが、数百万件のレコードを記録した後、メモリが不足します。これらの1670万レコードを処理するより良い方法があるかどうか疑問に思っていました。 import arcpy arcpy.TableToTable_conversion("combine_2013", "D:/mosaic.gdb", "combo_table") c_table = "D:/mosaic.gdb/combo_table" fields = ['dev_agg', 'herb_agg','forest_agg','wat_agg', 'cate_2'] start_time = time.time() print "Script Started" with arcpy.da.UpdateCursor(c_table, fields) as cursor: for row in cursor: # row's 0,1,2,3,4 = dev, herb, forest, water, category #classficiation water = 1; herb = 2; dev = 3; forest = …

13 arcpy  cursor  performance  arcgis-10.4  loop 

市街地のシェープファイルがありますが、残念ながらシェープファイルには1行しか含まれていません。シェープファイルは次のようになります。 これらの線をQGISの面積ポリゴンシェープファイルに変換したいと思います。「Lines to polygons」ツール（Vector> Geometry Tools> Lines to polygons）で試しましたが、正しく機能しませんでした。このツールを使用すると、シェープファイルは次のようになります。 QGISでラインシェープファイルをエリアポリゴンシェープファイルに変換する他の方法はありますか？ ところで、私はQGIS 2.14.3 Essenを使用しています。

13 qgis  polygon  line 

TileStacheを使用してベクタータイルを提供しています。すべての設定を思いどおりに行います。データはPostgresに保存されており、VecTilesプロバイダーを使用してGeoJSONタイルを提供しています。 すべてのタイルをキャッシュして、タイルの配信を高速化したい。私はtilestache-seed.pyを使用してキャッシュをシードしています。複数のマシンでtilestache-seedを実行しています。Tilestache-seedはズームレベル13まで本当にうまく機能しましたが、その後、タイルをキャッシュするには時間がかかりすぎています。ズームレベル16の場合、キャッシュする5023772タイルがあり、各マシンで1日あたり100k〜200kのタイルしか取得できません。 タイルキャッシュを高速化するにはどうすればよいですか？tilestache-seed.pyを微調整してシードを高速化する方法はありますか？ 更新：テーブルに空間インデックスを構築しようとしました（ジオメトリカラムとwhere句でデータをフィルタリングするために使用されるカラム）。それでも、タイルスピードの大幅な向上は見られませんでした。このレートでは、ズーム17のみで1か月かかりますが、この時間はズーム21に向かって指数関数的に増加するだけです。 更新2：マテリアライズドビューも作成しようとしましたが、パフォーマンスに目に見える変化はないため、データベースの最適化は機能していません。tilestache-seed.py自体を最適化するか、タイルをキャッシュする新しい方法を考案する必要があると思います。 ハードウェア情報 8つの異なるPCでキャッシングプロセスを実行しています。1つは32GB RAMを搭載したi7で、もう1つは4GB RAMを搭載したi3ですが、どちらもほぼ同じキャッシュ速度を提供します

13 python  postgis  cache  vector-tiles  tilestache 

ラスターをポリゴンに変換するオープンソースのPythonソリューションを探しています（ArcPyは使用しません）。 ラスターをポリゴンに変換するGDAL関数を知っていました。これはマニュアルです：http : //pcjericks.github.io/py-gdalogr-cookbook/raster_layers.html#polygonize-a-raster-band それにもかかわらず、出力はファイルとして保存されず、一時的にメモリ内にある、形の良いポリゴンまたは任意のオブジェクトになる可能性があります。この問題を処理するパッケージまたはコードはありますか？ ラスターがnumpy配列で処理されている場合、アプローチは以下のとおりです。

13 python-2.7  shapely  polygonize 

ArcGIS にCache / Mixed形式のラスターデータセットがあります。QGISのような他のデスクトップベースのGISソフトウェアでベースマップとして使用するために、これをジオリファレンスされたtiffまたはその他のラスターイメージ形式にエクスポートする必要があります。 これまでのところ、タイル形式を.tpkファイルまたは展開/コンパクトキャッシュ形式のいずれかにのみ変更できる、エクスポートタイルキャッシュ（データ管理）という名前のツールのみがArcGISで見つかりました。これらのタイルデータをイメージに変換するツールが見つかりませんでした。 ArcGISにあるデータのエクスポートオプションを使用すると、結果の画像は単なる黒の画像になります。 これらのタイルデータを画像にエクスポートする方法を知っている人はいますか？ 編集済み @felixIPの回答は解決策になる可能性がありますが、別の回避策を探しています。下の画像のように、タイルにはいくつかの設定ファイルがあります conf.cdiルックスは下記好き <?xml version="1.0" encoding="utf-8" ?> <EnvelopeN xsi:type='typens:EnvelopeN' xmlns:xsi='http://www.w3.org/2001/XMLSchema-instance' xmlns:xs='http://www.w3.org/2001/XMLSchema' xmlns:typens='http://www.esri.com/schemas/ArcGIS/10.1'> <XMin>8142366.0491449088</XMin> <YMin>4370513.4222595459</YMin> <XMax>8146042.4910550155</XMax> <YMax>4375009.1735663339</YMax> <SpatialReference xsi:type='typens:ProjectedCoordinateSystem'> <WKT>PROJCS["WGS_1984_Web_Mercator_Auxiliary_Sphere",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator_Auxiliary_Sphere"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",0.0],PARAMETER["Standard_Parallel_1",0.0],PARAMETER["Auxiliary_Sphere_Type",0.0],UNIT["Meter",1.0],AUTHORITY["EPSG",3857]]</WKT> <XOrigin>-20037700</XOrigin> <YOrigin>-30241100</YOrigin> <XYScale>148923141.92838538</XYScale> <ZOrigin>-100000</ZOrigin> <ZScale>10000</ZScale> <MOrigin>-100000</MOrigin> <MScale>10000</MScale> <XYTolerance>0.001</XYTolerance> <ZTolerance>0.001</ZTolerance> <MTolerance>0.001</MTolerance> <HighPrecision>true</HighPrecision> <WKID>102100</WKID> <LatestWKID>3857</LatestWKID> </SpatialReference> </EnvelopeN> 一方でconfig.xml、以下の情報を持っています <?xml version="1.0" encoding="utf-8" ?> <CacheInfo xsi:type='typens:CacheInfo' xmlns:xsi='http://www.w3.org/2001/XMLSchema-instance' xmlns:xs='http://www.w3.org/2001/XMLSchema' …

13 arcgis-desktop  arcgis-server  export  tiles  cache 

NASAが提供する降水量に関するTRMM 3B31ファイルを使用しています。他のシェイプファイルと一緒にWGS84投影を選択してQGISにロードすると、マップが垂直に表示されます。つまり、QGISは緯度と解釈され、代わりに経度と緯度が解釈されます。私はgdal_translateオプション-a_ullrを選択して良い座標を選択しようとしましたが、問題は解決しません：画像は私の要件を満たすために圧縮および伸張されるだけですが、これは明らかに私が望んでいるものではありません。 .csvファイルを読み込むときに、緯度と経度を選択できますが、ラスタに対しても同じことを行う方法はありますか？ 193個のラスターがあるので、繰り返し実行できるコマンドを提案できれば素晴らしいと思います。 この画像は、インドの地図と3B31ファイルの最初のバンド、つまり地表降水量を示しています。2本の黒い線は40および-40平行です。それらを水平にしたいので、雨のデータはインドと重複しています。

13 qgis  raster  latitude-longitude  trmm 

分類する必要がある植生の航空写真のセットがあります。手動の解釈を使用する代わりに、QGIS ProcessingのOrfeo Toolboxに基づくオブジェクトベースの分類を使用してみます。 ただし、QGISでそれを行う方法を説明するチュートリアルは見つかりません。Orfeoガイドでは、手順は次のとおりです。 画像のセグメンテーション（その全体または一部のみ）; イメージからLabelObjectMap（一種のstd :: map）変換。 最終的な再ラベル付け。 セグメンテーション前の画像を使用した領域の属性計算： オブジェクトのフィルタリング LabelObjectMapから画像への変換。 しかし、どのジオアルゴリズムを使用するかについて、特定のステップバイステップを見つけることができません。私はLSMC（大規模平均シフト分類）ワークフローに従い、セグメンテーション（平均シフト）も実行しようとしました。LSMCの結果は非常に良好でしたが、オブジェクトにラベルを付けて（またはトレーニングオブジェクトを作成して）分類するために次に何をすべきかわかりません。 TrainImageClassifier（SVM）、TrainORGLayerClassifier、ORGLayerClassifierを試しましたが、出力はありませんでした。チュートリアルを見つけようとしましたが、見つかりませんでした。ワークフローについて完全に混乱しています。

13 qgis  classification  orfeo-toolbox 

代理店からプロジェクトのシェープファイルを受け取り、常にUTMを使用していると言われました。ただし、値の桁数が多すぎてUTMに投影できません。 min max x 32507319 32624343 y 5925974 6074881 ファイルをQGISに読み取ろうとしましたが、どの投影を試しても、世界地図を使用していても、形状は常にマップの境界の外側に表示されます。それらが表示される領域は、ドイツの最北端の州であり、UTMで次の最小値と最大値を持つシュレスヴィヒホルシュタインです。 min max x 426167.5 650147.9 y 5913377.7 6101334.8 簡単にわかるように、上位データに何か問題があります。これを修正するために使用できる投影法は誰にもわかりますか？

13 coordinate-system  germany 

Postgisの空間GiSTインデックスの内部データ構造とストレージメカニズムについて興味があります。Oracleでは、SDOインデックスが「単なる別のテーブル」であり、レベルがインデックスの属性であり、インデックス付きジオメトリのbboxがblob（抽出可能）として格納されている方法を示すのが好きでした。postgisはインデックスをどのように保存しますか？ システムテーブルhttp://www.postgresql.org/docs/8.3/static/catalogs.htmlを使用してインデックスを識別できますが、実際のストレージにアクセスする方法がわかりません（実際にテーブルである場合）

13 postgis  postgresql  spatial-index