地理情報システム

編集/更新：ArcGIS Flex API、SQL Server Express、ArcGIS Server 10.1 を使用してフィーチャを編集してオンラインマップを作成するためのステップバイステップガイドを取得することを期待して、これに賞金を追加します。私は非常に多くの問題（Expressによるデータベース認証、機能の編集を可能にするエラー）に遭遇したので、最善の策は最初からやり直して各ステップを実行することだと思います。ある程度詳細なガイドをいただければ幸いです。 私は主にデスクトップGISプログラマーであり、複数のユーザーが機能を編集できるWebアプリケーションを開発する必要があります。ゼロから始める場合、このプロセスを実行するために必要な完全な手順の簡単な概要を探しています。 ソフトウェア：ArcGIS Flex APIおよびArcGIS Server 10.1を使用したSQL Server Express 私が今のところ理解している手順： 1）ArcGIS Serverをサーバー（実際には既存のサーバー上の仮想マシン）にインストールします 2）ArcGIS Web Adaptorのインストール 3）同じサーバーにArcSDEをインストールする 4）SQL Server Expressを同じサーバーにインストールする ここでいくつかの説明が必要です。理解したように、表示/編集したいフィーチャを格納するために新しいジオデータベースを作成する必要があります。また、フィーチャサービスとジオメトリサービスを作成する必要があることも理解しています。最後に、Webサーバーを使用して、開発したアプリケーションをホストする必要があります。 これまでに実行したことがないため、基本的にこのプロセスのウォークスルーを探しています。どんな情報/リソースも大歓迎です。

12 arcgis-server  web-mapping  enterprise-geodatabase 

バスルートを表す複数のラインレイヤーがあり、ルートの多くが重複しています（つまり、ターミナルの近く、転送場所など）。Qgisには、両方の機能を互いに/並列に表示する機能がありますか？例のリンクでは、赤と紫の線は2つの別々のルートであり、どちらも同じ道路を利用しています。紫色の線は赤色の線と重なり、赤色の線は見えません。それが不可能であるという悪いニュースでさえ、私はどんなフィードバックでも感謝します。

12 qgis  overlapping-features  display  parallel-lines 

次の例（http://en.wikipedia.org/wiki/Well-known_textから）で、シェープファイルを単純なWell Known Text（WKT）形式に変換しようとしています。 MULTIPOLYGON (((40 40, 20 45, 45 30, 40 40)), ((20 35, 45 20, 30 5, 10 10, 10 30, 20 35), (30 20, 20 25, 20 15, 30 20))) QGISを使用してシェープファイルをCSVとして保存し、ogr2ogrを使用してあらゆる変換を実行できることを知っています。実際のogr2ogr -f CSV out.wkt source.shp -lco GEOMETRY=AS_WKTところ、私はこれまで来たのと同じように私を近づけますが、そこまでは行きません。助言がありますか？ psマッパーズが述べたように、このスレッドはこれと非常によく似ています。このスレッドは、前述のogr2ogrアプローチを提供しましたが、この特定の課題を解決しませんでした。私が取得しているogr2ogrの出力は次のようになっています（行を切り捨てています）。どうやら、ogr2ogrでマルチポリゴンを使用する方法がわかりません。 WKT,AREA,PERIMETER,PINUPOND_,PINUPOND_I,CODE "POLYGON ((-120.630531650950061 50.838562484449184, . . . "POLYGON ((-123.206067372148397 51.038984074378327, …

12 shapefile  convert  wkt 

少し前に、属性テーブルをPython辞書に変換するための簡単なPython関数を作成しました。キーはユーザー指定の一意のIDフィールド（通常はOIDフィールド）から取得されます。さらに、デフォルトではすべてのフィールドが辞書にコピーされますが、サブセットのみを指定できるようにするパラメーターが含まれています。 def make_attribute_dict(fc, key_field, attr_list=['*']): dict = {} fc_field_objects = arcpy.ListFields(fc) fc_fields = [field.name for field in fc_field_objects if field.type != 'Geometry'] if attr_list == ['*']: valid_fields = fc_fields else: valid_fields = [field for field in attr_list if field in fc_fields] if key_field not in valid_fields: cursor_fields = valid_fields + …

12 arcpy  arcgis-10.1  attribute-table  cursor  dictionary 

シェープファイルをRにインポートし、テーブルに結合しました。シェープファイルには国勢調査IDがすべて含まれていますが、テーブルには選択した国勢調査IDのみが含まれています。私は今、すべての行を削除しようとしましたが、一致しませんでした。 これは私のデータセットのようです（NAですべての行を削除しようとしているので、最後の2つを削除する必要があります）。 CTUID Cluster Average 5350007.01 1 124.53 5350007.02 1 234.87 5350010.01 4 110.11 5350010.02 5 187.68 5350001 NA NA 5350002 NA NA 私はこの行コードを使用してみました： shape2[!(rowSums(is.na(shape2))==NCOL(shape2)),] これは私にこのエラーを与えました： Error in rowSums(is.na(shape2)) : 'x' must be an array of at least two dimensions In addition: Warning message: In is.na(shape2) : is.na() applied to …

