地理情報システム

GISシステムとGPUの新たな交差点に興味があります。これにより、特定のクラスのGIS問題が大幅に改善されます。この分野について議論する良いリソースを知っていますか？

11 geoprocessing  references  gpgpu 

列が次のようなスプレッドシートになるような何かを達成したいと思います。 列A-地名 列B-地名 列C-列Aと列Bの場所間の道路距離 列D-列Aと列Bの場所間の道路ごとの移動時間 Google Docsのスプレッドシートとその「Google」機能を使用してそれが可能かもしれないと思ったが、トリックを行う可能性のあるものは見つかりませんでした。 最終的には、Google Maps APIを使用してアプリをノックアップすることができますが、可能であれば回避します。

11 google-docs 

閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか？この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 2年前に閉店。 DEMのシンクを検出して埋めるためのオープンソースまたは安価なツールはありますか？ArcGIS Spatial Analystは私の価格範囲外です。

11 grass  dem  software-recommendations  open-source-gis 

4コアのdebianマシンでgdal2tilesを実行します。 sudo /usr/bin/gdal2tiles.py -r cubic -a 0,0,0 -z 10-15 /home/adm/topo/ekb.vrt /var/www/tiles ではhtopIを参照してください。 gdalが作業でコアを1つだけ使用するのはなぜですか？Gdalバージョンは2.1.2です。

11 gdal 

ArcToolboxでPythonスクリプトを実行することとスタンドアロンスクリプトとして実行することの違いを研究した人はいますか？バンド1を抽出して、RGBイメージのセットをシングルバンドに変換するための簡単なスクリプトを作成する必要がありました。スタンドアロンのスクリプトとしてPCに対して読み取りと書き込みを行うと、約350秒で1000の同じサイズのイメージを処理します。ArcToolboxから同じスクリプトを実行するには、約1250秒かかります。 import arcpy import csv from os import path arcpy.env.workspace = in_folder image_list = arcpy.ListRasters() #Create a CSV file for timing output with open(outfile, 'wb') as c: cw = csv.writer(c) cw.writerow(['tile_name', 'finish_time']) #Start the timer at 0 start_time = time.clock() for image in image_list: #Extract band 1 to create a …

11 arcpy  performance  python-script-tool  idle 

https://dwtkns.com/srtm30m/からいくつかのSRTM DEMラスタータイルをダウンロードし、タイルから陰影起伏バージョンを作成しました。ただし、詳しく調べると、陰影起伏のタイルには、非表示または削除する方法がわからない境界線があります（下の画像のように、4つのタイルの境界線が表示されています）。 [スタイル]オプションを確認しようとしましたが、線/境界線を非表示にする方法が見つかりません。タイルをシングルバンドの疑似カラーとしてレンダリングしようとすると、線が表示されないので、陰影起伏レンダリングを行ったときにのみ線が生成されるのではないかと思います。 これらの陰影起伏タイルの境界線を非表示/削除するにはどうすればよいですか？

11 qgis  raster  tiles  srtm  hillshade 

試験を採点するときに、QGIS 3.2と3.4 $areaがポリゴンの計算で異なる値を生成することを発見しました。 とを使用$areaしarea($geometry)て、平面計算と楕円体計算の違いを確認する実験を行いました。QGIS 2.18は平面領域をのデフォルトとして計算するようですが$area、QGIS 3.4は楕円領域を提供します$area（関数のヘルプセクションで説明されています）。 QGIS 3.2で、$areaのものとは異なる値（非常に高い）を生成する$areaのに対し、QGIS 3.4でのarea($geometry)2つのバージョンが同じ値を生成します。 コンテキストとして、私のデータはEPSG：32717（WGS84 UTM z17 S）に投影されます。 何が起こっているかについての洞察はありますか？これはバグですか？学生にこれらの注意事項を説明し、エリアから派生した変数を計算するときに注意したいと思います。

11 qgis  qgis-3  area 

緯度経度の値を修正する「式」を開発しようとしています。 私はVueリーフレットを使用していますが、「最初の」世界の外でパンすると、大きな数字が表示されます。+180以上または-180未満。 例：アメリカに向かって右（東方向）にパンすると、lng 215になります。私の考えでは、次のように修正するだけです。 215-360=-145 同じことは、東ロシアに向かって左（西方向）に画面移動した場合も同様で、たとえば-222になります。今計算する必要があります-222+360=138 ただし、世界は無期限であるため、ユーザーは8番目の世界に移動でき、値を調整する必要がありました。 正しい経度を計算することは可能ですか？（そして別の要件は、ユーザーが最初の世界にいるとき、24 lngはまだ24 lngであるべきです。

11 leaflet  javascript  latitude-longitude 

測地緯度は、基準楕円体の表面上の点の法線を基準に測定されることを知っています。しかし、表面の上下の点はどうですか？それらは双曲線経路をたどりますか？（私が作成したグラフィックを参照してください。）または、直線に従っていますか？ ウィキペディアは、「[楕円体]座標は、参照楕円体によって囲まれた質量の均一分布のための重力場のモデルにおける自然な選択です」と述べています。 緯度は、可能であれば重力に従うべきですよね。

