gdal_rasterizeで生成した画像がズレる。

Previous Topic Next Topic
 
classic Classic list List threaded Threaded
8 messages Options
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

gdal_rasterizeで生成した画像がズレる。

斎藤 直正
斎藤です。

奇妙な現象に悩まされています。
gdal_rasterizeを使って以下のコマンドを実行しました。

gdal_rasterize -a D_MAX -tr 0.0003125 0.000208333 -a_nodata 9999 in.shp
out.tif

手持ちのマシン5台で確認したところ、内3台でズレることなく生成した画像は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
[hidden email]
http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Re: gdal_rasterizeで生成した画像がズレる。

DaisukeTahara

田原@HCCです。

gdalに限りませんが、32bit/64bit版の違いで動作異なる事あったりしましたので
その違いだったりしないでしょうか。

とりあえずテスト比較したgdalのバージョン、マシンのOS/CPUなんかを
比較提示してみると原因特定しやすいかなと思います。

---
Centerが同じですけど四隅範囲がだいぶ違うので、ズレているほうの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
>
> 手持ちのマシン5台で確認したところ、内3台でズレることなく生成した画像は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
> [hidden email]
> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss

_______________________________________________
OSGeoJapan-discuss mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Re: gdal_rasterizeで生成した画像がズレる。

斎藤 直正
田原様

 斎藤です。
 ご意見いただきありがとうございます。
 
 異常な結果を示す環境も、正常な結果を示す環境も64Bit版でした。
 また、今回、正常に動作していたマシンでOSGeo4W64のバージョン更新を行いました
ところ、ズレてしまう現象が発生しました。

 現在、正常動作をしているマシンからOSGeo4Wのディレクトリ以下をコピーして、
障害が発生しているマシンへ部分コピーを行いながら検証してみることにしています。
 gdal_rasterize.exeだけコピーしても現象は変わらなかったので、何か他の物に影響されて
いると見ています。

 結果は別途報告させて頂きます。

 ちなみに、shapeファイルは全く同じファイルを使用して検証しています。

以上


On 2016/10/08 17:05, DaisukeTahara wrote:

>
> 田原@HCCです。
>
> gdalに限りませんが、32bit/64bit版の違いで動作異なる事あったりしましたので
> その違いだったりしないでしょうか。
>
> とりあえずテスト比較したgdalのバージョン、マシンのOS/CPUなんかを
> 比較提示してみると原因特定しやすいかなと思います。
>
> ---
> Centerが同じですけど四隅範囲がだいぶ違うので、ズレているほうの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
>>
>> 手持ちのマシン5台で確認したところ、内3台でズレることなく生成した画像は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
>> [hidden email]
>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>
> _______________________________________________
> OSGeoJapan-discuss mailing list
> [hidden email]
> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss


--
/////////////////////////////////////////////////////////
株式会社 エムエスケー
斎藤 直正
TEL:0598-51-7471 / FAX:0598-52-4849
E-mail:[hidden email]  URL:http://www.msk-web.co.jp
/////////////////////////////////////////////////////////

_______________________________________________
OSGeoJapan-discuss mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Re: gdal_rasterizeで生成した画像がズレる。

Nobusuke Iwasaki-3
齋藤様

いわさきです。

> 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 斎藤 <[hidden email]>:

