[OSGeoJapan-discuss] QGISにて横断ラインを作成する方法をご教示ください

[Minoru Akagi]

久保さま

MemoryErrorというメッセージから原因はメモリ不足ではないかと思います。対処方法としては，
・GdalToolsのClipperで必要な範囲のみを切り取ってサイズを小さくしておく。
・NoData Valueの処理部分でメモリ使用量が増えるのでNoData ValueなしのGeoTIFFにしておく。gdalwarp in.tif
out.tif とすれば，NoDataセルの値が0となり，NoData Valueが設定されないGeoTIFFが出力されます。

以上，ご参考までに。

赤木 実


2013年4月19日 10:55 Mayumi Kubo <mayumitt ＠ gmail.com>:

> 寒地土木研究所の久保と申します。
> 先日の宮園さんご質問の横断ライン作成方法について便乗させて下さい。
>
> 赤木さんからご提案のあったqProfプラグインをインストールして
> DEMからの横断ライン作成を試みました。
> QGIS起動→必要なレイヤ（DEMとラインシェープ）を読込→
> qProfでXYZ座標のcsvとポイントシェープファイルの出力成功したのですが、
> 二回目を行おうとすると下記のエラーが出ます。
> QGISを再起動すれば出力できるのですが二回目は出来ません。
>
> 使用環境はWin7(32bit) QGIS1.8.0(スタンドアロン) メモリ4Gです。
> qProfが試用版なので、プラグイン自体の問題なのかもしれませんが
> なにか対処法がありましたら、ご提案頂ければありがたいです。
>
> --------------
> Pythonプログラムを実行中にエラーが発生しました:
>
> Traceback (most recent call last):
>   File "C:\PROGRA~1\Quantum GIS
> Lisboa\apps\qgis\python\plugins\qProf\qProf_dialog.py", line 212, in
> calculate_profiles
>     dem_params, dem_array = read_raster_band( dem_fpath )
>   File "C:\PROGRA~1\Quantum GIS
> Lisboa\apps\qgis\python\plugins\qProf\gis_utils.py", line 152, in
> read_raster_band
>     data = np.where( abs( data - raster_params.noDataValue ) > 1e-05,
> data, np.NaN )
> MemoryError
>
> Pythonバージョン:
> 2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 bit (Intel)]
> ----------------
>
> 久保まゆみ
> 寒地土木研究所
>
>
> Date: Sun, 14 Apr 2013 21:57:51 +0900
>> From: Minoru Akagi <akaginch ＠ gmail.com>
>> To: 宮園正敏 <miyazono3725 ＠ gmail.com>
>> Cc: OSGeoJapan-discuss ＠ lists.osgeo.org
>> Subject: Re: [OSGeoJapan-discuss]
>>         QGISにて横断ラインを作成する方法をご教示ください
>> Message-ID:
>>         <CA+Bj0=o1=BnYJhuHVC=
>> T-q2nsg1m5Zg1aQJt0SCwqR_gFG8A5A ＠ mail.gmail.com>
>> Content-Type: text/plain; charset="iso-2022-jp"
>>
>>
>> 宮園様
>>
>> 赤木と申します。
>>
>>
>> プラグインリポジトリ[1]には断面形を計測するプラグインがいくつか公開されていますが，現在動作が確認できるもので，ご要望の機能を持ったプラグインは「qProf」[2]になるかと思います。このプラグインはまだ実験的なものとして公開されていますので，プラグインインストーラのオプションタブで「Show
>> all plugins, even those marked as experimental」オプションを選択しなければ一覧に表示されません。
>>
>>
>> ライン(シェープ)とDEMラスタ，測点間の最大距離を指定すれば計測が行え，結果はチャートに表示したり，CSV形式等で出力できます。測点での標高のサンプリング方法はバイリニア法です。現状ではドキュメントが十分ではないと思いますが試してみられてはいかがでしょう。
>>
>> [1] http://plugins.qgis.org/
>> [2] http://plugins.qgis.org/plugins/qProf/
>>
>> 赤木 実
>>
>
> _______________________________________________
> OSGeoJapan-discuss mailing list
> OSGeoJapan-discuss ＠ lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>
>
-------------- next part --------------
HTMLの添付ファイルを保管しました...
URL: <http://lists.osgeo.org/pipermail/osgeojapan-discuss/attachments/20130419/559c6009/attachment-0001.html>