地理情報システム

それはばかげた質問のように見えますが、それについて明確な答えを見つけることができません：どのユニットがジオパンダ/シェイプを使用してオブジェクト間の距離/面積を計算しますか？

13 python  shapely  geopandas 

私はVincentyのアルゴリズムの実装を使用しており、「位置」と呼ばれるクラスを持っています。このクラスは、「楕円体の表面の上の高度」で定義される緯度、経度、および高度を必要とします。 そのようなアプリケーションでの意味は何ですか？単純な「MSLを超える高度」ですか？

13 elevation  datum  gis-principle  ellipsoid  vincenty-formulae 

私が知っているOL3（OpenLayersを3）に動作していないラインダッシュ以前に頼まれたのが、提案された解決策は機能していません。 私のコードは次のとおりです。 var calculateStyle= new ol.style.Style({ fill: new ol.style.Fill({ color: [255, 255, 255, 0.6], opacity: 0.3 }), stroke: new ol.style.Stroke({ color: [0, 153, 255, 1], linedash: [40,40], width: 3 }) }); 私は両方を試してみました[4,4]し、[40,40]リンク質問で提案されています。しかし、それは機能していません。

13 javascript  openlayers 

これは私を怒らせています。レイヤーコントロールをマップに追加するためのLeafletチュートリアルを参照してください。 L.Icon.Default.imagePath = "Scripts/images"; var mapUrl = 'http://{s}.tile.openstreetmap.se/hydda/full/{z}/{x}/{y}.png'; var secondMap = 'http://{s}.tile.osm.org/{z}/{x}/{y}.png'; var littleton = L.marker([39.61, -105.02]).bindPopup('This is Littleton, CO.'), denver = L.marker([39.74, -104.99]).bindPopup('This is Denver, CO.'), aurora = L.marker([39.73, -104.8]).bindPopup('This is Aurora, CO.'), golden = L.marker([39.77, -105.23]).bindPopup('This is Golden, CO.'); var cities = L.layerGroup([littleton, denver, aurora, golden]); var grayscale …

13 javascript  leaflet  mapbox 

tifイメージファイルと、.tfwワールドファイル（OSMデータからMapertiveで作成）があります。（Linuxコマンドラインを使用して）どうすればそれをGeoTIFFファイルに変換できますか？（つまり、2つではなく1つのファイルが必要です）。私はそれがgdalまたは何かで可能だと確信しています。

13 gdal  geotiff-tiff 

このコードで属性名でフィーチャ属性を設定しようとしました pr = layer.dataProvider() pr.addAttributes([QgsField("test", QVariant.Int)]) layer.updateFields() for feature in layer.getFeatures(): attrName = 'test' feature[attrName] = 1 しかし、結果としてすべてのフィールドにNULLがあります。私が使用する場合 feature.setAttributes([1]) 正常に機能します。 私が間違っていることは何ですか？feature [attrName] = 1が機能しないのはなぜですか？ 更新：このソリューションを見つける pr = layer.dataProvider() pr.addAttributes([QgsField("test", QVariant.Int)]) layer.updateFields() for feature in layer.getFeatures(): attrName = 'test' feature[attrName] = 1 pr.changeAttributeValues({feature.id() : {pr.fieldNameMap()[attrName] : 1}})

13 pyqgis 

Qgisとデザイナーがどのように連携するかを理解しようとしています。Nathan Wの優れた投稿から、ユーザーフォームを作成する方法、qgisレイヤーにフックする方法、および検証コードを追加する方法を知っています。ここでの私の質問は非常に単純ですが、誰もが知っているものは、まさに印刷物では見つけることができないものです。 QGISに付属のQt Designerのコピーは、とにかくQGIS GUIと統合されていますか（ms-access内のフォームデザイナが残りのms-accessのように）？ QGIS内からQt Designerを起動できますか、それともWindowsエクスプローラーでQGISインストールを探し回って、ダブルクリックするQt Designer exeを見つける必要がありますか？ 統合されたQt Designerを使用するよりも、Qt Designer / Qt CreatorまたはpyQtの別個のバージョンをインストールした方が良いですか？ 新規ユーザー向けのパッケージ済みWindowsインストーラーには、フォームデザインに使用できるように構成されたQt Designerが付属していますか？ Python検証コードでフォームを作成する方法を探している場合、PyQtは良い解決策ですか？ 多くのデータ収集タスクは非空間的です。たとえば、道路標識の価格リストに追加します。GISは、どのサインを交換する必要があるかを記録できますが、コストの見積もりを提供するには価格表が必要です。とにかく非空間データテーブルのフォームの表示をトリガーする方法はありますか？

13 qgis  qt4  qt-designer 

