地理情報システム

私の質問は、ポイントの位置にあるラスターから値を抽出することです。関数抽出を使用すると、これは非常に簡単です。関数は、ポイント内のすべての変数の値を持つデータフレームを提供します。そのデータフレームに各ポイントの座標が必要です。どうすればそれを実現できますか？Rに、ラスタから値を抽出するときに、ロケーションポイントの列も追加すると言うことはできますか。 これは私の手順です： presencias=read.table("c:/SDM_R/presencias/P_lentiscus_pres.csv",header=TRUE,sep=";") lista_variables <-list.files(path="Variables_modelizacion/solo_ascii",pattern='*.asc',full.names=TRUE) variables <- stack(lista_variables) variables_presencia<-extract(variables, presencias) 結果は次のようになります。 > bio1 bio12 bio18 bio2 bio4 > 90 875 165 95 4886 > 115 1085 158 83 4075 > 135 1153 153 67 3402 > 85 1026 137 99 5203 > 96 667 128 108 5823 > 98 531 …

13 raster  r 

QGIS 1.8には、MSAccessデータベースからデータを読み取り、一連のメモリレイヤーに追加する作業プラグインがありました。暫定にはいくつかの処理が含まれているため、単にQGISを使用してデータベースから直接読み取ることは選択肢とは思いません。 QGIS 1.8から1.9に移行したいと思います（主に印刷コンポーザの品質が向上したため）。プラグインは新しいAPIでは機能しません。 Google検索で出てきたさまざまな方法を試しました。1つは、http： //www.qgis.org/pyqgis-cookbook/vector.html#memory-provider から以下のコードを変更します。つまり、ジオメトリと属性をデータプロバイダーに追加し、レイヤーを更新して、新しいAPIが機能するようにします少しですが、編集モードに入るまで属性は表示されませんでした（http://hub.qgis.org/issues/3713に似ています）。上記のリンクの返信＃1で詳しく説明されている別のアプローチでは、レイヤーと属性を正しく追加しましたが、レイヤーに機能を追加できませんでした。 それは非常に単純なタスクであるべきであるため、ここにいる誰かがこれがどのように行われるべきかの実例を提供できることを期待しています。（PS私はプロのプログラマーではなく、私のコーディングの大部分はかなり粗雑です-どんなガイダンスも歓迎しますが、私の側で無知を許してください） # Receivers = a list of lists returned from a database query # create layer vl = QgsVectorLayer("Point", item, "memory") pr = vl.dataProvider() # add fields pr.addAttributes( [ QgsField("Rec_No", QVariant.Int), QgsField("Include", QVariant.String), QgsField("Label", QVariant.String), QgsField("X", QVariant.Double), QgsField("Y", QVariant.Double), QgsField("Z", QVariant.Double), QgsField("Height", QVariant.Double), …

13 qgis  vector  qgis-plugins  pyqgis  qgis-1 

PythonおよびArcPy（非ArcPyモジュールなし）を使用したジオプロセシング操作中に、ArcGISレイヤーの単一のフィーチャ（または潜在的に多くのフィーチャ）をWell Known Text（WKT）形式に変換できるようにしたいと思います。その目的は、WKTをSQL Server Spatialに渡し、ArcGIS GPツールセットの外部で追加の処理を行うことです。ArcGISジオメトリフィーチャをArcPy経由でWKTに変換する方法はありますか？ 私は自分が何を求めているのかを見つけることなく、すでに以下を読んでいます： シェープファイルからテキスト（非推奨の「サンプル」ツールボックスのツールを参照） ESRIジオメトリとWKT間の変換（ArcObjectsまたは非ArcPyモジュールを使用） シェープファイルをWKTに変換する方法は？（QGISとOGRを使用） どうやら「と呼ばれるツールがありました書き込みは、テキストファイルに機能サンプルツールボックスにあった（Pythonのスクリプトのように見える）」が、そのツールボックスがされたバージョン10で非推奨と私はサンプルツールボックス（のコピーを見つけることができませんマシンにv10.0）があります。現在の解決策がない場合、誰かが私にサンプルのコピーを指し示すことができれば、そのツールを使用してもまったく問題ありません。

13 arcpy  convert  wkt 

PostGISデータベースにGISデータがあり、密度マップを生成してリーフレット/ Googleマップに表示したいと思います。 質問：マップの視点内のデータポイントをすべて（数百万）リーフレット/ Gmapsマップに転送せずに、これを行うための推奨/典型的な方法は何ですか？CLIなどを介して自動化が可能でない限り、ユーザーの介入（つまり、マップを生成するためのArcGISへのデータのインポート）は必要ありません。 現在、Postgresql / PostGISを空間データベースとして使用し、ブラウザでリーフレットを使用しています。

13 postgis  google-maps  leaflet 

