地理情報システム

DEM（デジタル標高モデル）とライン（x1、y1）->（x2、y2）（ペアの座標（x、y））が与えられた場合、そのラインに投影された標高プロファイルまたは断面をどのように計算できますか？ 私は自分のプロジェクトで使用できるオープンソースコード、またはガイドとして使用できる擬似コードを探しています。

16 raster  dem  open-source-gis  elevation  profile 

QGIS 2.12.2では、別のレイヤーのフィーチャが既に存在する場所にラベルを配置しないように、レイヤーラベルを設定するにはどうすればよいですか？ たとえば、湖の「中心線」を含むストリーム/川のポリラインレイヤーがあり、描画順序でその上に「湖」のポリゴンレイヤーを配置する場合、川のレイヤーに湖の内側にラベルを配置させたくない。代わりに、私はむしろ（必要に応じて）湖の外で川にラベルを付けたいと思います。そうすれば、lakesレイヤーからラベルを配置でき、ラベルの衝突に遭遇することはありません。 ここに例があります（視覚的な目的のために意図的に線を上に置いています）達成したいのは、湖のポリゴン内に川の中心線のラベルが表示されないことです：

16 qgis  labeling 

JavaScriptを使用してベクターレイヤーに手動でレイヤーを追加しようとしています。私はこれが失敗している理由を判断できないようです： http://jsfiddle.net/Kieveli/f4t6n6v1/4/ 私は、16,22のような健全な座標と、ビューのxy値に一致する大きな座標を試しました。ol3からjavascriptエラーが表示されます。TypeError：bQは関数ではありません。 HTML： <div id="map" class="map"></div> Javascript： var vectorSource = new ol.source.Vector({}); var map = new ol.Map({ layers: [ new ol.layer.Tile({ source: new ol.source.MapQuest({layer: 'sat'}) }), new ol.layer.Vector({ source: vectorSource }) ], target: 'map', view: new ol.View({ center: [-11000000, 4600000], zoom: 4 }) }); var thing = new ol.geom.Polygon( …

16 vector  openlayers 

追加したい既存のシェープファイルがあります。 現在、私のコードは次のようになっています。 self.iface.addVectorLayer(self.shpFilePath, "Track", "ogr") それはうまくいきますが、コード内でcrsを定義したいと思います。出来ますか？私はこのような似たような質問を見つけました：Pythonを使用してQGIS 2.4でプログラムでCRSを選択してください しかし、とにかく動作しません。

16 coordinate-system  pyqgis  vector 

経時的に同じ現象を表す一連の連続したラスターがあるため、それらすべてで同じシンボルを使用する必要があります。 そこで、最小/最大値からのカラーランプを使用して、最初のラスタのストレッチシンボルを設定します。 次に、2番目のラスタで、最初のラスタからシンボルをインポートします。残念ながら、ラベルは一致していますが、実際の最小/最大値は最初のレイヤーと同じではありません。 したがって、2番目のラスタの赤は、0.66ではなく0.47の値に対応します。 両方のラスターが一貫したシンボルを使用していることを確認するにはどうすればよいですか？

16 arcgis-desktop  arcmap  symbology 

経度緯度値をGEOGRAPHYデータ型として保存した後、個々の緯度経度値を取得するにはどうすればよいですか？ 失敗した試行： SELECT id, geog, ST_X(geog), ST_Y(geog) FROM locations; エラー： No function matches the given name and argument types. You might need to add explicit type casts.

16 postgis  coordinate-system  postgresql 

閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか？この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 4年前に閉鎖されました。 私は、新入社員向けの技術的なGISコンピテンシーマトリックスを開発中です。このマトリックスは、新しい従業員を評価するためだけでなく、従業員の成長を監視するためにも使用されます。マトリックスには、一般的なGISの概念から高度なサーバーサイドGISおよびWeb-GIS開発まで、すべてを含める必要があります。 誰もそのようなマトリックスの経験がありますか、またはこれがどのように見えるかについて光を当てることができますか？ 私の最初の試みはこれをこれを見ます。 GISの一般的な概念：CRSモデルと変換、形式、空間分析、地理参照、ヘッドアップデジタル化、シンボル体系、色 空間データベース：ESRI、Postgresql / PostGIS、Mysql WebGIS：オブジェクト指向プログラミング、Openlayers、サーバーサイドGIS、Gdal、Geoserver / Mapserver / Deegree 一般的なプログラミング：Python、JavaScript、HTML / CSS、JQuery、PHP Linux管理：シェル、スクリプト、コマンド、管理、監視 また、これらのさまざまな基準がどのような重み付けを行う必要があるかを調べることにも興味があります。私の個人的なバイアスは、技術的な側面よりもGIS-Workのコンテキストの性質を重視することです。つまり、最初に技術を達成するよりも、ツールを開発する前にデータを理解することが重要だと思います。データまたはツールが実際に表すものを理解しようとします。このため、後でGISスキルを開発したプログラマーよりも、GISスキルを学んだプランナーを採用します。明らかに、これはコンテキストの性質が強いGISプロジェクトにのみ適用されます。GIS管理者にLinux上のJAVAプロジェクトのみを監視させたい場合は、明らかにITプロフェッショナルを好むでしょう。

16 education 

サイトがあり、サイトの中心の5kmのバッファーに指示線を引きたいと思います。そこで、半径5kmの円を描きたいと思います。円、楕円、長方形、正方形、らせんを描くツールがあります。 これらのどれも中心点と半径を述べさせません。

