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

自然災害モデルのために、現在非常に厄介なベクター/ Pythonプロセスを改善しようとしています。現時点では、特定のポイントからの距離/方位線を生成して決定する長いスクリプトがあります。 交差するポリゴンのタイプ（森、草、沼など） そのポリゴンまでの距離 これらの線のどれだけがポリゴンと交差し、どれだけ「囲まれている」かを決定します。 より多くの関与がありますが、それはそれの要点です。私はこれを改善する方法を見つけようとしていて、現在パート3で困惑しています。アイデアは、ポイントがポリゴンで完全に囲まれているかどうか、たとえば200m以内かどうかを判断することです そのため、添付画像では、ポリゴンで完全に囲まれているため、ポイントAをポイントBよりもリスクが高いとマーク付けします。これは約1300万ポイント繰り返されるため、簡単な作業ではありません。スクリプトを実行するのではなく、値を導出するためのサーフェスを用意します。これを行うには、水文学のツールやコストパスのバリエーションが必要だと思いますが、頭を悩ませることはできません。 これについてどうすればいいですか？

9 arcgis-10.0  arcgis-desktop  raster  arcpy 

標高のラスターイメージがあります。地方の政治区分を含むベクトルシェープファイルを追加しました。標高の色は保持したいが、政治区分を表す黒い線も保持したい。 政治区分のシェープファイルの透明度を下げましたが、これが政治区分の線に影響を与えないようにする方法がわかりません。黒の境界線の色を完全に保ちながら、塗りの透明度を下げるためのアドバイスはありますか

9 qgis  raster  vector  elevation  transparency 

大きなECW（詳細は以下）をクリップしようとしていますが、ラスターファイルが大きすぎて完全に処理できません。 以下のECWの詳細 ドライバー：ECW / ERDAS圧縮ウェーブレット（SDK 5.0） ファイルサイズ：50 GBサイズは450000、565081ピクセルサイズ：0.15 0.15 COLORSPACE = RGB COMPRESSION_RATE_TARGET = 9 VERSION = 2バンド数：4 クリップしたい領域は、元のファイルの約1/5です。 ここに私が成功せずに試した方法があります： Arcgisを使用してecwをtiff /他の形式に保存しました...（すぐにあきらめました） Qgisとそのクリッパーツールを使用しました...ファイルの作成は約40％で止まっています。 Qgis以外のオプションでOSGeo4Wのgdal_translateを使用しました。（Qgisを使用しないで一部のメモリを解放する可能性があるという考えがトリックになるだろうと考えてみました） 画像をバラバラに切り、必要なものを取得すると考えてgdal_retileを使用しました。コマンド「gdal_retile -ps 10000 10000 -of ecw -tileIndex tile.shp -targetDir input.ecwこれはさらに速くクラッシュしました」 誰かがアイデアを持っていますか？ 参考までに、16GBのRAMを搭載したi5-3470 3.2GhzでWindows 7 64ビットを実行します。

9 raster  gdal  clip  ecw  big-data 

小さな庭の区画の単位面積あたりの平均太陽強度を表すラスターサーフェスの作成に興味があります。モデルへの入力は次のようになると思います。 プロットエリア-簡単にするために、庭のプロットが円形であると仮定しましょう。 周囲の森の端の平均木の高さ さまざまな太陽角度での太陽強度 時間の経過に伴う庭園の区画に関連する平均太陽角度。 最終的には、木々に囲まれた庭の区画の中で、植物の成長に最適な場所を示したいと思います。誰かがこの分析の解決策に私を向けることができますか？私は参照を見つけることができませんが、誰かがすでにこれを行っていると確信しています。 編集： 例として、添付の画像は、同じ年齢の針葉樹の林に囲まれた1エーカーの庭の区画を示しています。この分析をさらに単純化するために、6月から10月までの園芸シーズン中の太陽の天頂と強度のみに関心があると仮定します。伐採はオプションではありません！

9 arcgis-desktop  python  arcpy  raster  r 

特定の値を持つラスターの最小境界範囲を見つける方法があるかどうか疑問に思っています。グローバルイメージからラスターをクリップしましたが、範囲は、多くのNoData領域を持つグローバル範囲として設定されています。このラスターからNoData領域を削除し、特定の値のピクセルを含む領域のみを保持したいと思います。これどうやってするの？ これが私の例です。value= 1（青い領域）を抽出し、以降の処理には全世界ではなく青い領域の範囲を使用します。

9 arcgis-desktop  raster  extents 

