タグ付けされた質問 「clip」

あるデータセットの範囲を別のデータセットの範囲または境界で制限または削減すること。


2
Rでラスターをクリップする
米国北東部の地図を作成しています。地図の背景は、高度地図または平均年間気温地図である必要があります。Worldclim.orgからこれらの変数を提供する2つのラスターがありますが、関心のある州の範囲でそれらをクリップする必要があります。これを行う方法に関する提案。これは私がこれまでに持っているものです: #load libraries library (sp) library (rgdal) library (raster) library (maps) library (mapproj) #load data state<- data (stateMapEnv) elevation<-raster("alt.bil") meantemp<-raster ("bio_1.asc") #build the raw map nestates<- c("maine", "vermont", "massachusetts", "new hampshire" ,"connecticut", "rhode island","new york","pennsylvania", "new jersey", "maryland", "delaware", "virginia", "west virginia") map(database="state", regions = nestates, interior=T, lwd=2) map.axes() #add …
33 raster  r  clip  maps 

5
Rの多数のポリゴンによる切り抜き、マスク、およびラスターの抽出の速度を上げますか?
数千のポリゴン境界に基づいて、さまざまな土地利用タイプの面積と被覆率をラスタから抽出しています。個々のポリゴンを繰り返し処理し、特定のポリゴンのサイズまでラスターをマスクして切り取ると、抽出機能がはるかに高速に動作することがわかりました。それにもかかわらず、それはかなり遅いです、そして、誰かが私のコードの効率と速度を改善するための提案を持っているかどうか疑問に思っています。 私はこれに関連した見つけた唯一のものです。この応答使うことを提案ロジャーBivandによるGDAL.open()とGDAL.close()だけでなく、getRasterTable()とgetRasterData()。私はそれらを調べましたが、過去にgdalに問題があり、それを実装する方法を知るのに十分なほどよく知りません。 再現可能な例: library(maptools) ## For wrld_simpl library(raster) ## Example SpatialPolygonsDataFrame data(wrld_simpl) #polygon of world countries bound <- wrld_simpl[1:25,] #name it this to subset to 25 countries and because my loop is set up with that variable ## Example RasterLayer c <- raster(nrow=2e3, ncol=2e3, crs=proj4string(wrld_simpl), xmn=-180, xmx=180, ymn=-90, ymx=90) c[] …


2
Pythonでシェープファイルを読み取る方法は?
私の質問は、ポリゴンシェープファイルの垂直線の延長です。最初にその質問を参照してください。 表示されるのは、ユーザー定義の間隔で、バウンディングボックスに対して垂直線を生成する方法です。OGR、Fiona、Shapelyなどを使用して次のクリッピングのステップを実行できることは理解していますが、それらの使用法は理解していません。 ポリゴンシェープファイルの1行を読み取るにはどうすればよいですか?Shapelyを使用するすべてのアプリケーションは、LineString、Point、またはPolygonを生成する方法を示しますが、既存のシェープファイルを読み取ることはありません 少なくともスケルトン構造をサポートしてくれますので、その上に構築できます。

4
ラスターファイルのクリッピングエラー:カットラインの境界ボックスを計算できません
ラスターファイルがあり、それをクリップしたい場合、ベクターファイルがマスクレイヤーになるはずです。私は最も明白な方法を使用します(ラスター->抽出->クリッパー)。私が何をしても、このエラーが発生します(下図を参照)。 カットラインの境界ボックスを計算できません ラスタ形式の変更、投影の変更などで試してみました。もちろん、古い投稿/ Googleでも調査しました。この問題を解決する方法がわかりません。拡張してクリップすると、機能します。 奇妙なことに、3か月前に同じ方法を実行したことを覚えていますが、このエラーは発生しませんでした。 また、(このエラーの解決策の代わりに)このアクションを実行する別の方法を歓迎します。説明どおりにラスターファイルをクリップするだけです。誰か助けてもらえますか?
22 qgis  raster  clip 

