地理情報システム

地図製作者、地理学者、およびGIS専門家のためのQ&A

5
numpy配列をラスターファイルに書き込む
GISは初めてです。 火星の赤外線画像を熱慣性マップに変換するコードがあり、それを2D numpy配列として保存します。これらのマップをhdf5ファイルとして保存してきましたが、QGISで処理できるように、ラスターイメージとして保存したいのです。これを行う方法を見つけるために複数の検索を行ってきましたが、運はありません。http://www.gis.usu.edu/~chrisg/python/のチュートリアルの指示に従ってみましたが、彼のサンプルコードを使用して作成したファイルは、QGISにインポートするときにプレーンな灰色のボックスとして開きます。誰かが私がやりたいことの簡単な例に可能な限り簡単な手順を提案できれば、ある程度進歩することができるかもしれません。私はQGISとGDALを持っています。だれでもお勧めできる他のフレームワークをインストールできてとてもうれしいです。Mac OS 10.7を使用しています。 たとえば、次のような熱慣性のnumpy配列がある場合: TI = ( (0.1, 0.2, 0.3, 0.4), (0.2, 0.3, 0.4, 0.5), (0.3, 0.4, 0.5, 0.6), (0.4, 0.5, 0.6, 0.7) ) そして、ピクセルごとに緯度と経度があります: lat = ( (10.0, 10.0, 10.0, 10.0), ( 9.5, 9.5, 9.5, 9.5), ( 9.0, 9.0, 9.0, 9.0), ( 8.5, 8.5, 8.5, 8.5) ) …
30 python  raster  gdal  numpy 

4
QGISで列に連続番号を入力しますか?
属性テーブルに新しい列を作成し、デフォルト値(たとえば0)を設定します。(おそらくフィールド計算機を使用して)番号付きインデックスのように、行1の値が1、行2 = 2などになるという結果を得たいと思います。 フィールド計算機のコード例に感謝します(($rownum回線に切り替えるためにこれが必要だと思います)の使用を含む)。

9
ポリゴン内にあるポイントを生成する
ポリゴンフィーチャがあり、その中にポイントを生成できるようにします。これは1つの分類タスクに必要です。 ポリゴン内にランダムポイントを生成するのは、時間がかかるので予測できないため、機能しません。
30 python  qgis  polygon  ogr 

6
QGIS、Postgis:ジオメトリタイプが列タイプと一致しません
SPITを介してQGISのPostgisにいくつかのポリゴンシェープファイルをインポートしようとしています。そのうちの1つはインポートできず、このエラーを返します。 ERROR: Geometry type (Polygon) does not match column type (MultiPolygon) SPIT widnowのそのテーブルの属性「Feature Class」をMULTIPOLYGONからPOLYGONに変更しようとしましたが、何も起こりませんでした。 Postgisに両方のタイプ(ポリゴンとマルチポリゴン)を受け入れるようにする方法はありますか?または、どういうわけかシェイプファイルジオメトリをポリゴンからマルチポリゴンに変換できますか?何か案は? Qgis 2.0.1をWindowsで試し、Qgis 2.3をArchLinuxで試しました。
30 qgis  postgis  database 

3
中国からの出国が許可されていない地図資料
私は、ドイツのさまざまな国で自動運転車をテストする会社についての記事を読みました。 彼らは、ドイツで生成されたすべてのデータを分析すると述べています-1つの例外を除いて:中国。 「法的規制のため、地図資料は出国できません。」中国の地図を数秒間グーグルで検索すると、いくつかの結果が示されます。 誰でもこの矛盾を説明できますか?
29 license  china 

2
QGISの点から線を引く
ポイントレイヤーがあり、ポイントを通過する、またはポイントの近くにラインを描画する必要があります(ポイントとラインの間の距離は最大0.5m)。線のノードは、それらに最も近いポイントの属性を受け取ります。

4
GDALの学習リソース?
特定の問題を解決しようとしたり、構文を正しくしたりするのにしばしば動けなくなります。 GDAL Utilitiesのリンクで提供されている詳細は別として、サードパーティのリソース/マニュアルを知っている人はいますか?
29 gdal  references 

5
Rの多数のポリゴンによる切り抜き、マスク、およびラスターの抽出の速度を上げますか?
数千のポリゴン境界に基づいて、さまざまな土地利用タイプの面積と被覆率をラスタから抽出しています。個々のポリゴンを繰り返し処理し、特定のポリゴンのサイズまでラスターをマスクして切り取ると、抽出機能がはるかに高速に動作することがわかりました。それにもかかわらず、それはかなり遅いです、そして、誰かが私のコードの効率と速度を改善するための提案を持っているかどうか疑問に思っています。 私はこれに関連した見つけた唯一のものです。この応答使うことを提案ロジャーBivandによるGDAL.open()とGDAL.close()だけでなく、getRasterTable()とgetRasterData()。私はそれらを調べましたが、過去にgdalに問題があり、それを実装する方法を知るのに十分なほどよく知りません。 再現可能な例: library(maptools) ## For wrld_simpl library(raster) ## Example SpatialPolygonsDataFrame data(wrld_simpl) #polygon of world countries bound <- wrld_simpl[1:25,] #name it this to subset to 25 countries and because my loop is set up with that variable ## Example RasterLayer c <- raster(nrow=2e3, ncol=2e3, crs=proj4string(wrld_simpl), xmn=-180, xmx=180, ymn=-90, ymx=90) c[] …

