地理情報システム

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

4
RまたはPythonを使用して大きなラスターをポリラインに変換する最も速い方法は？
グローバルな水域（1ビット値0および1）の大きなラスターファイル（129600 x 64800ピクセル）があり、海と内陸の海岸線を抽出しようとしています。 ArcGISとQGISでラスターからポリラインに変換しようとしましたが、時間がかかります。 誰かがこのタスクのためのより良い/より速い方法（PythonまたはR）またはより良いツールを知っていますか？ 更新 R：rasterToContourは高速で正確かもしれませんが、私のような非常に大きなデータセット（8,398,080,000ピクセル）がある場合、非常に大きなRAM（16GB以上）が必要か、Rにハードドライブとそれ以上の処理を強制します年齢もかかります。 Python / GDAL：gdal_poligonizeはポリラインの代わりにポリゴンを作成します 更新2 R rasterToContour：rasterToContourは、必要な結果を提供しません。下の例に示すように、ArcGIS（ラスターからポリゴン、フィーチャからライン）と比較すると、正確なピクセルアウトラインは抽出されません。 rasterToContourの結果 ArcGISの結果 更新3 Python / GDAL：テストデータセットでコマンドラインからgdal_polygonizeをArcGISに対して実行しましたが、結果は非常に明確でした。 gdal：49秒 ArcGIS：1.84秒
14 python  r  performance 
4
.shpをPostgreSQLにインポートしますか？
shpファイルをpostgreSQLにインポートしたい。最初にsqlファイルを作成してから、PostgreSQLを実行します。SQLファイルを作成するには、Windows cmdで次のコマンドを実行します。 shp2pgsql -s 4326 worldCountries.shp worldcountries postgres > worldcountries.sql 次に実行します： psql -d postgres -U postgres -p 4321 -f worldcountries.sql しかし結果は： psql：worldcountries.sql：21：エラー：関数addGeometrycolumn（unknown、unknown、unknown、unknown、unknown、integer）は存在しません
5
線分からボロノイ図を作成
ポイントの代わりに線分に基づいてボロノイ図を作成する方法を探しています。以下の例を参照してください（このプレゼンテーションから引用）。 Pythonからスクリプトを作成できるものが欲しいのですが、ArcGISなどを使用したソリューションでも問題ありません。 私がこれまでに見つけた唯一のライブラリは、有望に見えるopenvoronoiです。他に何かありますか？
2
RのSpatialPointsDataFrameプロパティと演算子
R SpatialPointsDataFrameのspパッケージを使用してタイプのオブジェクトを作成しました。ただし、@, $, . and []演算子と、オブジェクトのさまざまなプロパティにアクセスするために演算子を使用するタイミングについて混乱しています。これが私のサンプルコードです。 library(sp) library(rgdal) #creating a SpatialPointsDataFrame with sample points in UTM x <- c(15.2, 15.3, 15.4, 15.5, 15.7) y <- c(50.4, 50.2, 50.3, 50.1, 50.4) v1 <- c(1.0, 2.0, 3.0, 4.0, 5.0) v2 <- c("a","b","b","c","a") attributes <- as.data.frame(cbind(v1,v2)) xy <- cbind(x,y) locationsDD <- SpatialPointsDataFrame(xy, attributes) …
14 r  syntax 
1
PostGIS：テーブル内のラインとポリゴンの頂点の数を数える
PostGISスキーマには、頂点の数をカウントする必要がある23個のテーブルがあります。テーブルには線と多角形が混在しているため、ST_NPoints（geom）を使用する必要があることに気付きました。 だから私は次のクエリを実行しました SELECT count(ST_NPoints(geom) FROM lines; 結果/カウントは、そのテーブル内のすべてのフィーチャの頂点の総数ではなく、そのテーブル内のフィーチャの数に等しくなります。 私は何かを見逃しているに違いありませんが、それを理解することはできません（月曜日の朝でなければなりません;））
14 postgis 
2
Python、GDAL、およびラスター属性テーブルの構築
PythonとGDALを使用してラスター属性テーブルを構築したい整数ラスターがあります。次のように、PythonでGDALラスター属性テーブルを作成できます。 >>> rat = gdal.RasterAttributeTable() ご覧のとおり、これは正常に機能します。 >>> rat <osgeo.gdal.RasterAttributeTable; proxy of <Swig Object of type 'GDALRasterAttributeTableShadow *' at 0x0000000002A53D50> > このようにして作成されたテーブルには行も列もありません。 >>> rat.GetRowCount() 0 >>> rat.GetColumnCount() 0 「値」という列を作成して、ラスターに一意の値を保存します。 >>> rat.CreateColumn("Value", gdalconst.GFT_Integer, gdalconst.GFU_MinMax) 0 これで問題ありません。列数が更新されます。 >>> rat.GetColumnCount() 1 ここで、列を使用するために値（レコード）を列に追加する必要があります。次のように、ラスターのバンドから一意の値のリストを取得できます。 >>> data = band.ReadAsArray(0, 0, dataset.RasterXSize, dataset.RasterYSize) >>> vals = list(numpy.unique(data)) >>> …
2
ラスターから大きなベクターポリゴンを滑らかにする方法
大面積の土地利用分類ラスターファイル（6クラス）があります。しかし、ベクターポリゴンで必要なので、ベクター化しました。これまでのところ問題はありませんが、結果は-もちろん-すべてのピクセルがすべてのegdesで正確に渦化されていることです:)しかし、輪郭のように滑らかにする必要があります。だから、大丈夫、たとえばボイルのアルゴリズムで一般化することはできますが、それを計算するにはあまりにも多すぎます... だから私は何ができますか？ありがとう！マーティン
2
QGISで不連続な機能を分割するにはどうすればよいですか？
連続していない機能を持つシェープファイルがあります。左側の領域と右側の領域を分割します。ただし、QGISの機能分割ツールでは、機能の上に線を引く必要があります。この機能の2つの領域の間に線を引くと、何も起こりません。 この機能をどのように分割すればよいですか？ これがスクリーンショットです。分割したい機能は黄色です。
14 qgis  shapefile 
2
画像の場所を特定していますか？
植生の多い都市の一部のみの衛星画像があります。その画像がどこで撮影されているかを知る必要があります。 Google Maps、Google Earth、または同様のものでどこにあるかを知る方法はありますか？ これは.jpg画像なので、検索を開始する方法がわかりません。
4
ArcGIS Serverからキャッシュサービス上のホットスポットを識別するHeatMap
ArcGIS Serverでは、キャッシュされたマップサービスを作成できます。これらはさまざまなアプリケーション（ArcGIS Desktop、Webアプリなど）で使用できます。キャッシュのどの領域が最もヒットしているのかを知りたいです。同じキャッシュマップにオーバーレイされたヒートマップを使用して、これらの結果を視覚化したいと思います。HeatMapsAPIとESRIが提供するAPIの1つを使用することを考えていました（Javascript APIで十分です） キャッシュされたマップはArcGIS Serverのv10でわずかに変化する可能性がありますが、ArcGIS ServerまたはWebサーバー（IIS）に大きな負荷をかけることなく、どのタイルが最も多くリクエストされるかを判断する最良の方法を考えています。 私はもともと、arcgiscacheフォルダーの下にある実際のキャッシュディレクトリを見て、おそらくDate Accessed属性にフックすることを考えていました（ただし、実際にこれを実装する方法についてはほとんど知りません）。 ArcGIS Serverの経験がある人から、これを行う方法について何か提案があれば教えてください。
2
QGIS印刷コンポーザーでヒストグラムの凡例を作成する方法は？
印刷コンポーザでグラフ（ヒストグラム）を動的に作成する方法を探しているときに、このブログ投稿を見つけました。 ブログから複製したい機能は、qgis print composerで静的な長方形を描画し、この場合のフィールド値に基づいてヒストグラムのバーとしてサイズを変更することです。 したがって、凡例/ヒストグラムを作成します このチュートリアルに従うことはできますが、a）matplotlibで必要なグラフを作成して画像をインポートする以外に簡単な方法があるか、b）この機能が将来のバージョンでqgisに組み込まれるかどうかに興味がありました。
3
SQL ServerジオメトリBLOBから他の何かに変換する方法は？
エクスポートされたSQL ServerテーブルであるサードパーティからのCSVデータファイルがあります。彼らは単にselect * from foo結果をテキストファイルに出力し、それを送信しました。 テーブルにはGeometryタイプの列があるので、生のテキストには「0xE610000010C47 ...」などがあります。現時点では、SQL Serverのテーブルにnvarcharとしてロードしています。 私は自分の側でGeometryフィールドに戻ることができると予想していましたが、それはそれほど簡単ではないようです。 STGeomFromWKB実際にはWKBではないため機能しません。文字列をGeometryとしてキャストできません。WKTではないと文句を言うからです。 それでは、この値を通常のGeometry BLOBのようにSQL Serverに取得する方法はありますか？SQL Serverにそのように処理するように指示できますか？ このリンクを見つけて、少なくともSQL Serverの内容に関する質問に答えることができましたが、そこまでは行きませんでした： SQLServer 2008のGeometryデータ型の形式は何ですか
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.