<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">If you can interpolate and organize the data as follows, I can incorporate it fairly painlessly;<div><br></div><div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">/* not real data */</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">sensor= Spot 2</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">band= 3 hvr2</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Lower edge wavelength = 0.715</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Upper edge wavelength = 0.9175</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Spectral Response =</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">.0000, .0001, .0001, .0002, .0004, .0005, .0007, .0009</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">.0012, .0014, .0017, .0021, .0024, .0028, .0032, .0036</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">.0040, .0045, .0054, .0065, .0080, .0112, .0170, .0245</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">.0330, .0450, .0740, .0990, .1520, .2140, .3050, .4200</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">.5970, .6750, .7900, .8270, .8400, .8408, .9010, .9050</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">.9100, .9140, .9050, .9260, .9470, .9670, .9760, .9840</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">.9710, .9790, .9880,1.0000, .9902, .9840, .9780, .9710</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">.9470, .9510, .9450, .9380, .9340, .9300, .9360, .9420</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">.9140, .8970, .7980, .7000, .5970, .4610, .3700, .2630</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">.1730, .1259, .0888, .0619, .0450, .0346, .0267, .0207</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">.0160, .0120, .0086, .0060, .0040, .0025, .0013, .0003</span></font></div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">.0000</span></font></div><div><br></div><div>Note that;</div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">(Upper edge wavelength - lower edge wavelength)/2.5nm</span></font> should equal the number of elements in the Spectral response.</div><div><font class="Apple-style-span" face="'Courier New'" size="3"><span class="Apple-style-span" style="font-size: 12px; ">(Lower edge wavelength - 250nm)/2.5nm</span></font> should be a integer. ie 250, 257.5 & 260 not 251,258.5 & 261.</div><div><br></div><div>Linear interpolation is probably the safest, but bilinear should be fine. I'd plot the "before" and "after" one on top of the other to make sure there are no major deviations between the two. If you're familiar with, I found using Octave/matlab and interp1 to greatly simplify things. In the mean time I'll try to figure out GeomCond.cpp so I can simply submit a patch.</div><div><br></div><div>Cheers,</div><div><br></div><div>Mike</div><div><br></div></div><div><div>On 28-Jul-09, at 5:23 PM, andrew haywood wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Hi Mike</div> <div> </div> <div>thanks for this supporting information. I will start the interpolation. Is there any prefered method for this? Can I use a bilinear approach.</div> <div>I think this is as far as I can take it though - as I have no understanding of lwave.cpp and the file construct. Will it be ok to post a txt file on the list with the interpretation results?</div> <div> </div> <div>Cheers</div> <div>Andrew<br><br> </div> <div><span class="gmail_quote">On 7/29/09, <b class="gmail_sendername">Michael Perdue</b> <<a href="mailto:michael_perdue@yahoo.ca">michael_perdue@yahoo.ca</a>> wrote:</span> <blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid"> <div style="WORD-WRAP: break-word">Hi Andy, <div>Yes, it looks like the right info for the SPOT sensor platforms. Iwave.cpp expects the data to be in 2.5nm increments so you will have to interpolate the data to the correct interval. It's a bit tedious, but is straight forward. I might be able to dedicate some time to it, but won't have any till the weekend now.</div> <div>Cheers,</div> <div>Mike</div> <div><span class="e" id="q_122c1e743df6d678_1"> <div><br> <div> <div>On 27-Jul-09, at 10:35 PM, andrew haywood wrote:</div><br> <blockquote type="cite"> <div>Hi All,</div> <div> </div> <div>Im just adding to this discusssion. I would be interested in using i.atcorr with SPOT 2, 3, 4 and 5. I have found a file on the <a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.spotimage.com/web/en/555-spot-images.php" target="_blank">http://www.spotimage.com/web/en/555-spot-images.php</a> website that contains the spectral sensitivity data. Is this the same as the spectral profile? It looks like the data is in 5nm increments. Do i need to create 2.5nm increments or can i tell the function its in 5nm or 10nm increments?</div> <div>Any help would be appreciated.</div> <div>thanks</div> <div>Andy</div> <div> </div> <div><br> </div> <div><span class="gmail_quote">On 7/28/09, <b class="gmail_sendername">Michael Perdue</b> <<a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:michael_perdue@yahoo.ca" target="_blank">michael_perdue@yahoo.ca</a>> wrote:</span> <blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Hi Markus,<br><br>Here is a text file with everything organized as I think it should be. I didn't add the this Iwave.cpp to submit a patch as I wasn't really sure how to modify GeomCond.cpp.<br> <br>Cheers,<br>Mike<br><span><br>On 26-Jul-09, at 1:34 PM, Markus Neteler wrote:<br></span> <div><span> <blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid"><br>Great - the header of Iwave.cpp contains instructions how to add<br>a new sensor like this.<br><br>thanks<br> Markus<br></blockquote></span></div><br>_______________________________________________<br>grass-user mailing list<br><a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br> <a onclick="return top.js.OpenExtLink(window,event,this)" href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br><br><br clear="all"></blockquote> </div><br></blockquote></div><br> </div></span></div></div></blockquote></div><br></blockquote></div><br></body></html>