[Gdal-dev] AAIGrid driver: small PROJCS detection problem

[Markus Neteler]

On Fri, Apr 23, 2004 at 09:45:43AM -0400, Frank Warmerdam wrote:
> Markus Neteler wrote:
> 
> > PROJCS["Monte_Mario_Italy_1_BZ",
> >    GEOGCS["GCS_Monte_Mario",
> >        DATUM["Monte_Mario",
> >            SPHEROID["International_1924",6378388,297]],
> >        PRIMEM["Greenwich",0],
> >        UNIT["Degree",0.017453292519943295]],
> >    PROJECTION["Transverse_Mercator"],
> >    PARAMETER["False_Easting",500000.0],
> >    PARAMETER["False_Northing",-5000000.0],
> >    PARAMETER["Central_Meridian",9],
> >    PARAMETER["Scale_Factor",0.9996],
> >    PARAMETER["Latitude_Of_Origin",0],
> >    UNIT["Meter",1]]
> 
> Markus,
> 
> I have committed a bunch of additional changes to GDAL to preserve the
> PROJCS and GEOGCS names as citations in the geotiff files, and to recover
> them when reading and forming WKT again.  I can't seem to find the
> exact message that triggered me into doing that but it was related to
> this thread.


Thanks, I have retried it:

gdalwarp -order 2 -t_srs '+proj=tmerc +lat_0=0 +lon_0=9 \
      +k=0.999600 +x_0=500000 +y_0=-5000000 +ellps=intl +units=m \
      +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68' \
      -tr 28.5 28.5 p192r028_7t20010826_z32_nn10.tif p192r028_7t20010826_z32_nn10_PBZ.tif

# creates:
gdalinfo p192r028_7t20010826_z32_nn10_PBZ.tif
Driver: GTiff/GeoTIFF
Size is 8765, 7797
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["International 1909 (Hayford)",
        DATUM["unknown",
            SPHEROID["unnamed",6378388,297.000000000005]],
        PRIMEM["Greenwich",0],
        UNIT[,0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",9],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",-5000000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (609787.107032,213299.997571)
Pixel Size = (28.50000000,-28.50000000)
Corner Coordinates:
Upper Left  (  609787.107,  213299.998) ( 10d26'44.63"E, 47d 3'47.46"N)
Lower Left  (  609787.107,   -8914.502) ( 10d23'40.03"E, 45d 3'49.97"N)
Upper Right (  859589.607,  213299.998) ( 13d43'39.92"E, 46d58'28.59"N)
Lower Right (  859589.607,   -8914.502) ( 13d33'38.62"E, 44d58'52.44"N)
Center      (  734688.357,  102192.748) ( 12d 1'57.18"E, 46d 1'55.89"N)
Band 1 Block=8765x1 Type=Byte, ColorInterp=Gray

#... which looks ok, but then, when cutting out a piece:

gdal_translate -of GTiff -projwin 610000 205000 767840 118240 \
       p192r028_7t20010826_z32_nn10_PBZ.tif \
       p192r028_7t20010826_z32_nn10_PBZ_small.tif
Input file size is 8765, 7797
Computed -srcwin 7 291 5538 3044 from projected window.
0.Segmentation fault

With gdb I got:
...
Program received signal SIGSEGV, Segmentation fault.
0x42074462 in _int_malloc () from /lib/tls/libc.so.6
(gdb) bt
#0  0x42074462 in _int_malloc () from /lib/tls/libc.so.6
#1  0x4207378d in malloc () from /lib/tls/libc.so.6
#2  0x40160957 in VSIMalloc () from /usr/local/lib/libgdal.so.1
#3  0x40143fe7 in _TIFFmalloc () from /usr/local/lib/libgdal.so.1
#4  0x4013e2d7 in TIFFClientOpen () from /usr/local/lib/libgdal.so.1
#5  0x40143f89 in TIFFOpen () from /usr/local/lib/libgdal.so.1
#6  0x4014f980 in XTIFFOpen () from /usr/local/lib/libgdal.so.1
#7  0x400dc542 in GTiffCreate(char const*, int, int, int, GDALDataType, char**) ()
   from /usr/local/lib/libgdal.so.1
#8  0x400dcc8e in GTiffCreateCopy(char const*, GDALDataset*, int, char**, int (*)(double, char const*, void*), void*) () from /usr/local/lib/libgdal.so.1
#9  0x4015181e in GDALDriver::CreateCopy(char const*, GDALDataset*, int, char**, int (*)(double, char const*, void*), void*) () from /usr/local/lib/libgdal.so.1
#10 0x40151c75 in GDALCreateCopy () from /usr/local/lib/libgdal.so.1
#11 0x0804a142 in main ()
#12 0x42015704 in __libc_start_main () from /lib/tls/libc.so.6

Maybe there is still an initialization problem left.

Kind regards

 Markus Neteler