16 qgis  circle 

フィールド計算機（ArcMap 10.0）を使用して、アドレスフィールド内の二重引用符を削除しようとしています（そこに入れませんでした、信じています）。私は多くの呪文を試してみましたが、計算で地面から降りることさえできません。はい、両方のフィールドは文字列であり、両方とも許容可能な長さです。 この計算： この結果： [結果]ダイアログに次のエラーメッセージが表示されます。 繰り返しますが、そこには入れませんでした。

16 arcgis-desktop  arcgis-10.0  arcmap  field-calculator  python-parser 

QGISで最も近い点までの距離を計算する方法はありますか？ArcGISには、ユークリッド距離と呼ばれる関数があり、入力ベクトルまたはラスターを取得し、距離のラスターを生成できます。 QGISに同様の機能はありますか？ドキュメント内で見つけることができませんでした。RasterQGIS の関数とMMGISHub Distanceからの計算（非ラスターソリューションの場合）を調査しましたが、どちらを使用してこれを行うかはわかりませんでした。ベクトルアプローチが唯一の解決策である場合、それでも問題はありません。

16 qgis  raster  distance  point 

現在、GISソフトウェアでサポートされていないデータ形式からGeoJSONデータを生成するために使用している小さなスクリプトがあります。 出力を検証したいのですが、ogr2ogrを使用してそれを読み取ると、どの側面、行、または文字に無効なデータが含まれているかについての有益な情報が得られません。 どこかにまともなGeoJSONバリデーターがありますか、またはOGRでGeoJSONデータを検証する良い方法はありますか？ 編集：JSONバリデーターまたはGeoJSON仕様を探していません。私は両方を見て、このポイントに到達するために使用しましたが、OGRはまだエラーを返します。 ありがとう！

16 geojson 

私は本当にPythonが初めてで、ArcGIS 10のPythonを使用して、ピクセルごとのラスターのセル値と座標（各ピクセルの中心のXY座標をマップ）を取得する簡単な方法があるかどうかを知りたいですか？ これをさらに説明するには、最初のピクセルのマップX、マップY、およびセル値を取得し、それらの3つの値を3つの変数に割り当て、他のピクセルの残りに対してこの手順を繰り返す必要があります（ラスター全体をループします）。 私の質問をもっと説明する必要があると思います。問題は、最初のラスタのピクセルのXY位置を取得し、そのXY位置に対応するいくつかの他のラスタのセル値を取得する必要があることです。このプロセスは、中間ポイントシェープファイルを作成せずに、最初のラスターのすべてのピクセルをループ処理する必要があります。80億ピクセル近くのラスターを処理する必要があるため、本当に時間がかかります。また、ArcGIS 10でPythonを使用してこれを行う必要があります。 @JamesS：ご提案ありがとうございます。はい、これは1つのラスターで機能しますが、他のいくつかのラスターのセル値も収集する必要があります。問題は、最初のラスタの最初のピクセルのX座標とY座標を取得した後、最初のラスタのX、Y位置、3番目のラスタなどに対応する2番目のラスタのセル値を取得する必要があることです。だから、最初のラスターをループするとき、ピクセルのXとYの位置を取得し、その位置に対応する他のラスターのセル値を取得することは同時に行われるべきだと思いますが、わかりません。これは、最初のラスターをポイントシェープファイルに変換し、ArcGIS 10でポイント値への複数値の抽出を実行することで実行できますが、 @hmfly：ありがとう、はい、このメソッド（RastertoNumpyarray）は、配列の既知の行と列の値の座標を取得できれば機能します。 @whuber：計算を実行したくないので、テキストファイルにXY座標とセル値を書き込むだけで十分です。

16 arcgis-10.0  python  raster 

すべてのデータの80％が空間コンポーネントをどのように持っているかについてのこの引用は、至る所に広まっています。GISクラス、製品の説明、講演などの紹介。誰が最初に言ったのか誰にもわかりますか？または、元の記事へのリンクがいいでしょう。

16 history 

次の質問がやや馬鹿げている場合はおologiesび申し上げますが、私はこのGISのことをまったく知りません。 Pythonでgdalを使用して、投影されたgeoTiff画像をWGS84に変換しようとしています。次のようなものを使用して、投影されたGeoTiffs内のポイントを変換するプロセスの概要を示す投稿を見つけました。 from osgeo import osr, gdal # get the existing coordinate system ds = gdal.Open('path/to/file') old_cs= osr.SpatialReference() old_cs.ImportFromWkt(ds.GetProjectionRef()) # create the new coordinate system wgs84_wkt = """ GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0, AUTHORITY["EPSG","8901"]], UNIT["degree",0.01745329251994328, AUTHORITY["EPSG","9122"]], AUTHORITY["EPSG","4326"]]""" new_cs = osr.SpatialReference() new_cs .ImportFromWkt(wgs84_wkt) # create a transform object …

16 python  gdal  coordinate-system  wgs84 

ここでスタック交換する初心者。 現在、ここにある自動タイルカッターを使用して、マップのカスタムマップタイプを作成するために必要な標準の256x256タイルを作成しています。マップタイルカッタースクリプトを使用するのは非常に不格好だと感じ、マシンの4GBのRAMをすべて使い切って、Photoshopを開いて使用します。それは、画像を正しく配置するために、画像の上下に空きスペースを埋めることは非常に面倒です。 Googleマップのタイルを作成するために使用できる他のツールを知っている人はいますか？

16 google-maps  tiles