地理情報システム

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


2
GeoJSONをTopoJSONに変換
質問 GeoJSONデータをTopoJSONに変換するにはどうすればよいですか? なぜ Topojsonファイルは多くの場合geojsonよりも小さく、必要に応じて変換によりデータセットを段階的に簡素化できます。

7
土地被覆分類のための機械学習アルゴリズム
機械学習アルゴリズム(k-NN、ランダムフォレスト、決定木など)を使用して土地分類にどのソフトウェアが存在するかを知りたい。RのrandomForestパッケージとPythonのMILKとSPyを知っている。 土地被覆分類に適したオープンソースまたは商用の機械学習アルゴリズムは何ですか?


1
QGISで新しいカスタム投影を定義しますか?
QGISで投影を定義しようとしています。北米Albers Equal Area Conicを追加します。空間参照サイトから次の行をコピーしましたproj4。 +proj=aea +lat_1=20 +lat_2=60 +lat_0=40 +lon_0=-96 +x_0=0 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs 次にQGIS -> Settings -> Custom CRS、次のウィンドウに貼り付けました。 しかし、QGIS内で投影を表示すると、まったく異なるパラメーターセットが表示されます。 これでいい? 生成されたパラメーターが正しいものであるかどうか、または他に何かが起こっているかどうかを判断するための投影については、十分に知りません。また、指定した名前も無視したため、パラメータが間違っている可能性があります。 Windows 7でQGIS 1.7.3を使用しています。


1
オブジェクトのジオメトリの中心を見つけますか?
2Dまたは3Dポイントのセットが与えられた場合: オブジェクトのジオメトリの中心を見つける方法は? 次の図によると、幾何学の中心は、最も単純な形式、つまり均質な質量密度で計算される場合、質量の中心とは異なります。実際、問題はそれらの計算に現れます。一般的に、1つのアプローチは、X座標とY座標を個別に平均することです。つまり、指定されたポイント(ここでは2D)の平均位置を見つけます。これは、オブジェクトを表す一連のポイントの重心として使用できます。示されているように、下端に沿って余分な頂点があるため、単純な長方形の場合、結果の重心は(0.5,0.4)であり、正解は(0.5,0.5 )です。 指定された例は単純すぎることに注意してください。ただし、2Dの複雑な形状と、頂点の座標のみが使用可能な3Dのオブジェクトの場合に問題になります。 ところで、効率的な計算方法が興味深いです。 WikipediaのようないくつかのWebリンクをチェックしましたが、私の現在の問題は、それらの代表としてポイントを見つけたい2Dおよび3Dポイントのグループがあることです。したがって、重心が関心を集めました。ポイントはトポロジ情報なしで与えられます。それらを点群と見なすことができます。ここでのデモンストレーションは、座標の一般的に知られている平均化(たとえば、このStack Overflow Q&Aを参照)が例に示すように間違っている可能性があることを明確にするために提供されました。 比較のための実装を次に示します。 aa =以下の受け入れられた回答 chull =点の凸包、つまり黄金の多角形 cent =ウィキペディアで提案された重心で、ポリゴン重心としてaaで議論されています centl = aaで説明されているポリラインの重心 視覚centl的には、に比べて与えられたジオメトリをよりよく表しcentます。他の2人はここで有望に見えますが、通常、ポイントの分散が不均一である場合、通常は偏りすぎます。 また、凸包は問題を合理的に単純化しますが、空間内の対称位置付けなしでは長すぎるエッジと短すぎるエッジを生成する可能性があることを考慮してください。 :ポイント全体(緑)または凸包ポリゴン頂点(青)。 1つのアプリケーションは、特定のポイントの最小面積長方形を見つけることで見つけることができますか?。

4
Pythonを使用したシェープファイルの簡単なテーマ別マッピング?
ArcGIS / ArcPyを使用せずにPythonで地理データを視覚化し、地図を作成したいと思います。 インターネットで、Pythonを使用して主題図を作成する方法を見つけました。 私が試したコードは次のとおりです。 import shpUtils import matplotlib.pyplot as plt shpRecords = shpUtils.loadShapefile('C:\\Users\\shapefile.shp') for i in range(0,len(shpRecords)): x = [] y = [] for j in range(0,len(shpRecords[i]['shp_data']['parts'][0]['points'])): tempx = float(shpRecords[i]['shp_data']['parts'][0]['points'][j]['x']) tempy = float(shpRecords[i]['shp_data']['parts'][0]['points'][j]['y']) x.append(tempx) y.append(tempy) plt.fill(x,y) plt.axis('equal') plt.title("Testing") plt.show() ただし、これを実行すると、ランダムな色が表示されます。 シェープファイルの特定の列を視覚化する場合、同様のコードを使用してこれを実装するにはどうすればよいですか? これは、彼が色の使用法についてのみ議論している上記のリンクでは非常に不明確です... デカルトやPySALなど、これを実現するためにおそらく追加のモジュールが必要でしょうか?

