<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<DEFANGED_META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=us-ascii">
<DEFANGED_META NAME="Generator" CONTENT="MS Exchange Server version 6.0.6603.0">
<TITLE>[GRASSLIST:9858] Re: GeoTiffs shifted by 360 degrees after processing in GRASS</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/rtf format -->
<P><FONT SIZE=2 FACE="Courier New">Markus, thanks for your reply.</FONT>
</P>
<P><FONT SIZE=2 FACE="Courier New">I currently use Grass 6.0.0 and it seems 6.1 is still under development. Is it possible to get the i.in.spotvgt script only and to implement it in the old version or do we have to compile the whole source code of 6.1.</FONT></P>
<P><FONT SIZE=2 FACE="Courier New">Dorit</FONT>
</P>
<BR>
<BR>
<BR>
<P><FONT SIZE=2 FACE="Courier New">On Mon, Jan 16, 2006 at 04:44:57PM +0100, Gross, Dorit (SDRN) wrote:<BR>
><I> Dear all,<BR>
</I>><I><BR>
</I>><I> I have a set of GeoTiffs that I import into GRASS to do some processing and<BR>
</I>><I> export them afterwards as GeoTiffs again. The problem now is that GRASS<BR>
</I>><I> exports the images shifted by 360 degrees to the east.<BR>
</I>><I><BR>
</I>><I><BR>
</I>><I> Corner Coordinates (before the import into GRASS):<BR>
</I>><I> Upper Left (-180.0044643, 75.0044643) (180d 0'16.07"W, 75d 0'16.07"N)<BR>
</I>><I> Lower Left (-180.0044643, -56.0044643) (180d 0'16.07"W, 56d 0'16.07"S)<BR>
</I>><I> Upper Right ( 179.9955353, 75.0044643) (179d59'43.93"E, 75d 0'16.07"N)<BR>
</I>><I> Lower Right ( 179.9955353, -56.0044643) (179d59'43.93"E, 56d 0'16.07"S)<BR>
</I>><I> Center ( -0.0044645, 9.5000000) ( 0d 0'16.07"W, 9d30'0.00"N)<BR>
</I>><I><BR>
</I>><I> Corner Coordinates (after the export from GRASS):<BR>
</I>><I> Upper Left ( 179.9955357, 75.0044643) (179d59'43.93"E, 75d 0'16.07"N)<BR>
</I>><I> Lower Left ( 179.9955357, -56.0044643) (179d59'43.93"E, 56d 0'16.07"S)<BR>
</I>><I> Upper Right ( 539.996, 75.004) (539d59'43.93"E, 75d 0'16.07"N)<BR>
</I>><I> Lower Right ( 539.996, -56.004) (539d59'43.93"E, 56d 0'16.07"S)<BR>
</I>><I> Center ( 359.996, 9.500) (359d59'43.93"E, 9d30'0.00"N)<BR>
</I>><I><BR>
</I>><I><BR>
</I>><I> Has someone of you faced a similar problem? How do I receive a proper<BR>
</I>><I> georeference when exporting from GRASS?<BR>
</I>><I><BR>
</I>><I> The input GeoTiffs are a set of converted Spot-Veg HDF files using GDAL.<BR>
</I><BR>
Did you try to use the dedicated script i.in.spotvgt:<BR>
<BR>
</FONT><A HREF="http://grass.itc.it/grass61/manuals/html61_user/i.in.spotvgt.html"><U><FONT COLOR="#0000FF" SIZE=2 FACE="Courier New">http://grass.itc.it/grass61/manuals/html61_user/i.in.spotvgt.html</FONT></U></A><BR>
<BR>
<FONT SIZE=2 FACE="Courier New">Maybe this solves the problem of the global wrap-around.<BR>
Please let us know.<BR>
<BR>
Markus<BR>
<BR>
></FONT><I><FONT SIZE=2 FACE="Courier New"> The -180.0044... degrees for the upper left X seems a bit odd, but GDAL imports<BR>
</FONT></I><FONT SIZE=2 FACE="Courier New">></FONT><I><FONT SIZE=2 FACE="Courier New"> the Spot-Veg properly with these coordinates. I checked the images in Erdas<BR>
</FONT></I><FONT SIZE=2 FACE="Courier New">></FONT><I><FONT SIZE=2 FACE="Courier New"> and they are correct. I also checked the corner coordinates of the imported<BR>
</FONT></I><FONT SIZE=2 FACE="Courier New">></FONT><I><FONT SIZE=2 FACE="Courier New"> file insight GRASS and also here the coordinates are correct. So apparently<BR>
</FONT></I><FONT SIZE=2 FACE="Courier New">></FONT><I><FONT SIZE=2 FACE="Courier New"> the shift occurs during the export from GRASS into GeoTiff. At the moment I<BR>
</FONT></I><FONT SIZE=2 FACE="Courier New">></FONT><I><FONT SIZE=2 FACE="Courier New"> solve the problem by reassigning the georeference to the GeoTiff after the<BR>
</FONT></I><FONT SIZE=2 FACE="Courier New">></FONT><I><FONT SIZE=2 FACE="Courier New"> export from GRASS, but I would be more happy with properly exported images<BR>
</FONT></I><FONT SIZE=2 FACE="Courier New">></FONT><I><FONT SIZE=2 FACE="Courier New"> right away.<BR>
</FONT></I><FONT SIZE=2 FACE="Courier New">></FONT><I><FONT SIZE=2 FACE="Courier New"><BR>
</FONT></I><FONT SIZE=2 FACE="Courier New">></FONT><I><FONT SIZE=2 FACE="Courier New"> Can someone help please?<BR>
</FONT></I><FONT SIZE=2 FACE="Courier New">></FONT><I><FONT SIZE=2 FACE="Courier New"><BR>
</FONT></I><FONT SIZE=2 FACE="Courier New">></FONT><I><FONT SIZE=2 FACE="Courier New"> Thank you,<BR>
</FONT></I><FONT SIZE=2 FACE="Courier New">></FONT><I><FONT SIZE=2 FACE="Courier New"> Dorit<BR>
</FONT></I><FONT SIZE=2 FACE="Courier New">></FONT><I><FONT SIZE=2 FACE="Courier New"><BR>
</FONT></I><BR>
<FONT SIZE=2 FACE="Courier New">--<BR>
Markus Neteler <neteler itc it> </FONT><A HREF="http://mpa.itc.it"><U><FONT COLOR="#0000FF" SIZE=2 FACE="Courier New">http://mpa.itc.it</FONT></U></A><BR>
<FONT SIZE=2 FACE="Courier New">ITC-irst - Centro per la Ricerca Scientifica e Tecnologica<BR>
MPBA - Predictive Models for Biol. & Environ. Data Analysis<BR>
Via Sommarive, 18 - 38050 Povo (Trento), Italy</FONT>
</P>
</BODY>
</HTML>