1
GDALを使用したナノスケールDEMの作成
少し奇妙な質問かもしれませんが、実際の質問の前に、背景について簡単に説明しましょう。 原子間力顕微鏡（AFM）は、要するに（そして、私の限られた知識で）研究者がマイクロおよびナノスケールで領域をスキャンできる方法です。ある種のプローブを使用して領域を「スキャン」することで機能します。それ以上のことを説明するのは難しいです。私はそれを本当に理解していないからです。私が知っていること、そして私の好奇心を引き起こしたのは、結果が実際には「高さ」値の「グリッド」であるということでした（そのポイントでのプローブの高さを表す512x512の値のマトリックス）。 それから私は考えました：さて、スケールとは別に、これは実際にはデジタル標高モデルです！また、これは、GISツールで理解できるようにDEMファイルを作成できれば、GIS分析を適用できることを意味します。 現状では、AFMマシンを備えたラボで私の重要な仕事をしており、彼女のプロジェクトの1つでそれを使用しています。私は彼女からいくつかのスキャンファイルを取得し、Python（structおよびnumpy）を使用してこれらのバイナリファイルを解析しました。 次に計画しているのは、「適切なDEMへのマッピング」パートです。私はDEMSについてある程度の知識を持っていますが、実際にそれらを生成することになると、私はかなり新しいです。 私が考えているのは、何らかの方法でデータをジオリファレンスする必要があるということです。そのためには、カスタム（平面）座標系が必要です。私の座標系では、マイクロメートルまたはナノメートルを単位として使用することを想定しています。次に、AFMでスキャンされた領域のサイズを見つけるだけです（これはバイナリファイルのどこかにあると考えられます。これは既知であると仮定してください）。 更新：解像度は異なるが同じ領域のスキャンもいくつか行っています。たとえば、2つのスキャンに関する次の情報があります。 大きな画像： Scan Size: 51443.5 nm X Offset: 0 nm Y Offset: 0 nm 小さい（詳細）画像： Scan Size: 5907.44 nm X Offset: 8776.47 nm Y Offset: 1486.78 nm 私が考えているのは、カスタム座標系の原点は0,0である必要があり、大きな画像ではピクセル0,0に座標値（0,0）を割り当て、ピクセル512,512に座標値（51443.5、51443.5）を割り当てることです）（必要な他のポイントの写真を取得すると思います）。 次に、大きい画像はピクセル（0,0）を（8776.47、1486.78）に、（512,512）を（8776.47 + 5907.44、1486.78 + 5907.44）にマッピングします 1番目の質問は次のとおりです。このような座標系のproj4 defを作成するにはどうすればよいですか？すなわち、これらの「現実世界の座標」をカスタム座標系に割り当てるには（または、whubersの提案に従い、ローカル座標系を使用してユニットについて横になっている場合（つまり、ナノメートルをキロメートルとして扱う） 次に、numpy 2次元配列をGeoreferenced DEMファイル形式に転送する必要があります。私はGDAL（または、むしろPythonバインディング）の使用を考えていました。 2番目の質問は次のとおりです。私のような「任意の」データからジオリファレンスされたDEMを作成するにはどうすればよいですか。できればPythonで、オープンソースライブラリを使用してください。 そうすれば、適切な分析ツールを使用するだけで、残りはかなり簡単になります。問題は、このタスクは私自身の好奇心によって駆動されるため、ナノスケールDEMで実際に何をするべきかはよくわかりません。これは頼みます 3番目の質問：ナノスケールDEMをどうするか？どのような分析を行うことができ、DEM分析に適したツールは何ですか？最後に、このデータから陰影起伏と等高線を使用してマップを作成することは可能ですか？ 私はすべての提案とポインタを歓迎しますが、予算や資金がない厳密な趣味ベースのプロジェクトであり、ライセンスされたGISアプリケーションにアクセスできないため、無料の代替案を探していることに注意してください。さらに、これらのAFMマシンを販売しているBrukerがソフトウェアを出荷していることは知っていますが、それを使用するのは面白くないでしょう。