地理情報システム

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

1
WMSを使用したジオサーバーでのポイントフィーチャの消失
世界中に約6,500ポイントのシェープファイルがあり、WMSを使用してGeoserver 2.2.1で提供しようとしています。リーフレットを使用するクライアントアプリケーションにフィルタリング機能を実装するまでは、すべて問題ないようです。CQL_FILTER(空間ではなく属性フィルター)をWMSクエリに追加すると、ズームアウトすると機能が欠落していることに気付きました。ズームインすると、ときどき再表示されますが、常に表示されるとは限りません。下の画像を参照してください- 左側のズームレベルでは、アトランタは表示されません。ズームインするとそうなります。ただし、左側のズームレベルにタンパのポイントさえ表示されない場合があります。さらに3レベル縮小すると、ポイントがまったく表示されません。問題がCQL_FILTERパラメーターであるかどうかはわかりません。6,500ポイントでは、グローバルスケールでいくつかの欠落ポイントに気づくことは困難ですが、例としてここに示している特定のフィルターは、3つの機能のみにフィルターをかけます。特に目立つズームレベルによっては、1〜3個が欠けています。 高度なマップオプションのCQLフィルターを使用して、ジオサーバーの独自のレイヤープレビューでこれと同じ動作を再現できるため、クライアントアプリケーションではないことを確信しています。 この問題をデバッグするための私自身の手順については、デフォルトのポイントスタイル/ sldファイルを試して、自分のレイヤースタイルを除外しました。私が知っているすべてのキャッシュをオフにしました。プロジェクションが正しいことを再確認しました-プロジェクションとしてWGS_1984_Web_Mercator_Auxiliary_Sphereを使用してArcGIS 10でシェープファイルを作成し、ジオサーバーでレイヤーがEPSG:3857に設定されています。また、geoserver 2.2から2.2.1にアップグレードしましたが、両方で同じ問題が発生しました。ジオサーバーの空間インデックスファイル(.qix)も削除し、空間インデックスが破損しているArcで同様の問題が発生したので再作成しましたが、明らかにそれも機能しませんでした。 これは、CQLフィルターをオンにして、上に示したのと同じ領域に拡大した、Geoserverの独自のレイヤープレビューからのスナップショットです。赤い円は、ほぼ別の点(アトランタ)が見える場所です。 考えられる他のすべての設定を微調整してみましたが、うまくいきません。また、ジオサーバーログを調べて詳細ログを有効にしましたが、エラーや例外は表示されません。また、Chromeの開発者ツールに悪いリクエストはありません。 重要な情報が不足している場合は、できる限り提供しますが、これは内部/非公開アプリケーションです。

3
米国の場所の過去24時間の合計降雨量のソースを知っていますか?
特定の場所(米国マサチューセッツ州ボストン)の降雨イベントの規模に関する履歴情報を取得したいと考えています。たとえば、特定の場所の降水量ごとに降雨イベントのヒストグラムを作成します。これまでのところ、「過去3日間」または「先月」タイプのデータセットのみがNOAAなどから見つかりました。できれば、何年分分のデータが欲しいです。 表形式のデータ(データベース、Excel、CSVなど)ほど、シェープファイルを探していません。データは任意の緯度/経度で利用できる必要はありませんが、空港などの近くのPOIで利用できます。

2
ポリゴン内にあるポイントでデータベース行を更新する
2つのテーブルを含むPostGIS / Postgresql DBがあります。1つはポイントジオメトリを持ち、もう1つは国の境界をポリゴンとして表します。各ポイントが交差する国名を、ポイントテーブルの各行に追加したいと思います。多分1つの大きな更新クエリとして。これは単純なSQLを使用して実行できると思いますが、どこから始めればよいかわかりません。これについてのアドバイスは大歓迎です...


1
ラスターフィーチャをポリゴンにスムージングする
都市化されたエリアを表すラスターデータがあります。これをトレースして、ポリゴンに変換します。 ArcMapのRaster to Polygonツールを使用して、それに近いものを取得できます。 単純化: 簡略化: ただし、コーナーを共有するポリゴンのチェーンを作成する多くの対角線(道路)があることがわかります。 これにより、ポリゴン数が大幅に増加します。これらのラスターデータセットは州レベルで処理されています。これらの線を滑らかにしてポリゴン数を減らすためのArcまたはGeoToolsツールはありますか?結果の精度がいくらか失われても問題ありません。最後に、単純化されたポリゴンで得られるものに似たものを探しています。隣接するポリゴンが結合されているだけです。