12 shapefile  r  data-management 

タイプポイントの単一レイヤーを持つシェープファイルがあります。 Google CRS（EPSG：900913）にあります。 [フィーチャの識別]ツールを使用すると、クリックされたポイントに関する多くの情報が表示されますが、座標が表示されません。ポイントの緯度/経度をインタラクティブに調べる方法はありますか？

12 qgis 

ちょうど私が正しい軌道に乗っていることを確認するために： 球上のすべての大円と正距円筒図法（緯度、経度のペア）は次のいずれかです。 経線（つまり、極から極へ） フォームの tan latitude = sin360(longitude + rotation) * amplitude + offset （オフセット/振幅の組み合わせに追加の制限があります-明らかに、振幅が0の大円経路もすべてオフセット0-赤道です）。 または、このスキームに適合しない大円のパスがありますか（再度、他の地図投影ではなく、経度と緯度の座標系でのみ）。 注：tan質問を投稿した後、whubersの優れた返信に対する反応として、上記を追加しました。これは、ことが判明しoffset、その後は常に0です。

12 coordinate-system  latitude-longitude  spherical-geometry  great-circle 

この質問の答えはコミュニティの努力です。この投稿を改善するには、既存の回答を編集してください。現在、新しい回答やインタラクションを受け付けていません。 最近、GIS Analyst / DeveloperがArcGISのジオプロセシングとArcPyサイトパッケージでPythonを使用して自動化しようとしている操作の種類について質問がありました。作成するものが業界に関連し、後で職場で再利用される可能性があることを確認するために、いくつかの演習を検索するときに知っておくと役立つと思います。 最も簡単な答えは「Esriヘルプを読んでサンプルを調べる」ことでしょうが、実装するのに最も一般的なより具体的なシナリオを探していました。したがって、共有ワークフローは、「シェープファイルを含む.zipファイルを取得し、Pythonを使用して解凍し、それらをすべてX座標系に投影し、ArcSDEジオデータベースにロードし、これらのデータへのアクセスをユーザーに許可する」ことを歓迎します。ワークフローを簡潔に説明してください。極端な詳細は必要ありません。

12 arcpy  references 

QGIS 1.7.4-Wroclawを使用し、WGS84でSHPレイヤーを使用しています。ftoolsバッファーツールを使用して、このレイヤーから特定のキロメートル数までバッファーアウトします。 バッファツールは常にレイヤー単位を使用することを理解しています。これはWGS84の場合は10進度です。度がメートルに一貫して変換されないことは誰もが知っているので、バッファを作成するにはどうすればよいですか？ シェープファイルを、メーターをネイティブで使用する別のCRSに変換する必要がありますか？ その場合、どのように選択すればよいですか？

12 qgis  buffer  units 

ファイルジオデータベース（fgdb）には、多数のファイルジオデータベーステーブルが含まれています。私の知る限り、それらはdbfファイルとして存在しますが、内にありDatabase.gdbます。 ArcCatalogでは、ファイルパスウェイはに似ていC:\Users\...\Database.gdb\Stats_AAます。 これらすべてのdbfファイルをR（統計ソフトウェア）に読み込む方法は？供給する正しい経路は何ですか？使用される関数はread.dbf（外部パッケージ内）です。 のバリアント test<-read.dbf(file="C:/Users..Database.gdb/Stats_AA.dbf") そして test<-read.dbf(file="C:/Users..Database/Stats_AA.dbf") 動作しません。使用するファイル名の正しい「フォーム」とは何ですか、またはすべてのファイルジオデータベーステーブルを他のフォームまたは場所にエクスポートする必要がありますか？

12 r  file-geodatabase 

私はGeoserverと協力して、米国の下位48郡をオープンレイヤーに提供しています（3109ポリゴン-さらに多くの頂点）。郡はpostgisデータベースにロードされます。その量の頂点をクライアントにプッシュしようとするときの開発者の経験に興味があります。 どのWFS形式で最高の結果を達成しましたか？Geoserverへの追加のチューニングが使用されましたか？ タイル化されたWMSの方が高速であることを認識していますが、openLayersを使用してコロプレスマップの動的な変更を許可したいのです。ユーザーがフォームを送信し、Pythonスクリプトが呼び出され、マップdivを再ロードするために新しいレイヤーがopenlayersに返されます。また、オープンレイヤーのポリゴンの複雑さを軽減する前に、これをフル解像度の形式で試してみたいと思います。