> 田原様
>
>  斎藤です。
>  ご意見いただきありがとうございます。
>
>  異常な結果を示す環境も、正常な結果を示す環境も64Bit版でした。
>  また、今回、正常に動作していたマシンでOSGeo4W64のバージョン更新を行いました
> ところ、ズレてしまう現象が発生しました。
>
>  現在、正常動作をしているマシンからOSGeo4Wのディレクトリ以下をコピーして、
> 障害が発生しているマシンへ部分コピーを行いながら検証してみることにしています。
>  gdal_rasterize.exeだけコピーしても現象は変わらなかったので、何か他の物に影響されて
> いると見ています。
>
>  結果は別途報告させて頂きます。
>
>  ちなみに、shapeファイルは全く同じファイルを使用して検証しています。
>
> 以上
>
>
>
> On 2016/10/08 17:05, DaisukeTahara wrote:
>>
>>
>> 田原@HCCです。
>>
>> gdalに限りませんが、32bit/64bit版の違いで動作異なる事あったりしましたので
>> その違いだったりしないでしょうか。
>>
>> とりあえずテスト比較したgdalのバージョン、マシンのOS/CPUなんかを
>> 比較提示してみると原因特定しやすいかなと思います。
>>
>> ---
>> Centerが同じですけど四隅範囲がだいぶ違うので、ズレているほうの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
>>>
>>> 手持ちのマシン5台で確認したところ、内3台でズレることなく生成した画像は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
>>> [hidden email]
>>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>>
>>
>> _______________________________________________
>> OSGeoJapan-discuss mailing list
>> [hidden email]
>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>
>
>
> --
> /////////////////////////////////////////////////////////
> 株式会社 エムエスケー
> 斎藤 直正
> TEL:0598-51-7471 / FAX:0598-52-4849
> E-mail:[hidden email]  URL:http://www.msk-web.co.jp
> /////////////////////////////////////////////////////////
>
>
> _______________________________________________
> OSGeoJapan-discuss mailing list
> [hidden email]
> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss



--
岩崎 亘典
_______________________________________________
OSGeoJapan-discuss mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Re: gdal_rasterizeで生成した画像がズレる。

斎藤 直正
いわさき様

 さいとうです。
 ご意見いただきありがとうございます。

 結果的に、gdal_rasterize.exeを変更しました。
 具体的には、障害が発生していたgdal_rasterize.exeを古いバージョンのgdal_rasterize.exeに戻しました。
 障害が発生したバージョン:2.1.1
 正常に動作したバージョン:1.11.3

 障害が確認されたのは、弊方における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 斎藤 <[hidden email]>:
>> 田原様
>>
>>  斎藤です。
>>  ご意見いただきありがとうございます。
>>
>>  異常な結果を示す環境も、正常な結果を示す環境も64Bit版でした。
>>  また、今回、正常に動作していたマシンでOSGeo4W64のバージョン更新を行いました
>> ところ、ズレてしまう現象が発生しました。
>>
>>  現在、正常動作をしているマシンからOSGeo4Wのディレクトリ以下をコピーして、
>> 障害が発生しているマシンへ部分コピーを行いながら検証してみることにしています。
>>  gdal_rasterize.exeだけコピーしても現象は変わらなかったので、何か他の物に影響されて
>> いると見ています。
>>
>>  結果は別途報告させて頂きます。
>>
>>  ちなみに、shapeファイルは全く同じファイルを使用して検証しています。
>>
>> 以上
>>
>>
>>
>> On 2016/10/08 17:05, DaisukeTahara wrote:
>>>
>>> 田原@HCCです。
>>>
>>> gdalに限りませんが、32bit/64bit版の違いで動作異なる事あったりしましたので
>>> その違いだったりしないでしょうか。
>>>
>>> とりあえずテスト比較したgdalのバージョン、マシンのOS/CPUなんかを
>>> 比較提示してみると原因特定しやすいかなと思います。
>>>
>>> ---
>>> Centerが同じですけど四隅範囲がだいぶ違うので、ズレているほうの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
>>>>
>>>> 手持ちのマシン5台で確認したところ、内3台でズレることなく生成した画像は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
>>>> [hidden email]
>>>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>>>
>>> _______________________________________________
>>> OSGeoJapan-discuss mailing list
>>> [hidden email]
>>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>>
>>
>> --
>> /////////////////////////////////////////////////////////
>> 株式会社 エムエスケー
>> 斎藤 直正
>> TEL:0598-51-7471 / FAX:0598-52-4849
>> E-mail:[hidden email]  URL:http://www.msk-web.co.jp
>> /////////////////////////////////////////////////////////
>>
>>
>> _______________________________________________
>> OSGeoJapan-discuss mailing list
>> [hidden email]
>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>
>


--
/////////////////////////////////////////////////////////
株式会社 エムエスケー
斎藤 直正
TEL:0598-51-7471 / FAX:0598-52-4849
E-mail:[hidden email]  URL:http://www.msk-web.co.jp
/////////////////////////////////////////////////////////

