Hi,
I have a PCI pix raster, and it's pretty big - 2.4 GB. I'm not familiar with this image format as I've never worked with it before. Gdalinfo seems to recognize it, but importing it with r.in.gdal crashes with the following error:
$
$ gdalinfo DEMvalley.pix
Driver: PCIDSK/PCIDSK Database File
Size is 65362, 41926
Coordinate System is:
PROJCS["UTM Zone 20, Northern Hemisphere",
GEOGCS["Unknown datum based upon the WGS 84 ellipsoid",
DATUM["Not specified (based on WGS 84 spheroid)",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-63],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["Meter",1]]
Origin = (275082.000000,5021966.000000)
Pixel Size = (2.00000000,-2.00000000)
Metadata:
SOFTWARE=PACEV10.
FILE_ID=S_\Peter\Lidar_DEM_MAIN\DEMvalley.pix
GENERATING_FACILITY=PCI Inc., Richmond Hill, Canada
DESCRIPTION1=
DESCRIPTION2=
DATE_OF_CREATION=10_55 15Nov2006
DATE_OF_UPDATE=10_55 15Nov2006
Corner Coordinates:
Upper Left ( 275082.000, 5021966.000) ( 65d52'10.28"W, 45d18'54.52"N)
Lower Left ( 275082.000, 4938114.000) ( 65d49'55.86"W, 44d33'40.47"N)
Upper Right ( 405806.000, 5021966.000) ( 64d12'8.49"W, 45d20'41.55"N)
Lower Right ( 405806.000, 4938114.000) ( 64d11'12.10"W, 44d35'24.74"N)
Center ( 340444.000, 4980040.000) ( 65d 1'22.35"W, 44d57'21.17"N)
Band 1 Block=128x128 Type=Float32, ColorInterp=Undefined
Description = DEMfull.pix_1. focus _Exported from 9 on D_\valley_lidar_DEM\Li
Metadata:
DATE_OF_CREATION=10_55 15Nov2006
DATE_OF_UPDATE=10_55 15Nov2006
UNITS=
LAYER_CRC=0x2a8e
$ r.in.gdal -e input=DEMvalley.pix output=DEMvalley
WARNING: Datum 'Not_specified_based_on_WGS_84_spheroid' not recognised by
GRASS and no parameters found. Datum transformation will not be
possible using this projection information.
Projection of input dataset and current location appear to match.
Proceeding with import...
*** glibc detected *** double free or corruption (!prev): 0x08484740 ***
Aborted
I'm assuming that the reason this fails is because there's too many rows and columns? Alternatively, is there some command line parameter that can be passed to gdal_translate to import it in sections or tiles?
Thanks,
~ Eric.