12 geoserver  kml  wfs  geojson  gml 

QGIS 1.8のマップウィンドウと同じように、Print Composerでマップをレイヤー/オブジェクト/マップの中心にどのように中央揃えできるかわかりません。地図の中央に配置したいレイヤーがあるとしましょう。マップウィンドウでは、非常に簡単です。[ レイヤー範囲にズーム ]をクリックするだけです。 残念ながら、Map Composerでマップフレームを作成すると、レイヤーは幅の中央に配置されますが（良好）、高さの中央には配置されません（下に揃えられていると思います）。 [ 設定 ]ボタンを使用してPrint Composer /アイテムプロパティ/エクステントのキャンバス範囲をマップすると、再び幅の中央に配置されますが、今回は上に揃えられます。 どちらの場合でも、中央に配置する方がはるかに良いと思います。マップをPrint Composerの中央に配置し、マップウィンドウのマップと同じ中央に配置する方法はありますか？

12 qgis  printing  layouts  print-composer 

リアルタイムのマップマッチングを必要とするプロジェクトに取り組んでいます。現在ArcGISを使用していますが、動作するかどうかはわかりません。 これが私のシステムで起こることです： 携帯電話はサーバーに5秒ごとにGPS修正を送信し、マップのリンクの1つにポイントをスナップして一致したポイントを保存するには、リアルタイムのマップマッチングを行う必要があります。 ArcGISには、ポイントを最も近いリンクにスナップするための「ニア」ツールがあります。ただし、ツールへの入力は、ポイントのシェープファイルである必要があります。 ArcGISで取得したリアルタイムデータをどのように処理できるのでしょうか。各ポイントの形状ファイルを作成しますか？ この作業を達成するための他の提案はありますか？

12 arcgis-desktop  gps  snapping  real-time 

Windows 7 64ビットにインストールされたNumpy 1.3およびGdal 1.9.1でPython 2.6.5（32ビット）を使用しています。800 MBのImagine（.img）ラスターデータセットをNumpy配列に読み取ってラスター代数を実行しようとしていますが、次のコードを実行するとすぐにPython.exeがクラッシュします。 from osgeo import gdal g = gdal.Open(r'path\to\dataset', gdal.GA_Readonly) b = g.GetRasterBand(1) data = b.ReadAsArray() Python.exeはb.ReadAsArray()呼び出し時にクラッシュします。Google検索をいくつか行って、Gdal 1.6からWindows 7 64bitでこの問題に言及している記事を見つけましたが、当時の最新の開発バージョンで修正されたと述べました。 他の誰かがこの問題を抱えていましたか？回避策はありますか？ 更新： PyDevでコードをデバッグして、どこでエラーが発生しているかを特定することにしました。私が伝えることができる（まだエラーメッセージが表示されない）ことから、g​​dal_array.pyの22行目で失敗しています。 _mod = imp.load_module('_gdal_array', fp, pathname, description) 上記のコード行に足を踏み入れると、numpyのinit .pyモジュールに移動します。numpyの最後に到達したとき。__ init __ .pyモジュール、上記のコード行に戻ります。次に、ステップインボタンを押すと、gdal_array.py内の次の行に移動しますが、スクリプトはエラーメッセージも何も表示されずに終了します。 更新＃2： GDAL 1.9.1をアンインストールし、Python CheeseshopおよびOSGeoのWindows BinariesからGDAL 1.6.1をインストールしました。まだ同じ問題がありました。

12 gdal  numpy  python-2.6  array 

私の最初の修士研究テーマは空間インデックスについてでした。RツリーとR +ツリー、R *ツリー、その他の非常によく似たトピックについて学びました。これは公開された情報ではありませんが、ESRIはどの空間インデックスタイプを使用しますか？二次的な質問として、このタイプはどのように変更できますか？異なるタイプのインデックスがパフォーマンスを向上させる状況がいくつかあります。 Oracleの場合、Rツリーまたはクワッドツリーのいずれか、または両方にインデックスを付けることができるST_Geometry空間データタイプを使用できます。http：//docs.oracle.com/html/A88805_01/sdo_intr.htm 他のデータベースは、デフォルトでは非空間であり、利用可能な空間インデックスを使用していないように見えるバイナリストレージメカニズムを使用します。SQL 2005以前では、空間データ型は存在しませんでしたが、ESRIはSQLデータベースを使用できました。このケーススタディSQL 2012：http : //www.microsoft.com/casestudies/Case_Study_Detail.aspx?CaseStudyID=710000000547によると、自動化された空間パフォーマンスタイプを使用しますが、R-Treeの種類はリストしません。 この質問に答える資料を読むことに対するボーナスポイント。 ありがとう。

12 arcgis-server