$B4d:j$5$s!"$3$s$K$A$O(B<br><br>$B$d$jJ}$r8+$D$1$^$7$?!#$"$^$jBg$-$$2hA|$K$O8~$$$F$$$J$$$+$b$7$l$^$;$s$,!"(BR$B$N(Braster$B%i%$%V%i%j!<$r;H$C$F(Bstack$B$7$^$7$?!#0J2<$K(BR$B$N%3!<%I$rIU$1$^$7$?!#$b$7B>$KNI$5$=$&$JJ}K!$,$"$l$P65$($F2<$5$$!#$o$?$7$,(BGDAL$B$r;H$&$H$$$C$?$N$O(BGDAL$B$G$+$/%l%$%d!<$r%P%s%I%*%V%8%'%/%H$H$7$FFI$_9~$_!"3F%l%$%d!<$N%5%$%:$rB7$($F$+$i=q$-$@$9$H$$$&<j=g$G$7$?$,BgJQ$=$&$J$N$G$9$G$K$"$k%D!<%k$rC5$7$F$$$^$7$?!#(B<br><br>
$B0J2<$,%3!<%I$G$9!#(B<br><br># This code stack 3 band tif image and 1 band dem together and create 4 band fusion image.<br># 20110906 Hiroo Imaki (<a href="mailto:hiroo@angeli.org">hiroo@angeli.org</a>)<br><br>library(raster)<br>library(rgdal)<br>
<br># check available file format<br>writeFormats()<br><br># read vege data (3 bands)<br>veg <- "G:\\Precision\\imakihi\\project\\snohomish\\gisdata\\vegetation_mapping\\qwuloolt_veg_3.tif"<br>veg.r <- stack(veg)<br>
nlayers(veg.r)<br>veg.r<br><br># read elev data (one band)<br>elev <- "G:\\Precision\\imakihi\\project\\snohomish\\gisdata\\vegetation_mapping\\hh_be_3.tif"<br>elev.r <- raster(elev)<br>nlayers(elev.r)<br>
elev.r<br><br># crop image to fit each other<br>veg.crop <- crop(veg.r, elev.r)<br>elev.crop <- crop(elev.r, veg.crop)<br>veg.crop<br>elev.crop<br><br># resample elevation image (if resolutions are different each other)<br>
#elev.resample <- resample(elev.crop, veg.r, method='bilinear')<br><br># strech values from 0 to 255 (you don't need this if you are not interested in streaching values)<br>elev.min <- cellStats(elev.crop, stat="min")<br>
elev.max <- cellStats(elev.crop, stat="max")<br><a href="http://elev.int">elev.int</a> <- ceiling((elev.crop - elev.min) * 1000 / (elev.max - elev.min))<br><a href="http://elev.int">elev.int</a><br><br><br>
# stack elevation layer on top of vegetation layers.<br>veg.stack <- stack(<a href="http://elev.int">elev.int</a>, veg.crop)<br>veg.stack<br><br># output the stacked image<br>writeRaster(veg.stack, "G:\\Precision\\imakihi\\project\\snohomish\\gisdata\\vegetation_mapping\\wetland_stack3.tif", overwrite=TRUE, datatype="INT2U")<br>
#writeRaster(elev.crop, "c:\\temp\\stacked\\elev.tif", overwrite=TRUE)<br><br># plot elevation and vegetation layers<br>plot(<a href="http://elev.int">elev.int</a>)<br>plotRGB(veg.r, r=1, g=2, b=3)<br><br>$B$$$^$-(B<br>
<br><br><div class="gmail_quote">2011/9/6 Nobusuke Iwasaki <span dir="ltr"><<a href="mailto:wata909@gmail.com">wata909@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
$B$$$^$-$5$s!$$_$J$5$s(B<br>
<br>
$B$$$o$5$-$G$9!#(B<br>
$B2<5-$N7o!$$J$s$+$$$$%"%$%G%#%"$O$"$j$^$7$?$G$7$g$&$+!)(B<br>
$B;d$b6=L#$,$"$j$^$9!#(B<br>
<br>
$B$A$J$_$K!$(BGDAL$B$r;H$&$H$$$&$N$O!$%P!<%A%c%k%i%9%?!<!J(BVRT$B!K$r;H$&$H$$$&$3$H$G$7$g$&$+!)(B<br>
<br>
<br>
2011$BG/(B8$B7n(B30$BF|(B11:13 Hiroo Imaki <<a href="mailto:hiroo@angeli.org">hiroo@angeli.org</a>>:<br>
<div><div></div><div class="h5">> $B%3%_%e%K%F%#!<$N3'MM!"(B<br>
><br>
> $B8=:_!"%P!<%/%l!<$G:n$C$?%$%a!<%8%;%0%a%s%F!<%7%g%s%=%U%H%&%'%"$G2hA|$N2r@O$r$7$F$$$k$N$G$9$,!"(BTIFF$B$GJ]B8$5$l$F$$$k9R6u<L??!J#3%P%s%I!K$H(BLIDAR$B$NI89b%G!<%?!J#1%P%s%I!K$r$R$H$D$N2hA|%G!<%?!J#4%P%s%I$H$7$F!KJ]B8$7$?$$$N$G$9$,$I$N$h$&$JJ}K!$,$$$$$N$+$$$^$$$A;W$$$D$+$J$$$N$G$I$J$?$+$4B8$8$NJ}$,$$$l$P65$($FM_$7$$$N$G$9$,!#(B<br>
><br>
> $B;THN$N%=%U%H%&%'%"!<$G$N$d$jJ}$b$$$$$N$G$9$,!"%*!<%W%s%=!<%9$N%D!<%k$r;H$C$F=PMh$?$i$H;W$C$F$$$^$9!#:#$N$H$3$m(BGDAL$B$H%Q%$%=%s$r;H$C$FCOL#$K$d$kJ}K!$O9M$($i$l$k$N$G$9$,!"$b$&$A$g$C$H8zN($N$h$$J}K!$,$"$k$N$G$O$H;W$$$^$7$F!#(B<br>
><br>
> $B$h$m$7$/$*4j$$$7$^$9!#(B<br>
><br>
> $B:#LZ(B<br>
><br>
> --<br>
> Hiroo Imaki<br>
> <a href="mailto:hiroo@angeli.org">hiroo@angeli.org</a><br>
> <a href="http://www.geopacific.org" target="_blank">http://www.geopacific.org</a><br>
><br>
</div></div>> _______________________________________________<br>
> OSGeoJapan-discuss mailing list<br>
> <a href="mailto:OSGeoJapan-discuss@lists.osgeo.org">OSGeoJapan-discuss@lists.osgeo.org</a><br>
> <a href="http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss" target="_blank">http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss</a><br>
><br>
><br>
<font color="#888888"><br>
<br>
<br>
--<br>
$B4d:j(B $BOKE5(B<br>
</font></blockquote></div><br><br clear="all"><br>-- <br>Hiroo Imaki<br><a href="mailto:hiroo@angeli.org">hiroo@angeli.org</a><br><a href="http://www.geopacific.org">http://www.geopacific.org</a><br>