テンプレートジオデータベースがあり、それを地域のデータクリエーターに送信します。これらのジオデータベースを更新し、事前定義されたフィーチャクラスの地域のすべてのデータをデジタル化します。 ファイルジオデータベースを取得したので、多くのフィーチャクラスが変更されていることがよくあります。フィールドが追加または削除される場合があります。場合によっては、フィーチャクラス自体が削除されるか、新しいフィーチャクラスが作成されます。 どのフィーチャクラスとテーブルが変更されたか、および共通のテーブルとフィーチャクラスの異なるフィールドを示すレポートを取得したいと思います。 私はすでに次のようないくつかの質問を見てきました。 地理空間データのバージョン管理されていないコピーを比較するツール SDEデータベースの古いバージョンと新しいバージョンの差分を取得する良い方法はありますか？ しかし、これらの質問で与えられた答えは役に立たなかった。

13 arcgis-10.0  esri-geodatabase  schema-architecture 

私の組織では、ユーザーがデータを表示してmxdsを構築する手段として* .lyrファイルを一貫して使用しています。 私たちの目的にはうまく機能しますが、.lyrファイルが何をするのか、なぜ使用されるのかをユーザーが実際に理解していないと思われる問題に遭遇しました。定義/用語。 .lyrファイルとは何ですか？ユーザーがGISの軽度/中程度の知識を持っているとしますか？

13 arcgis-desktop  layer-file  terminology 

LeafletやOpenLayersのようなオープンソースのJavaScriptライブラリがあります。これらのライブラリは、Webアプリケーションやモバイルアプリケーションで使用でき、画像の代わりにベクトルベースのマップを読み込みますか？道路情報といくつかのPOIだけが必要です。 私が見つけた唯一のライブラリはポリマップでしたが、それでもベクターSVGデータ（OpenStreetMapから推奨）を提供するサーバーまたはサービスが必要になります。また、Mapsforgeと呼ばれる素敵なAndroidプロジェクトもありますが、Webアプリケーションや他のモバイルデバイスに必要です。

13 openstreetmap  web-mapping  vector  web-service  vector-tiles 

Tilemillをインストールしてマップをmbtilesとしてエクスポートし、mbutilを使用してPNG画像とUTFGridとしてエクスポートしました。など、MSS（Carto）ファイルを含むTilemillプロジェクトがあります。データソースは毎週変更されます。このアクションを自動化します。 それで、Tilemillはコマンドラインスクリプトまたは特定のプロパティでプロジェクトをエクスポートする他の何かをサポートしていますか？ 現在のシステム：Windows 7 x64 @Sunilに応じて編集： TileMillには、Cartoスタイルシートを自動的に操作する機能はありませんが、bashスクリプトを使用してCartoファイルを直接操作できます。TileMillには、bashから呼び出すことができるエクスポートコマンドもあります。 Mac OS Xでは、エクスポートコマンドの使用の詳細については、ターミナルで次のようなコマンドを実行できます。 /Applications/TileMill.app/Contents/Resources/node /Applications/TileMill.app/Contents/Resources/index.js export --help 正確なコマンドは、TileMillのインストール方法によって異なります 残念ながら、Mac OS Xの同じコマンドは現在のシステムでは機能しません。最後の文はどういう意味ですか、Tilemillの異なるインストールはありますか？

13 tilemill 

私が見る歩行者のルーティングを提供するアプリケーションのほとんどは、多くのニーズには十分ですが、最適ではありません。街の広場のようなオープンスペースを横断することはできません、ルーティングは通りの真ん中にあります。 歩行者専用のルーティングを備えたWebサイトまたはアプリケーションの例はありますか？

13 algorithm  routing  transportation 

QGISとそのプラグインでバグをよく見つけます。 インターネットでこれについて尋ねたいのですが、この問題に対処する適切な場所を見つけられないようです。 指示はありますか？

13 qgis  qgis-plugins 

コロプレスマップのプロットに関して質問があります。カスタムブレークを使用してマップをプロットする場合の推奨事項を教えてください。そして、それを少し魅力的に見せてください。 私のスクリプトでは、plotおよびspplotコマンドを使用して一部のデータを視覚化しましたが、結果に完全には満足していません。理由は簡単にわかります。spplotを使用して、異なる値の色を繰り返すカラーランプを取得しますが、これは少し奇妙です。カスタムブレークでプロットを使用すると、たとえばゼロ値がマップ上で白になるように、値の間隔を改善できませんでした。凡例とマップの間の歪んだ比率は、おそらくRstudioの結果だと思います。 だから私は誰かがカスタムブレークを使用してRでコロプレスマップをプロットするためのより良い提案があるのだろうか？ 編集20-02-2013 示唆されたように、コロプレスマップを作成するためにggplot2を使用しました。以下のコードを使用しました。2つの問題が発生しました。1つは、データを強化した後、シェープファイルのすべてのデータがデータフレームに転送されるわけではないことです。したがって、データをマージできるように「id」変数の名前を変更する必要があります。軽微な問題。 より大きな問題は、データフレーム内のデータが正しいという事実にもかかわらず、ggplotはデータが正しくないマップを表示することです。凡例に応じた色の値は、データ内の変数の値と一致しません。ここで何かを見落としていますか？ # Plotting polygon shapefiles (try-out) # Load the packages require(rgdal) require(maptools) require(ggplot2) gpclibPermit() # required for the fortify method # Load the shapefile africa=readOGR("/home/GIS",layer="africaII") names(africa) # Load the data on conflict conflict<-read.csv("africa_conflict.csv", header=TRUE) names(conflict) # Merge the data together africadat<-merge(africa, conflict, by="ISO3") # Changing …

