地理情報システム

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

4
R spplotチュートリアル
R spplotのspパッケージの関数を適切に使用する方法を扱ったチュートリアル/本/その他のトレーニング資料が必要です。 これらの機能の多くと同様に、それは非常に強力なように見えますが、何をしたいかを正確に知らない限り、ドキュメントはそれほど役に立ちません。さまざまな例のページを見つけましたが、それらはまだそれほど明確ではありません。 spplot便利なマップを作成するために使用できるリソースがありますか（たとえば、ラスターデータ、ポイント上のポリゴンのオーバーレイ、さまざまなカラースケール、凡例、北矢印など）。 ？ 理想的には、Rのラスターデータと地図の作成に焦点を当てたUseR！スタイルの本が欲しいのですが、その種の唯一の空間の本は、主に地球統計学とポイントデータに焦点を当てているようです。 何か案は？
11 cartography  r 
4
CSVへの表形式データのエクスポートをサポートしているArcGISジオプロセシングツールは何ですか？
表形式の出力を許可するすべてのジオプロセシングツールではないにしても、ほとんどの場合、サポートされている出力形式のジオデータベーステーブルとDBFファイルに限定されることがわかりました。CSV出力もサポートしているものはありますか？Table to TableやCopy Rowsなど、多数のツールがCSV入力をサポートしています。 私が見た唯一の例外は、空間属性ツールボックス内のExport Feature Attribute to ASCIIというスクリプトツールですが、これはテーブルではなくフィーチャクラスのみをサポートし、非ASCII文字でチョークします。 <type 'exceptions.UnicodeEncodeError'>: 'ascii' codec can't encode character u'\xa0' in position 213: ordinal not in range(128)` 編集：これまでの回答の選択肢に感謝します。私は実際にこれらのいくつかを実際に使用していますが、明確にするために、ArcGIS 10の組み込みジオプロセシングツールに関する質問を明確にしています。 ArcMapが、どうやらGPフレームワークを介していない（Pythonスクリプトなし）。 ちなみに、ArcGIS Ideasサイトを少し調べてみたところ、いくつかの既存のアイデアに出くわしました。 テーブルにテキストオプションを追加 モデルビルダー専用ツール-CSVへのエクスポート テーブルをMS Excelファイルに直接エクスポートします（特にCSVについてではありませんが、関連する多くの議論があります）
1
QGISからのリクエストをキャプチャするためにフィドラーをセットアップすることは可能ですか？
QGISからWMS / WFSリクエストをキャプチャするセットアップフィドラーが可能かどうか、そしてその方法を知りたいです。 FiddlerはMapInfoからのリクエストをキャプチャしていますが、それを実現するためにfiddlerで設定を行ったことを覚えていません。 QGISからフィドラーキャプチャリクエストを作成できますか？ バージョン情報：Fiddler 2.3.9.5、QGIS 1.7.0
11 qgis  wms  wfs  request  fiddler 
2
ArcGIS ModelBuilderのCollect Values出力からデータをエクスポートしますか？
現在、ModelBuilderで遊んでいます。現在持っているモデルの写真を追加しました。 このモデルは現在6つのシェープファイルを反復しているため、「出力値」のそれぞれにはリストに6つの数字が含まれています。これらのリストからテーブル/テキストファイルなどに値を抽出する方法が見つかりません。 これを行う方法はありますか？ これは、モデルの実行後に「出力値」を開いたときに表示されるものです。 どういうわけかこの6つの数字を取得したいだけです...
4
QGISを使用して多基準分析を実行していますか？
「開発するのに最適なロット」という質問に答えるために、多基準分析を行う必要があります。 基準のいくつかは次のとおりです。 最寄りのバス停までの距離（バス停があるポイントレイヤー） 最寄りの店舗の距離（店舗のあるポイントレイヤー） 洪水の危険性（ポリゴンレイヤー、危険度属性が1〜4） 自然保護地域（ポリゴンレイヤー）のロット 所有者は既に自分のロットで何かを計画しています（ロットの属性に手動で入力された情報）など... QGISで試してみたいと思ったので、次のようにしました。 ロットレイヤー属性テーブルに次の列を追加します。 「analysis_BUS」 「analysis_SHOPS」 「analysis_FLOOD」 「analysis_PROJECT」 「...」 「analysis_MEAN」 「ポリゴンを重心に」を使用して、ロットレイヤーをポイントに変換します 「距離行列」ツールを実行する CSVを開いてExcelで操作を実行します（バス停のグレードは200mに近い場合は1.0、750mを超える場合は0.0ですが、QGISでMIN（）関数を見つけることはできません） 結果のCSVをQGISに戻します ショップについても同じことを繰り返します 「ポリゴン内のポイント」ツールを実行して、自然保護エリア内のすべてのポイントを選択します 選択したすべてのポイントに0.0を設定します 他の「領域内」基準について繰り返します 「空間結合」ツールを実行して、洪水危険区域情報をマージします 列計算機を使用して計算を実行し、平均評点を取得します（各基準に対して決定された係数を使用） すべて完了したら、各条件ごとにBUILDING LOTSシェープファイルを1回追加し直します 各基準について、変換されたレイヤー（重心を持つレイヤー）をLOT IDで結合します 対応する基準属性と平均成績属性に従って、表示を赤から緑へのグラデーションに設定します 今では、2日間の良い仕事を経て、構築に適している場合はすべての基準が緑色で表示され、悪い場合は赤色で表示され、すべての基準を1つの美しい赤緑マップに集約する合成があります。（そして、「shapefiles」フォルダに巨大な混乱があります） 今問題。 仮に ： 別のバスネットワークシナリオで同じ分析を試してみたいですか？ 更新されたロットシェープファイルを受け取ります（たとえば、13000ロットすべてで13の変更があります） 基準に応じて異なる重みをテストしたいですか？ 最初からやり直す必要がありますか？ 間違ったツールを正しく使用していますか、それとも間違ったツールを使用していますか？ 市販のGISソフトウェアの方が簡単でしょうか？ 回答者/コメント者の意味がわかりますが、実際にはラスターの使用を考えていませんでした。 ただし、主な問題は、すべてのプロセスを最初から再起動することなく、さまざまなシナリオを試す機能やベースデータを更新する機能に関するものでした。 ニュースのステップがあるので、あなたの提案は私が提案したものよりもはるかに柔軟ではないようです（さらに複雑かもしれません）：-（基準ごとに）ラスター化。-（最終的に）サンプリング（部分的なオーバーラップを含める場合は非常に複雑です） そのSextanteモデルビルダーは素晴らしいようです。実際、最後のコメントを投稿するとき、私はまさにそのようなことを考えていました。 Rhino3Dモデラーの優れたプラグインであり、ノードグラフワークフロー構築の同じ概念を使用するGrasshopper3Dをかなり多く使用しました（GISソフトウェアとは関係ありません）。（例：http : //designreform.net/2009/07/rhino-grasshopper-parametric-truss） これは、多くのGISデータ分析に非常によく適合しているため、このようなノードグラフツールを中心に構築されたGISソフトウェアをぜひご覧ください。 Sextante Modelerを試して、どのように機能するかをお知らせします。私はそれをグーグルで自分で見つけたかったのですが、キーワード「モデルビルダー」を知りませんでした。
11 qgis  raster  workflow 
3
ArcCatalogのフィールドフィルターを使用して、ArcSDEフィーチャクラスの多くのフィーチャデータセットを検索しますか？
最大30個のフィーチャデータセットと、最大100個のフィーチャクラスがフィーチャデータセット全体に分散しています。フィーチャデータセットの中から名前でフィーチャクラスを検索したいと思います。また、すべてのフィーチャクラスから、someDateからsomeOtherDateの「Date_Mod」を持つレコードをクエリしたいと思います。 これはすべて、SQL Server 2008上のArcGIS SDE 10.0にあります。
3
ArcGIS ModelBuilderに任意のコードを追加しますか？
ArcGISモデルに任意のコードスニペットを追加する方法はありますか？ 高度なラベル付けやフィールド計算の式ビルダーの「コードブロックの表示」に似たものを思い描いています。 たとえば、今の私のプロジェクトでは、ツールを実行している人のユーザー名を取得することで大きなメリットが得られます。私はPythonでこれが簡単であることを知っていimport os; os.environ.get('USERNAME')ます。そのスニペットをget_username.pyに記述してArccatalogのツールボックスに追加し、そのツールをモデルに追加して、現在のツールの必要なモデルパラメーターとして使用できると思います。それはこのような小さなものには多くの作業であり、次回はインライン変数置換や複合変数など、何か他のものが必要になるかもしれません。 コードスニペットツールなどは、一般的で広く役立つ方法でこれらのニーズすべてに対処できると思います。
3
複雑なシェープファイルをグリッドに分割する
ポリゴン/マルチポリゴン機能を備えたきちんとした詳細なシェープファイルがあります（ファイルは約500MBです）。実際には、全世界のシェープファイルであり、フィーチャは海岸線を表します。グリッドを使用してこのデータを分割する必要があります。明確にするために、データを「ソート」したくはありませんが、実際にはポリゴンをタイルに分割します。この質問は以前に聞かれたことに気づきましたが、見つけた解決策はうまくいきませんでした。 私はもう試した： QGISを使用し、シェイプファイルのコンテンツをベクターグリッドと交差させます-結果はひどいです。主要な陸地の大部分は魔法のように消えますが、土地の小さな塊が時々それを作るようです。この方法は、はるかに単純なデータ（つまりポイントが少ない）で非常にうまく機能することに注意してください。 OGRの交差ツールを使用します。私はogr2ogrを介して、さらには独自のC ++ツールを使用して試してみました。どちらもQGISと同じ問題を抱えています。また、単純なファイルではこの問題は発生しませんが、より複雑なファイルでは失敗します。参考までに、サイズが20 MB未満のオーストラリアとニュージーランドのシェープファイルを使用していますが、QGISとOGRの両方が「グリッド化」できません。 PostGISには交差機能があるため、ある時点でPostGISを使用することを提案しましたが、PostGISのST_IntersectはOGRと同じGEOSバックエンドを使用します。実際、私が知る限り、両方とも同じ関数を呼び出しているため、PostGISが異なる結果をもたらすとは思いません。 私は他に何ができるかについての提案を探していました。非常に詳細なシェープファイルをタイルに分割できる堅牢なアプリケーションまたはツールキットが必要です。 編集：いくつかの情報を追加する シンバマングへの応答： シェープファイルは、基本的にOpenStreetMapの海岸線データです。これは、「processed_p」ファイルのマージされたバージョンであるため（タイルに分割されません）、開発者リストをメールで送信しました。タイルの分割（重複した100 km x 100 kmチャンク）は必ずしも必要なものではないことに注意してください-重複したくないので、グリッドサイズを自由に選択したい、または単にデフォルトのprocessed_p。 デフォルトでは、海岸線データにはQGISによって報告されたジオメトリエラーがあります。これらのエラーは、この問題に特に対処するために設計されたいくつかのコード（海岸線データのジオメトリエラーの修復：https : //github.com/tudelft-gist/prepair）を使用して作成した小さなツールで修正します。このツールを使用してファイルを実行すると、QGISで検出されるほぼすべてのエラーが修正されます。ファイルをクリーニングした後にのみ、交差を試みます。 QGISを使用して行ったこととまったく同じ：データを開いて、QGISで正常に表示されることを確認します。指定した間隔でベクターグリッドを使用してタイルのレイヤーを作成し、2つのレイヤーを交差させてタイルに分割してみてください。より小さなデータセットを使用してみてください-オセアニア（オーストラリア、ニュージーランド）の機能を選択して、より小さなデータセットを試します-この形状ファイルのサイズは20 MB未満です。もう一度分割してみてください、動作しません。 OGRで行ったこと：ogr2ogrは、spat_extentで '-spat'および '-clipsrc'オプションを直接使用しました。また、WKTで動作する小さなC ++ツールを作成したため、ogr2ogrを使用してシェープファイルをWKTに変換し、テキストファイルをアプリケーションにフィードします。ファイル全体を実行し、http：//www.gdal.org/ogr/classOGRGeometry.htmlに記載されているIntersection（）メソッドを呼び出します。最終的にはogr2ogrを直接使用するのとまったく同じことを行うと思います。 ブレントへの応答： します。すべてがWGS84 Lat / Lonにあります グリッドタイルの特定のセットでは、各タイルにより空間的にローカライズされる可能性のある断片化されたフィーチャの束ではなく、1つの巨大なマルチポリゴンと交差するのに時間がかかると思いました興味深い提案-私はそれを試して報告します。 プロセス中に属性フィールドは保持されません。ジオメトリにのみ興味があります。 よくわかりませんが、特定のグリッドタイルと重なるポリゴンを選択してから交差を実行する必要があると言っていると思います。これはQGISでは手動で面倒です。私のツールは、バウンディングボックスチェックを使用して、すでにある程度これを行っています。少しスピードアップしますが、最終結果はまだ悪く、目立った違いはありません。 これはオプションではありません。現在、データを1 deg lat x 1 deg lonに分割しようとしています。すべてのケースで機能する一般的で堅牢な方法論を探しています。グリッドサイズ（10x10など）を増やして、より良い結果が得られるかどうかを確認しましたが、グリッドサイズと出力の品質との間に相関関係は見られません。 編集＃2： 私はこれをもっと試してみましたが、一般的には、GEOSとQGIS（fToolsを使用していますが、GEOSを再び使用するかどうかはわかりません）の両方で結果が信頼できないようです。グリッドのサイズが結果とは何の関係もないと述べるのは間違っていました-グリッドが大きいほど、結果は良くなります（それは知っておくのは良いことですが、解決策ではありません）。これは、ほとんど機能しているが、1つのタイルで部分的に失敗した、実際に間隔をあけたグリッドのスクリーンショットです。 ジオメトリはきれいです-QGISは「妥当性のチェック」ツールでエラーを表示しません。私はこの問題に一歩ずつ取り組むつもりはありません。特定のフィーチャが視覚的に明らかではない場合に、この大きなデータセットの交差に失敗したかどうかを確認することは（そして、小さなタイルでは不可能です）実用的ではありません。
2
縮尺地形図を表示するときの私の距離はどれくらいですか？
私の会社の親友が地形縮尺地図（1：25.000）について質問しました。 what was the altitude to a location in km (as you are viewing Earth from space) when he was viewing the map from 1 meter away? これは本当に興味深い質問です。私は長い間GISの世界で縮尺された地図を扱ってきましたが、地形図から見た高度が何であるか疑問に思いませんでした... 数式を知っている場合は、説明と共有してください...
5
失敗している大規模なモザイクプロセスを回避するにはどうすればよいですか？
約550Gbのtif画像をモザイク化する必要があり、試したソフトウェアが失敗し続けます。領域はゾーンに分割されているため、最小のタイルは約200タイルです。 3.30ギガヘルツIntel Xeon E31245、DELL、16GB RAM、64ビットWin 7 ProfessionalでERDAS（ImagineおよびMapper）、ArcINFO、およびGlobal Mapperの最新バージョンを使用しました。マルチコア（合計4）、ハイパースレッド（合計8）マシン。私のCには700GBの空き容量があり、Dには1.5TBがあります。 Grassの使用を検討しています（以前はありません）が、i.image.mosaicは4つのファイルしか処理していないようです。他に試してみるオプションやオープンソースソフトウェアはありますか？ 申し訳ありませんが、モザイクデータセット（または他のソフトウェアの同等物）は使用できません。ecwとして定義されたデータなしの領域を持つゾーンを作成し、任意のGISソフトウェアで開き、より低い解像度/古い新しいデータがシームレスに存在しない場合のデータ。 モザイク化されたファイルが異なるソフトウェアでどのように見えるかの例。Global Mapper / ERDASは問題ありませんが、arcgisでは正しくありません。 ---旧情報--- ラフな描画でごめんなさい。したがって、5つのゾーンとして色付きの領域を使用すると、より大きなAOIでデータ領域がなくなります。 arcgisのコードは次のとおりです（これはモデルとして実行され、PythonではなくtifList入力を取得できないため）。 arcpy.MosaicToNewRaster_management(tifList+";" +mask,RootOutput,"Tile1.tif","PROJCS['GDA_1994_MGA_Zone_55',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',147.0],PARAMETER['Scale_Factor',0.9996],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]","16_BIT_UNSIGNED","0.5","3","MAXIMUM","#") # Replace a layer/table view name with a path to a dataset (which can be a layer file) or create the layer/table view within the script # The following inputs are …
1
データドリブンページを使用してテーブルを作成するPythonスクリプト
特定のデータドリブンページ内にある機能のテーブル（dbfに基づく）を表示するために、Pythonスクリプトを変換しようとしています。これまでのところ、特定のテーブルへのマップを正常に更新できるスクリプトがありますが、テーブルは更新されません。 ユーザーがArcToolboxからスクリプトを実行するときに3つの特定のフィールドで更新する必要がある3つのテキストボックスとして設定しています。 テーブルが更新されない理由に関する提案はありますか？ import arcpy, sys, os #Reference current MXD mxd = arcpy.mapping.MapDocument("current") #Get input parameter Name = arcpy.GetParameterAsText(0) #Reference data frames mapatlasDF = arcpy.mapping.ListDataFrames(mxd, "Layers")[0] locatorDF = arcpy.mapping.ListDataFrames(mxd, "Locator Map")[0] #Reference appropriate layers atlasLyr = arcpy.mapping.ListLayers(mxd, "PinalCreekMapAtlas_HalfMile", mapatlasDF)[0] locatorLyr = arcpy.mapping.ListLayers(mxd, "Locator Map", locatorDF)[0] atlasoutlineLyr = arcpy.mapping.ListLayers(mxd, "Map Atlas …
1
IPジオコーディングはどのように機能しますか？
IPアドレスの各部分の使用がジオコーディングの精度にどのように影響するかを理解しようとしています。たとえば、 15.201.xxx.xxx の代わりに使用されます 15.201.123.xxx KMのおおよその距離は？IP2Location Liteデータベースを使用して、誰がWebサイトにアクセスしたかを調べています。README次のことを述べています： 1. IP2Location Lite Edition is free package with accuracy up to Class C (192.168.1.X) only. It is restricted for non-commercial use. ただし、これの正確性が実際に意味するものに関するドキュメントを見つけるのは困難です。この精度は時間とともに変化しますか？
11 geocoding 
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.