3
QGISにラベルを手動で配置しますか?
QGISにラベルを手動で配置する方法はありますか? 一部の場所(E 1st Stと1st Stなど)にほぼ重複したラベルがあるマップがあり、それらのいくつかを削除したいです。 また、いくつかのラベルを移動し、現在ではない機能にラベルを付けたいと思います。 ラベルをフィーチャに沿って湾曲させた新しいラベルを使用しています。

2
天体座標の地理データを保存するには、どの座標系を使用する必要がありますか?
私は天文学プロジェクトをしています。画像に関する情報を空間的に有効化されたデータベースに保存したい。空は完全に球体として扱うことができ、地球の表面のような楕円形の処理を必要としないため、これはGIS機能の非常に簡単な特別なケースになるはずです。残念ながら、私はまだこれを行う方法を見つけていませんし、楕円地球を使用する空間関数で地雷をかわしています。(度ではなくメートルを返すほとんどの関数は、楕円計算を使用している可能性があります。幸いなことに、必要なPostGIS関数の多くは、返された結果が球ではなく楕円体。しかし、それは将来のバージョンで変わる可能性があり、懸念の原因です) 背景:現在、PostGISおよびWGS 84座標(SRID = 4326)でPostgreSQLを使用しています。これはかなりうまくいきます。画像の四隅の赤経と赤緯から閉じたポリゴンを作成しています。たくさんの画像(10k以上)があり、空の広い領域をカバーしています。各画像は約1度の正方形です。これらの画像のセットから、15〜30個の画像の小さなサブセットからモザイクを作成しています。各モザイクは約1.5度の正方形です。 現在、モザイクの地理を、モザイクに入った各画像に対応するすべてのポリゴンで構成されるマルチポリゴンとして保存しています。[より良い解決策は、すべての個々のポリゴンの結合の周囲を記述する単一のPOLYGONを作成することです。これが球面座標(つまり、地理タイプ)で実行できるかどうかはわかりません。これも私にとって興味深い答えです。]日付変更線と天の極がデータセットの画像に含まれている可能性があるため、可能な限り平面座標への投影を避けています。 PostGIS関数を使用して天体座標に使用する必要がある座標系は何ですか? 私はhttp://spatialreference.org/を見まし たが、これまでのところ何も見つかりませんでした。Googleはほとんど姿を見せていない。私は困惑しています。基本的に、関数が距離としてメートルを返す場合、それは球上の大円に沿ったメートルであることを確認したいです。 より一般的には、空間データベースで天体座標を使用するためのアドバイスも歓迎します。 PostGISを選択することでエラーが発生しましたか? はるかに優れた商業的選択肢はありますか? FOSSの選択? PostGIS 1.5.2を使用しています。PostGIS 2.0はまだ試していません。ST_CoveredBy関数がタイプgeographyのPOLYGONおよびMULTIPOLYGONで機能するかどうか興味があります。誰かが2.0を実行している場合、これと同じエラーが発生するかどうか教えてください: mydb=# select ST_CoveredBy(ST_GeographyFromText('MULTIPOLYGON(( (10.37795 -69.57926,8.9498 -69.54875,9.0178 -69.21643,10.4242 -69.24648,10.37795 -69.57926),(10.42436 -69.24618,9.01774 -69.2162, 9.08363 -68.88389,10.46914 -68.91344,10.42436 -69.24618)))'),ST_GeographyFromText('POLYGON((10.46915 -68.91315,9.08371 -68.88364,9.14755 -68.5513,10.5125 -68.58038,10.46915 -68.91315))')); ERROR: geography_covers: only POLYGON and POINT types are currently supported CONTEXT: SQL function "st_coveredby" …


12
html5地理空間アプリケーションの最新技術とは何ですか?
canvas地図のhtml5 要素に非常に興味があります。たとえば、MIT Media Labの研究者のCartagenのような仕事は非常に有望に見えます。carto.netには、SVGベースの興味深い作品がいくつかあります。 WebSocketライブ地理空間データ用の素晴らしいAPIのようです。FlickrとPicasaのデータを使用して、いくつかのhtml5実験を行いました。たとえば、ここにあります。これらの新しいテクノロジーに人々は何に慣れていますか、それとも何を試しましたか?

15
無料ツールまたはArcGIS for Desktopコア製品を使用してMapInfo TABおよび/またはMIFをSHPに変換しますか?
無料のツールまたはArcGISのみを使用して、MapInfo TABやMIFファイルをSHPに変換するにはどうすればよいですか? 過去にこれをやろうとしてあまり運がありませんでしたし、ArcGIS Data Interoperabilityエクステンションを探し出す必要はありません。

21
ArcSDEを使用せずに表示とクエリのためにArcMapからPostGISデータベースに接続しますか?
ArcGIS Desktop 9.3以降を使用して、ArcMapからPostGISデータベースに接続するにはどうすればよいですか? テーブルの内容を単にダンプするのではなく、空間的に有効化されたクエリを実行し、結果(空間および非空間結合、フィルタリングなど)を受信できるようにしたいと思います。 ArcSDE空間拡張を使用したくない、ArcGIS DesktopでPostGIS空間拡張を使用したい。


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