簡単なことのようですが、うまく機能しないようです。Linuxでcurlを使用して、REST APIを介してラスターデータファイル（GeoTIFF）をジオサーバーに追加しようとしています。REST API は、それを使用して新しいワークスペースを作成したため、機能しています。 curl -u admin:geoserver -v -XPOST -H 'Content-type: text/xml' \ -d '<workspace><name>restProba</name></workspace>' \ http://localhost:8080/geoserver/rest/workspaces Geoserverの管理UIに示されているように正しく作成されていますが、次のコマンドを実行してそのワークスペースにカバレッジストアを追加しようとすると curl -u admin:geoserver -v -XPOST -H 'Content-type: text/xml' -d '<coverageStore><name>int_dec</name><enabled>true</enabled> <type>GeoTIFF</type> <url>$home/int_dec.tif</url></coverageStore>' "http://localhost:8080/geoserver/rest/workspaces/restProba/coveragestores?configure=all" 正しく読み込まれていても、デフォルトのワークスペースに直接移動します。これは、希望するワークスペースではありません。希望するワークスペースを定義する方法はありますか？restProbaのワークスペースURIを指すURLで言及することで十分だと思いましたが、省略されているようです。 ご協力いただきありがとうございます ：-） pd：geoserverの初心者、私はこれを検索してみましたが、正しい検索基準を使用していないか、愚か者であるか、または私がそれを尋ねるのがだまされていると誰もが尋ねるのは簡単すぎます：p

9 geoserver  raster  rest 

次のコマンドを使用して、ラスタをpostgresデータベースにロードしています。 raster2pgsql -I -C -e -Y -F -s 3086 -I -C -M myraster.tif myraster -F -t 30x30 | psql -U postgres -d database -h localhost -p 5432 myraster.tif155Mです。QGISでラスタープラグインを使用しようとすると、ロードに非常に長い時間がかかります（15〜20分）。これをスピードアップするために私にできることがあるか、またはラスターを表示する別のアプローチがあるかどうかはわかりません。私ができることは、ラスターをすばやく表示して、データベースに格納されているベクターデータと正しく一致することを確認することだけです。（対照的に、ベクトルデータは非常に速く表示されます。）現時点では、QGISはフリーズし、ラスターを表示します。黒い四角形が表示されるだけで、相互作用によってプログラムがフリーズするため、カラースケールを変更することはできません。 私はラスターの読み込みに関連する微妙な点に精通していないので、おそらく使用する必要のあるさまざまなフラグと引数があるでしょう。これは大きなラスターであることを理解しているので、表示のためだけに再サンプリングしたり、同様に処理したりできます（どうすればよいかわかりません）。 wktrasterQGISでプラグインの最新バージョンを使用しています。リモートサーバーに接続してデータを表示しています。QGIS 1.8でOSXを使用しています。任意のプログラムを使用して、任意のOSでラスターを表示できます。

9 qgis  postgis  raster  postgresql 

ArcGIS DesktopにGeoTIFFラスターがあります。これは、次の2つのリンクからダウンロードできます。 http://www.rtwilson.com/downloads/MODIS_MOD04_16-06-2006_Subset.tfw http://www.rtwilson.com/downloads/MODIS_MOD04_16-06-2006_Subset.tif これらをArcGIS Desktopに読み込むと、大きな均一な領域が表示されます。情報ツールでこれらの領域をクリックすると、ピクセル値がデータなしと表示されます。ただし、レイヤーの[シンボル]設定で[データなし]の値を透明（または任意の色）として表示することを選択した場合、[データなし]の値はそのルールに従いません。 興味深いことに、画像のメイン領域外のデータなしの値はルールに従います（色を赤に設定してそれを確認します）が、画像内の値はそうではありません。 私はこの問題を解決するのに本当に苦労しました-そして、それを行う方法をまったく見つけることができません。 誰かアイデアはありますか？

9 arcgis-desktop  raster  arcmap  arcgis-10.0  spatial-analyst 

