地理情報システム

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

1
QGISエクスポートコンポーザーをPythonでPDFとしてエクスポート
作曲家を含むQGISプロジェクトがあります。もちろん、GUIに移動すると、pdfなどとしてエクスポートできます。makeスクリプトの場合は、Pythonスクリプトを使用してこれを実行します。ここでオンラインで読むことから私の試みです： #!/usr/bin/env python import sys from qgis.gui import QgsMapCanvas, QgsLayerTreeMapCanvasBridge from PyQt4.QtCore import QFileInfo from PyQt4.QtXml import QDomDocument from qgis.core import * qgs = QgsApplication(sys.argv, True) QgsApplication.setPrefixPath("/usr", True) qgs.initQgis() def make_pdf(): canvas = QgsMapCanvas() bridge = QgsLayerTreeMapCanvasBridge( QgsProject.instance().layerTreeRoot(), canvas) bridge.setCanvasLayers() QgsProject.instance().read(QFileInfo('../board.qgs')) composition = QgsComposition(canvas.mapSettings()) map_item = composition.getComposerItemById('board36x48') map_item.setMapCanvas(canvas) map_item.zoomToExtent(canvas.extent()) composition.refreshItems() …
11 qgis  python  pyqgis  pdf 
1
GDALは並列処理をサポートしていますか？
このコマンドに対応するプロセスを高速化したい： gdalwarp -t_srs EPSG:4326 -overwrite input.ntf output.tif GDALで並列処理を使用する方法はありますか？セミオフトピック：そうでない場合、gdalwarpのようなプロセスを高速化するために非GDALソリューションをお勧めしますか？ 私が見たウェブページは次のとおりです。 GDALとQGISで並列処理を行うことはできますか？ [gdal-dev] GDALラスター処理：並列計算 編集：gdalwarpの処理中にCPU％が100％を下回っていると思ったため、この質問をしました。ただし、再確認すると、555.5％に達しました。
2
QGISでレイヤーパネルの幅を設定しますか？
QGISを開くと、レイヤーパネルのデフォルトの幅が[レイヤー/グループの削除]ボタンを表示しません。したがって、QGISを開くたびにボタンを表示するようにパネルのサイズを変更する必要があります。 このパネルのデフォルトの幅を設定して、サイズを変更する必要はありませんか？ Windows 10でQGIS 2.18.5を使用しています。
3
Pythonで複数のポリゴンの交差点を効率的に取得する
複数のポリゴンの交差点を取得したいです。Pythonのshapelyパッケージを使用すると、intersection関数を使用して2つのポリゴンの交差を見つけることができます。複数のポリゴンの交差を取得するための同様の効率的な機能はありますか？ これが私の意味を理解するためのコードスニペットです。 from shapely.geometry import Point coord1 = ( 0,0 ) point1 = Point(coord1) circle1 = point1.buffer(1) coord2 = ( 1,1 ) point2 = Point(coord2) circle2 = point2.buffer(1) coord3 = ( 1,0 ) point3 = Point(coord3) circle3 = point3.buffer(1) 2つの円の交点はで見つけることができますcircle1.intersection(circle2)。で3つすべての円の交差点を見つけることができますcircle1.intersection(circle2).intersection(circle3)。ただし、このアプローチは、より多くのコードを必要とするため、多数のポリゴンに対応できません。任意の数のポリゴンを取り、それらの交差を返す関数が欲しいです。
1
r.mapcalcはQGISで出力を提供しません
しばらくの間、私はQGIS Processing Toolboxからr.mapcalcを使用しようとしています。ただし、このアルゴリズムでは何も出力されません。いくつかのラスターレイヤーと異なる式を試しましたが、アルゴリズムは「結果ディレクトリ」に何も出力しません。以下は、アルゴリズムのセットアップの簡単な例を示す画像です。たぶん私は何かを見逃していますが、似たようなオンラインは見つかりませんでした。
2
パンバンドを使用したパンシャープンランドサット反射率製品は科学的に正しいですか？
Landsat 8表面反射率製品とその各帯域のパンバンドを融合/パンすることが科学的に有効かどうかを知る必要がありますか？Landsat反射率製品の詳細については、こちらをご覧ください。この製品を入手するには、表面反射率製品を個別に注文する必要があることに注意する必要があります。この製品には、IRおよびパンバンドではなく、7バンド（30m）しか含まれていません。それで、私の質問は、表面反射率積の7バンド（30m）を通常の（表面反射率ではない）パンバンド（15m）と融合することは妥当です。私はこのパンシャープン画像をセグメンテーションと土地被覆マッピングに使用したいので、学界でこのタイプのパンシャープニングの実践が確立されていることを参照する必要があります。
3
GDALを使用してシェープファイル内のポリゴンの総面積を計算しますか？
British National Gridプロジェクションにシェープファイルがあります。 Geometry: 3D Polygon Feature Count: 5378 Extent: (9247.520209, 14785.170099) - (638149.173223, 1217788.569952) Layer SRS WKT: PROJCS["British_National_Grid", GEOGCS["GCS_airy", DATUM["OSGB_1936", SPHEROID["Airy_1830",6377563.396,299.3249646]], PRIMEM["Greenwich",0], UNIT["Degree",0.017453292519943295]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",49], PARAMETER["central_meridian",-2], PARAMETER["scale_factor",0.9996012717], PARAMETER["false_easting",400000], PARAMETER["false_northing",-100000], UNIT["Meter",1]] cat: Integer (9.0) GDAL / OGRを使用して、シェープファイル内のすべてのポリゴンの総面積をヘクタール単位で取得できますか？ 私はこれが可能かどうか疑問に思っています-sql： ogrinfo -sql "SELECT SUM(ST_Area(geom::geography)) FROM mytable" myshapefile.shp しかし、私はそれをしようとしていますERROR 1: Undefined function 'ST_Area' …
11 gdal  ogr  area  ogrinfo 
2
特定のLong / Lat座標のSentinelタイルを見つける
「Sentinel-2 on AWS」から画像をダウンロードするツールを構築しています Long / Lat座標を特定のSentinelタイルに変換するのに問題があります。Sentinel 2はMGRSを使用していることを知っており、Long / LatをMGRSに単純に変換しようとしましたが、低精度でS2AタイルIDを取得しました。そして、それは動作しますが、常にではありません。 たとえば、Long / Lat 34.665,31.625はMGRSの「36R XA」に解決されますが、Sentinelタイル36RXAは存在しません。 Long / Lat座標を使用してSentinelタイルを決定する正しい方法は何ですか？
1
境界ボックス内のすべてのポイントを選択します
この質問は以前に聞いたことがありますが、私が試みている答えは機能していません。 境界ボックス内にあるすべてのOSMウェイを照会したいと考えています。OSMデータは、デフォルトの球状メルカトールでインポートされました。LAT / LONを使用してクエリを実行しているため、変換 SELECT * FROM planet_osm_ways WHERE ST_Transform( ST_MakeEnvelope(37.808264, 37.804339, -122.275244, -122.271189, 4326),3857 ); これを実行すると、次のエラーが表示されます。 エラー：WHEREの引数は、タイプジオメトリではなくブール型でなければなりませんLINE 3：WHERE ST_Transform（ST_MakeEnvelope（37.808264、37.804339、-12 ...
2
巨大なXYZ CSVをGeoTIFFに変換する
私はUTM座標を含むCSV形式のデータの膨大な量を持っているXとYように標高値Z情報。さらに分析するには、これらのデータをGeoTIFFとしてDEMに変換する必要があります。この場合、膨大な量は16 mを意味します。行、1つのポイントX、YおよびZ行ごと。点は均等に分布しているため、補間は必要ありません。各ポイントをラスタセルに変換する必要があります。 元のデータにはセパレータがなく、固定列幅が使用されていました。ストリームテキストエディターsedを使用して、固定幅の代わりにセパレーターを使用し、すべてのスペース文字を削除するようにファイル構文を変換する方法を既に考えました。ここからは、通常、私のワークフローからフィーチャクラスを作成することにより、ArcGISのデータをインポートするだろうX、YとZ使用して、データと第2段階としてのGeoTIFFにポイントシェープファイルを変換ポイントラスタにツールを。しかし、私が現在持っているファイルは、このプロセスには大きすぎます。 上記のワークフローの代わりに、効率的な代替手段を探していて、GDALを発見しました。ただし、ではgdal_translate、サポートされているファイルタイプリストで最も近いサポートされている形式はASCIIグリッドですが、カンマ区切りのXYZはありません。もう1つの難点は、UTM座標を持っていることですが、ほとんどの例では10進度座標を使用しているようです。ただし、UTMシステム内にとどまる必要があります（または少なくとも、出力GeoTIFFはUTM座標系にある必要があります）。 したがって、GDALを使用してCSV XYZをGeoTIFFに変換する方法を探していますが、これまでのところ、この正確な問題に対処する例を見つけることができませんでした。いくつかのヒントやコード例があればとても嬉しいです。
2
ogr2ogrは複数のシェープファイルをマージします：-nlnタグの目的は何ですか？
サブフォルダーを再帰的に繰り返し、すべてのシェープファイルを単一のシェープファイルにマージするための基本的なスクリプトは次のとおりです。 #!/bin/bash consolidated_file="./consolidated.shp" for i in $(find . -name '*.shp'); do if [ ! -f "$consolidated_file" ]; then # first file - create the consolidated output file ogr2ogr -f "ESRI Shapefile" $consolidated_file $i else # update the output file with new file content ogr2ogr -f "ESRI Shapefile" -update -append $consolidated_file …
1
これらはどのような緯度/経度座標ですか？
私は受け取った座標セットをマップしようとしています： N 39,21.440 W 120,48.076 しかし、私はフォーマットに精通していません。DMSのlat / longsと標準（？）形式（39.123234、-121.348348など）に精通している これらをDMSまたは標準に変換する方法を教えてください。
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.