4
Rのシンプルフィーチャオブジェクトを切り取る
maptools::pruneMap(lines, xlim= c(4, 10), ylim= c(10, 15))SpatialPolygonまたはSpatialLine に使用されるような、sfマップオブジェクトを切り取る機能はありますか? 私は検討してst_intersection()いますが、適切な方法があります。
20 r  clip  sf 

2
QGISを使用してラスターレイヤーをトリミングしますか?
QGISでラスターレイヤーをトリミングするにはどうすればよいですか? レイヤーのセクションから輪郭を作成したいのですが、レイヤー全体を含めると、必要以上のデータが作成され、レンダリングが遅くなります。 ラスター計算機を使用して新しいレイヤー範囲を指定することでそれができることがわかりましたが、行と列に入力する値を知るにはどうすればよいですか? Photoshopでジオティフのトリミングも試みましたが、座標データが失われたようです。 最初にトリミングせずに、ラスターレイヤーの一部から等高線を作成する方法はありますか?

2
Rのポリゴンレイヤーを使用して、ポリゴンレイヤーの真のGISクリップを実行する方法は?
一連の単一境界ポリゴンを使用して土壌ポリゴンのRで真のGISクリップを実行したいのですが、適切に実行するR関数が見つかりません。clipESRIのArcMapの関数と同じように機能するはずです。パッケージover内のメソッドを試しましたが、sppolys over polysでは機能しないようです。 1つの提案は、次のコードを使用して、gIntersectionin rgeosパッケージをクリップとして使用することです。 #------------------------------------ library(rgeos) library(maptools) #Read layers as SpatialPolygonsDataFrame (both the same Albers projection) Soils_poly = readShapePoly("Soils_polygons") #Note - Has 400 polygons clipper_poly = readShapePoly("clipper_polygon") #Note - Has 1 polygon #Try gintersection as clip Clipped_polys = gIntersection(Clipper_Tile_poly, Soils_poly) #----------------------------------- これには、実行に5分かかり(速度が遅すぎる)、次のエラーが発生します。 RGEOSBinTopoFunc(spgeom1、spgeom2、byid、id、drop_not_poly、 "rgeos_intersection")のエラー:TopologyException:-721459.77681285271 2009506.5980877089で発信dirEdgeが見つかりません また、重複をチェックするためにこのコードを試しました。 gIntersects(Clipper_Tile_poly, Soils_poly) そして結果は真でした。clipESRI …
16 r  polygon  layers  vector  clip 

3
PyQGISでのGIS操作の並列化?
GISの一般的な要件は、処理ツールを多数のファイルに適用するか、1つのファイルの多数のフィーチャのプロセスを別のファイルに適用することです。 これらの演算の多くは、計算結果がループ内の他の演算に影響を与えないという点で、恥ずかしいほど並行しています。それだけでなく、多くの場合、入力ファイルはそれぞれ異なります。 典型的な例は、シェイプファイルを、ポリゴンを含むファイルに対してクリップします。 これは、QGISのPythonスクリプトでこれを達成するための(テスト済みの)古典的な手続き型の方法です。(テストファイルを処理する時間を半分にしたよりも、実際のファイルへの一時メモリファイルの出力をfyi) import processing import os input_file="/path/to/input_file.shp" clip_polygons_file="/path/to/polygon_file.shp" output_folder="/tmp/test/" input_layer = QgsVectorLayer(input_file, "input file", "ogr") QgsMapLayerRegistry.instance().addMapLayer(input_layer) tile_layer = QgsVectorLayer(clip_polygons_file, "clip_polys", "ogr") QgsMapLayerRegistry.instance().addMapLayer(tile_layer) tile_layer_dp=input_layer.dataProvider() EPSG_code=int(tile_layer_dp.crs().authid().split(":")[1]) tile_no=0 clipping_polygons = tile_layer.getFeatures() for clipping_polygon in clipping_polygons: print "Tile no: "+str(tile_no) tile_no+=1 geom = clipping_polygon.geometry() clip_layer=QgsVectorLayer("Polygon?crs=epsg:"+str(EPSG_code)+\ "&field=id:integer&index=yes","clip_polygon", "memory") clip_layer_dp = clip_layer.dataProvider() clip_layer.startEditing() clip_layer_feature …