5
地形の曲率を計算する方法は?
標高ラスターがあり、勾配、頑丈さ、アスペクト、曲率などの基本的な地形変数を取得したい。 問題は、テレインアナリストで曲率を計算する可能性がないように見えることです。 このインジケータがリストにない理由を誰にも教えてもらえますか?それ以外の場合はどうすれば計算できますか?
29 qgis  dem  terrain  curvature 

6
2つの円の交差を計算しますか?
中心の緯度/経度と各点の半径を指定して、地球上の2つの交差する円の共通点を数学的に導き出す方法を見つけようとしています。 たとえば、次の場合: 緯度/経度(37.673442、-90.234036)半径107.5 NM 緯度/経度(36.109997、-90.953669)半径145 NM 2つの交点を見つけて、そのうちの1つ(36.948、-088.158)を見つけます。 これを平面で簡単に解決することは簡単ですが、地球の表面などの不完全な球で方程式を解く経験はありません。

4
Linux Ubuntu ServerにGDALをインストールしますか?
apt-get install python-gdalバージョン1.6.0をインストールしますが、これはすでにインストールされています。1.9.0バージョンをインストールしたいと思います。 最初に1.6.0をアンインストールしてから手動でコンパイルする必要がありapt-getますか、それとも見逃したパッケージがありますか?

5
GeoServerレイヤーを認証の背後に配置する方法は?
OpenLayersとGeoServerを使用してマップレイヤーを公開したい。これは簡単に設定できますが、認証の可能性に関しては少し迷っています。私が欲しいのは、独自のユーザー名とパスワードでログインすることで、独自のマップウィンドウにリダイレクトされるさまざまなユーザーを作成することです。どうすればこれを安全に実行して、ユーザーがお互いのレイヤーにアクセスできないようにすることができますか? たとえば、利用可能なすべてのレイヤーを一覧表示するように要求するGeoServerに対してクエリを実行することにより、セットアップを「ハッキング」することは非常に簡単であることがわかりました。 私のオプションは何ですか?私のubuntuサーバーでApacheをセットアップする際にApacheが役立つと思いますか?このテーマに関するチュートリアルはありますか?

4
空間インデックスを使用すべきでない場合
私は主にOracleで働いていたのでこれを求めていますが、過去1年間はPostGISとSQLServer 2008で倍増しています。Oracleのほとんどの空間関数はORA-13226エラーを返す空間インデックスなしでは機能しません。 13226、00000、「インタフェースは空間索引なしではサポートされていません」原因:ジオメトリ表に空間索引がありません。// *アクション:空間演算子で参照されるジオメトリテーブルに空間インデックスがあることを確認します。 これは理にかなっています。空間クエリを実行します=空間インデックスが必要です。しかし、私が理解している限りでは、SQL ServeではなくPostGISもこれを必要としません。PostGISには、空間インデックスを明示的に使用しない関数(_ *例_STContains)があるようです。 問題は、空間インデックスを使用してはならないケースがあるかどうかです。「取る」か「残す」かは必ずしも違いません。つまり、違いはありませんが、空間インデックスを使用しないとパフォーマンスが低下します。私にとって、最後の文は用語の矛盾ですが、そうでなければPostGISがこれらの機能を提供するのはなぜですか?

5
QGISの最も近いポイントから属性を結合しますか?
ポイントを持つ2つのレイヤーがあり、それぞれ異なる属性を持つAとBと呼ばれます。 AIの各ポイントについて、Bの最も近いポイントからのすべての属性を結合します。 理想的には、これには両方のポイント間の距離を測定するフィールドも含まれます。 ArcGISでこれを行う方法を知っています:レイヤーA> joins and relateds> joinを右クリックし、空間位置に基づいてデータを結合し、最も近いポイントから属性を結合することを選択します。 QGISで同様の機能を見つけることができませんでした。ポインタはありますか?

8
QGISを使用してコントロールポイントでベクターレイヤーをジオリファレンスしますか?
ジオリファレンスする必要のある非ジオリファレンスベクターレイヤーがあります。ラスターレイヤーを使用すると、タスクは簡単で簡単ですが、ベクターレイヤーで何をすべきかわかりません。既知の座標を持ついくつかのコントロールポイントがあり、それらは変換の基盤を提供するはずです。したがって、1、2、および3のid-sのポイントがx1、y1の座標を持つ必要があることを知っているとしましょう。x2、y2; x3、y3。単純なシフトに加えて、回転およびスケール変換が行われる場合があります。 何か案は?

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.