_______________________________________________
OSGeoJapan-discuss mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Re: gdal_rasterizeで生成した画像がズレる。

斎藤 直正
さいとうです。

 連投で失礼します。

 仮想でWindows
10を真っ新な状態で立ち上げて、OSGeo4W64をインストールしました。
 OSGeo4WのShellで動作を確認しましたが、同様の現象を確認しました。

 25mメッシュのSHPを使用したのですが、ギリギリの分解能だとズレが生じる様でした。
 分解能を細かくすれば現象は収まるのですが・・・処理に時間が掛かってしまうのが厄介です(^^ゞ

 細かなメッシュデータでの確認までは行っていないので、原因は定かではありませんが、
桁落ちとか?従来と少し変わってしまっているところがあるのかな?

 取り敢えず、しばらくは様子を見て、改善されたら?改善方法が見つかったら、改めて
投稿させて頂きたいと思います・・・・(覚えてるかな?w)

以上

On 2016/10/13 15:18, 斎藤 wrote:

> いわさき様
>
>  さいとうです。
>  ご意見いただきありがとうございます。
>
>  結果的に、gdal_rasterize.exeを変更しました。
>  具体的には、障害が発生していたgdal_rasterize.exeを古いバージョンのgdal_rasterize.exeに戻しました。
>
>  障害が発生したバージョン:2.1.1
>  正常に動作したバージョン:1.11.3
>
>  障害が確認されたのは、弊方における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 斎藤 <[hidden email]>:
>>> 田原様
>>>
>>>  斎藤です。
>>>  ご意見いただきありがとうございます。
>>>
>>>  異常な結果を示す環境も、正常な結果を示す環境も64Bit版でした。
>>>  また、今回、正常に動作していたマシンでOSGeo4W64のバージョン更新を行いました
>>>
>>> ところ、ズレてしまう現象が発生しました。
>>>
>>>  現在、正常動作をしているマシンからOSGeo4Wのディレクトリ以下をコピーして、
>>>
>>> 障害が発生しているマシンへ部分コピーを行いながら検証してみることにしています。
>>>
>>>  gdal_rasterize.exeだけコピーしても現象は変わらなかったので、何か他の物に影響されて
>>>
>>> いると見ています。
>>>
>>>  結果は別途報告させて頂きます。
>>>
>>>  ちなみに、shapeファイルは全く同じファイルを使用して検証しています。
>>>
>>> 以上
>>>
>>>
>>>
>>> On 2016/10/08 17:05, DaisukeTahara wrote:
>>>>
>>>> 田原@HCCです。
>>>>
>>>> gdalに限りませんが、32bit/64bit版の違いで動作異なる事あったりしましたので
>>>>
>>>> その違いだったりしないでしょうか。
>>>>
>>>> とりあえずテスト比較したgdalのバージョン、マシンのOS/CPUなんかを
>>>> 比較提示してみると原因特定しやすいかなと思います。
>>>>
>>>> ---
>>>> Centerが同じですけど四隅範囲がだいぶ違うので、ズレているほうの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
>>>>>
>>>>> 手持ちのマシン5台で確認したところ、内3台でズレることなく生成した画像は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
>>>>> [hidden email]
>>>>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>>>>
>>>> _______________________________________________
>>>> OSGeoJapan-discuss mailing list
>>>> [hidden email]
>>>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>>>
>>>
>>> --
>>> /////////////////////////////////////////////////////////
>>> 株式会社 エムエスケー
>>> 斎藤 直正
>>> TEL:0598-51-7471 / FAX:0598-52-4849
>>> E-mail:[hidden email]  URL:http://www.msk-web.co.jp
>>> /////////////////////////////////////////////////////////
>>>
>>>
>>> _______________________________________________
>>> OSGeoJapan-discuss mailing list
>>> [hidden email]
>>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>>
>>
>
>


--
/////////////////////////////////////////////////////////
株式会社 エムエスケー
斎藤 直正
TEL:0598-51-7471 / FAX:0598-52-4849
E-mail:[hidden email]  URL:http://www.msk-web.co.jp
/////////////////////////////////////////////////////////

_______________________________________________
OSGeoJapan-discuss mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Re: gdal_rasterizeで生成した画像がズレる。

Nobusuke Iwasaki-3
さいとうさま

いわさきです。追加情報、ありがとうございます。

ちなみに、ずれているというのは、具体的にどんな感じでしょうか?
私もためにしにやってみたところ、以下のリンク画像のように端がずれている状態でした。ただ、地図内部にはずれはありませんでした。

https://dl.dropboxusercontent.com/u/1876391/img001697.png

GDALのバージョンは同じく2.1.1になります。


2016年10月13日 18:00 斎藤 <[hidden email]>:

> さいとうです。
>
>  連投で失礼します。
>
>  仮想でWindows 10を真っ新な状態で立ち上げて、OSGeo4W64をインストールしました。
>  OSGeo4WのShellで動作を確認しましたが、同様の現象を確認しました。
>
>  25mメッシュのSHPを使用したのですが、ギリギリの分解能だとズレが生じる様でした。
>  分解能を細かくすれば現象は収まるのですが・・・処理に時間が掛かってしまうのが厄介です(^^ゞ
>
>  細かなメッシュデータでの確認までは行っていないので、原因は定かではありませんが、
> 桁落ちとか?従来と少し変わってしまっているところがあるのかな?
>
>  取り敢えず、しばらくは様子を見て、改善されたら?改善方法が見つかったら、改めて
> 投稿させて頂きたいと思います・・・・(覚えてるかな?w)
>
> 以上
>
>
> On 2016/10/13 15:18, 斎藤 wrote:
>>
>> いわさき様
>>
>>  さいとうです。
>>  ご意見いただきありがとうございます。
>>
>>  結果的に、gdal_rasterize.exeを変更しました。
>>  具体的には、障害が発生していたgdal_rasterize.exeを古いバージョンのgdal_rasterize.exeに戻しました。
>>  障害が発生したバージョン:2.1.1
>>  正常に動作したバージョン:1.11.3
>>
>>  障害が確認されたのは、弊方における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 斎藤 <[hidden email]>:
>>>>
>>>> 田原様
>>>>
>>>>  斎藤です。
>>>>  ご意見いただきありがとうございます。
>>>>
>>>>  異常な結果を示す環境も、正常な結果を示す環境も64Bit版でした。
>>>>  また、今回、正常に動作していたマシンでOSGeo4W64のバージョン更新を行いました
>>>> ところ、ズレてしまう現象が発生しました。
>>>>
>>>>  現在、正常動作をしているマシンからOSGeo4Wのディレクトリ以下をコピーして、
>>>> 障害が発生しているマシンへ部分コピーを行いながら検証してみることにしています。
>>>>  gdal_rasterize.exeだけコピーしても現象は変わらなかったので、何か他の物に影響されて
>>>> いると見ています。
>>>>
>>>>  結果は別途報告させて頂きます。
>>>>
>>>>  ちなみに、shapeファイルは全く同じファイルを使用して検証しています。
>>>>
>>>> 以上
>>>>
>>>>
>>>>
>>>> On 2016/10/08 17:05, DaisukeTahara wrote:
>>>>>
>>>>>
>>>>> 田原@HCCです。
>>>>>
>>>>> gdalに限りませんが、32bit/64bit版の違いで動作異なる事あったりしましたので
>>>>> その違いだったりしないでしょうか。
>>>>>
>>>>> とりあえずテスト比較したgdalのバージョン、マシンのOS/CPUなんかを
>>>>> 比較提示してみると原因特定しやすいかなと思います。
>>>>>
>>>>> ---
>>>>> Centerが同じですけど四隅範囲がだいぶ違うので、ズレているほうの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
>>>>>>
>>>>>> 手持ちのマシン5台で確認したところ、内3台でズレることなく生成した画像は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
>>>>>> [hidden email]
>>>>>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>>>>>
>>>>>
>>>>> _______________________________________________
>>>>> OSGeoJapan-discuss mailing list
>>>>> [hidden email]
>>>>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>>>>
>>>>
>>>>
>>>> --
>>>> /////////////////////////////////////////////////////////
>>>> 株式会社 エムエスケー
>>>> 斎藤 直正
>>>> TEL:0598-51-7471 / FAX:0598-52-4849
>>>> E-mail:[hidden email]  URL:http://www.msk-web.co.jp
>>>> /////////////////////////////////////////////////////////
>>>>
>>>>
>>>> _______________________________________________
>>>> OSGeoJapan-discuss mailing list
>>>> [hidden email]
>>>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>>>
>>>
>>>
>>
>>
>
>
> --
> /////////////////////////////////////////////////////////
> 株式会社 エムエスケー
> 斎藤 直正
> TEL:0598-51-7471 / FAX:0598-52-4849
> E-mail:[hidden email]  URL:http://www.msk-web.co.jp
> /////////////////////////////////////////////////////////
>
> _______________________________________________
> OSGeoJapan-discuss mailing list
> [hidden email]
> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss



--
岩崎 亘典
_______________________________________________
OSGeoJapan-discuss mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Re: gdal_rasterizeで生成した画像がズレる。

斎藤 直正
いわさき様

 こちらで、新たにデータを作って確認してみました。
 ズレている現象の画像を2パターン以下の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 斎藤 [hidden email]:
さいとうです。

 連投で失礼します。

 仮想でWindows 10を真っ新な状態で立ち上げて、OSGeo4W64をインストールしました。
 OSGeo4WのShellで動作を確認しましたが、同様の現象を確認しました。

 25mメッシュのSHPを使用したのですが、ギリギリの分解能だとズレが生じる様でした。
 分解能を細かくすれば現象は収まるのですが・・・処理に時間が掛かってしまうのが厄介です(^^ゞ

 細かなメッシュデータでの確認までは行っていないので、原因は定かではありませんが、
桁落ちとか?従来と少し変わってしまっているところがあるのかな?

 取り敢えず、しばらくは様子を見て、改善されたら?改善方法が見つかったら、改めて
投稿させて頂きたいと思います・・・・(覚えてるかな?w)

以上


On 2016/10/13 15:18, 斎藤 wrote:
いわさき様

 さいとうです。
 ご意見いただきありがとうございます。

 結果的に、gdal_rasterize.exeを変更しました。
 具体的には、障害が発生していたgdal_rasterize.exeを古いバージョンのgdal_rasterize.exeに戻しました。
 障害が発生したバージョン:2.1.1
 正常に動作したバージョン:1.11.3

 障害が確認されたのは、弊方における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 斎藤 [hidden email]:
田原様

 斎藤です。
 ご意見いただきありがとうございます。

 異常な結果を示す環境も、正常な結果を示す環境も64Bit版でした。
 また、今回、正常に動作していたマシンでOSGeo4W64のバージョン更新を行いました
ところ、ズレてしまう現象が発生しました。

 現在、正常動作をしているマシンからOSGeo4Wのディレクトリ以下をコピーして、
障害が発生しているマシンへ部分コピーを行いながら検証してみることにしています。
 gdal_rasterize.exeだけコピーしても現象は変わらなかったので、何か他の物に影響されて
いると見ています。

 結果は別途報告させて頂きます。

 ちなみに、shapeファイルは全く同じファイルを使用して検証しています。

以上



On 2016/10/08 17:05, DaisukeTahara wrote:
田原@HCCです。

gdalに限りませんが、32bit/64bit版の違いで動作異なる事あったりしましたので
その違いだったりしないでしょうか。

とりあえずテスト比較したgdalのバージョン、マシンのOS/CPUなんかを
比較提示してみると原因特定しやすいかなと思います。

---
Centerが同じですけど四隅範囲がだいぶ違うので、ズレているほうの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

手持ちのマシン5台で確認したところ、内3台でズレることなく生成した画像は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
[hidden email]
http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
_______________________________________________
OSGeoJapan-discuss mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss


_______________________________________________
OSGeoJapan-discuss mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss



_______________________________________________
OSGeoJapan-discuss mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
Loading...