3
DEM上にベクトルを適切にドレープする方法は?
調査地域のDEMとその中の小川があります。両方をArcSceneに追加し、DEMの上にベクターをドレープすると、一部のベクターが表面から消え、表面の下に表示されます。(下の画像を参照してください)。 DEM上にベクトルを適切にドレープするための特定の設定はありますか?DEMを介してベクトルをドレープする必要がある場合、(視覚的に)より良い結果を生成する他のオプションまたはソフトウェアはありますか?

2
スキャンされたマップは、元のソースの投影で地理参照する必要がありますか?
スキャンされたインド亜大陸の歴史的な紙の地図があります。地図はアルバースの等面積円錐図法であり、その上に緯度と経度のグリッドがあります。 ソースの投影で地理参照を試みましたが、この投影に必要なすべてのパラメーターがありません。地図には2つの標準緯線のみが記載されています。中央子午線または起源の緯度、またはその問題のデータムは知りません。 このような状況で何をすべきですか?これは4000km×4000kmの地図なので、高精度な出力は得られないと思います。私はそれをwgs84地理座標系で地理参照することを考えていました。 これは良い考えですか?それともあなたが提案する他の何かがありますか?

1
層別OpenLayersクラスタリング
ポイントクラスタリングを備えたOpenLayersマップがありますが、マップで表されるデータにさらに細かい粒度を適用したいと考えています。具体的には、機能の属性に従ってクラスターを層別化したいと思います(これを「カテゴリ」と呼びます)。 つまり、5つのカテゴリがあるとします。事実上、同じカテゴリの他のポイントとのみポイントをクラスタリングする必要があります。これは5つの異なるレイヤーで達成できると思いますが、カテゴリーの数が増えるとすぐに面倒になります(現時点では50以上あります)。 OpenLayersでこのようなものを見たことがありますか?

2
ArcMapの内部で実行されるPythonスクリプトと外部で実行されるPythonスクリプト
私は仕事のためにPythonスクリプトを使い始めたばかりです。 現在、プロセスを自動化するスクリプトを作成しています。 基本的に、ユーザーにクライアント名を尋ね、可能な場合はプロジェクションを取得し、C:ドライブにクライアント用のディレクトリを作成し、クライアント固有のファイルジオデータベースを作成し、必要なデータセットを作成し、クライアントデータ固有のフィーチャクラスを作成します。最終的には、各フィーチャクラスに必要なフィールドを追加し、おそらく他の項目も追加します。 私は、ArcMapのPythonスクリプトの適切なエチケットを本当に知らずにこれを始めました。しかし、これまでに作成したものは、ArcMapの外でしか実行されません。 これは受け入れられますか? たった今見つけたarcpy.getparamaterastext()を介してユーザー入力を取得する代わりに、raw_input()を使用しています。 よろしいですか? これは機能しますが、これがスクリプトを作成する適切な方法であるかどうかはわかりません。 ここに私がこれまでに持っているコードがあります。 import sys import arcpy import os #Records name of the client client = raw_input("Enter the name of the client: (letters and underscores only) \n") #Records filepath of client to be created clientpath = "C:/" + client #Inquires if projection file exists …
10 arcpy 

2
タイルごとに繰り返されるポリゴンラベル
スタンドアロンのGeoWebCacheを使用して、リモートのGeoServerからタイルを提供しています。私の問題は、ポリゴンの重心に1回だけではなく、ポリゴンラベルがタイルのそれぞれに追加されることです。 この問題を説明する投稿を見つけました:http : //osgeo-org.1560.n6.nabble.com/polygon-label-repeated-for-each-tile-td4995203.html 最初の返答は可能な解決策に言及しました: 「全体として、私はとにかく、GeoServerから小さな画像を要求する代わりに、タイルレンダリングエンジン(GeoWebCache、MapProxy、TileCache)を使用して、タイルレンダリングエンジンにタイルスライスを実行させることをお勧めします。GeoServer(複数の小さい画像の代わりに1つの大きい画像)、これにより、全体的なタイルキャッシュの作成時間が短縮されます。」 問題は、GeoWebCacheのドキュメントを参照してもその方法を見つけることができず、上記の投稿ではそれを実装する方法が説明されていないことです。 また、GeoWebCacheの「Tiled」ドキュメントにリンクする回答を含む投稿を見つけましたが、私のコードはすでに必要なすべての属性を使用していますが、ラベルは何度も表示されます。 var Layer_1874 = new OpenLayers.Layer.WMS( 'Grundkort', '/wms10.ashx' , { format: 'image/png', srs: 'EPSG:25832', layers: 'ballerupkommune_grundkort_bk', tiled: true, tilesOrigin: '698804,6173460' } , { displayInLayerSwitcher: true, isBaseLayer: true, transitionEffect: 'resize', displayOutsideMaxExtent: true, visibility: false } ); 誰かがアイデアを持っていますか?

