地理情報システム

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

4
マップポイントの周りに運転時間ポリゴンを作成しますか?[閉まっている]
閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか?この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 3年前に閉店しました。 Googleマップにプロットした400個の個別のポイントの周りに30分の運転時間ポリゴンを作成したいと思います。 http://www.freemaptools.com/how-far-can-i-travel.htmを発見しましたが、海に近づくとこれはおかしいようです(私のポイントの多くは)。 また、http: //cartoo.dyndns.org/も機能しますが、.KML形式でポリゴンを抽出できません。 誰かが私を助けることができるツールに出会いますか?

5
GPSログの地理データベースを構築する方法は?
航空野生生物調査では、通常、数十のGPXファイルを生成し、通常、それらのほとんどをGISアプリケーションで表示するためのシェープファイルに変換します。GPXファイルは、GPSデータのすべて(またはほぼすべて)を1か所に保管する素晴らしい方法です。これらは、トラックログ、ウェイポイント、ルート情報を1つの便利な形式で保存するxmlファイルです。それらは複数のトラックを保存し、Garminから保存された古いトラックも保存します。 同じ追跡からのポイント観測データを持つ他のテーブルと並行して、このすべての追跡データを空間データベースに保存したいと思います。 ogr2ogrを使用すると、トラックを空間データベースに変換できます。 ogr2ogr -append -f "SQLite" -dsco SPATIALITE=yes # add to table with Spatialite filetype -dsco INIT_WITH_EPSG=yes -t_srs epsg:4326 # using WGS84 SURVEYDATA.sqlite GPXFILE.gpx tracks -nln tracktable; # adding tracks to 'tracktable' ただし、これにより、保存された各トラックが線としてインポートされ、各トラックポイントの時間/位置情報が失われます。 その情報を保存するトラックログデータベースを構築する賢い方法を知っている人はいますか? 編集: トラックログは単なる線ではなく、それぞれが標高と時間の属性を持つ連続したポイントのコレクションです。各ポイントは、特定のトラックに属するトラックセグメントに属し、各トラックセグメント内のポイントは中断することなく取得され、リンクされていると見なされます。例えば: <trk><name>ACTIVE LOG</name> <trkseg> <trkpt lat="-13.471403" lon="31.382982"><ele>467.818725585938</ele><time>2009-09-09T03:53:38Z</time></trkpt> <trkpt lat="-13.471403" lon="31.382982"><ele>468.780029296875</ele><time>2009-09-09T03:53:50Z</time></trkpt> <trkpt lat="-13.471403" lon="31.382982"><ele>465.896118164063</ele><time>2009-09-09T03:54:43Z</time></trkpt> </trkseg> …
13 gps  convert  ogr  spatialite  gpx 

3
Pythonリストにある場合、属性によって機能を選択しますか?
Pythonで属性による選択を完了しようとしていますが、属性がリストに存在するかどうかのクエリに基づいています。 そのような最も単純なクエリは、次のようになります。 qry = " \"OBJECTID\" in oid_list" arcpy.SelectLayersByAttribute_management(inft, "NEW_SELECTION", qry) しかし、そのアプローチは無効な式エラーを返します。 過去には、このタイプのクエリには、次のようなより複雑な構文を使用する必要がありました。 sqlQuery2 = "nid in (" + ','.join(["'"+x+"'" for x in delta_list]) +")" しかし、このスニペットの適応も私にとってはうまくいかないようです。 "OBJECTID_1 in (" + ','.join(["'"+str(x)+"'" for x in oid_list]) +")" ここで何が欠けていますか?

2
GDALを取得してPythonでGTiffの統計を作成する方法
私は定期的にPythonでGDALを使用して自分のGeoTIFFラスターを作成します。例: from osgeo import gdal from numpy import random data = random.uniform(0, 10, (300, 200)) driver = gdal.GetDriverByName('GTiff') ds = driver.Create('MyRaster.tif', 200, 300) band = ds.GetRasterBand(1) band.WriteArray(data) ds = band = None # save, close ただし、結果をArcCatalog / ArcGISで表示すると、統計がないため、黒または灰色に見えます。これは、ArcCatalogでラスターを右クリックして[統計の計算...]を選択するか、コマンドプロンプトでgdalinfoを使用することで解決できます。 gdalinfo -stats MyRaster.tif MyRaster.tif.aux.xmlArcGISがラスターを適切にスケーリングするために使用するを生成します。PAM(Persistent Auxiliary Metadata)ファイルには、統計、特に最小値と最大値が含まれています。 <PAMDataset> <PAMRasterBand band="1"> <Metadata> <MDI key="STATISTICS_MINIMUM">0</MDI> …

