タグ付けされた質問 「vector-grid」

どのようにRを使用して シェープファイルを200メートルの正方形/サブポリゴンに分割し、 以下の元のマップ上にこのグリッド（各正方形のID番号を含む）をプロットし、 特定の地理座標がどの正方形にあるかを評価します。 私はGISの初心者であり、これはおそらく基本的な質問ですが、Rでこれを行う方法に関するチュートリアルは見つかりませんでした。 これまでに行ったことは、NYCのシェープファイルを読み込んで、いくつかの典型的な地理座標をプロットすることです。 以下のデータを使用してこれを行う例（Rコード）を探しています。 # Load packages library(maptools) # Download shapefile for NYC # OLD URL (no longer working) # shpurl <- "http://www.nyc.gov/html/dcp/download/bytes/nybb_13a.zip" shpurl <- "https://www1.nyc.gov/assets/planning/download/zip/data-maps/open-data/nybb_13a.zip" tmp <- tempfile(fileext=".zip") download.file(shpurl, destfile=tmp) files <- unzip(tmp, exdir=getwd()) # Load & plot shapefile shp <- readShapePoly(files[grep(".shp$", files)]) plot(shp) # Define …

32 r  vector-grid  point-in-polygon 

QGISを使用してフィッシュネットシェープファイルを作成するにはどうすればよいですか？ ArcGISの「フィッシュネットの作成」と同等のタスク。

23 qgis  vector-grid 

不規則な間隔のポイントデータベースからラスターイメージを作成しようとしていました。データは次のようになります > head(s100_ras) x y z 1 267573.9 2633781 213.29545 2 262224.4 2633781 69.78261 3 263742.7 2633781 51.21951 4 259328.4 2633781 301.98413 5 264109.8 2633781 141.72414 6 255094.8 2633781 88.90244 作成したメッシュ内にこれらの「z」値が必要です # Create a fine mesh grid my_mesh=expand.grid(seq(min(s100_ras$Y),max(s100_ras$Y),l=100), seq(min(s100_ras$X),max(s100_ras$X),l=100)) また、データポイントの外側にあるメッシュポイントのz値を「NA」として割り当てる必要があります。メッシュ上のポイントは次のようになります。https : //drive.google.com/file/d/0B6GUNg-8d30vYzlwTkhvaHBFTnc/edit?usp=sharing when plot plot(my_mesh) points(s100_ras$Y, s100_ras$X, pch="*", col='blue') 問題は、これに基づいて構築する方法がわからないということです。メッシュグリッドとデータポイントは同じスケールではないため、次の手順は機能しません!! …

13 raster  r  vector-grid 

調査ツール「Regular Points」を使用してポリゴン内に特定の数のポイントを作成しましたが、ポリゴンの形状内だけでなく、長方形のポイントグリッドを形成します。ポリゴン内に特定の数のポイントが必要なので、試行錯誤せずにそれを行う方法はありますか？ 例えば、私が望む25ポイントのうち8ポイントのみがポリゴン内にありますが、25ポイントが欲しいので、50ポイントまで増やして、ポリゴン内にいくつ表示されるかを確認できます

12 qgis  polygon  vector-grid 

ヒートマップを保存および表示するために使用されるベクトルグリッドを作成する必要があるアプリケーションを作成しています。次の要件があります。 地球全体をカバーできます。 グリッドの正方形の大部分には値がありません。 グリッド自体を保存する必要はありません。その場で計算したいです。 グリッドで使用されるデータの規模は大きく異なる可能性があります。 グリッドは、1 kmから100 kmの範囲の正方形にしたいと思っています。（これがいくつになるかは承知しています（1kmで〜510million、100kmで〜51,000））。 値は各グリッドの正方形に対して累積/集約されます。 理想的には、大きなグリッドセル値を保存するのではなく、小さなグリッドセルを使用して、大きなグリッドセルの値を簡単に計算できるようにします。 私はOpenLayersを使用してOpenStreetMapの上にドレープします。 SpatiaLiteまたはSQLiteに保存するので、できればネイティブでサポートすることが望ましい（つまり、SpatiaLite =サポートされているCRS、またはSQLite =純粋な数値ベースのシステム）。 だから私の質問は次のとおりです。このグリッドにはどの投影法を使用する必要がありますか？ また、これを設計する良い方法はありますか？誰もがこの問題の潜在的な解決策を知っているか、以前に同様の問題を解決したことがありますか？または、有用な方向に私を向けることができます。 ユースケースの編集 -基本的に、さまざまな形状とサイズの境界ボックスを集約しています。それらは、数ヘクタールから数千平方キロメートルまでのサイズのものです。彼らはまた、異なる予測にあるかもしれません。 以下は、私が大規模に自動的に達成することを目指している種類の特注バージョンです。 どうもありがとう。

