地理情報システム

rasterioを使用してラスター内のポイントで単一のピクセル値を取得するには、https：//github.com/mapbox/rasterio/pull/275に例があります。 しかし、ラスタ内の単一のポイントで値を抽出するために使用できる、ラスタリオ（CLIではなく）内に直接APIがありますか？ -編集 with rasterio.drivers(): # Read raster bands directly to Numpy arrays. # with rasterio.open('C:\\Users\\rit\\38ERP.tif') as src: x = (src.bounds.left + src.bounds.right) / 2.0 y = (src.bounds.bottom + src.bounds.top) / 2.0 vals = src.sample((x, y)) for val in vals: print list(val)

14 python  numpy  rasterio 

DBマネージャーでテーブルを編集した後（列を追加するなど）、レイヤーの属性テーブルを更新する方法はありますか？

14 qgis  attribute-table  database 

Shapely / Geopandasを使用して2本の線を互いにスナップしようとしていますが、スナップの結果は非常に奇妙です。私は試した ： import geopandas as gpd from shapely.geometry import * from shapely.ops import snap lines1 = gpd.GeoDataFrame.from_file('lines1.shp') lines1 = lines1.to_crs({'init': 'epsg:2227'}) lines2 = gpd.GeoDataFrame.from_file('lines2.shp') lines2 = lines2.to_crs({'init': 'epsg:2227'}) res = lines1 lines2_union = lines2.geometry.unary_union res.geometry = res.geometry.apply(lambda x: snap(x, lines2_union, 14)) res.to_file('result.shp', driver="ESRI Shapefile") そして、この結果を得ました： lines1 = 赤い線 …

14 shapely  snapping  geopandas 

私はリモートセンシングの分野では新しいです。画像の前処理について読んだときに時々混乱しました。誰かが私に違いと放射線測定補正と放射線測定較正を行ういくつかの方法例を詳しく説明してもらえますか？DN値を、補正またはキャリブレーションとしてカウントされた異なる時間からの2組の画像（Landsat 7など）のTop of Atmospheric（TOA）反射率値に変換していますか？DOS（Dark Object Subtraction）とサンアングル補正は放射補正の一部ですか？

14 remote-sensing 

サーバーからWMSタイルを取得する前に、HTTPリクエスト（一部のマップ変換を処理する）を実行する必要があるアプリケーションがあります。 マップをドラッグした後にTileWMSの自動更新を一時停止する方法はありますか？ 非表示の設定またはレイヤーの削除と追加は機能しますが、HTTPリクエストの実行中にレイヤーが消えないようにします。私も試しましたsetUrl('')が、レイヤーも消えます。 また、OpenLayersのソースコードを読み取って、リクエストを実行する前にチェックするプロパティを定義しようとしましたが、リクエストが行われたコードで見つけることができませんでした。

14 wms  openlayers 

フロントエンドにリーフレット付きのdjangoアプリがあり、座標を渡すことができるようにgeojsonオブジェクトのbboxを計算する必要があります map.fitBounds() 私はいくつかのライブラリを試し、getBoundsリーフレットの特定のオブジェクトのメソッドを試しましたが、などFeatureGroup、境界が定義されていないことを訴えました。 誰も私にこれの簡単な解決策を教えてくれますか？

14 python  leaflet  geojson  django 

パンやズームなどのためにできるだけ多くの画像をGPUに保存すると大いに役立つようです。 非常に大きな画像を多数取り扱っています。多くの場合、8 GBから数百GBのサイズです。自宅には、40インチの4kモニターと12 GBのビデオメモリを備えたGTX Titanがあり、投げたものすべてに対して驚くほど優れたパフォーマンスを発揮します。作業？ うQGIS、ArcMapのとPCI Geomaticaは数百ギガバイトの範囲内でマルチギガバイトのイメージとイメージモザイクを扱うときに豊富なグラフィックスメモリの利点を取ることができますか？ 大きな画像を処理するときに、豊富なGPUメモリから最も恩恵を受けるのはどのGISアプリケーションですか？

