<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
On 22-02-2012 00:54, Jay L. wrote:
<blockquote
cite="mid:CA+bfmPuWDbw1C2H3474zaw7_gRUd5uC2D6GeHprNZ+cdH-gggQ@mail.gmail.com"
type="cite">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>
<br>
Jay,<br>
<br>
I crossed what seams to be the same problem a couple of days ago and
opened this ticket<br>
<a class="moz-txt-link-freetext" href="http://trac.osgeo.org/gdal/ticket/4523">http://trac.osgeo.org/gdal/ticket/4523</a><br>
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.<br>
if you are interested, I have one Mirone solution to warp from [0
360] to [-180 180].<br>
<br>
Joaquim<br>
<br>
<blockquote
cite="mid:CA+bfmPuWDbw1C2H3474zaw7_gRUd5uC2D6GeHprNZ+cdH-gggQ@mail.gmail.com"
type="cite">
<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 style="background-color:rgb(255,255,255)">Any insight
would be greatly appreciated.</span></div>
<div><span style="background-color:rgb(255,255,255)">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>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
gdal-dev mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="http://lists.osgeo.org/mailman/listinfo/gdal-dev">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a></pre>
</blockquote>
<br>
</body>
</html>