2
マップを仕上げるためのオープンソースソリューション
Adobe Illustratorの600ドルの値札を見た後、GISで作成されたマップを仕上げるためのオープンソースソリューションの学習に特に興味があります。Arc、QGIS、または選択したGISソフトウェアで作成されたマップに最後の仕上げ(矢印、テクスチャ、色、グラフィック、テキストなど)を追加するために、どのオープンソースまたは一般に利用可能なソフトウェアを使用しますか?例は歓迎します。



2
ArcSDE(SQL Server)からPostGIS(PostgreSQL)への移行
ArcSDEデータベース内に(SQL Server経由で)弊社のすべての地理空間情報があります。 ポイント(アンテナの位置)、 線(標高、河川)、 ポリゴン(保護されたエリア、家)。 これらのデータは、ArcGIS DesktopおよびArcGIS Serverで使用します。 すべての情報を(PostgreSQLを介して)PostGISに移行することを計画しており、いくつかの手順と留意すべき考慮事項に従うべきドキュメント/ウェブページがあるかどうか疑問に思っています。データモデルをエクスポートするのは非常に難しいですか? 私は調査し、GDAL-Pythonスクリプトを使用してSDEから読み取り、PostGISに書き込むいくつかのテストを行ってきましたが、経験はありがたいです。 ところで、ArcGIS Desktop(QGISとともに)およびGeoServerでPostGISからの情報を使用する予定です。 事前にどうもありがとうございました。

1
多角形内の線の長さを取得する方法
赤いポリゴン内の「Road 7」の長さ(茶色でマークされた部分)を取得するにはどうすればよいですか? 私のスタート: SELECT ST_Length(way) FROM lines WHERE road = "Road 7" AND ST_Intersects(line, polygon)? それで、私は全長だけを取得します:(誰かがアイデアを持っていますか?

4
米国のDEMデータをダウンロードしていますか?
USGSは、The National Map Viewerを通じてデータを提供する方向に動いています。Webサイトには、標高データを一般公開する(陰影起伏&DEM)いくつかのWMServiceがあります。問題は、ArcToolboxのジオプロセシングツールがこれらのレイヤーを入力として受け入れないことです。 分析に使用できる(無料の)DEM gisデータはどこでダウンロードできますか?特に興味があるのは、1/3秒(〜10m)または1秒(〜30m)のデータです。 現在、USGS NED Tileダウンロードツールを使用しています。このツールには、郡ごとに最大4つのファイルがあります。関心のある分野は約50の郡をカバーしています-これらのファイルを1つずつ手動でダウンロードすることはありません。

4
Googleマップで住所を一度ジオコーディングし、データベースに保存し、そこからジオコーディングされた結果を読み取ります
住所をジオコーディングし、結果をデータベースに保存し、そこから結果を読み取ることは可能ですか? ジオコードの制限に明らかに達する数百の住所をジオコーディングしているGoogleマップを使用しています。ジオコーディングされた住所は再び変更されることはほとんどないため、ジオコードの結果をデータベースに保存することは理にかなっています。 プラットフォームとしてWordpressを使用しているため、可能であれば、特定の会社に関連するデータベース内のフィールドとしてジオコーディングされた結果を保存したいと思います。 クライアントは、ヤフーやその他のあらゆるものを排除するために、プロジェクト全体でGoogleを使用することを特に求めています。 UPDATE:みんなの言って、私が求めているものについては、そのTOS逆行が、同じことをやって非常に独自の例ではありません- developers.google.com/maps/articles/phpsqlgeocode?

2
ArcGISバッファー精度の問題?
もう一度、ArcGISで非常に奇妙な何かに遭遇します。ポリラインの片側でバッファー分析を実行します。深刻な精度の問題があるようで、原因がわからない。だから私はそこにいくつかの専門知識を見つけたいと思っています...私が成功せずに試したこと: 距離設定を「フィールド」から「線形単位」に変更します 距離設定を「フィールド」から「線形単位」に変更します 終了タイプの設定を「ラウンド」から「フラット」に変更します ラインフィーチャでジオメトリチェックツールを実行します 回線のコピーを作成し、それを何らかの影響を与えたが問題を完全に解決しなかった関心のある領域に短絡させました。特に、鋭角の周りでは、バッファはまだショートカットを使用します! バッファーは、反対側で最大4メートルまで傾いています!問題のより良い印象を得るために私が添付した写真を見てください。誰が私がこの行動を引き起こしているのか考えていますか?ArcGIS 9.3.1を使用しています。WinXPで。 更新 以下の私のコメントも参照してください。 説明した一連のアクション(分割、交差など)を実行して、ソースデータの精度レベルを上げ、ラインに沿って3メートルごとに頂点を配置しました。依然として不正確さが残っていますが、特に大まかなターンの周りのバッファリングは特にトラブルを引き起こしています。2枚目の写真をご覧ください。

4
シェープファイル内の長い線は、測地線と見なされるべきですか、それとも2D緯度経度空間内の直線と見なされますか?
線がシェープファイル形式で頂点を接続する方法についての定義はありますか? 最も単純な場合、標準のWGS84地理座標系を使用して、40、-118から40、-112までの2つのポイント(米国ではどこかランダム)のラインを想像してください。.prjファイルの内容は次のとおりです。 GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]] ポイントは、ラインの40.1、-116北または南ですか? 線が緯度の長い空間で直線的に補間されると考えると、それは40度の平行(小さな円)に従い、点は線の北になります。 線が地球表面上の最短経路であると考える場合、それは線の中央で最大緯度が40.1度を超える測地線(大円)です。その後、ポイントはラインの南になります。 それとも、単に未定義ですか?シェープファイル形式には曲線の概念はないため、直線を結ぶ直線セグメントのみがあります。この答えを明確にするには、ラインを高密度化する必要があります(ラインに沿ってポイントを追加)。 QGISでこのようなシナリオを作成すると、線は40度の平行線をたどり、答えは1であるとわかります。しかし、これを明確な答えとして受け取らず、より堅実な答えを聞きたいと思います。

