[OSGeoJapan-discuss] gdal_rasterizeで生成した画像がズレる。
Nobusuke Iwasaki
wata909 ＠ gmail.com
2016年 10月 12日 (水) 00:16:25 PDT
齋藤様
いわさきです。
> 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
>
>
>
> --
> /////////////////////////////////////////////////////////
> 株式会社 エムエスケー
> 斎藤 直正
> TEL:0598-51-7471 / FAX:0598-52-4849
> E-mail:nsaito ＠ msk-web.co.jp URL:http://www.msk-web.co.jp
> /////////////////////////////////////////////////////////
>
>
> _______________________________________________
> OSGeoJapan-discuss mailing list
> OSGeoJapan-discuss ＠ lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
--
岩崎 亘典
OSGeoJapan-discuss メーリングリストの案内