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

ラスターは、通常は画像のような形式で格納される、値の規則的なグリッドで構成されるデータ形式です。

5
DEMの色を均一にする方法は?
これはばかげた質問ですが、私がしばらくの間理解しようとしている質問です。私は2つのAster DEMタイルを並べています。両方の色を均一にするにはどうすればよいですか?
8 qgis  raster  dem  color 

2
Pythonを使用してラスターの範囲を取得する方法
保存して、抽出された範囲で別のラスターをクリップするためにそれを使用できるように、特定のラスターの範囲を取得する必要があります。今のところ、私は範囲の部分で立ち往生しています。 import arcpy from arcpy import env desc= arcpy.env.extent("D:\GIS @ UTD\Sem 1\GISC 6317\Lab\Lab 10\temp\Tahoe\Tahoe\Emer\erelev.grid") print desc 以下は、出力を示します0 0 0 0 NaN NaN NaN NaN。 助言がありますか。
8 raster  arcpy  clip  extents 

3
マスク出力でラスターを抽出すると、マスクが完全にオーバーラップしませんか?
流域集積グリッドと集水域を含むシェープファイルがあります。次に、それらに対していくつかの計算を行いたいので、特定の集水域に隣接するラスター内のセルからデータを抽出する必要があります。 そのプロセスでは、集水域ポリゴンの1つをマスクとして使用して、ラスターで[マスクによる抽出]を使用しようとします(後でループしますが、今は1回だけ作業する必要があります)。 問題は、元のラスターとマスキングポリゴンの両方と比較して、抽出されたラスターがセルの半分だけ左にシフトすることです。ラスターをポイントに変換し、後で他のポリゴンと交差させようとするので、これは私のプロセスの残りを台無しにしています。そして、抽出されたラスターが元のラスターと重なっていない限り、正しいポリゴンが得られません。 ポリゴンはラスターから作成されるため、元のラスターとポリゴンは完全に整列しています。どういうわけか投影図になんらかの変化があったかどうかを確認しようとしましたが、見つけることができません。また、新しいArcMapドキュメントを試し、すべてを数回再起動し(これにより、驚くほど頻繁に問題が解決します)、スクリプトを呼び出すのではなく、Pythonプロンプトからツールを実行しました。ツールボックスからマスクによる抽出ツールを実行しようとしても、同じ結果になりました。 これは、ArcMapのpythonプロンプトで使用する(非常に単純な)コードです。 inpgs = r'C:\GIS\Catchm_ex.shp' flowacc = r'C:\GIS\flowacc' arcpy.MakeFeatureLayer_management(inpgs, 'currentMask', '"HydroID" = 269634') outRaster = arcpy.sa.ExtractByMask(flowacc, 'currentMask') ここにサンプルファイルがあります Windows 7、ArcInfo 10.0、Python 2.6を使用しています。 私は今、Extract by Polygonを試しましたが、それはうまくいくようです。ExtractByMaskでマシンが誤動作しているだけですか?唯一の欠点は、マスキングポリゴンの座標を取得する方法がわからないため、スクリプトでその関数を使用するときに問題が発生することです(9つのセルのみをカバーし、中心にポイントがある場合)。 次のコードを使用して、ポリゴンによる抽出で動作するようになりました(Esriのサポートページから変更)。残念ながら、これは私の問題を解決しません。出力(このツールとマスクによる抽出)がセルの半分だけ左にシフトする場合とシフトしない場合があります。ほとんどの場合はそうですが、実際には本来あるべき場所にあることもあります。 だから今、私は本当に混乱しています。これは私のコンピューターで起こっているのですか?他の誰かが同じことを経験している場合は、コメントを書き込んでください。私だけではないことを私に知らせてください:)私にとって、この問題はバグとしては大きすぎる問題であり、これが問題となっています途中で自分自身を作成し​​たこと。 それを避ける方法に関するアイデアはありますか? def ReadPoints(infc): try: import arcpy, numpy except ImportError: sys.exit('Unable to import arcpy') # Identify the geometry field desc = …

