<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
color:black;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
pre
{mso-style-priority:99;
mso-style-link:"HTML Preformatted Char";
margin:0in;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New";
color:black;}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
mso-margin-top-alt:auto;
margin-right:0in;
mso-margin-bottom-alt:auto;
margin-left:0in;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
color:black;}
span.HTMLPreformattedChar
{mso-style-name:"HTML Preformatted Char";
mso-style-priority:99;
mso-style-link:"HTML Preformatted";
font-family:Consolas;
color:black;}
span.EmailStyle21
{mso-style-type:personal;
font-family:"Calibri",sans-serif;
color:windowtext;}
span.EmailStyle22
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
/* List Definitions */
@list l0
{mso-list-id:129564708;
mso-list-template-ids:567861862;}
@list l1
{mso-list-id:168519960;
mso-list-template-ids:739307374;}
@list l1:level1
{mso-level-tab-stop:.5in;
mso-level-number-position:left;
text-indent:-.25in;}
@list l1:level2
{mso-level-tab-stop:1.0in;
mso-level-number-position:left;
text-indent:-.25in;}
@list l1:level3
{mso-level-tab-stop:1.5in;
mso-level-number-position:left;
text-indent:-.25in;}
@list l1:level4
{mso-level-tab-stop:2.0in;
mso-level-number-position:left;
text-indent:-.25in;}
@list l1:level5
{mso-level-tab-stop:2.5in;
mso-level-number-position:left;
text-indent:-.25in;}
@list l1:level6
{mso-level-tab-stop:3.0in;
mso-level-number-position:left;
text-indent:-.25in;}
@list l1:level7
{mso-level-tab-stop:3.5in;
mso-level-number-position:left;
text-indent:-.25in;}
@list l1:level8
{mso-level-tab-stop:4.0in;
mso-level-number-position:left;
text-indent:-.25in;}
@list l1:level9
{mso-level-tab-stop:4.5in;
mso-level-number-position:left;
text-indent:-.25in;}
ol
{margin-bottom:0in;}
ul
{margin-bottom:0in;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body bgcolor=white lang=EN-AU link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='color:windowtext;mso-fareast-language:EN-US'>Mike,<o:p></o:p></span></p><p class=MsoNormal><span style='color:windowtext;mso-fareast-language:EN-US'>So everything in the capture below of the Raster calculator should work? I cannot find any output.<o:p></o:p></span></p><p class=MsoNormal><span style='color:windowtext;mso-fareast-language:EN-US'>I’m in Perth, Australia so home now and away from the office and QGIS, but would like to get it working tomorrow.<o:p></o:p></span></p><p class=MsoNormal><span style='color:windowtext;mso-fareast-language:EN-US'>How do apply a grayscale to a raster? As I said, I’m new to QGIS<o:p></o:p></span></p><div><p class=MsoNormal><span style='color:windowtext'>Chas Sheen<o:p></o:p></span></p></div><p class=MsoNormal><span style='color:windowtext;mso-fareast-language:EN-US'><o:p> </o:p></span></p><div><div style='border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0in 0in 0in'><p class=MsoNormal><b><span lang=EN-US style='color:windowtext'>From:</span></b><span lang=EN-US style='color:windowtext'> Micha Silver [mailto:tsvibar@gmail.com] <br><b>Sent:</b> Wednesday, January 17, 2018 8:39 PM<br><b>To:</b> Chas Sheen <chassheengeo@gmail.com>; qgis-user@lists.osgeo.org<br><b>Subject:</b> Re: [Qgis-user] Generating black white from RGB landsat using Raster Map calculator?<o:p></o:p></span></p></div></div><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal style='margin-bottom:12.0pt'>Hello Chas<o:p></o:p></p><div><p class=MsoNormal>On 01/17/2018 01:52 PM, Chas Sheen wrote:<o:p></o:p></p></div><blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><p class=MsoNormal><span style='color:windowtext;mso-fareast-language:EN-US'>Richard,</span><o:p></o:p></p><p class=MsoNormal><span style='color:windowtext;mso-fareast-language:EN-US'> </span><o:p></o:p></p><p class=MsoNormal><span style='color:windowtext;mso-fareast-language:EN-US'>I should have made my initial query clearer.</span><o:p></o:p></p><p class=MsoNormal><span style='color:windowtext;mso-fareast-language:EN-US'>I am trying to output a black and white geotiff.</span><o:p></o:p></p></blockquote><p class=MsoNormal>It's still not clear what you are asking. AFAIK, there's no such thing as a "black and white geotiff". You can apply any color table you want to any raster. You can apply a gray scale to any of the landsat bands.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>If you are trying to create a new raster from the three RGB bands using some expression, that should work. You can use either of three tools: the built in QGIS raster calculator, or thru the Processing toolbox there's SAGA Raster calculus and the GRASS module r.mapcalculator. <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Again, none of this has anything to do with coloring. After running the raster calculator, you can color the result with a gray scale color table, a red color table, or any other color ramp you choose.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><p class=MsoNormal><span style='color:windowtext;mso-fareast-language:EN-US'>Have split the original (landsat) image into RGB images.</span><o:p></o:p></p><p class=MsoNormal><span style='color:windowtext;mso-fareast-language:EN-US'>When I apply the expression in the capture below, and hit run, nothing is saved to the results directory.</span><o:p></o:p></p><p class=MsoNormal><span style='color:windowtext;mso-fareast-language:EN-US'>I’m assuming that my expression or something else in the raster calculator is wrong?</span><o:p></o:p></p><p class=MsoNormal><span style='color:windowtext;mso-fareast-language:EN-US'>The final product that I’m after is a xyz grid where z is the black-white amplitude…which I think should be relatively easy (?) once I have the black-white image.</span><o:p></o:p></p></blockquote><p class=MsoNormal>After you running the calculator expression, then you can save the result in many formats, including gridded XYZ, ArcINFO ASCII, or many others, using the gdal_translate tool. This is under the Raster->Conversion menu.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><p class=MsoNormal><span style='color:windowtext;mso-fareast-language:EN-US'> </span><o:p></o:p></p><div><p class=MsoNormal><span style='color:windowtext'>Chas Sheen</span><o:p></o:p></p><p class=MsoNormal><span style='color:windowtext'>0448017282</span><o:p></o:p></p></div><p class=MsoNormal><span style='color:windowtext;mso-fareast-language:EN-US'> </span><o:p></o:p></p><div><div style='border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0in 0in 0in'><p class=MsoNormal><b><span lang=EN-US style='color:windowtext'>From:</span></b><span lang=EN-US style='color:windowtext'> Richard McDonnell [<a href="mailto:richard.mcdonnell@opw.ie">mailto:richard.mcdonnell@opw.ie</a>] <br><b>Sent:</b> Wednesday, January 17, 2018 7:13 PM<br><b>To:</b> Chas Sheen <a href="mailto:chassheengeo@gmail.com"><chassheengeo@gmail.com></a>; <a href="mailto:qgis-user@lists.osgeo.org">qgis-user@lists.osgeo.org</a><br><b>Subject:</b> Re: [Qgis-user] Generating black white from RGB landsat using Raster Map calculator?</span><o:p></o:p></p></div></div><p class=MsoNormal> <o:p></o:p></p><p class=MsoNormal>Chas,<br>Two options for splitting RGB into separate images in QGIS.<o:p></o:p></p><ol start=1 type=1><li class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;mso-list:l1 level1 lfo3'>Semi-Automatic Classification Plugin, Preprocessing, Split Raster Bands<o:p></o:p></li><li class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;mso-list:l1 level1 lfo3'>Grass r-rgb (Processing Toolbox)<o:p></o:p></li></ol><p class=MsoNormal style='margin-bottom:12.0pt'>Any one of these tools will split the Raster for you easy enough<br>Regards,<br><br>Richard<o:p></o:p></p><div><p class=MsoNormal>On 17/01/2018 00:48, Chas Sheen wrote:<o:p></o:p></p></div><blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><div><div><p class=MsoNormal><span style='font-size:9.5pt'>I split the image into R, G, B</span><o:p></o:p></p></div><div><p class=MsoNormal><span style='font-size:9.5pt'>Worked out in QGIS wher the Raster map calculator is.</span><o:p></o:p></p></div><div><p class=MsoNormal><span style='font-size:9.5pt'>Run...several times</span><o:p></o:p></p></div><div><p class=MsoNormal><span style='font-size:9.5pt'>No result foud in Results Directory. It runs quickly then refreshes inputs and logs to blanks.</span><o:p></o:p></p></div><div><p class=MsoNormal><span style='font-size:9.5pt'>I'm new to QGIS....so have probably left out something basic.</span><o:p></o:p></p></div><div><p class=MsoNormal><span style='font-size:9.5pt'>Any clues</span><o:p></o:p></p></div><div><p class=MsoNormal><span style='font-size:9.5pt'> </span><o:p></o:p></p></div><div><p class=MsoNormal><span style='font-size:9.5pt'><img border=0 width=616 height=720 style='width:6.4166in;height:7.5in' id="_x0000_i1025" src="cid:image001.png@01D38FDA.C7A79ED0" alt="Displaying image.png"></span><o:p></o:p></p></div></div><p class=MsoNormal><br><br><br><br><o:p></o:p></p><pre>_______________________________________________<o:p></o:p></pre><pre>Qgis-user mailing list<o:p></o:p></pre><pre><a href="mailto:Qgis-user@lists.osgeo.org">Qgis-user@lists.osgeo.org</a><o:p></o:p></pre><pre>List info: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user">https://lists.osgeo.org/mailman/listinfo/qgis-user</a><o:p></o:p></pre><pre>Unsubscribe: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user">https://lists.osgeo.org/mailman/listinfo/qgis-user</a><o:p></o:p></pre></blockquote><p class=MsoNormal> <o:p></o:p></p><div><p class=MsoNormal>-- <br><b><u>Richard McDonnell</u></b><br><b><span style='font-size:7.5pt'>GIS Specialist PgD GIS AssocSCSI</span></b><br><b><span style='font-size:7.5pt'>OPW FRM Data Management</span></b><br><b><span style='font-size:7.5pt'>52 Stephens Green, Dublin 2.</span></b><br><span style='font-size:7.5pt'>TEL: 01 6476543</span><o:p></o:p></p></div><p class=MsoNormal><a href="http://www.opw.ie/"><span style='text-decoration:none'><img border=0 width=405 height=225 style='width:4.2187in;height:2.3437in' id="_x0000_i1026" src="http://www.opw.ie/en/media/opw-logo-december-2016.gif"></span></a><o:p></o:p></p><p class=MsoNormal><br><br>Email Disclaimer: <a href="http://www.opw.ie/en/disclaimer/">http://www.opw.ie/en/disclaimer/</a> <o:p></o:p></p><p class=MsoNormal><br><br><br><o:p></o:p></p><pre>_______________________________________________<o:p></o:p></pre><pre>Qgis-user mailing list<o:p></o:p></pre><pre><a href="mailto:Qgis-user@lists.osgeo.org">Qgis-user@lists.osgeo.org</a><o:p></o:p></pre><pre>List info: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user">https://lists.osgeo.org/mailman/listinfo/qgis-user</a><o:p></o:p></pre><pre>Unsubscribe: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user">https://lists.osgeo.org/mailman/listinfo/qgis-user</a><o:p></o:p></pre></blockquote><p class=MsoNormal><br><br><o:p></o:p></p><pre>-- <o:p></o:p></pre><pre>Micha Silver<o:p></o:p></pre><pre>Ben Gurion Univ.<o:p></o:p></pre><pre>Sde Boker, Remote Sensing Lab<o:p></o:p></pre><pre>cell: +972-523-665918<o:p></o:p></pre></div></body></html>