11 coordinate-system  spatialite  heat-map  sqlite  vector-grid 

ポリゴン/マルチポリゴン機能を備えたきちんとした詳細なシェープファイルがあります（ファイルは約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は「妥当性のチェック」ツールでエラーを表示しません。私はこの問題に一歩ずつ取り組むつもりはありません。特定のフィーチャが視覚的に明らかではない場合に、この大きなデータセットの交差に失敗したかどうかを確認することは（そして、小さなタイルでは不可能です）実用的ではありません。

11 qgis  shapefile  ogr  qgis-processing  vector-grid 

1x1 kmセルのフィッシュネットを受け取る陸上プロトコルがあります。一部のセルはランダムに選択されます。各セルに4つのポイントを配置する必要があり、これらのポイントも道路上になければなりません。ポイント間の最小距離は、すべての細胞のすべてのポイントのために500メートルでなければならない、可能な場合、またはそれがない場合は、我々は可能な最大距離を求めています。 最初の試行では、ST_CreateFishnetを使用してすべてのセルを4つの500x500 mセルに分割し、次にサブセルの重心にポイントを配置してから、最も近い道路（ST_ClosestPoint）に配置します。いくつかの良い結果が得られますが、以下の例では、ポイント5が6から近すぎるため、左側の道路に移動できることがわかります。 WITH r1 AS ( -- only sub-cells which intersects random cells SELECT id_maille, ROW_NUMBER() OVER() AS id_grille, fishnet_500.geomgrille FROM fishnet_500 JOIN t_mailles ON ST_Intersects(ST_Buffer(t_mailles.geom,-200), fishnet_500.geomgrille) -- buffer < 0 to not select neightbours ) , r2 AS ( -- cut roads in every cells SELECT id_maille, …

10 postgis  geoprocessing  qgis-processing  algorithm  vector-grid 

QGISを使用しています。 グリッドのサイズが緯度0.5度、経度1度のグリッドで構成されるシェープファイルを作成するにはどうすればよいですか？

10 qgis  vector-grid 

パラメーター「templateExtent」をシェープファイルで定義すると、ツール「pypy.CreateFishnet_management」を使用できません。パラメーター「originCoordinate」と「yAxisCoordinate」は自動的に入力されません。 import arcpy from arcpy import env env.overwriteOutput = True env.workspace = r"D:\Users\julia\erste_aufg" #Process: Create Fishnet outFeatureClass = r"D:\Users\julia\erste_aufg\at001l_wien\at001l_wien\wien.shp" cellSizeWidth = '200' cellSizeHeight = '200' templateExtent = r"D:\Users\julia\erste_aufg\at001l_wien\at001l_wien\at001l_wien.shp" arcpy.CreateFishnet_management(outFeatureClass, "", "", cellSizeWidth, cellSizeHeight, '0', '0', "", "NO_LABELS", templateExtent, "POLYGON") これはModelBuliderで機能しているため、ModelBuliderのバックグラウンドで何かが実行されており、「templateExtent」があるときに「originCoordinate」と「yAxisCoordinate」のパラメータを作成できます。パラメータ「templateExtent」のみを指定して、ArcPyでこのツールを実行するにはどうすればよいですか？ 誰かが解決策を持っていると本当に嬉しいです。スクリプトツールでフィッシュネットが必要ですが、最終的にループがあり、範囲の値が常に異なるため、解決策がないと解決できません。 スクリプト全体の最初の部分

9 arcpy  vector-grid 

次の座標があります minx, maxx, miny ,maxy = 448262.080078, 450360.750122, 6262492.020081, 6262938.950073 Pythonを使用してサイズ1 mの正方形グリッドを作成したいと思います。 import math minx,maxx,miny,maxy = 448262.080078, 450360.750122, 6262492.020081, 6262938.950073 size = 1 def set_bbox(minx, maxx, miny, maxy, distx, disty): nx = int(math.ceil(abs(maxx - minx)/distx)) ny = int(math.ceil(abs(maxy - miny)/disty)) new_maxx = minx + (nx*distx) new_miny = maxy - (ny*disty) …

9 python  polygon  vector-grid 

地図上のポイントとして多くの建物を含むレイヤーがあります。各建物の属性の1つは"download_speed"です。 地図の上に正方形のグリッド（100 x 100メートル）を配置したいと思います。正方形は次のように動作します。 広場は、広場に少なくとも1つの建物がある場合にのみ表示されます。 正方形内のどの建物も"download_speed"10（Mbit / s）を超えていない場合、正方形は赤になります。 正方形の一部の建物が"download_speed"10（Mbit / s）を超える場合、正方形は灰色になります。 正方形内のすべての建物が"download_speed"10（Mbit / s）を超える場合、正方形は黒である必要があります 私はQGIS（および一般にGISソフトウェア）の完全な初心者ですが、データサイエンスにおけるPythonの経験豊富なユーザーです。 最終結果は、次の図のようになります。

8 qgis  pyqgis  style  vector-grid  categorized-renderer 

DEMから抽出された高度の値を持つポイント（添付を参照）で構成されるグリッドがあります。これらの値をリストとしてではなく、SciLabで読み取るサイズ50x50のマトリックスとしてエクスポートする必要があります。グリッドはすでに50x50ポイントで構成されているため、サイズは実際には問題ありません。問題は、データをエクスポートするために見つけたすべてのメソッドが、探している行列の形ではなく、出力として値のリストを持っていることです。行列の原点は、添付の画像（強調表示）の最上部のポイントでなければなりません。この行列は、空間的に明示的な環境で水流速度を計算するモデルで使用されます。

8 arcgis-desktop  arcgis-10.0  point  export  vector-grid 

QGISで作成した次の10x10度の緯度/経度グリッドのシェープファイルがあります。rgdalパッケージを使用してRに読み込むことができます。 Grid<-readOGR(".","GridShapeFile") 次の属性と構造があります。 summary(Grid) Object of class SpatialPolygonsDataFrame Coordinates: min max x -180 190 y -100 90 Is projected: FALSE proj4string : [+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0] Data attributes: ID XMIN XMAX YMIN YMAX Min. : 0.0 Min. :-180 Min. :-170 Min. :-100 Min. :-90 1st Qu.:175.5 1st Qu.: -90 …

8 r  overlay  vector-grid 

私はこのプロジェクトに取り組んでおり、たとえば60度程度の特定の回転のグリッドが必要です。 それから始めて、20mx50mのブロックを作成する必要があります。最後に、各ブロック内にサブブロックの小さなセットが必要です。 結果はこれです： もちろん手作業で行っていますが、このように30ブロック、多分もっと多くの時間を費やして、自動で作る方法を考えていました。

8 arcgis-desktop  arcgis-10.2  vector-grid 

GCS_NAD83に255,000ポリゴンのアルバータ州（カナダ）セクションポリゴンシェープファイルがあります。私は各セクションを4つの等しい象限（NW、NE、SW、SE）に分割したい アルバータ州にはいくつかのUTMゾーンがあるため、セクションスキャンが投影単位でのみ分割される場合は、ランベルト正角円錐またはアルバース均等面積CSに変換できると思います。この情報を販売している人がいますが、これは私が持っているデータを使用して作成できるものかどうか疑問に思いましたか？私にお知らせください。ありがとうございました。

8 arcgis-desktop  polygon-creation  vector-grid  township-range-section