13 r  choropleth 

以下の画像は、4つのデータフレームを持つArcGIS 10.1のレイアウトを示しています。各データフレームは同じ地理的領域ですが、画像は異なります。たとえば、左上の画像は1989 DOQ、右上の画像は1989 DOQの反転、左下の画像は2012 DOQ、右下の画像は2012 DOQの反転です。 データフレームをリンクして、あるデータフレームで画像をパンすると、他のすべてのデータフレームで画像がパンするようにすることは可能ですか（つまり、Erdas Imagineで「ビューをリンク」する方法と同様） これはQGISで可能ですか？

13 qgis  arcgis-desktop  arcgis-10.1 

私はこれが完全なnoobの質問であることを知っていますが、それからPostgreSQLに関して言えば、私は完全なnoobです... Windows 7、64ビットを実行しているコンピューターにOpenGeoスイートをインストールしました。チュートリアルデータをダウンロードし、作業を進めています。ArcGISで作成したデータベースにアクセスして、それがどのように機能するかを学び始めたいと思います。ただし、「データベース接続の追加...」を使用してデータベースに接続しようとすると、「インスタンス」に何を入力するかわからないようです-私が試したすべて（localhost、localhost、54321、ポートがある場合とない場合のIPアドレス（私のコンピューター名）は機能しないようです。 OpenGeoチュートリアルページのすべての指示に従ってnycテストデータベースをインストールして作成しましたが、ここでは完全に失われています。これは簡単に対処できなければならないことはわかっていますが、頭を悩ませることはできません。もっと経験のある人がこれに骨を投げつけることはできますか？ イライラするほど、QGISはほとんど労力をかけずにデータベースに接続し、すべてが正常に動作しますが、私の組織は非常にESRIに重点を置いており、QGISへの切り替えは長いショットです... 編集 この件についてご意見をお寄せいただきありがとうございます。私はこのアドバイスに本当に感謝しています。それがこのサイトを素晴らしいものにしているのです。インストールしました 拡張なしのOpenGeo Suite 3.0.1 ArcGIS Desktop 10.1、サーバーのライセンスなし Opengeo Webサイトにアクセスして、Windowsインストーラーをダウンロードし、デフォルト構成をインストールしました。コロンとコンマ、マシン名、localhost、postgresql、サーバー、IPアドレス、ランダムな呪いの言葉、考えられるものを使用して、考えられるインスタンス名のすべての順列を試しました。OpenGeoワークショップが、使用する「nyc」練習データベースを設定するよう指示したため、5432ではなくポート54321を使用していました。 localhostが動作するようです。pgAdminのサーバープロパティのショットを次に示します。

13 postgresql  arcgis-10.1 

arcgisフォーラムに投稿しましたが、回答がありませんでした。誰もこれを解決する方法を教えてもらえますか？ 10月1日から12月31日までの風速情報を含む車両追跡記録があるとします。すべてのデータはジオデータベース（sde-sqlserver）に保存され、フィーチャレイヤーとしてArcGIS Server 10.1に保存されます。タイムスライダーは、車の位置を風速で表示します。 ユーザーが時間範囲（たとえば、10月2日から10月4日）を変更すると、最初のクエリタスク（カウント用）が範囲内のフィーチャの数を計算します。通常、2日（例：1750）でも1000を超える結果があります（ただし、この制限を変更したくありません）。 別のクエリタスク（executeforIds）を使用して、すべてのレコードを手元に保持しましたが、風速の全体的な傾向を示すグラフを作成するのに十分な数であるモジュラス（1/10）で量を減らしました。ただし、csvでデータセット全体（この場合は1750行）をダウンロードするオプションも提供したい ここでは、findtaskを使用して、時間範囲内の属性データセットを取得しました。 try { //console.log(app.objIDs); var attribs; csvdata = ""; var find = new esri.tasks.FindTask("http://xyz.com/arcgis/rest/services/Realtime/Car_Wind_WM/MapServer"); var params = new esri.tasks.FindParameters(); params.layerIds = [0]; params.searchFields = ["OBJECTID"]; for (var i = 0; i < app.objIDs.length; i++) { params.searchText = app.objIDs[i]; find.execute(params, function (results) { attribs = …

13 arcgis-10.1  arcgis-server  arcgis-javascript-api  query 

QGISで小数点以下2桁までの行の長さを計算する方法を試しています。フィールド計算機を使用して新しいフィールドを作成する場合、フィールドタイプ、フィールド幅、精度にどのような変更を加えても、長さは常に小数点以下13桁になります。

13 qgis