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>
<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>