1
ArcPyを使用してcon()スクリプトを使用して新しいラスターを作成しますか?
私はPythonが初めてで、ラスター計算機スタイルのスクリプトをコーディングしようとしています。基本的に、値が0〜10,000のラスターを、特定の値の間に1が設定されたラスターに変換します。 私がPythonでそれをしたい理由は、最終的にこのプロセスを自動化して、さまざまな値(つまり、0から-10、-5から-25など)を選択し、最終的にはそれぞれ1のラスターになるようにするためです。 。私の目的は、これらをポリゴンに変換して、約800のシェープファイルで多数のベクタークリップ操作を実行できるようにすることです。 最後に新しいラスターを生成するスクリプトを記述しましたが、このグリッドが1でいっぱいではなく、値が0〜256であるため、コードのどこかにバグがあります。どなたでも見つけて頂ければありがたいです!誰かがコードを高速化する方法、または私がこれをすべて行う方法についてのより良いアイデアを持っているなら、それもまた興味深いことです。 これが私のコードです: import arcpy from arcpy import env from arcpy.sa import * arcpy.env.overwriteOutput = 1 #this also doesn't work by the way, I have also tried #TRUE, no errors come up but it still doesn't let me overwrite the output file #Set environment settings env.workspace = "C:/folderlocation …

1
シェープファイルを投影してから、シェープファイルをラスタライズ&ラスタライズしてから投影するのに違いはありますか?
現在、シェープファイルを単純な存在/不在(1&0)ラスタにラスタライズし、各ラスタを積算合計に追加して、最後に合計リッチネスラスタを作成するプロジェクトがあります。シェープファイルはWGS84投影にありますが、最終的なラスターはBehrmann Equal Area投影にある必要があります。 私が知りたいのは、すべてのシェープファイルをベールマンに投影し、それらをそれぞれラスタライズおよびマージして、最後にラスター化およびマージしてから最終ラスターをベーアマンに投影した場合に、最終ラスターに違いがあるかどうかです。


3
raster2pgsqlを使用してPostGISにジオティフをロードする
いくつかのOrdnance SurveyラスターデータをPostGISに読み込むときに問題が発生します。次のraster2pgsqlコマンドを使用して、C:\ Program Files \ PostgreSQL \ 9.1 \ binフォルダに配置したtifをPostGISにロードしています。 raster2pgsql -s 27700 -t 5000x5000 -I -C -M *.tif > os_strview.sql 次に、次のpsqlコマンドを実行します。 psql -h localhost -U postgres -d raster_test -f os_strview.sql これは正常に機能しますが、QGISでデータを表示しようとすると、PostGISレイヤーを追加しようとしたときに、選択ボックスにraster_columnsテーブルしか表示されませんでした。 代わりに、プラグインLoad Post raster to PostGISを使用してみました。これは機能しますが、ラスターイメージに色がありませんでした。その後、この投稿を見ました。「ラスター」>「変換」>「PCT to RGB」テクニックを使用して、単一のtifをgeotifに変換しました。 残念ながら、上記のpsqlコマンドを実行すると、次のエラーメッセージが表示されます。 C:\ Program Files \ PostgreSQL \ 9.1 \ bin> psql …

1
ArcPyまたはArcGIS Desktopでのラスター処理の範囲を削減しますか?
Pythonで何らかの自動手順を実行して、ラスターの処理範囲を減らす必要があります。現在、私が使用しているデータセットの境界に異常があります(ランドサット走査線の分類のため)。 上記の情報をSetNull関数で取得したときに境界線が処理されないように、ラスター処理の範囲を縮小したいと思います。何かのようなもの これは、マスク領域の手動での選択を必要としますが、マスク関数による抽出(つまり、figure2)で簡単に実行できます。エクステント領域(nodata値を無視)から各行/列のx方向とy方向の境界幅を引いたものを計算arcpy.env.extentしてPythonで利用する方法はありますか?

