問題タブ [satellite-image]

場所の衛星画像を無料で入手できる画像を知っていますか?

高さマップを視覚化するために GIS データを使用する地形エディタがありますが、それらを衛星画像でテクスチャリングしたいと考えています。

いくつかの情報源を提供していただけないでしょうか。

人工衛星の構造を特定するために、衛星画像の分析アプリケーションの作業を開始する必要があります。これにはCまたはJavaを使用したいと思います。

衛星については、Googleマップのデータを使用する予定です。

ここに3つの質問があります：

1. グーグルマップ/地球以外のGISデータの最良の情報源は何ですか。
2. サードパーティのAPIを使用する必要があることを考慮して、このようなアプリケーションを作成するのに最適な言語
3. 人工構造物を識別するオープンな画像処理エンジンはありますか？

それはたくさんの質問ですが、ここの賢い人たちがここで私を助けてくれることを願っています。

特定のズーム レベルで衛星画像が利用できない場合、Google マップ API (javascript、v2) を使用してプログラムで検出することは可能ですか? いくつかのマップを自動的に作成し、それらを特定のズーム レベルに設定していますが、そのレベルで使用できる衛星画像がない場合がいくつかあります。その場合、ズーム レベルを自動的にバックアップしたいと考えています。API は、この画像の欠如を判断する何らかの方法を公開していますか?

これまでのところ、私が考えることができる唯一のことは、APIが行う画像リクエストが失敗しているかどうかをjavascriptで調べ、失敗した画像リクエストに基づいて反応することですが、実際にはそれを得ることができません仕事...そしてそれは少しエレガントではないようです。ご協力いただきありがとうございます！

バックグラウンド

私は合成開口レーダー衛星からの非常に大きなデータセットを扱っています。これらは、一辺が10kピクセル程度の高ダイナミックレンジのグレースケール画像と考えることができます。

最近、私は、SAR画像の線形特徴を検出するためのリンドバーグのスケールスペースリッジ検出アルゴリズム法のシングルスケールバリアントのアプリケーションを開発しています。これは、方向フィルターの使用またはハフ変換の使用の改善です。これは、どちらよりも計算コストが低いため、以前は両方とも使用されていた方法です。（4月のJURSE 2011で最近の結果をいくつか発表します。それが役に立ったら、プレプリントをアップロードできます）。

私が現在使用しているコードは、ピクセルごとに1つずつ、レコードの配列を生成します。各レコードは、ピクセルの右下にある長方形の尾根セグメントを記述し、隣接するピクセルで囲まれています。

のエントリには、、、の2つだけが、0〜128の範囲の値を持つ場合、ridgesリッジセグメントが含まれます。topleftrightbottom

次に、尾根セグメントの開始（x1、y1）と終了（x2、y2）を見つけることができます。

これらの個々の尾根セグメントがレンダリングされると、次のような画像が得られます（はるかに大きな画像の非常に小さなコーナー）。

これらの長い曲線はそれぞれ、一連の小さな尾根セグメントからレンダリングされます。

尾根セグメントを含む2つの隣接する場所が接続されているかどうかを判断するのは簡単です。ridge1at（x、y）とat ridge2（x + 1、y）がある場合、0 <= ridge1.right<= 128および ridge2.left=の場合、これらは同じ行の一部になりridge1.rightます。

問題

理想的には、すべての尾根セグメントを線につなぎ合わせて、画像で見つかった各線を繰り返し処理して、さらに計算を適用できるようにします。残念ながら、これを行うためのアルゴリズムを見つけるのは難しいと感じています。これは、複雑さが低く、メモリ効率が高く、マルチプロセッシングに適しています（非常に大きな画像を処理する場合のすべての重要な考慮事項です！）

私が検討したアプローチの1つは、リンクされた尾根セグメントが1つしかない尾根が見つかるまで画像をスキャンし、結果の線を歩いて、線内の尾根に訪問済みのフラグを立てることです。ただし、これはマルチプロセッシングには適していません。これは、高価なロックなしで、同じ行を他の方向から（たとえば）歩いている別のスレッドがないかどうかを判断する方法がないためです。

読者は可能なアプローチとして何を提案しますか？誰かが過去に効率的な方法を考え出したようなもののようです...