11 latitude-longitude  elevation  depth  geodetic 

私は2日間問題に取り組んでいますが、トピックがPostGIS（v2.5）の交差部分である場合、多くの人が行き詰まることにも気づきました。そのため、より詳細で一般的な一般的な質問をすることにしました。 次の表があります。 DROP TABLE IF EXISTS tbl_foo; CREATE TABLE tbl_foo ( id bigint NOT NULL, geom public.geometry(MultiPolygon, 4326), att_category character varying(15), att_value integer ); INSERT INTO tbl_foo (id, geom, att_category, att_value) VALUES (1, ST_SetSRID('MULTIPOLYGON (((0 6, 0 12, 8 9, 0 6)))'::geometry,4326) , 'cat1', 2 ); INSERT INTO tbl_foo (id, …

11 postgis  sql 

私がやろうとしているのは、一連のポイントフィーチャにラベルを付けることですが、ポイントフィーチャの密度が高いため、ラベルをポイントフィーチャから遠ざける必要があります。次に、（矢印ではなく）単純な直線を使用して、どのラベルがどのポイントを参照しているかを示します。 プラグイン「Easy Custom Labels」を使用してラベルを生成し、希望する場所に移動しました。しかし、私ができないことは、ポイントとテーブルの間に単純な細い直線を作成することです。 私が達成した最高の方法は、「シンボルレイヤータイプ」（「ラベルレイヤー」>「プロパティ」>「スタイル」）を「矢印」に変更し、「曲線矢印」のチェックを外し、「各セグメントの矢印を繰り返す」のチェックを外すことです。そして、頭の長さと太さを0にし、これをできるだけ「矢印の幅」と「開始時の矢印の幅」にします。しかし、これはまだ分厚いラインを作成します。画像Aを参照してください。 「シンボルの後でタイプ」を「単純な線」に設定してみましたが、これにより、ねじれ/頂点/コーナーが線で表示されます。画像Bを参照してください。 理想的には、画像Bと同じ幅のラベルと対象物（画像Aのような）の間の直線が欲しいです。

11 qgis  qgis-plugins  labeling 

皆さん、こんにちは。PyQGISは初めての方で、QGIS 3.0でスクリプトを実行しようとしましたが、成功せず（完全に2.18で実行されました）、常に同じエラーが発生します。 AttributeError：モジュール 'processing'に属性 'runalg'がありません上記の例外の処理中に、別の例外が発生しました：トレースバック（最後の最新の呼び出し）：ファイル "C：\ PROGRA〜1 \ QGIS3〜1.0 \ apps \ Python36 \ lib \ code.py "、91行目、runco​​de exec（code、self.locals）ファイル" "、8行目、TypeError：BaseExceptionから継承しないクラスのキャッチは許可されていません これはPYTHONPATHと関係がありますか？次のように入力すると、アルゴリズムのリストも表示されません：processing.alglist（）

11 pyqgis  qgis-processing  qgis-3 

先ほど、2017年にリリースされる予定のNASADEMという新しいグローバルDEMについて聞いたことがあります。これは、これまでで最高のグローバルDEMになることを約束しました。NASADEMとは正確には何で、いつ利用可能になりますか？

11 dem  srtm  nasa  nasadem 

EPSG：4326にgeopandasデータフレームがあり、特定の境界ボックス内にあるすべての行で構成される新しいデータフレームを作成します。 最初に、気になるバウンディングボックスを取得します（これは実際には別のデータフレームのバウンディングボックスです）。 print df_sussex.total_bounds [ -1.57239292 50.57467674 0.14528384 51.27465152] 次に、その境界ボックスのみで構成されるデータフレームを作成します。 pts = gpd.GeoDataFrame(df_sussex.total_bounds) そして最後に、その境界ボックスと交差するすべての機能を取得しようとします。 sac_sussex = gpd.overlay(pts, df_sac, how='intersection') しかし、これは私に与えますAttributeError: No geometry data set yet (expected in column 'geometry'。 何が悪いのですか？

11 python  geopandas 

既存のレイヤーから一時レイヤーを作成し、選択されている既存のレイヤーからそれらの機能のみを取得したいと考えています。ArcGISには、選択したフィーチャからレイヤーを作成できるこの機能があることを知っています。 しかし、QGISでこれに対する有効な解決策を見つけることができませんでした。私がここでの目的ではない選択を新しいファイルに保存する方法を説明するこの質問を知っています。 「名前を付けて保存」のコンテキストメニューを確認しても、結果は得られませんでした。プラグインセクションで「selection」という検索語を使用しても、意味のある結果は得られませんでした。 いろいろ試してみたところ、pyqgisを使用した解決策が思い付きましたが、これを同僚に示す必要があるため、コーディングを含まない解決策を選びます。 * pyqgisの方法： layer = iface.activeLayer() # layer which has a selection feat = layer.selectedFeatures() fields = layer.fields() selection = QgsVectorLayer('Point', 'temp', 'memory') dp = selection.dataProvider() dp.addAttributes(fields) dp.addFeatures(feat) selection.commitChanges() selection.updateExtents() QgsMapLayerRegistry.instance().addMapLayer(selection)

11 qgis  qgis-2