<html>
<head>
<meta content="text/html; charset=utf-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">いわさき様<br>
<br>
こちらで、新たにデータを作って確認してみました。<br>
ズレている現象の画像を２パターン以下のURLに示します。<br>
<br>
＜10m解像度でのサンプル画像＞<br>
<a class="moz-txt-link-freetext"
href="http://www.msk-web.co.jp/gdal_rasterize_problem/10m.PNG">http://www.msk-web.co.jp/gdal_rasterize_problem/10m.PNG</a><br>
10mの解像度でラスタ画像を生成してみた例です。<br>
<br>
＜50m解像度でのサンプル画像＞<br>
<a class="moz-txt-link-freetext"
href="http://www.msk-web.co.jp/gdal_rasterize_problem/50m.PNG">http://www.msk-web.co.jp/gdal_rasterize_problem/50m.PNG</a><br>
50mの解像度でラスタ画像を生成してみた例です。<br>
<br>
黄色のメッシュは、SHPファイルをそのまま表示させたものです。<br>
左上のメッシュは50mメッシュ<br>
右下のメッシュは10mメッシュ<br>
になります。<br>
<br>
憶測に過ぎませんが、gdal_rasterize.exe内で、<br>
①実数の扱いミスで、桁落ちが生じている？<br>
②間違って、メッシュの左上ではなく中心座標を起点にしている？<br>
などの原因を想定しています。。。<br>
<br>
いずれにしても、プログラムソースを追いかける必要がありますが、<br>
体力と時間が乏しいので、しばし放置します。<br>
＃バージョン更新で修正が成されることを、少し期待しながら(^^ゞ<br>
<br>
しかし、同現象の発生が私だけではなかったことに安堵しました。<br>
現象再現を行って頂きありがとうございました！<br>
<br>
以上<br>
<br>
<br>
On 2016/10/14 9:13, Nobusuke Iwasaki wrote:<br>
</div>
<blockquote
cite="mid:CAFvDKTbSxWrmFfND+v3dYE37v0fLxoAQX5ukAUxe--=Ec2-QPA@mail.gmail.com"
type="cite">
<pre wrap="">さいとうさま
いわさきです。追加情報、ありがとうございます。
ちなみに、ずれているというのは、具体的にどんな感じでしょうか？
私もためにしにやってみたところ、以下のリンク画像のように端がずれている状態でした。ただ、地図内部にはずれはありませんでした。
<a class="moz-txt-link-freetext" href="https://dl.dropboxusercontent.com/u/1876391/img001697.png">https://dl.dropboxusercontent.com/u/1876391/img001697.png</a>
GDALのバージョンは同じく2.1.1になります。
2016年10月13日 18:00 斎藤 <a class="moz-txt-link-rfc2396E" href="mailto:nsaito@msk-web.co.jp"><nsaito@msk-web.co.jp></a>:
</pre>
<blockquote type="cite">
<pre wrap="">さいとうです。
連投で失礼します。
仮想でWindows 10を真っ新な状態で立ち上げて、OSGeo4W64をインストールしました。
OSGeo4WのShellで動作を確認しましたが、同様の現象を確認しました。
25mメッシュのSHPを使用したのですが、ギリギリの分解能だとズレが生じる様でした。
分解能を細かくすれば現象は収まるのですが・・・処理に時間が掛かってしまうのが厄介です(^^ゞ
細かなメッシュデータでの確認までは行っていないので、原因は定かではありませんが、
桁落ちとか？従来と少し変わってしまっているところがあるのかな？
取り敢えず、しばらくは様子を見て、改善されたら？改善方法が見つかったら、改めて
投稿させて頂きたいと思います・・・・（覚えてるかな？ｗ）
以上
On 2016/10/13 15:18, 斎藤 wrote:
</pre>
<blockquote type="cite">
<pre wrap="">いわさき様
さいとうです。
ご意見いただきありがとうございます。
結果的に、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:
</pre>
<blockquote type="cite">
<pre wrap="">齋藤様
いわさきです。
</pre>
<blockquote type="cite">
<pre wrap="">gdal_rasterize -a D_MAX -tr 0.0003125 0.000208333 -a_nodata 9999 in.shp
out.tif
</pre>
</blockquote>
<pre wrap="">のコマンドで、-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 斎藤 <a class="moz-txt-link-rfc2396E" href="mailto:nsaito@msk-web.co.jp"><nsaito@msk-web.co.jp></a>:
</pre>
<blockquote type="cite">
<pre wrap="">田原様
斎藤です。
ご意見いただきありがとうございます。
異常な結果を示す環境も、正常な結果を示す環境も64Bit版でした。
また、今回、正常に動作していたマシンでOSGeo4W64のバージョン更新を行いました
ところ、ズレてしまう現象が発生しました。
現在、正常動作をしているマシンからOSGeo4Wのディレクトリ以下をコピーして、
障害が発生しているマシンへ部分コピーを行いながら検証してみることにしています。
gdal_rasterize.exeだけコピーしても現象は変わらなかったので、何か他の物に影響されて
いると見ています。
結果は別途報告させて頂きます。
ちなみに、shapeファイルは全く同じファイルを使用して検証しています。
以上
On 2016/10/08 17:05, DaisukeTahara wrote:
</pre>
<blockquote type="cite">
<pre wrap="">
田原＠HCCです。
gdalに限りませんが、32bit/64bit版の違いで動作異なる事あったりしましたので
その違いだったりしないでしょうか。
とりあえずテスト比較したgdalのバージョン、マシンのＯＳ/CPUなんかを
比較提示してみると原因特定しやすいかなと思います。
---
Ｃｅｎｔｅｒが同じですけど四隅範囲がだいぶ違うので、ズレているほうのshapeは
ファイル名が同じで違う範囲のファイルを変換していたりしてないでしょうか。
On 2016/10/06 18:49, 斎藤 wrote:
</pre>
<blockquote type="cite">
<pre wrap="">斎藤です。
奇妙な現象に悩まされています。
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
<a class="moz-txt-link-abbreviated" href="mailto:OSGeoJapan-discuss@lists.osgeo.org">OSGeoJapan-discuss@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss">http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss</a>
</pre>
</blockquote>
<pre wrap="">
_______________________________________________
OSGeoJapan-discuss mailing list
<a class="moz-txt-link-abbreviated" href="mailto:OSGeoJapan-discuss@lists.osgeo.org">OSGeoJapan-discuss@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss">http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss</a>
</pre>
</blockquote>
<pre wrap="">
_______________________________________________
OSGeoJapan-discuss mailing list
<a class="moz-txt-link-abbreviated" href="mailto:OSGeoJapan-discuss@lists.osgeo.org">OSGeoJapan-discuss@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss">http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss</a>
</pre>
</blockquote>
<pre wrap="">
</pre>
</blockquote>
</blockquote>
</blockquote>
</blockquote>
<br>
</body>
</html>