1
ArcGIS 10の「フィールドマッピング」-ArcPy
空間結合といくつかの簡単な計算を行うPythonスクリプトを作成しました。私の問題は、ある特定のフィールドにマージルールを設定し、残りのフィールドをそのままにしておくことです。たとえば、空間の場所で結合されたときに、最初に発生した人口カウントを取得するマージルール「最初」を使用する人口フィールドがあります。別の多角形の空間範囲で見つかったすべての多角形間の人口値を合計するために、マージルールを「合計」に設定できるようにします。 フィールドマップとフィールドマッピングオブジェクトを徹底的にいじくり回しましたが、適切に機能していないようです。具体的には、メソッドpopFieldMap.mergeRule = 'Sum'を使用してmergeRuleを設定しましたが、常に「最初」に戻ります。 空間結合の1つのフィールドのマージルールをプログラムで変更する方法はありますか? ありがとう! ここに私のコードがあります(それは私のデータに非常に固有のものであり、スクリプトの特定の段階をテストする行が含まれていることに注意してください): import arcpy,sys,os #Get the Files involved, set some variables. SectorTable = sys.argv[1] SectorShape = sys.argv[2] MaxDev = sys.argv[3] PopulationFC = sys.argv[4] OutputFC = sys.argv[5] DeviationField ="Angle_Deviation" ID = "SectorID" newID = "BP_ID" mxd = arcpy.mapping.MapDocument('CURRENT') df = arcpy.mapping.ListDataFrames(mxd)[0] #Check to see if ID …

7
無料/オープンソースのGIS教育ソース
高度なGISユーザー向けの無料/オープンソースのGIS教育ソースを探しています。具体的には、プログラミング、空間分析、モデリング、空間データベース/データ管理をカバーするソース。私は主に連邦政府の従業員としてArcGISを使用していますが、概念の概念に触れる情報源も高く評価しています。 私はGISプログラムのペンシルベニア州立マスターズの大ファンです。なぜなら、彼らはMITライセンスの下でフルクラスの多くを提供しているからです。ダブリン工科大学の講義も読んでいます。私はまだ出会っていない他の高品質のリソースがあることを望んでいます。 ある程度のインタラクティブ性を備えた実際のコース(クイズやレッスンに関連するタスクなど)を好む一方で、本が高品質で無料のリソースである場合は本を検討する必要があります。

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