GDALを使用して、ラスタ上のすべての値が0のピクセルをnodataに設定する方法はありますか？ fillnodataを使用していくつかの穴を埋めたいのですが、穴が0に設定されているため、最初にそれらを削除する必要があります

13 raster  gdal 

地震を表示する便利な方法を見つけようとしているときに、管理がそれほど簡単ではないことに気付きました。 ポイントのスタイルから「マグニチュード」、「日付」、および「深度」を認識できるようにしたいと思います。何かのようなもの： サイズ：大きさに基づく 深度：基本色（例：<15赤、15-30黄色、> 30青） 日付：100％の色から白/ ...への陰影としての色の強度 何十ものルールを定義せずにそうする方法はありますか？

13 qgis  cartography  style 

私はpythonスクリプトを書いています。mxdがあり、そのワークスペースは必要なすべてのファイルを含むジオデータベースファイルです。コードのテストと作成で直面している問題は、コマンドにもかかわらず、 arcpy.env.overwriteOutput=True ジオデータベースに含まれるすべてのファイルは上書きできません。これを回避する方法はありますか？

13 arcpy 

私の問題は次のとおりです。ラインシェップ（ターゲットレイヤー）とポイントシェップ（結合レイヤー）があります。2番目（ポイントレイヤー）は、ラインに最も近い距離に基づいてラインレイヤーと結合する属性を保持します。点と線は交差していません ArcGISと空間結合（一致オプション：最も近い）でこれを行うことができますが、QGISでこれをどのように達成できるかわかりません。

13 qgis  attribute-joins  attribute-table 

背景として、シェープファイルの使用から移行することを検討しています。デスクトップ作業用のジオパッケージ（QGIS / ArcGIS）の使用を開始しましたが、多くの場合、PostgreSQLデータベースにデータをアップロードする必要があります。 QGISを使用してGeopackageコンテンツをPostgresにアップロードできることを知っていますが、ジオパッケージまたはジオパッケージ内のテーブルをPostgresにアップロードするためのshp2pgsqlのようなコマンドラインユーティリティがあるかどうか疑問に思いますか？おそらくogr2ogrを想定していますが、ドキュメントが見つかりません。

13 postgis  postgresql  ogr  ogr2ogr  geopackage 

私は持っています GeoPandas 0.2.1 py27_0 fiona 1.7.0 np110py27_1 Python 2.7.10 Anaconda 2-4.1.1-Windows-x86環境にインストールされます。GeoDataFrame入力データセットを読み取ってデータを操作して構築できますが、出力データセットを保存しても座標系は保持されません。 import geopandas as gpd world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres')) world.to_file(driver='ESRI Shapefile',filename=r'C:\GIS\Temp\world_out.shp') world_out.shp特徴およびソースシェープファイルのスキーマを含むが、座標ないシステムが定義されていません（<Unknown>）。world_out.prjファイルが0キロバイトで、何が含まれていません。 out = gpd.read_file(r'C:\GIS\Temp\world_out.shp') print out.crs >>> {} 純粋なフィオナを使用してこの操作を実行するようにテストしましたが、正常に機能し、crs情報を保持します。 import fiona datasrc_path = r'C:\Program Files (x86)\ArcGIS\Desktop10.4\TemplateData\TemplateData.gdb' with fiona.drivers(): with fiona.open(datasrc_path,layer='cities') as source: meta = source.meta meta['driver'] = 'ESRI Shapefile' meta['layer'] …

13 python  fiona  geopandas  pandas  anaconda 

シェープファイル、フィオナ、ogrなどのライブラリを使用してコードのいくつかの例を試して、ポイント（x、y）がArcMapで作成された（したがって、シェープファイル形式で）マルチポリゴンの境界内にあるかどうかを確認しました。ただし、通常の単一ポリゴンシェープファイルでは問題なく機能しますが、マルチポリゴンではうまく機能しない例もあります。私が試したスニペットの一部を以下に示します。 # First example using shapefile and shapely: from shapely.geometry import Polygon, Point, MultiPolygon import shapefile polygon = shapefile.Reader('shapefile.shp') polygon = polygon.shapes() shpfilePoints = [] for shape in polygon: shpfilePoints = shape.points polygon = shpfilePoints poly = Polygon(poly) point = Point(x, y) # point in polygon test if polygon.contains(point): print …

13 shapefile  ogr  shapely  point-in-polygon 

緯度と経度の列を含む、80万件を超えるレコードを持つテーブルがあります。データソースはジオメトリを作成しなかったため、これを行う方法を探していました。ST_MakePoint（x、y、SRID）があることは知っていますが、すべてのドキュメントと他のスレッドは、単一のポイントを作成する方法を示しています。これらの列を使用してジオメトリタイプを作成するにはどうすればよいですか？

13 postgis  sql