I can confirm that gdal_retile.py also exhibits this behavior. I am testing the shift to negative longitude now to see if that "fixes" the issue.<div><br></div><div>Should I also post to the ticket concerning the persistance of this error in gdal_retile?</div>
<div><br></div><div>Jay<br><br><div class="gmail_quote">On Wed, Feb 22, 2012 at 9:48 AM, J Luis <span dir="ltr"><<a href="mailto:jmfluis@gmail.com">jmfluis@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br><br><div class="gmail_quote"><div class="im">On Wed, Feb 22, 2012 at 12:54 AM, Jay L. <span dir="ltr"><<a href="mailto:jzl5325@psu.edu" target="_blank">jzl5325@psu.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
List, <div><br></div><div>I am attempting to reproject 8 gtiffs, each in a local equirectangular projection, back to gcs for use in geoserver. I immediately went to gdalwarp to perform this task and am having an issue. The input files are stored in positive longitude, while the output SRS is in negative longitude. The "shift" to negative longitude is a primary reason for wishing to perform the warp.</div>
<div><br></div><div>The input images which cover 0-90 north and 0-90 south are warping correctly. The others are not.</div><div><br></div><div>The other warped output images have an extent of -180 180, when the should cover 90 degrees. The images which remain in positive longitude are drawing correctly, but have the remained of the extent filled in black. Meaning the image which spans 90-180n and 90-180s renders, but has an extent of -180 to 180. The input image which are shifted to negative longitude have a global extent and are completely black.</div>
</blockquote><div><br></div></div><div>Jay,</div><div><br></div><div>I crossed what seams to be the same problem a couple of days ago and opened this ticket</div><div><a href="http://trac.osgeo.org/gdal/ticket/4523" target="_blank">http://trac.osgeo.org/gdal/ticket/4523</a>
</div><div>With a bit more of experimenting after that report I found that the, at least with global grids, the problem is when they span [0 360]. If they are in [-180 180] gdalwarp works fine.</div><div>if you are interested, I have one Mirone solution to warp from [0 360] to [-180 180].</div>
<div><br></div><div>Joaquim</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
<div><br></div><div>I have tried the following without success:</div><div><br></div><div><div><span style="background-color:rgb(204,204,204)">gdalwarp -t_srs moon2000.prf input output_warped.tif</span></div></div><div><br>
</div><div><span style="background-color:rgb(204,204,204)">gdalwarp -s_srs <copy of the gdalinfo srs information in a prf> -t_srs moon2000.prf input.tif output.tif</span></div><div><br></div><div><span style="background-color:rgb(204,204,204)">gdalwarp -t_srs moon2000.prf -te <desired extent as xmin ymin xmax ymax> input.tif output.tif</span></div>
<div><span style="background-color:rgb(204,204,204)"><br></span></div><div><div>Where moon2000.prf is:</div><div><br></div><div><div><span style="background-color:rgb(204,204,204)">GEOGCS["Moon 2000",</span></div>
<div><span style="background-color:rgb(204,204,204)"> DATUM["D_Moon_2000",</span></div><div><span style="background-color:rgb(204,204,204)"> SPHEROID["Moon_2000_IAU_IAG",1737400.0,0.0]],</span></div>
<div><span style="background-color:rgb(204,204,204)"> PRIMEM["Greenwich",0],</span></div><div><span style="background-color:rgb(204,204,204)"> UNIT["Decimal_Degree",0.0174532925199433]]</span></div>
</div></div><div><span style="background-color:rgb(204,204,204)"><br></span></div><div><span>Any insight would be greatly appreciated.</span></div><div><span>Jay</span></div>
<div><br></div><div>Here is the projection info as output by gdalinfo:</div><div><br></div><div><div><span style="background-color:rgb(204,204,204)">Coordinate System is:</span></div><div><span style="background-color:rgb(204,204,204)">PROJCS["EQUIRECTANGULAR MOON",</span></div>
<div><span style="background-color:rgb(204,204,204)"> GEOGCS["GCS_MOON",</span></div><div><span style="background-color:rgb(204,204,204)"> DATUM["D_MOON",</span></div><div><span style="background-color:rgb(204,204,204)"> SPHEROID["MOON_localRadius",1737400,0]],</span></div>
<div><span style="background-color:rgb(204,204,204)"> PRIMEM["Reference_Meridian",0],</span></div><div><span style="background-color:rgb(204,204,204)"> UNIT["degree",0.0174532925199433]],</span></div>
<div><span style="background-color:rgb(204,204,204)"> PROJECTION["Equirectangular"],</span></div><div><span style="background-color:rgb(204,204,204)"> PARAMETER["latitude_of_origin",0],</span></div>
<div><span style="background-color:rgb(204,204,204)"> PARAMETER["central_meridian",0],</span></div><div><span style="background-color:rgb(204,204,204)"> PARAMETER["standard_parallel_1",0],</span></div>
<div><span style="background-color:rgb(204,204,204)"> PARAMETER["false_easting",0],</span></div><div><span style="background-color:rgb(204,204,204)"> PARAMETER["false_northing",0],</span></div><div>
<span style="background-color:rgb(204,204,204)"> UNIT["metre",1,</span></div><div><span style="background-color:rgb(204,204,204)"> AUTHORITY["EPSG","9001"]]]</span></div><div><span style="background-color:rgb(204,204,204)">Origin = (100.000000000000000,1819400.000000000000000)</span></div>
<div><span style="background-color:rgb(204,204,204)">Pixel Size = (100.000000000000000,-100.000000000000000)</span></div></div><div><br></div><div>The coordinate information:</div><div><div><span style="background-color:rgb(204,204,204)">Corner Coordinates:</span></div>
<div><span style="background-color:rgb(204,204,204)">Upper Left ( 100.000, 1819400.000) ( 0d 0'11.87"E, 59d59'59.88"N)</span></div><div><span style="background-color:rgb(204,204,204)">Lower Left ( 100.000, -100.000) ( 0d 0'11.87"E, 0d 0'11.87"S)</span></div>
<div><span style="background-color:rgb(204,204,204)">Upper Right ( 2729200.000, 1819400.000) ( 90d 0'11.69"E, 59d59'59.88"N)</span></div><div><span style="background-color:rgb(204,204,204)">Lower Right ( 2729200.000, -100.000) ( 90d 0'11.69"E, 0d 0'11.87"S)</span></div>
<div><span style="background-color:rgb(204,204,204)">Center ( 1364650.000, 909650.000) ( 45d 0'11.78"E, 29d59'54.00"N)</span></div></div><div><br></div><div>Oddly the central meridian and standard_parallel_1 parameters do not change with each input file. Since I know that gdal can see the s_srs from within the files I have tried using:</div>
<div><br></div><div><div><br></div></div>
<br></div></div><div class="im">_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br></div></blockquote></div><br>
</blockquote></div><br></div>