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

2013年 4月 23日 (火) 02:03:35 PDT

赤木 様

ご丁寧な返信ありがとうございます。
風邪をひいて返信が遅くなってしまいました。
これでQGISでも必要な横断データが取れそうです。
ありがとうございました。

久保まゆみ

2013年4月21日日曜日 Minoru Akagi akaginch ＠ gmail.com:

> 久保さま，皆様
>
> スタンドアロン版のQGIS 1.8でProfile from
>
> lineプラグインの動作が確認できました。OSGeo4W版ではPythonのShapelyパッケージを別途インストールする必要があったのでパスしました。
>
>
> プラグインのインストール後，メニューアイテムやツールバーのアイコンが淡色表示になっていますが，ラインレイヤを開いた後にDEMレイヤを開けばクリックを受け付ける状態になります。計測するレイヤを選び測点間隔を入力したあとOKボタンを押すと，久保さんが示されたのと同様のエラーメッセージが表示されました。このエラーに対しては，エラーメッセージに示されているProfileFromLine.pyの252行目のstr()で囲まれたbandNameのstr()を取ってやり，QGISを再起動すれば動作するようになります。
>
>
> 出力されたポイントレイヤを保存する際にカンマ区切りファイルを選び，レイヤ作成オプションにGEOMETRY=AS_XYを入力すればXYZ座標値を出力できます。このプラグインやProfile
> toolプラグイン[1]では測点での標高の計算方法は最近傍が用いられています。
>
> [1] http://plugins.qgis.org/plugins/profiletool/
>
>
> 赤木 実
>
> 2013年4月19日 16:21 Mayumi Kubo <mayumitt ＠ gmail.com <javascript:;>>:
> >
> > 赤木さま、みなさま
> >
> > ご教授ありがとうございます。
> > TIFFファイルの容量を下げてみたら問題なく出力できるようになりました。
> >
> > しかし、「qProf」は線が折れているとそこからまたサンプリングを始めるので
> > 折れ線ではサンプリング間隔が一定にはならないようでした。
> > 一定間隔でサンプリングするために下記のサイトで紹介されていた
> >
> http://gis.stackexchange.com/questions/27102/how-to-create-equidistant-points-in-qgis
> > 「Profile from line」というプラグインを試してみましたが、
> > こちらはUnicodeEncodeErrorが出て出力できませんでした。
> > 質問続きで申し訳ありませんが、、対処法ご存じの方いらっしゃいますでしょうか。
> >
> > ----------------------------------
> >
> > Traceback (most recent call last):
> >   File
> "C:/Users/kan-pc8_user/.qgis//python/plugins\profilefromline\ProfileFromLineDialog.py",
> line 54, in accept
> >     self.ui.progressBar)
> >   File
> "C:/Users/kan-pc8_user/.qgis//python/plugins\profilefromline\ProfileFromLine.py",
> line 166, in create_profile_layer
> >     rasterLayers, rasterFieldsStart)
> >   File
> "C:/Users/kan-pc8_user/.qgis//python/plugins\profilefromline\ProfileFromLine.py",
> line 184, in create_points_from_line
> >     rasterFieldsStart)
> >   File
> "C:/Users/kan-pc8_user/.qgis//python/plugins\profilefromline\ProfileFromLine.py",
> line 252, in create_new_point
> >     if str(bandName) == rLayer.bandName(1):
> > UnicodeEncodeError: 'ascii' codec can't encode characters in position
> 0-2: ordinal not in range(128)
> >
> >
> > Pythonバージョン:
> > 2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 bit (Intel)]
> > -------------------------------
> >
> > 久保まゆみ
> >
> >
> > 2013年4月19日 14:33 Minoru Akagi <akaginch ＠ gmail.com>:
> >
> >> 久保さま
> >>
> >> 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:
> >>>>
-------------- next part --------------
HTMLの添付ファイルを保管しました...
URL: <http://lists.osgeo.org/pipermail/osgeojapan-discuss/attachments/20130423/972e4e27/attachment.html>