14 qgis  arcmap  performance  hardware 

このソースから複数のDEMラスターを1つに結合する必要があります。 http //srtm.csi.cgiar.org/SELECTION/inputCoord.aspなので、Mosaic to New Rasterツールを使用します。 ツールをセットアップしました： ラスターを入力します（すべて同じソース、同じサイズ、投影なし）。 ピクセルタイプを設定します（元のラスターのように32ビット浮動小数点）。 セルサイズを元のラスターとして設定します。 元のラスターとしてのバンド数= 1。モザイクオペレーター、私はBLENDとMEANでそれを行いました（同じ結果が得られます）。 私が抱えている問題は、結果のモザイクが、個々のラスターの最大値と最小値とは異なる最大値と最小値の範囲を示すことです。たとえば、ラスター1（-5123.8、23.25）、ラスター2（-5974.6、40.09）、ラスター3 （-57770.2、38）、ラスター4（-2534.3、23.55）、および最終モザイクラスター（-5975.8、81.1）。 この解決策は正しくないと思います。少なくともそれを期待していませんでした。これが大丈夫かどうか、そしてそうでない場合、誰でもそれを解決し、正しい最大値と最小値を持つ適切なモザイクラスターを取得する方法を知っていますか？ ArcGIS 10.2.2 for Desktopを使用しています。

14 arcgis-desktop  raster  arcgis-10.2  dem  mosaic 

私は地理データの絶対的な初心者なので、質問が適切でない場合はご容赦ください。 NCDC NARRからデータをダウンロードし、rasterパッケージを使用してRにロードすることに成功しました。緯度、経度、値のリストを取得したいのですが。私はそれrasterToPoints()が私が望むことを正確に行うべきであることを理解していますが、私の緯度と経度の値は奇妙に見えます： r <- raster(myfile) data_matrix <- rasterToPoints(r) head(data_matrix) x y value [1,] -5405401 4347242 70 [2,] -5372938 4347242 88 [3,] -5340475 4347242 76 [4,] -5308012 4347242 85 [5,] -5275549 4347242 87 [6,] -5243086 4347242 88 私は、現在Lambert Conformal Conic（LCC）である投影法で何かをすべきだと思います。ラスターに関する詳細情報は次のとおりです。 > r class : RasterLayer dimensions : 277, 349, 96673 …

14 raster  r 

QGIS Desktop（2.7.0）のトップレベルメニューに新しいメニュー項目を追加する必要があります。つまり、「プロジェクト」、「編集」、「表示」、「レイヤー」などと同じレベルになります。 私が知っている限り、QgisInterface（http://qgis.org/api/classQgisInterface.html）を使用してメニュー項目を非表示にする方法（設定->カスタマイズ）があり、QGISで既に定義されているメニューに新しい項目を追加できます（プラグインコードからaddLayerMenu、addPluginToDatabaseMenuなど）。 しかし、最上層に新しいメニュー項目が必要です（これは顧客の要求です）。

14 qgis  pyqgis  qt  gui  context-menu 

特定の地理位置（緯度、経度）およびズームレベルのWMTSタイルのインデックス（x、y）を取得する方法を知りたいと思いました。 たとえば、（48.675、2.7）にPOIがあります。ズーム10に対応するオープンストリートマップタイルを取得します。 数学はできますか？Webサービスが必要ですか？精度：これをプログラムで行う必要があります。

14 convert  development  wmts 

FMEでは、複数のリーダーを単一のトランスフォーマーに接続する方法はありますか？私はFME 2015を使用しており、多数のフィーチャクラスを1つのフィーチャにクリップしたいのですが、これを行う簡単な方法が見つからないようです。 編集：入力データセットから複数のフィーチャタイプを選択し、それらをすべてフィーチャタイプの個々の出力矢印をクリックしてトランスフォーマの入力矢印にリンクすることなく単一のトランスフォーマにリンクしようとします。 編集2：これはスクリーンショットです、うまくいけばこれが役立つでしょう：