私のバックグラウンド：経験豊富なプログラマー、アマチュア歴史家、GISの完全な初心者（したがって、QGISを使用しています）。私はQGISジオリファレンスツールの使い方を学びました。 入力：このWikipediaのような画像があるとします：File：VoyagesOfRabbanBarSauma.jpg。それは地理的メタデータのない「裸の」JPGファイルであり、作成者はそのような情報を提供できないと仮定します。入力画像は、GNU FDLの下でライセンスされています。これは、派生物をリリースできることを嬉しく思います。 目標： Wikipediaマップに表示されたパスをQGISにインポートして、投影の変更、標高の重ね合わせ、物理的な地理のラベルの追加などを行うことができます。 注意点：回避策の1つは、マップ内の緯度/経度の場所を新しいGISプロジェクトに入力することです。この特定の作者による地理情報の推定の記録が欲しいので、このルートは使いたくありません。 私が試したこと：大陸のエッジの周りに一連のポイントを含むラスターJPGファイルにQGIS Georeferencerを適用しました。大陸内部の点の配置はまあまあです。 知っておきたいこと：この特定の地図のゆがみは、ユーラシアの北端付近でかなり目立ち、かなり標準的な投影のようです。このマップを生成した投影と投影パラメータを特定し、画像ピクセルと緯度/経度間の正確な非補間変換を復元するアルゴリズム的な方法があるはずだと私には思われますか？ これについて2つの質問が表示されます。 不明なときにシェープファイルの座標系を識別しますか？ 古い図面の不明な座標系 これらは、GDALまたはQGISを含むブルートフォース方式について話しているようですが、GIS初心者にとっては、やや手間がかかるようです（基本的に地理参照で、既知の投影法に適合させます）。それらを自動化します。 これらのGIS.SEの両方での説明は、私の経験で許可されているよりも少し高く、このブルートフォースメソッドを正確に実行する方法についてのヒントは非常に高く評価されます。たとえば、上記の質問2は、gdaltransformを使用して座標のリストを緯度/経度に変換する方法を示しています。入力座標はxyピクセルの位置ですか？

9 qgis  raster  coordinate-system  georeferencing 

ブール分類（1または2）の2つのラスターイメージがあります。それらは特定の領域で重なり合っています。 重なり合う部分のみのラスターが必要です。 Arcinfo 10の使用。Answerはラスターメソッドのみを使用でき、ポリゴン、ポイントなどへの変換はできません。

9 arcgis-10.0  raster 

私はpython gdalライブラリを使用してラスターを作成しようとしていますが、データが出力されているポイントに到達しましたが、出力データは原点のx軸で反転しています。私は何かを見落としている必要があることを知っていますが、どこが間違っているのか理解できません。何か案は？ ラスターを作成するときに、左上のx / y値を設定すると、配列には左上からインデックスが付けられ、右下に向かってインデックスが付けられているように見えます。以下のコードでは、配列に行の値を入れています。 配列を出力すると、次のようになります。 [[ 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. …

9 python  raster  gdal 

前の質問を参照して、Rを使用して、生成されたラスターイメージをマップ上のレイヤーとして表示できる形式に「ジオコード」/変換するにはどうすればよいですか。要するに、私は商用のSpatialKey（スクリーンショット）ソフトウェアの結果を再現することを目指しています。 地図のデータをホストするためにGeoServerインスタンスを使用していますが、必要な投影と座標情報を画像に組み込んで、実際の地理的位置に正しく位置合わせする方法を教えてください。密度ヒートマップは、空間データを使用して作成されました。

9 raster  r  georeferencing 

（等高線ツールからの）ポリラインがオーバーレイされたクリップされたラスターデータセットがあります（添付を参照）。ポリラインをラスターデータセットの範囲までクリップしようとしています。これは非常に単純なように見えますが、複数のことを試しましたが、機能しないようです。最終的には、arcpyを使用してこれをコーディングする必要がありますが、今のところ、どのツールがこれを達成するかについて正しい方向に少し進むと、今のところ役立ちます。 ありがとう。 PS ArcInfo 10.0を使用しています

9 arcgis-10.0  raster  arcpy  vector  clip 

イタリア北部の土地利用に関するラスターファイルがあります（作物/都市部など）。 ベクターレイヤーの特定のポイントの周りにいくつかのバッファーゾーンを定義しました。選択した領域の平均使用量を推定するために、ラスターファイルの各値が占めるこれらの領域の割合を計算したいと思います（ 70％作物、30％都市など）。 ArcGISには、「Tabulate Area」と呼ばれるSpatial Analyst関数があることを知っていますが、QGISでそのような機能を見つけることができません。 誰かが私を助けることができれば私はとても感謝しています。

9 qgis  raster 

山岳地帯のDEMがあります。合計DEM領域の60％を超えるDEMの領域を抽出したいと思います。だから、本質的には山の上部を表す領域を抽出したいと思います。 私はこれをhypometric線（H60線）と呼んでいると聞いたことがありますが、それを達成する方法に関するドキュメントを見つけることができませんでした。 私の本能は、何とかDEMの統計を計算して60％レベルがどこにあるかを教えてから、この標高でDEMの輪郭を描くことです。次に、この輪郭によってDEMをクリップし、出力DEMサブセットに対して計算を実行できます。 この60％のレベルはどのように決定しますか？

9 arcgis-desktop  raster  arcgis-10.0  dem