地理情報システム

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

1
Mapboxエラー:スタイルの読み込みが完了していません
私は、suncalcの結果に基づいてgeojsonを使用して、マップボックスマップに描画しようとしています。最初に、描画しようとしていた各行に1つずつ、2つの関数を作成しようとしました。しかし、それを行うと、最後に呼び出された関数のみが表示されます。だから私は、mapboxとリーフレットの構文は初めてなので、レイヤーの管理方法がわからないと結論付けました。 リーフレットの例を使用して、次のコードに到達しました。 ジャバスクリプト function drawOnMap(sAz){ //sun calculation stuff, doesn't matter for the problem var lt = parseFloat(document.getElementById('lat').value); var ln = parseFloat(document.getElementById('long').value); var R = 6371000; //Earth's radius var d = 10000* Math.sin(1); // Distance var brngSR = deg2rad(sAz); var ltr = deg2rad(lt); var latSR = rad2deg(Math.asin( Math.sin(ltr)*Math.cos(d/R) + Math.cos(ltr)*Math.sin(d/R)*Math.cos(brngSR) )); …

1
QGIS、自己更新式(スプレッドシートなど)を使用してフィールドに入力する
fieldの値として数式を入力する方法(プラグインなど)はありますか?基づいている他のフィールドが調整されると、自己更新しますか?フィールド計算機がフィールドの値を更新する方法とまったく同じですが、毎回フィールド計算機に戻る必要はなく、依存するレイヤーが更新されるとスプレッドシートのように機能します。 たとえば、ポリゴン内の住宅ユニットの数を調べたい場合、フィールド計算機でのようなことができますround("Density"*"Area")が、ポリゴンを調整するたびにそれをやり直す必要があります。計算がより複雑になると、手動で実行して手動で再調整するのが非常に面倒になります。

2
ArcCatalogでISO 19139メタデータプロファイルをカスタマイズしますか?
ISO 19115/19139メタデータ標準を使用しているため、一部のオプション要素を必須要素に変換する必要があります。ArcCatalogをメタデータエディターとして使用してこれを実現することは可能ですか? メタデータの規格とプロファイルについて「と言うプロファイルは、さらにそれが必須ところ、それはオプションでした前になって、メタデータ要素の選択性を制限することができます。」 ArcCatalogでは、メタデータ標準としてISO 19139を設定できます。 しかし、ISO 19139のインスタンスを編集して(カスタムプロファイルを作成する)ため、(たとえば)系統がオプションではなく必須になりますか? 下でのArcGIS \ Desktop10.4 \メタデータ\翻訳\ルールファイルがあるISO19139_min_schema.txtが。 このファイルを編集できますか?その場合、系統を必須にする構文は何ですか?

3
隣接するラスターセルへの値の条件付き割り当て?
値ラスターがあります: m <- matrix(c(2,4,5,5,2,8,7,3,1,6, 5,7,5,7,1,6,7,2,6,3, 4,7,3,4,5,3,7,9,3,8, 9,3,6,8,3,4,7,3,7,8, 3,3,7,7,5,3,2,8,9,8, 7,6,2,6,5,2,2,7,7,7, 4,7,2,5,7,7,7,3,3,5, 7,6,7,5,9,6,5,2,3,2, 4,9,2,5,5,8,3,3,1,2, 5,2,6,5,1,5,3,7,7,2),nrow=10, ncol=10, byrow = T) r <- raster(m) extent(r) <- matrix(c(0, 0, 10, 10), nrow=2) plot(r) text(r) このラスターから、この図に従って現在のセルの8つの隣接するセルに値を割り当てる(または値を変更する)にはどうすればよいですか?このコード行から現在のセル内に赤い点を配置しました。 points(xFromCol(r, col=5), yFromRow(r, row=5),col="red",pch=16) ここで、期待される結果は次のようになります。 現在のセルの値(つまり、値ラスタの5)は0に置き換えられます。 全体として、8つの隣接セルの新しい値は次のように計算する必要があります。 新しい値=赤い長方形に含まれるセル値の平均*現在のセル(赤い点)と隣接するセル間の距離(つまり、斜めに隣接するセルの場合はsqrt(2)、そうでない場合は1) 更新 隣接セルの境界がラスターの制限を超えている場合、条件を尊重する隣接セルの新しい値を計算する必要があります。条件を尊重しない隣接セルは「NA」に等しくなります。 たとえば、[row、col]表記を使用して参照位置がc(5,5)ではなくc(1,1)である場合、右下隅の新しい値のみを計算できます。したがって、期待される結果は次のようになります。 [,1] [,2] [,3] [1,] NA NA NA [2,] NA …
12 raster  r  adjacency 


1
ogr2ogrを使用してpostgisにシェープファイルをインポートすると、次のようになります。データソースを開くことができません
postgisデータベースにシェープファイルをインポートするには、ogr2ogrを使用します。ogr2ogrを正常にインストールし、pgsqlから次のコマンドを実行します。 ogr2ogr -f "PostgreSQL" PG:"host=localhost user=user_1 password=***** dbname=imports" world_boundaries.shp 私が返すのはエラーメッセージです: Unable to open datasource `world_boundaries.shp' with the following drivers: --a list of drivers follows (ESRI Shape File etc.) また、シェープファイルの完全なパスを定義しようとしましたが、同じメッセージが表示されました。 また、私は実行しようとしました: ogrinfo world_boundaries.shp 同じこと。 ファイルの権限に関する問題を修正した後、次のエラーが表示されます。 ERROR 1: AddGeometryColumn failed for layer world_boundaries, layer creation has failed. ERROR 1: Terminating translation prematurely …

4
リーフレットマーカーのポップアップ位置を変更しますか?
Leafletのマーカーアイコンの下部にポップアップを開きます。 私のコード: var mymap = L.map('mapid').setView([51.505, -0.09], 13); L.tileLayer('https://api.tiles.mapbox.com/v4/{id}/{z}/{x}/{y}.png?access_token=pk.eyJ1IjoibWFwYm94IiwiYSI6ImNpandmbXliNDBjZWd2M2x6bDk3c2ZtOTkifQ._QA7i5Mpkd_m30IGElHziw', { maxZoom: 18, attribution: 'Map data © <a href="http://openstreetmap.org">OpenStreetMap</a> contributors, ' + '<a href="http://creativecommons.org/licenses/by-sa/2.0/">CC-BY-SA</a>, ' + 'Imagery © <a href="http://mapbox.com">Mapbox</a>', id: 'mapbox.streets' }).addTo(mymap); var lati = 51.51; var longi = -0.09; var popupLocation1 = new L.LatLng(lati, longi); var popupContent1 = 'This …
12 leaflet  popup 

2
QGISフィールドにデフォルト値を設定する
レイヤーのフィールドにデフォルト値を追加したいと思います。つまり、新しい機能が作成されるたびに、そのフィールドにはデフォルト値が自動的に入力されます。 私の場合、デフォルト値は@project_filename変数(プロジェクトのファイル名)です。 この機能はどこにもありません。 よろしく、


1
LeafletでlayerGroupの境界を取得するにはどうすればよいですか?
複数のgeoJSONレイヤーを最初にLayerGroupに追加してから、layerGroupをマップに追加して、マップに追加しました。 今、私はすべてのレイヤーの結合された境界でマップ上のfitBoundsを呼び出したいです、問題はlayerGroupにメソッドgetBoundsがないことです 私に何ができる?
12 leaflet 

1
Googleマップで2つのポリゴンが重なっているかどうかを判断しますか?
Googleマップとポリゴンを使用していますが、ポリゴンが互いに重なり合っていないことを確認する必要があります。 ポイント(lat、lng)がポリゴン内にあるときにTRUEを返す関数がありますが、ポリゴンのすべてのポイントが別のポリゴン内にあるかどうかを判断するだけでは不十分です。 なにか提案を?

5
分散写真からオブジェクトの写真測量ポイントクラウドを作成する、ここ:UAV画像を使用して建物の正確な3Dモデルを作成しますか?
オブジェクトの密集して分布した写真のセットから写真測量の点群を作成する方法はありますか? 応用: UAV(Phantom 3 Pro)の画像から建物の外側の3Dモデル(点群)を作成することが、地上雲スキャナー(TLS)を使用して点群を作成するための実行可能な代替手段になるかどうかを調査しようとしています。UAVには12メガピクセルのカメラとGPSがあります。 目的は、比較的近距離で画像を撮影し、高密度画像マッチングソフトウェア(できればPix4D)からポイントクラウドを作成して、建物の周りにUAVを飛ばすことです。計画では、建物を見下ろす通常の直下と斜めの画像を撮影することですが、低高度で飛行し、窓や軒先などを見上げてまっすぐ斜めの写真を撮影して、できるだけ多くの詳細をキャプチャしようとします。

2
衛星データを使用した氷の種類の区別
氷河の領域を雪(存在する場合)と氷のカテゴリで分類したいと思いますが、最も重要なことは、古い氷と新鮮な氷の間です。フィールドで認識できるさまざまなプロパティがありますが、衛星データを使用してこれを実行できますか?(30 / 15mの空間解像度のため、Landsatが望ましい)

1
同じマシンでArcGIS ProとArcGIS DesktopでArcPyを使用していますか?
ArcGIS Proを初めて試したので、PYTHONPATHをpython34に変更する必要がありました 現在、ArcGIS Proは[解析]タブを開くたびにクラッシュすることなく動作しています。しかし、ArcMap 10.3.1でプロセスを実行しようとすると、さまざまなPythonエラーが発生します。同じ正確なデータ入力で30分前に機能したツールについて。 ArcGIS 10.3.1は2.7を使用しているため、これはPythonバージョンの競合ですか? 同じマシン上でArcGIS ProとArcGIS 10.3.1を同時に実行する方法はありますか?


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