14 geoprocessing  fme  clip 

特定の日の旅行時間の等時線を表すポリゴンのテーブルがあります。各原点には、5つのアイソクロネジオメトリがあります（別々の行に格納されます）。各原点について、5つのアイソクロナス（バイナリNULLまたは1）をラスタライズし、それらを1つのラスタレイヤーに結合します。このラスターレイヤーには単純なマップ代数sum / 5が必要です。そのため、各原点は、最終的に[NULL、0.2、0.4、0.6、0.8、1]の値を持つ単一のラスターレイヤーに関連付けられます。構成レイヤーはオーバーラップします。それは確率曲面です。 私のデータはすべてPostgres 9.3（PostGISを使用）に保存されます。私の問題は、PostGISラスタの使用方法を学びたいが、非常に急な学習曲線を持っているように思え、すべての例が単一のラスタレイヤに対処できることです。例では、このレイヤーはポリゴンオーバーレイの一部として使用され、各ポリゴンのラスターの値を平均化しています。以下を組み合わせた複製可能な例は見つかりませんでした。a）ベクトル->ラスターb）マップ代数。およびc）最初の段落ごとのGROUP BY属性。 このタスクを実行するために必要な場合、GDALまたはGRASSを使用しても問題ありませんが、これはPostGISが処理できるもののようです。入力データがすでにPostGISジオメトリである場合、そうするのが便利です。そして、私は本当に PostGISラスターに同意したいと思っています。 サンプルデータ構造： areaid time date isogeom (polygon) 1000 07:15:00 2014-05-05 xxx 1000 07:15:00 2014-05-06 xxy ... 1006 07:15:00 2014-05-05 zzz ラスタライズし、areaidでグループ化し、マップ代数演算を実行して次のことを行います。 areaid isorast (raster) 1000 aaa 1006 bbb 私はこれをPostGISに含めることに成功していません。私のアプローチは、ベクターをラスターに変換し、ラスターを配列にダンプし、psycopg2を介してnumpy配列との組み合わせを実行してから、GeoTIFFに書き込む（PostGISに戻す）ことです。理想的ではありませんが、実行可能です。

14 postgis  raster  map-algebra  postgis-raster 

閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか？この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 閉まっている 3年前にました。 私は、Google Earthプラグインを彼のWebサイトに統合するように依頼しました。 Google Earthプラグインサイトからこれを読みました。 Google EarthプラグインAPIは、2014年12月12日をもって廃止されました。APIは、2015年12月12日までサポートされているブラウザで引き続き動作し、その日までシャットダウンします。 どのような選択肢がありますか？ Google Maps API、Google Maps Engine API、Google Street View APIがあることは知っていますが、それらとGoogle Earthプラグインを「代替」できるものについて混乱しています。 必要な機能は次のとおりです。重要度の高い順に示します。 サテライトビューでKMZファイルを読み込む KMZファイルエンティティに存在するHTMLデータを含むバルーンを表示します。htmlにはストリートビューへのリンクが含まれます。バルーン内でストリートビューを有効にしたいと思います。Google Earthでなんとかできました。 衛星画像履歴画像セレクター 測定ツール クリックされたポイントのアドレスを取得します（geocooding） 3Dの建物

14 google-earth  google-maps-api  api  software-recommendations  google-maps-engine 

Rを使用してtwitterをデータマイニングすることに関するいくつかの研究を行ってきましたが、私の質問に対する答えや適切なチュートリアルは実際には見つかりませんでした。 特定の時間枠内で、特定のハッシュタグを使用してTwitterからツイートを取得し、QGISまたはArcMapのいずれかのマップにそれらのツイートの場所をプロットすることに興味があります。 ツイートにはジオロケーションが関連付けられている可能性があることは知っていますが、そもそもこの情報を抽出するにはどうすればよいですか？

14 python  point  twitter