[gdal-dev] Conversion from ASCII to GeoTIFF with gdalwarp: incorrect output

Jamie Adams jaadfoo at gmail.com
Tue Oct 14 15:52:27 PDT 2014


I would use the CENTER_LONG config argument.  That will extend the raster
past the -180/180 boundary depending on which side you want to extend.

--config CENTER_LONG 180
Upper Left  ( 176.0010322, -16.1523859) (176d 0' 3.72"E, 16d 9' 8.59"S)
Lower Left  ( 176.0010322, -18.2685614) (176d 0' 3.72"E, 18d16' 6.82"S)
Upper Right (     181.182,     -16.152) (181d10'55.25"E, 16d 9' 8.59"S)
Lower Right (     181.182,     -18.269) (181d10'55.25"E, 18d16' 6.82"S)

--config CENTER_LONG -180
Upper Left  (    -183.999,     -16.152) (183d59'56.28"W, 16d 9' 8.59"S)
Lower Left  (    -183.999,     -18.269) (183d59'56.28"W, 18d16' 6.82"S)
Upper Right (-178.8179863, -16.1523859) (178d49' 4.75"W, 16d 9' 8.59"S)
Lower Right (-178.8179863, -18.2685614) (178d49' 4.75"W, 18d16' 6.82"S)


On Sun, Oct 12, 2014 at 1:43 PM, Lukasz Tracewski <
lukasz.tracewski at outlook.com> wrote:

> Hi,
>
> I have a few thousands of ESRI ASCII rasters in World_Mollweide
> (EPSG:54009) projection, which I need to convert to GeoTIFF. Usually it
> goes fine, but there are certain special cases that resist peaceful
> conversion. They are all placed on Fiji, which is special in a sense being
> on very border of a world's map. Problem is that in those cases somehow
> output is completely messed up, with incorrect resolution.
>
> Here is how I do the conversion:
> gdalwarp -s_srs '+proj=moll +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84
> +datum=WGS84 +units=m no_defs' -t_srs EPSG:4326 -co NBITS=1 -co
> COMPRESS=CCITTRLE -co PHOTOMETRIC=MINISWHITE -ot Byte $ascfile $newname
>
> My ASCII files are bilevel, with 0 and 1 only.
>
> Here is an example of a problematic ASCII:
>
> https://drive.google.com/file/d/0B46Cy3WOKgyJYVdMLWV0OGFJQ2M/view?usp=sharing
>
> That will produce GeoTIFF with only 3 pixels in latitude:
>
> https://drive.google.com/file/d/0B46Cy3WOKgyJTnoxSHdyaUR5cWs/view?usp=sharing
>
> Thing is if I crop the ASCII in QGIS to include only left side of the
> image, then it's fine. If I take the right side, then again I get incorrect
> output.
>
> Any ideas?
>
> Cheers,
> Lucas
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20141014/3eaabd6a/attachment.html>


More information about the gdal-dev mailing list