[OSGeoJapan-discuss] gdal_rasterizeで生成した画像がズレる。

斎藤 nsaito ＠ msk-web.co.jp
2016年 10月 13日 (木) 22:23:05 PDT 

いわさき様

　こちらで、新たにデータを作って確認してみました。
　ズレている現象の画像を２パターン以下のURLに示します。

　＜10m解像度でのサンプル画像＞
　http://www.msk-web.co.jp/gdal_rasterize_problem/10m.PNG
　10mの解像度でラスタ画像を生成してみた例です。

　＜50m解像度でのサンプル画像＞
　http://www.msk-web.co.jp/gdal_rasterize_problem/50m.PNG
　50mの解像度でラスタ画像を生成してみた例です。

　黄色のメッシュは、SHPファイルをそのまま表示させたものです。
　左上のメッシュは50mメッシュ
　右下のメッシュは10mメッシュ
　になります。

　憶測に過ぎませんが、gdal_rasterize.exe内で、
　①実数の扱いミスで、桁落ちが生じている？
　②間違って、メッシュの左上ではなく中心座標を起点にしている？
　などの原因を想定しています。。。

　いずれにしても、プログラムソースを追いかける必要がありますが、
体力と時間が乏しいので、しばし放置します。
＃バージョン更新で修正が成されることを、少し期待しながら(^^ゞ

　しかし、同現象の発生が私だけではなかったことに安堵しました。
　現象再現を行って頂きありがとうございました！

以上
　

On 2016/10/14 9:13, Nobusuke Iwasaki wrote:
> さいとうさま
>
> いわさきです。追加情報、ありがとうございます。
>
> ちなみに、ずれているというのは、具体的にどんな感じでしょうか？
> 私もためにしにやってみたところ、以下のリンク画像のように端がずれている状態でした。ただ、地図内部にはずれはありませんでした。
>
> https://dl.dropboxusercontent.com/u/1876391/img001697.png
>
> GDALのバージョンは同じく2.1.1になります。
>
>
> 2016年10月13日 18:00 斎藤<nsaito ＠ msk-web.co.jp>:
>> さいとうです。
>>
>> 　連投で失礼します。
>>
>> 　仮想でWindows 10を真っ新な状態で立ち上げて、OSGeo4W64をインストールしました。
>> 　OSGeo4WのShellで動作を確認しましたが、同様の現象を確認しました。
>>
>> 　25mメッシュのSHPを使用したのですが、ギリギリの分解能だとズレが生じる様でした。
>> 　分解能を細かくすれば現象は収まるのですが・・・処理に時間が掛かってしまうのが厄介です(^^ゞ
>>
>> 　細かなメッシュデータでの確認までは行っていないので、原因は定かではありませんが、
>> 桁落ちとか？従来と少し変わってしまっているところがあるのかな？
>>
>> 　取り敢えず、しばらくは様子を見て、改善されたら？改善方法が見つかったら、改めて
>> 投稿させて頂きたいと思います・・・・（覚えてるかな？ｗ）
>>
>> 以上
>>
>>
>> On 2016/10/13 15:18, 斎藤 wrote:
>>> いわさき様
>>>
>>> 　さいとうです。
>>> 　ご意見いただきありがとうございます。
>>>
>>> 　結果的に、gdal_rasterize.exeを変更しました。
>>> 　具体的には、障害が発生していたgdal_rasterize.exeを古いバージョンのgdal_rasterize.exeに戻しました。
>>> 　障害が発生したバージョン：2.1.1
>>> 　正常に動作したバージョン：1.11.3
>>>
>>> 　障害が確認されたのは、弊方における３台のWindows PCです。
>>> 　OS：Windows 10 x 2台
>>> 　　　  Windows 8.1 x 1台
>>>
>>> 　使用しているパッケージ：OSGeo4W64
>>> ＜障害発生ファイルの詳細＞
>>> 　対象ファイル：gdal_rasterize.exe
>>> 　ファイルサイズ：49,152バイト
>>> 　sumの値：26269    48
>>> 　バージョン　2.1.1
>>> 　（OSGeo4W64のsetupでは、2.1.1-3で発生。2.1.1-1でも同様の現象を確認）
>>>
>>> 　一応、真っ新なマシンにOSGeo4Wを入れて確認してみようとは思っていますが、
>>> 一旦ここまでで止めて、本来の画像生成作業に戻ろうと思います。
>>>
>>> 　ご協力ありがとうございました。
>>>
>>> 以上
>>>
>>> On 2016/10/12 16:16, Nobusuke Iwasaki wrote:
>>>> 齋藤様
>>>>
>>>> いわさきです。
>>>>
>>>>> gdal_rasterize -a D_MAX -tr 0.0003125 0.000208333 -a_nodata 9999 in.shp
>>>>> out.tif
>>>> のコマンドで、-s_srs で座標系を明示的に指定してみてはいかがでしょうか？
>>>> あまり、関係ないかもしれませんが･･･。
>>>>
>>>> また、ずれている画像では、gdalinfoの結果のうち、
>>>>
>>>>     Min=0.000 Max=3.950
>>>>     Minimum=0.000, Maximum=3.950, Mean=0.807, StdDev=0.760
>>>>     NoData Value=9999（ここは同じ）
>>>>     Metadata:
>>>>       STATISTICS_MAXIMUM=3.95
>>>>       STATISTICS_MEAN=0.80744435917114
>>>>       STATISTICS_MINIMUM=0
>>>>       STATISTICS_STDDEV=0.76004957512707
>>>>
>>>> という列が、上手くいった場合にくらべ多いようです。
>>>> ここの部分は、通常のgdalinfoでは出力されずに、gdalinfo に-statsオプションを付けると生成されるようです。
>>>>
>>>> デフォルトでついているとすると、なにかおかしくなっているかもしれません。
>>>>
>>>> 以上、参考までに。
>>>>
>>>> 2016年10月12日 15:09 斎藤<nsaito ＠ msk-web.co.jp>:
>>>>> 田原様
>>>>>
>>>>> 　斎藤です。
>>>>> 　ご意見いただきありがとうございます。
>>>>>
>>>>> 　異常な結果を示す環境も、正常な結果を示す環境も64Bit版でした。
>>>>> 　また、今回、正常に動作していたマシンでOSGeo4W64のバージョン更新を行いました
>>>>> ところ、ズレてしまう現象が発生しました。
>>>>>
>>>>> 　現在、正常動作をしているマシンからOSGeo4Wのディレクトリ以下をコピーして、
>>>>> 障害が発生しているマシンへ部分コピーを行いながら検証してみることにしています。
>>>>> 　gdal_rasterize.exeだけコピーしても現象は変わらなかったので、何か他の物に影響されて
>>>>> いると見ています。
>>>>>
>>>>> 　結果は別途報告させて頂きます。
>>>>>
>>>>> 　ちなみに、shapeファイルは全く同じファイルを使用して検証しています。
>>>>>
>>>>> 以上
>>>>>
>>>>>
>>>>>
>>>>> On 2016/10/08 17:05, DaisukeTahara wrote:
>>>>>> 田原＠HCCです。
>>>>>>
>>>>>> gdalに限りませんが、32bit/64bit版の違いで動作異なる事あったりしましたので
>>>>>> その違いだったりしないでしょうか。
>>>>>>
>>>>>> とりあえずテスト比較したgdalのバージョン、マシンのＯＳ/CPUなんかを
>>>>>> 比較提示してみると原因特定しやすいかなと思います。
>>>>>>
>>>>>> ---
>>>>>> Ｃｅｎｔｅｒが同じですけど四隅範囲がだいぶ違うので、ズレているほうのshapeは
>>>>>> ファイル名が同じで違う範囲のファイルを変換していたりしてないでしょうか。
>>>>>>
>>>>>>
>>>>>>
>>>>>> On 2016/10/06 18:49, 斎藤 wrote:
>>>>>>> 斎藤です。
>>>>>>>
>>>>>>> 奇妙な現象に悩まされています。
>>>>>>> gdal_rasterizeを使って以下のコマンドを実行しました。
>>>>>>>
>>>>>>> gdal_rasterize -a D_MAX -tr 0.0003125 0.000208333 -a_nodata 9999
>>>>>>> in.shp
>>>>>>> out.tif
>>>>>>>
>>>>>>> 手持ちのマシン５台で確認したところ、内３台でズレることなく生成した画像はSHPの画像と
>>>>>>> 一致するのですが、内2台で座標系を間違ってしまったかの様に画像がズレます。
>>>>>>>
>>>>>>> osgeo4wのパッケージを再インストールしてみましたが、現象は変わりません。
>>>>>>> また、古いパッケージでも同様の現象が発生する様なので、バージョンに依存している
>>>>>>> とも思われません。
>>>>>>>
>>>>>>> 原因が特定できず、お知恵を拝借できないでしょうか？
>>>>>>>
>>>>>>> 正常に生成された画像とズレた画像で、gdalinfoを確認してみました結果は
>>>>>>> 以下の通りです。
>>>>>>>
>>>>>>> ＜＜ズレた画像のgdalinfo＞＞
>>>>>>> Driver: GTiff/GeoTIFF
>>>>>>> Files: aaa.tif
>>>>>>> Size is 676, 478
>>>>>>> Coordinate System is:
>>>>>>> GEOGCS["WGS 84",
>>>>>>>       DATUM["WGS_1984",
>>>>>>>           SPHEROID["WGS 84",6378137,298.257223563,
>>>>>>>               AUTHORITY["EPSG","7030"]],
>>>>>>>           AUTHORITY["EPSG","6326"]],
>>>>>>>       PRIMEM["Greenwich",0],
>>>>>>>       UNIT["degree",0.0174532925199433],
>>>>>>>       AUTHORITY["EPSG","4326"]]
>>>>>>> Origin = (136.377968750000180,34.944687499833471)
>>>>>>> Pixel Size = (0.000312500000000,-0.000208333000000)
>>>>>>> Metadata:
>>>>>>>     AREA_OR_POINT=Area
>>>>>>> Image Structure Metadata:
>>>>>>>     INTERLEAVE=BAND
>>>>>>> Corner Coordinates:
>>>>>>> Upper Left  ( 136.3779688,  34.9446875) (136d22'40.69"E,
>>>>>>> 34d56'40.87"N)
>>>>>>> Lower Left  ( 136.3779688,  34.8451043) (136d22'40.69"E,
>>>>>>> 34d50'42.38"N)
>>>>>>> Upper Right ( 136.5892188,  34.9446875) (136d35'21.19"E,
>>>>>>> 34d56'40.87"N)
>>>>>>> Lower Right ( 136.5892188,  34.8451043) (136d35'21.19"E,
>>>>>>> 34d50'42.38"N)
>>>>>>> Center      ( 136.4835938,  34.8948959) (136d29' 0.94"E,
>>>>>>> 34d53'41.63"N)
>>>>>>> Band 1 Block=676x1 Type=Float64, ColorInterp=Gray
>>>>>>>     Min=0.000 Max=3.950
>>>>>>>     Minimum=0.000, Maximum=3.950, Mean=0.807, StdDev=0.760
>>>>>>>     NoData Value=9999
>>>>>>>     Metadata:
>>>>>>>       STATISTICS_MAXIMUM=3.95
>>>>>>>       STATISTICS_MEAN=0.80744435917114
>>>>>>>       STATISTICS_MINIMUM=0
>>>>>>>       STATISTICS_STDDEV=0.76004957512707
>>>>>>>
>>>>>>> ＜＜正常な時にgdalinfo＞＞
>>>>>>> Driver: GTiff/GeoTIFF
>>>>>>> Files: bbb.tif
>>>>>>> Size is 675, 477
>>>>>>> Coordinate System is:
>>>>>>> GEOGCS["WGS 84",
>>>>>>>       DATUM["WGS_1984",
>>>>>>>           SPHEROID["WGS 84",6378137,298.257223563,
>>>>>>>               AUTHORITY["EPSG","7030"]],
>>>>>>>           AUTHORITY["EPSG","6326"]],
>>>>>>>       PRIMEM["Greenwich",0],
>>>>>>>       UNIT["degree",0.0174532925199433],
>>>>>>>       AUTHORITY["EPSG","4326"]]
>>>>>>> Origin = (136.378125000000180,34.944583333333469)
>>>>>>> Pixel Size = (0.000312500000000,-0.000208333000000)
>>>>>>> Metadata:
>>>>>>>     AREA_OR_POINT=Area
>>>>>>> Image Structure Metadata:
>>>>>>>     INTERLEAVE=BAND
>>>>>>> Corner Coordinates:
>>>>>>> Upper Left  ( 136.3781250,  34.9445833) (136d22'41.25"E,
>>>>>>> 34d56'40.50"N)
>>>>>>> Lower Left  ( 136.3781250,  34.8452085) (136d22'41.25"E,
>>>>>>> 34d50'42.75"N)
>>>>>>> Upper Right ( 136.5890625,  34.9445833) (136d35'20.63"E,
>>>>>>> 34d56'40.50"N)
>>>>>>> Lower Right ( 136.5890625,  34.8452085) (136d35'20.63"E,
>>>>>>> 34d50'42.75"N)
>>>>>>> Center      ( 136.4835938,  34.8948959) (136d29' 0.94"E,
>>>>>>> 34d53'41.63"N)
>>>>>>> Band 1 Block=675x1 Type=Float64, ColorInterp=Gray
>>>>>>>     NoData Value=9999
>>>>>>>
>>>>>>>
>>>>>>> _______________________________________________
>>>>>>> OSGeoJapan-discuss mailing list
>>>>>>> OSGeoJapan-discuss ＠ lists.osgeo.org
>>>>>>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>>>>>> _______________________________________________
>>>>>> OSGeoJapan-discuss mailing list
>>>>>> OSGeoJapan-discuss ＠ lists.osgeo.org
>>>>>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>>>>>
>>>>>
>>>>> _______________________________________________
>>>>> 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/20161014/67b62f45/attachment-0001.html>

OSGeoJapan-discuss メーリングリストの案内