2
pyshpを使用して.csvファイルを.shpに変換しますか?
Pythonでcsvモジュールを使用して、pythonスクリプトと同じフォルダーにあるcsvファイルを開き、シェープファイルモジュールpyshpを使用してシェープファイルを作成する方法を理解しようとしています。 csvファイルは次のようになりますが、数千行のレコードを持つことができます。 id_nr;date;target;start_lat;start_lon 1;2012-05-21;navpoint 25x;55.123654;13.456954 1;2012-05-23;navpoint 11f;55.143654;12.456954
10 python  shapefile  csv  pyshp 

1
Postgis 2.0-public.geometry_columnsビューにビューの正しい情報が含まれていない
最近、データベースをPostgis 2.0に更新しました。 私たちのシステムでは、geometry列を含むテーブルのビューを作成するために使用していました。ビューの作成後、geometry_columnsテーブルを必要な情報で手動で更新しました。ビューは、Geoserverレイヤーの作成とogr2ogrのデータのフィードに使用されます。 Postgis 2.0 の問題は、geometry_columnsテーブルが、システムカタログからデータを取得するビューに置き換えられたことです。私たちの場合、作成されたビューでは、「coord_dimension」、「srid」、および「type」列の正しい情報を取得していません。 情報:いくつかの調査の結果、geometry_columnsビューがpg_constraintテーブルを使用してこれらの列を埋めていることがわかりましたが、そのビューには制約が定義されていません。 誰かがこの問題をどのように修正できるか考えていますか?gemetry_columnsビューに必要な情報を取得させるにはどうすればよいですか?これを行うには回避策または他の方法がありますか?

2
Nominatimを使用して通りの交差点をジオコーディングすることは可能ですか?
Nominatimを使用して道路交差点の位置をジオコーディングすることは可能ですか? ジオコードクエリの例 Commonwealth Ave at Washington Street, Boston, MA Googleによる成功したジオコード https://maps.google.com/maps?q=Commonwealth+Avenue+at+Washington+street,+Boston,+MA&hl=en&sll=42.351351,-71.118573&sspn=0.054868,0.132093&t=h&hnear=Commonwealth+Ave+%26+Washington+St,+Boston,+Suffolk,+Massachusetts+02135&z=17 Nominatimによる失敗したジオコード http://nominatim.openstreetmap.org/search.php?q=Commonwealth+Ave+at+Washington+Street%2C+Boston%2C+MA&viewbox=-337.5%2C85.34%2C337.5%2C-80.5

2
ArcGIS Desktopを使用してSQLiteファイルを処理しますか?
仕事にはqgisとgrassを使用しています。100mbsを超えるシェイプファイルを作成しました(統計分析の結果です)。 それをsqlite形式(空間的に有効)に変換すると、データは15 mbsに縮小されます。メールやドロップボックスと共有することは許容されます。 しかし、私の同僚はarcgis 9.3を使用しています-彼は問題なくデータを読み取ることができますか?

5
OSM道路を交差点で個々のセグメントに分割する方法は?
OpenStreetMapデータを使用してpgRoutingで使用する道路ネットワークを作成したいと思います。GeoFabrikからシェープファイルをPostgresテーブルにロードしました(PostGISを有効にした状態)。しかし、私が抱えていた問題の1つは、道路が常に交差点で終わらないということでした。そのため、すべての交差点または交差点で道路をすべて分割することにしました。 道路が交差または交差するすべての交差点を特定するために、私は次のものを使用しましたSQL(前の質問と同様)。 CREATE TABLE split_points as SELECT DISTINCT ST_GeometryN(ST_Intersection(a.geom, b.geom),1) as geom FROM roads as a, roads as b WHERE ST_Touches(a.geom, b.geom) OR ST_Crosses(a.geom, b.geom) AND a.gid != b.gid GROUP BY ST_Intersection(a.geom, b.geom); これらのポイントを使用して道路を分割したいと思います。私は次のアプローチを使用しました: CREATE TABLE split_roads as SELECT ST_GeomFromEWKB((ST_Dump(ST_Split(g.geom, blade.geom))).geom) As geom, generate_series(1,ST_NumGeometries((ST_Split(g.geom, blade.geom)))) as gid FROM split_points as …

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