地理情報システム

各頂点/ポイントに関連付けられたいくつかのデータを持つ一連のラインストリングのジオパンダデータフレームがあります。 Point_x = (Lat, Lon, Time, ID, Data1, Data2, Data3) ポイントは、IDに基づいてラインストリングに変換され、時間順に並べられます。 ある条件が満たされた時点で、ラインストリングを分割したい。現在、ポイント間の距離が特定の値よりも大きい場合です。将来的には、データフィールドの機能が何らかの値になる可能性があります。たとえば、速度が5 kphを超えるとラインストリングを分割します。 現在の問題は、トラックの一部が重複したIDを持つポイントから形成されているため、ラインストリングが非常に長い距離で前後にジャンプし、これらのラインを分割するためのしきい値が必要なことです。 これまたはライブラリ/メソッドを構造化する正しい方法に関するアイデアはありますか？ データフレームには15,000を超えるトラックがあり、トラックごとに多くのポイントが含まれているため、効率が良いでしょう。 トラックDFの例を次に示します。 ID geometry 204235000 LINESTRING (37.62001 -28.99535, 37.62015 -28.9... 205400000 LINESTRING (3.807816666666666 -18.083181666666... 207138000 LINESTRING (22.73206 -34.97915833333333, 22.73... 209016000 LINESTRING (8.447673333333331 -23.522783333333... これがポイントDFのサンプルです。Datetime、Point（Lon、Lat）、Speed、Sizeなどを含む18列があります。 Index Heading Latitude Longitude ID 20 92.8 -35.946802 13.089695 210725000 21 …

11 python  shapely  geopandas 

GDAL（OSGeo4W）を使用して、16ビットの符号なしTIFFイメージを8ビットイメージにスケーリングしようとしています。ただし、8ビット画像の範囲内に収まるようにピクセル値を制限します。つまり、16ビット画像の視覚的な割合が8ビット画像（0〜255）で表現されるようにします。gdal_translateを使用して16ビットから8ビットに変換すると、ピクセル値がカットされ、RGBイメージは生成されません。gdal_translate内からどの関数と設定を呼び出すか、またはより良いオプションがあるかどうかはわかりません。 GDALINFO C:\>gdalinfo C:\Projects\Stormwater\ForPCI\images\1537TORO_6cm_04761_RGBI.tif Driver: GTiff/GeoTIFF Files: C:\Projects\Stormwater\ForPCI\images\1537TORO_6cm_04761_RGBI.tif Size is 17310, 11310 Coordinate System is `' Metadata: TIFFTAG_DATETIME=2015:05:03 12:27:13 TIFFTAG_IMAGEDESCRIPTION=UltraCam-Lvl03 -------------- CAM_ID: UC-SXp-1-50215465 [3] IMG_NO: 4761 RECORD_GUID: 8fc3d7a9-d2e5-40e5-babf-939f1f803dcf IMG_GUID: C850D0A4-39CA-4BC8-ACD7-81A363D78C6E FILE_GUID: 86C4F8E4-A2FD-43FB-B3B1-4C4ED8DF0E42 LICENSE_ID: 820342819 SOFTWARE: UltraCam Aerial Radiometry Core 12.6.1408.2501 PIXEL_SIZE_WIDTH: 6 [micron] PIXEL_SIZE_HEIGHT: 6 [micron] APERTURE: F_8 EXPOSURE_TIME: 0.002000000000 …

11 gdal  imagery  gdal-translate 

大きなGTFSフィードから小さな領域を抽出することをサポートするツールはありますか？GoogleのTransitfeedライブラリには、マージツールは含まれていますが、抽出/クリップツールは含まれていません。私が見つけた他のライブラリは、データの読み取りとプロットのみが可能ですが、操作はサポートしていません。 背景情報：Transitfeed Feedvalidatorには、GTFSファイルに関する多くのエラーが表示されます。関心のある地域を抽出して、地域の品質を評価したい。 Twitterの推奨事項は、GTFSリーダーとライターの両方を備えたFMEを示しています。私の簡単なワークフローは、GTFSリーダー（左下）、Shapefileリーダー（左上の関心領域のポリゴン）、クリッパー、および出力用のGTFSライター（現在のところストップとシェープファイルのみを書き込みますが、最終的にはすべてのファイルが必要です）： これまでの問題は、出力シェイプのコンテンツ（オレンジポイント）に、関心のある黄色の領域内にあるすべての入力ポイント（緑）が含まれていないことです。 （データソース：https : //transitfeeds.com/p/helsinki-regional-transport/735/latest/download） [集計の作成] = [いいえ]を使用すると、クリップされた図形は完全に見えます しかし、停止は間違いなくありません

11 fme  general-transit-feed-specification 

QGISで平行線の複数の頂点を同時にドラッグすることは可能ですか？ 平行線のラインレイヤーがあります。一部の地域では、すべての平行線の頂点を移動して、互いの距離を維持する必要があります。 ノードツールを使用すると、ドラッグしたいラインの1つの頂点のみを選択できますが、同時に複数のラインの頂点は選択できません。 Multiline-Layersでは、すべての行を同時にドラッグできますが、単一行のすべての属性を保持する必要があるため、これはオプションではありません。 行をすべて同時にドラッグする他のオプションはありますか？ EDIT1：だけが見つかりました。CadNodeTool（https://github.com/wonder-sk/CadNodeTool）を。このプラグインを使用すると、複数の行の頂点を同時に選択することが少なくともすでに可能です。おそらく、このプラグインを拡張して、平行線の頂点を移動することができます。

11 qgis  pyqgis  editing  vertices 

QGIS 2.16.2を使用すると、次の式を使用してシェープファイル属性テーブルからレコードを正常に選択できます。 left("start_time", 10) = '2015-08-01' start_timeは文字列フィールドであることに注意してください。 ただし、同じ式をシェープファイルのレイヤープロパティ>一般>機能フィルター>クエリビルダーにコピーして貼り付けると、次のエラーが表示されます。 同じ表現がQGISのどこでも機能すると期待しています。私は何を間違えていますか？

11 qgis  query  filter  expression 

Amazon S3からSentinel-2タイルをダウンロードしていますが、100％整列しているようには見えません。右の画像は、ArcGis Desktopに追加されたバンド08 Sentinel-2（08.jp2をマップに追加するだけです）、左のesriベースマップ写真です。Sentinelのオフセットは約35〜40メートルです。それらを修正するにはどうすればよいですか？ これがダウンロードされたタイルです。レベル1C、センチネル2： 解決： HDunnの答えから、私はこれをしました： 道路ベクトルレイヤーiを使用して、5つのポイントペアを追加し、x座標とy座標の差を計算しました。平均差は、Shiftツールでxvalueおよびy値として使用されます。結果：

11 arcgis-desktop  remote-sensing  sentinel-2 

最近のトレーニングコースで、QGISがアトラスジェネレーターを使用して作成されたマップブックの次/前および上/下のページ番号を自動的に計算できるかどうかを尋ねられました。グリッドの幅と高さがわかっていれば、通常のグリッドに対してかなり合理的なラベル表現を作成できました。 しかし、その後、私たちの興味のある地区を含まないページを描画したくない現実的な例について考え始めました。たとえば、私の本拠地郡の1つです。 それで、今日の午後、私はpythonスクリプトを使って、各グリッドセルに興味のある4人の隣人を見つけ出し、それらの値をグリッドに追加しました（これは、Ujaval Gandhiのチュートリアルに基づいています）。 for f in feature_dict.values(): print 'Working on %s' % f[_NAME_FIELD] geom = f.geometry() # Find all features that intersect the bounding box of the current feature. # We use spatial index to find the features intersecting the bounding box # of the current feature. This will …

11 qgis  pyqgis  print-composer  atlas 

細い線を使用して水域を表示した古い地形図の外観をどのように再現できますか？塗りつぶしオプションを見て、見つけられる最も近いものはシェイプバースト塗りつぶしですが、線ではなく色のグラデーションを使用します。 この画像は、私が作成しようとしている効果を示しています。 画像ソースウィキペディアコモンズ このスタイルの名前は「ウォーターライニング」と呼ばれ、平版印刷版では色のない部分を印刷できないために使用されたと思います。 QGISでこのスタイルを作成する方法はありますか？

11 qgis  polygon  style 

とが NE lat lngありSW lat lngます。私の目標は単純です。RStudioを使用しており、上記の2つの緯度と経度から境界ボックスを作成し、ポリゴンを生成して、正方形ポリゴンにランダムポイントを追加します。data.frameにインポートできるcsvファイルにNELat / LngとSW-Lat / Lngがあります。 始めるには少しガイダンスが必要です。R言語は非常に新しいもので、これまで慣れていたものとはまったく異なります。 これまでのところ私はこれを持っています： coords = cbind(78.46801, 19.53407) coordsmax = cbind(78.83157, 19.74557 ) sp = SpatialPoints(coordsmax) sp2 = SpatialPoints(coords) r1 = rbind(coords, coordsmax[1, ]) # join P1 = Polygon(r1) Ps1 = Polygons(list(P1), ID = "a") plot(Ps1) 他の場所でいくつかの例をみましたが、4つの角のポリゴンをプロットできません。

11 r  rgdal 

背景を透明に設定して、ArcGIS Pro 1.2.0からPNG画像ファイルをエクスポートしようとしています。ArcMap 10.3ではこれは簡単で、エクスポート中に透明色を白に設定するためのボタンがありました。 ArcGIS Proで画像を「共有」するときに、背景または透明色を設定するための同じオプションが表示されません。このボタンは他の場所に移動しましたか？または、もうオプションではありませんか？ 目標は、上記のスクリーンショットと同様の設定を見つけることです。このプロジェクトでは、スタンドアロンPDFファイルとしてベースマップがあります。コールアウトボックスをそのファイルの上にオーバーレイする必要があります。これには、コールアウトボックスが白いボックスではなく、周囲全体に透明性を持っている必要があります（ある程度クリップしたときのように）。

11 arcgis-pro  transparency  png 

Ubuntu 14.4を実行していて、Python仮想環境をインストールしました。その後、いくつかのPythonライブラリ（gdal、lxml、psycopg2など）をpipインストールしました 私のシステムにもQGIS用のモジュールがあります。したがって、pythonと入力すると、このモジュールを次のようにインポートできます。 import qgis.core しかし、virtualenvで同じことを行おうとすると、「モジュールが見つかりません」というエラーが表示されます。 virtualenvを作成したときに、次のオプションを使用しました。 virtualenv --no-site-packages env これはシステムの環境から私の環境を制限します。問題は、virtualenvでQGISモジュールにアクセスする方法ですか？

11 qgis  python 

私は、左から、左から、右から、右へのフィールドフィールドを持つ通りのshpファイルでQGISに取り組んでおり、ジオコーダーを構築しています。 これらの機能がどの方向を向いているかを知る必要があります。これは、線を作る2つ（またはそれ以上）のポイントのどちらが最初に来るかを知るのと同じだと思います。 マーカーラインスタイルの三角形マーカーを自動的に回転して、ラインの向きを示す方法はありますか？

11 qgis  style  line  markers 

閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか？この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 3年前に閉店しました。 多くの場合、地図を見て、いつ作成されたのだろうと思います。 日付のない世界地図の年齢を知るにはどうすればよいですか？ 世界の政治的境界のベクターファイルを持っているが、日付がわからない場合、どのリソースを使用して日付を決定できますか？ 日付を入力して、その日付の世界の境界の地図を取得できる地図サービスがウェブ上にありますか？ あると仮定すると、どのプロセスを使用して一致の程度を判断しますか？ 地図の日付に固有のルールベースのシステムはありますか？

11 time 

次のコードを使用すると、ラスター値0に100％の透明度を設定できます。 map=None for lyr in QgsMapLayerRegistry.instance().mapLayers().values(): if lyr.name() == "some_layer": map= lyr tr=None tr=QgsRasterTransparency() tr.initializeTransparentPixelList(0) map.renderer().setRasterTransparency(tr) map.triggerRepaint() 2つの値（0と8）を100％透明に定義するにはどうすればよいですか？ ラスター値8で最後のブロックのコマンドを繰り返すと、レイヤーの透明度から最初の値（0）が削除されます。

11 qgis  python  pyqgis  transparency 

QGIS（2.14.3-Essen）プロジェクトがあり、職場でも自宅でも使用しています。プロジェクトには2つのレイヤー以外のすべてを保持します。 残りの2つのレイヤーは非常に大きく、合計20 GBを超えるDEMと陰影起伏があり、職場と自宅では異なるパスに保存されます。 Handle Bad Layersダイアログを無視し、それらのレイヤーの壊れたパスを保持したいと思います。私はしばしば彼らが悪いことを気にしません。プロジェクトの他の部分に変更を加えたい。 不正なレイヤーの処理ダイアログを無視またはスキップして、更新しない場合にレイヤーへの参照が失われないようにする方法はありますか？

11 qgis  layers