ClipperでTIFFファイルの切出し後ズレが発生します

292 views
Skip to first unread message

TJ

unread,
Jan 22, 2016, 5:42:08 AM1/22/16
to FMEユーザーフォーラム
Clipperを使って、メッシュのポリゴンでTIFF画像(tfw付)を切出して、GEOTIFFライダーでGEOTIFF形式で出力しています。


切出し後のGEOTIFF画像と切出し前のTIFF画像をArcGISで確認すると、ズレが発生しています。
しかし、ArcGISのClip機能を使って、同じメッシュとTIFF画像を使ってで切出してみた結果ズレが発生しませんでした。

このズレはどうすれば解消できるのでしょうか?
Auto Generated Inline Image 1

Takashi Iijima

unread,
Jan 22, 2016, 6:54:20 AM1/22/16
to FMEユーザーフォーラム
メッシュの境界とラスターのセルの境界がピッタリ一致していない場合にずれることはあると思います。
どのようなケースで、どの程度のズレが生じていますか?

TJ

unread,
Jan 22, 2016, 10:48:46 AM1/22/16
to fm...@googlegroups.com
TIFF画像に長方形の枠線があります。その4隅に合わせてメッシュを作成してあります(添付図の紫色の枠線)。
両方の座標系とも旧日本測地系に定義してあります。

クリップ後、左上の隅だけずれます。ArcGISの計測ツールで測ると約420mあります。


>メッシュの境界とラスターのセルの境界がピッタリ一致していない場合にずれることはあると思います。
この場合、画像全体ほぼ均等にずれると思いますが、今回の場合、左上だけずれが発生しますので、何か不思議です。


Takashi Iijima

unread,
Jan 22, 2016, 8:34:04 PM1/22/16
to fm...@googlegroups.com
再現できませんでした。添付のデモ(サンプルデータ含む)の結果は次のとおりです。

基盤地図情報10mメッシュDEM(三宅島)をラスターに変換してGeoTIFF形式で保存。
それをShape形式の矩形ポリゴンデータによってクリップした結果をGeoTIFF形式で保存したうえで、ArcMap で表示。


差し支えなければ、状況が再現できる最小限のデータとワークスペースを見せていただけませんか?
raster-clipping-example.fmwt

TJ

unread,
Jan 22, 2016, 9:56:21 PM1/22/16
to FMEユーザーフォーラム
DEMデータで検証していただいてありがとうございました。

データとワークスペースをメールで送りいたしましたので、ご検証をよろしくお願いいたします。

Takashi Iijima

unread,
Jan 24, 2016, 7:41:20 PM1/24/16
to FMEユーザーフォーラム
お送りいただいたデータを FME Data Inspector で調べたところ、画像中の枠線がメッシュポリゴンの境界とずれていました。
GeoTIFFファイルと同じ名前で、*.tif.aux.xml, *.tif.ovr という補助ファイルが含まれています。
これらはArcGISで何らかの作業をしたときに作成されるもので、GeoTIFFデータ自体が修正されていなくても、ArcGISはこれらの補助ファイルを参照して、幾何補正や座標変換した後の状態での表示や処理を行う場合があります。
しかし、FMEはデータをあるがままに取り扱うことを大原則としており、これらの補助ファイルを参照して自動的に幾何補正などを行うことはありません。
おそらくこれが、ArcMapで表示したときには画像中の枠線とポリゴンの境界が一致しているのに、FMEでクリップするとずれてしまう、ということの理由だと思います。
データを修正しなくても、幾何補正などをした状態での表示や処理ができるのはArcGISの優れた機能ですが、反面、実際のデータの状態を調べるのには不向きであると言えます。

TJ

unread,
Jan 24, 2016, 8:56:07 PM1/24/16
to FMEユーザーフォーラム
>データを修正しなくても、幾何補正などをした状態での表示や処理ができるのはArcGISの優れた機能ですが、反面、実際のデータの状態を調べるのには不向きであると言えます。
ご指摘のように、今回はArcGISで表示した状態の画像を正としていたのが問題の原因だと思います。今回の経験はいい教訓だと思います。今後、データ検証する場合、なるべくFME Data Inspectorを使って検証するようにします。

教えていただいた情報を参考し、以下の2つの解決法を考えました(両方検証済み)。
  1. ArcMapおよびArcCatologのデータ管理ツール ⇒ ラスタ ⇒ ラスタ プロセシング ⇒ クリップツールで処理する
  2. 画像を一旦ArcMap上に追加し、追加した画像レイヤを右クリック ⇒ データ ⇒ データのエクスポート(レクティファイと相当する?)を実行する ⇒ FMEのClipperで書き出した画像に対してクリップ処理をする

いろいろ検証していただき、本当にありがとうございました。

Reply all
Reply to author
Forward
0 new messages