I'm trying to reproject the Natural Earth I 1:50m raster into an Azimuthal Equidistant projection centered on my hometown. As you may know, the source raster covers from -180° to 180° latitude, and -90° to 90° Here's the command I'm using:
gdalwarp -multi -dstnodata 0 -ot Byte -t_srs "+proj=aeqd +lat_0=48.037 +lon_0=-122.408 +ellps=sphere +units=m +no_defs" -r cubic -of GTiff -ts 10800 10800 NE1_50M_SR_W.tif natural_earth_10800.tif
This works pretty OK, except that it leaves a big, ugly gap west of 180°. Also, it cuts off some of each side of the output raster.
I've experimented with setting a radius instead of the
+ellps parameter, with no change. Using WGS84 gives only a small slice of the Earth as output.
Here's what I got when I tried
+R_A instead of
The result with WGS84 was much the same. Also tried reprojecting to Spherical Mercator as an intermediate projection before reprojecting to my custom AEqD, same artifacts appear.
So, how do I fill in the gap? The pixels are definitely there in the input raster, they're just not making it into the output.