USDA Crop Data Layer (データはこちらから入手できます。大きなファイルに注意してください: http://www.nass.usda.gov/research/Cropland/SARS1a.htm ) をダウンロードし、分析のために R にインポートしたいと考えています。他の大きな GTIFF ファイル (LANDSAT データなど) をインポートしましたが、関数を呼び出すと次のエラーが発生します。

ここで、「ファイル」は USDA が提供する .tif ファイルのいずれかです (中西部地域のコンパイルを使用しています)。

問題は、GTiff が真の序数/スペクトル データではなく、カテゴリ データ (土地利用カテゴリ) であることだと思いますが、画像エディター (プレビューなど) で正常にレンダリングされ、R / GDAL でボークする理由がわかりません。 .

助けてくれてありがとう。

編集: これが GIS フォーラムにある場合は、私に知らせてください。一般的に、空間問題に対するプログラミングのアプローチはここにあるようです。

個別のヘッダー ファイルと共に、バンド インターリーブ ピクセル (BIP) 形式で保存されたマルチバンド衛星画像があります。ヘッダー ファイルは、イメージの行数と列数、バンド数 (標準の 3 より多い場合もあります) などの詳細を提供します。

画像自体は次のように保存されます (5 バンドの画像を想定)。

[B1][B2][B3][B4][B5][B1][B2][B3][B4][B5] ... など (基本的に 5 バイト - 各バンドに 1 つ - 各ピクセルの開始)画像の左上隅から）。

これらの各バンドを Python 3.2 (Windows 7 64 ビット) で PIL イメージとして分離する必要がありますが、現在、問題に間違って取り組んでいると思います。私の現在のコードは次のとおりです。

このコードは BIP ファイルを開くのに時間がかかりすぎます。これを行うためのより良い方法があるはずです。numpy および scipy ライブラリも持っていますが、それらをどのように使用できるか、またはそれらが何らかの形で役立つかどうかはわかりません。

画像内のバンド数も可変であるため、ファイルをすばやく読み取り、画像をコンポーネント バンドに分離する方法を見つけるのが難しいと感じています。

記録のために、ループ内でリストメソッドをいじってみました（スライスを使用する、スライスを使用しない、追加のみを使用する、拡張のみを使用するなど）が、主要な時間が失われるため、特に違いはありません。関連する反復回数 - (幅 * 高さ * numberOfBands)。

提案やアドバイスは本当に役に立ちます。ありがとう。

Pythonを使用して、衛星からの画像データを含むhdfファイル（Hierarchical Data Format）を読み取りました（名前はMODISである必要があります）。Python コードは、画像の経度、緯度、および放射輝度を取得できます。次に、データを使用して画像を生成します。

先生は、目的は画像を地図に載せることだと私に言いましたが、これの意味がわかりません。しかし、彼女は画像を地図に「配置」すると言いました。おそらく「近似最近傍」ライブラリ（ANNライブラリ）を使用する必要があります。

画像がマップに「配置」されると、次のようになります。

では、地図上に画像を「載せる」とはどういう意味でしょうか? ありがとうございました。

私は Mort Canty のhttp://mcanty.homepage.t-online.de/バイテンポラル RapidEye マルチスペクトル画像で Python iMAD 実装を実行しようとしています。これは基本的に、2 つの画像の正準相関を計算してから、それらを減算します。私が抱えている問題は、画像が 5000 x 5000 x 5 (バンド) ピクセルであることです。イメージ全体でこれを実行しようとすると、メモリ エラーが発生します。

pyTables のようなものを使用すると、これに役立ちますか?

Mort Canty のコードは、gdal を使用して画像をロードし、10 x 25,000,000 配列に格納しようとしています。

float の 10 x 25,000,000 numpy 配列を作成するだけでも、メモリ エラーがスローされます。これを回避する方法を知っている人はいますか？初めての投稿なので、投稿の仕方などアドバイスいただければ幸いです。

ご挨拶

MRT-MODIS Reprojection Tools を使用してモザイク MODIS MCD15A2 データをバッチ処理したいのですが、私のスクリプトは次のとおりです。

それは機能せず、「警告: CheckMosaicArgs : 一般処理: モザイク ツールの (MOSAICINPUT.TXT) からの input_filenames の処理中にエラーが発生しました」というメッセージが表示されました。

どうすればこれを修正できますか? 私はwindows7を使用しています。ありがとうございます。

重複の可能性:
Big Satellite Image Processing

バイテンポラル RapidEye Multispectral 画像で Mort Canty の Python iMAD 実装を実行しようとしています。これは基本的に、2 つの画像の正準相関を計算してから、それらを減算します。私が抱えている問題は、画像が 5000 x 5000 x 5 (バンド) ピクセルであることです。イメージ全体でこれを実行すると、コンピューターがひどくクラッシュし、オフにする必要があります。

Pythonがこのようにコンピューターをクラッシュさせる原因を知っている人はいますか? たとえば、バンドごとに 2999x2999 ピクセルを選択すると、すべて正常に動作します。

8 GB RAM、I7-2617M 1.5 1.5 GHz、Windows7 64 ビット。すべての64ビットバージョンを使用しています：python（2.7）、numpy、scipy、およびgdal。

前もって感謝します！

if name == ' main ': main()