4
FMEでは、リーダーから複数の入力フィーチャタイプを接続し、単一のトランスフォーマーにルーティングする方法はありますか?
FMEでは、複数のリーダーを単一のトランスフォーマーに接続する方法はありますか?私はFME 2015を使用しており、多数のフィーチャクラスを1つのフィーチャにクリップしたいのですが、これを行う簡単な方法が見つからないようです。 編集:入力データセットから複数のフィーチャタイプを選択し、それらをすべてフィーチャタイプの個々の出力矢印をクリックしてトランスフォーマの入力矢印にリンクすることなく単一のトランスフォーマにリンクしようとします。 編集2:これはスクリーンショットです、うまくいけばこれが役立つでしょう:

3
Rの逆クリッピング(消去)?
他のオブジェクトの内側にある 部分を保存する通常のクリップとは対照的に、リバースクリップは、別のオブジェクトの境界の外側にある空間オブジェクトの部分のみを保存します。 ArcMapで逆クリップを実行していますか?ArcMapでそれを行う方法を示します。 Rでこれを行うにはどうすればよいですか? 再現可能な例(Linuxマシンの場合): system("wget 'https://github.com/Robinlovelace/Creating-maps-in-R/archive/master.zip' -P /tmp/") unzip("/tmp/master.zip", exdir = "/tmp/master") uk <- readOGR("/tmp/master/Creating-maps-in-R-master/data/", "ukbord") lnd <- readOGR("/tmp/master/Creating-maps-in-R-master/data/", "LondonBoroughs") plot(uk) plot(lnd, add = T, col = "black") ここでやりたいことは、ロンドンを除くすべての英国を救うことです。視覚的には、結果の画像の黒い形状が穴になるようにします。
14 r  clip  erase 

3
Gdal:別のラスターでラスターをクリップ
私は、マルチバンドジオティフラスターファイルのバッチを同じ(より小さい)領域にトリミングする簡単なユーティリティを書いています。gdalwarpを使用すると、単一ポリゴンのクリッピングシェープファイルを使用してファイルを簡単にトリミングできます。 gdalwarp -cutline clipper.shp -crop_to_cutline input.tif output.tif ただし、クリップする実際の領域は、最初は常にシェイプファイルではなく、別のジオティフラスターファイルによって定義されます。そのラスターの範囲をクリッピングファイルとして使用できればいいのですが、どうすればよいかわかりません。当然のことながら、以下は機能しません(エラーは発生せず、何も生成されません)。 gdalwarp -cutline clipper.tif-crop_to_cutline input.tif output.tif だから、私の質問は、ラスターを提供する方法はありgdalwarp -cutlineますか?あるいは、別のラスターを使用してラスターをクリップできる別のgdal関数がありますか?これらのどちらも不可能な場合、ラスタの範囲によって定義された単一のポリゴンでシェープファイルを作成する非常に簡単な方法はありますか? このコードは、より広範なpythonスクリプトにラップされるため、コマンドラインのgdalユーティリティまたはgdalのPythonバインディングを使用できます。 補足として、QGISのラスターの範囲をカバーするクリッピングシェープファイルを簡単に作成できることを知っています。簡単な解決策が見つからない場合はそれを行うかもしれませんが、最終的には大規模な自動分析の一部として数百ではないにしても数十の領域でこのユーティリティを使用することになりますので、面倒なことはしたくないと思います非常に簡単な場合でも手動ステップ。
14 python  raster  gdal  clip  gdalwarp 

2
QGISでポリゴンに基づいてラインレイヤーをクリッピングしますか?
ラインレイヤーとポリゴンレイヤーがあります。両方のファイルに同じSRSがあります。新しいシェイプファイルには、ポリゴンで覆われた線のみが必要です。 QGISでVector -> Geoprocessing -> Clipツールを試しましたが、空のシェープファイルしか取得できません。 QGISでこれを行うにはどうすればよいですか?
13 qgis  clip 

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