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.tifThis 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 +ellps=sphere:
The result with WGS84 was much the same.
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.
أكثر...
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.tifThis 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 +ellps=sphere:
The result with WGS84 was much the same.
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.
أكثر...