1
PostGIS out-dbラスターの使用に制限はありますか?
大きなラスターをPostGIS 2.0データベースにインポートして、いくつかの障害にぶつかっています。 最初にraster2pgsqlを使用してラスターをインポートしようとすると、次のエラーが表示されます 。rt_band_set_pixel_line:座標が範囲外です 見つけたら、このページを、そしてgdalinfoを使用して、私は、ラスタが65535 X 65535の最大許容サイズを超えたことがわかりました。 gdalinfoの出力: $ gdalinfo ari100.tif ドライバー:GTiff / GeoTIFF ...サイズは42971、77138 座標系は: GEOGCS ["WGS 84"、 ... AUTHORITY ["EPSG"、 "4326"]] Origin =(152.741676637167842、-26.215302802012008 ) ピクセルサイズ=(0.000009239757419、-0.000009239757419) メタデータ: AREA_OR_POINT =エリア イメージ構造メタデータ: INTERLEAVE = BAND コーナー座標: 左上(152.7416766、-26.2153028)(152d44'30.04 "E、26d12'55.09" S) 左下(152.7416766 、-26.9280392)(152d44'30.04 "E、26d55'40.94" S) 右上(153.1387183、-26.2153028)(153d 8'19.39 "E、26d12'55。09 "S) 右下(153.1387183、-26.9280392)(153d 8'19.39 "E、26d55'40.94" S) …

1
カラー画像の表示を改善するには?
ラスターを表現するためのいくつかの本当に素晴らしいカラーランプを見つけましたが、ArcGISとQGISでは、複雑さが原因で滑らかなグラデーションではなく明白なステップが発生するため、完成品には簡単に使用できないようです。ラスターレンダラーは、24ビット標準RGBまたはディープ48ビットRGBを含む数学的な変換ではなく、8ビットパレット(256レベルのみ)で動かなくなっているように見えます。 おそらくラスターを複数のレイヤーに分割して個別にシンボル表示するという厄介な回避策がありますが、このようなストレッチ/カラー化されたグラデーションをより高い精度のフォーマットでネイティブに表示できるデスクトップおよびWeb GISシステムを知りたいのですが。

4
ArcGIS Desktopを使用して選択したラスターグラフィックをエクスポートしますか?
ArcGIS 10を使用してラスターグラフィックスをエクスポートする際に問題が発生しました。非常に単純なように見えますが、何らかの理由でArcMapを使用できません。これが私のステップです: ArcMapにラスターを追加する ラスターの属性テーブルを開き、いくつかの値を選択します ArcMapのコンテンツウィンドウでラスターを右クリック->データ->データのエクスポート データ範囲には、(a)データフレーム、(b)ラスターデータセット、(c)選択したグラフィックスの3つのオプションがあります。最後のオプション(必要なオプション)はグレー表示されています。 [ラスターデータのエクスポート]ダイアログで[選択したグラフィック]オプションをアクティブにできませんでした。奇妙なことに、私の同僚のArcMapのバージョンはまったく同じで、上記の手順を繰り返すと、選択したグラフィックスをエクスポートするオプションがグレー表示されません。調整が必要な設定があると想定しています。誰かがこれを引き起こしている可能性があることを知っていますか?

1
ファイルジオデータベース間で大きなラスターデータセットをすばやくコピーする方法はありますか?
大きなラスターデータセットを1つのファイルジオデータベースから別のジオデータベースにコピーするのに、なぜそれほど時間がかかるのですか?Windowsエクスプローラーでデータベース全体を数分でコピーできますが、ArcCatalogで単一のラスターをコピーすると時間がかかります。Windowsエクスプローラーを使用してファイルジオデータベース間でデータをコピーすることはお勧めできないと常に思っていましたが、これを行う安全な方法はありますか?もしそうなら、どうですか?


2
数十億のポリゴン(バウンディングボックス)の効率的なラスターサンプリング
何十億ものバウンディングボックス(ファイルから順番に読み込まれる)で構成されるセットが与えられ、各セルのラスター値が重複するバウンディングボックスの数を与える必要がある場合、ラスターを効率的に(Pythonで)計算するにはどうすればよいですか? 4000 * 4000ラスターの場合 私はnumpy行列の作成に時間をかけました: $ python -m timeit 'import numpy' 'a = numpy.zeros(shape=(4000,4000))' 10 loops, best of 3: 51.7 msec per loop 標準のPython行列の作成: $ python -m timeit 'a = 4000*[0]' 'for i in range(4000):' ' a[i]=4000*[0]' 10 loops, best of 3: 218 msec per loop したがって、numpyはより高速ですが、ループあたり50ミリ秒、10億回の反復により、実行時間は約1年(0.05ミリ秒* 1000000000/60/60/24/365 = 1.5年)になります。 …

2
オープンソースのGISを使用してラスター時系列をアニメーション化する方法は?
約10分間隔で取得した一連のラスターがあり、それらのかなり十分なマップを生成する素晴らしいQGISプロジェクトがあります。オープンソースのGISを使用して、このラスターシリーズをカラーマップで(静的)ベクターデータでオーバーレイしてアニメーション化する方法があるかどうか疑問に思っています。私はすでにプロジェクトをセットアップしているので、QGISソリューションが理想的ですが、他の提案も歓迎します。

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