[OSGeoJapan-discuss] Geotiffについて

["臼井"]

朝日様、今木様

お返事がおそくなってすいません

臼井です

教えていただいた方法でやってみましたところ

上手くいきそうです、また

QGISのラスタ計算機というものがあること自体しりませんでした

非常に、参考になりました

ありがとうございます

また、よろしくお願いします


----- Original Message -----
Subject: Re: [OSGeoJapan-discuss] Geotiffについて
Date: Fri, 31 Oct 2014 05:59:42 -0400
From: Hiroo Imaki <hiroo ＠ angeli.org>
To: "OSGeoJapan-discuss ＠ lists.osgeo.org" <osgeojapan-discuss ＠ lists.osgeo.org>

臼井さん、朝日さん、こんにちは、

もしかしてよくわかってないかもしれませんが、以下のポスティングが役立つかもしれません。
http://gis.stackexchange.com/questions/2919/how-do-i-change-geotiff-no-data-color-to-white-transparent-rather-than-black

一旦vrt（バーチャルラスタ）に書き出して、カラーテーブルをテキストでいじくれるようにしてから、カラーテーブルを書き直し、それからgdal_mergeで結合するという方法です。

今木

2014-10-31 3:23 GMT-04:00 ASAHI Kosuke <asahi ＠ hcc.co.jp>:

> 臼井さん
>
> 朝日です。
>
> gdal_calc.py辺りを上手く使えないでしょうか？
>  http://www.gdal.org/gdal_calc.html
>
> 例えば、A画像の0: 0,0,0,255を抜き出しておく(OSGeo4wのバッチで実行したと
> して)、
>  gdal_calc.bat -A A.tif --outfile=A_mask.tif --calc=" (A<1)"
> B画像の1: 0,0,0,255を抜き出しておく、
>  gdal_calc.bat -A B.tif --outfile=B_mask.tif --calc=" (A>0)"
> としてからマージするとか。
>
> これでいいのであれば、gdal_calcの--calcオプションを上手く使うと一度で出
> 来そうな気もしますし、QGISのラスタ計算機でも出来そうな気もしますね。
>
> (2014/10/31 8:50), "臼井" wrote:
> > すいません、臼井ともうします
> >
> > 今回、２回目の投稿となります
> >
> > よろしくお願いします
> >
> >
> > gdalinfo にて A,BふたつのGeotiffファイルがあります
> >
> > A,B二つのラスタを結合して一つのラスタにしたいのですが
> >
> > 背景色が違うため、結合すると背景が黒と白色の部分ができてしまい
> >
> > 困っています。
> >
> > 作業はgdalwarlpもしくはQGIS2.4のラスタ処理の中の結合（gdal_mergebat）
> >
> > で行いました
> >
> > windows環境で作業をしています
> >
> > GDALのバージョンは1.11.0 QGIS2.4です
> >
> >
> > gdalInfoでファイルの情報をみてみると
> > Color Table (RGB with 2 entries)
> > が原因なのは分かったので、この値をどちらかに統一すれば
> > いいと思うのですが、方法がわかりません
> > こちらについて、何かいい方法を教えていただけますでしょうか？
> >
> >
> > また、できたgeotiffが未圧縮のため非常に容量がおおきいので
> > gdal_translateで小さくしていますが
> > やはり 元 737kb、468kb にくらべ 作成されたものが 9906kbと
> > 非常に大きいです
> > こちらも、他に何かいい方法があれば教えていただけますでしょうか？
> >
> >
> >
> >
> >
> >
> > 一つ目のファイル（A）
> >
> > Band 1 Block=8936x11 Type=Byte, ColorInterp=Palette
> >    Image Structure Metadata:
> >      NBITS=1
> >    Color Table (RGB with 2 entries)
> >      0: 0,0,0,255
> >      1: 255,255,255,255
> >
> > 二つ目のファイル（B）
> >
> > Band 1 Block=8809x14 Type=Byte, ColorInterp=Palette
> >    Image Structure Metadata:
> >      NBITS=1
> >    Color Table (RGB with 2 entries)
> >      0: 255,255,255,255
> >      1: 0,0,0,255
> >
> >
> >
> > _______________________________________________
> > OSGeoJapan-discuss mailing list
> > OSGeoJapan-discuss ＠ lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
> >
> >
>
>
> --
> 北海道地図株式会社 研究開発課
> 朝日 孝輔
> mailto:asahi ＠ hcc.co.jp
> Tel : 0166-62-4711
>
> _______________________________________________
> OSGeoJapan-discuss mailing list
> OSGeoJapan-discuss ＠ lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>



-- 
Hiroo Imaki
Pacific Spatial Solutions, LLC
1523 Chatham Colony Ct.
Reston, VA  20190
hiroo.imaki ＠ pacificspatial.com
http://www.pacificspatial.com
http://www.geopacific.org (GIS info site)
hiroo ＠ angeli.org (private)