Reference

pkascii2img(1)

program to create raster image based on ascii file

NAME

pkascii2img - program to create raster image based on ascii file

SYNOPSIS

FCpkascii2imgF[0] [FC-iF[0] input] [FC-oF[0] output] [FC-otF[0] type] [FC-ofF[0] GDALformat] [FC-coF[0] NAME=VALUE] [FC-dxF[0] value] [FC-dyF[0] value] [FC-ulxF[0] value] [FC-ulyF[0] value] [FC-ctF[0] filename] [FC-a_srsF[0] EPSG:number] [FC-dF[0] description]

DESCRIPTION

pkascii2img creates a raster dataset from an ASCII textfile. The textfile is in matrix format (rows and columns). The dimensions in x and y are defined by the number of columns and rows respectively. The georeferencing can be defined by providing the options for cell size (FC-dxF[0] FC-dyF[0]), upper left position (FC-ulxF[0] FC-ulyF[0]) and the projection (FC-a_srsF[0]). Some dataset formats can also store a description (FC-dF[0]) and a color table (FC-ctF[0]).

OPTIONS

FC-iF[0] filename, FC--inputF[0] filename

input ASCII file

FC-oF[0] filename, FC--outputF[0] filename

Output image file

FC-otF[0] type, FC--otypeF[0] type

Data type for output image ({Byte / Int16 / UInt16 / UInt32 / Int32 / Float32 / Float64 / CInt16 / CInt32 / CFloat32 / CFloat64}). Empty string: inherit type from input image

FC-ofF[0] GDALformat, FC--oformatF[0] GDALformat

Output image format (see also gdal_translate(1)). Empty string: inherit from input image

FC-coF[0] NAME=VALUE, FC--coF[0] NAME=VALUE

Creation option for output file. Multiple options can be specified.

FC-ulxF[0] value, FC--ulxF[0] value

Upper left x value bounding box (in geocoordinates if georef is true)

FC-ulyF[0] value, FC--ulyF[0] value

Upper left y value bounding box (in geocoordinates if georef is true)

FC-dxF[0] value, FC--dxF[0] value

Output resolution in x (in meter)

FC-dyF[0] value, FC--dyF[0] value

Output resolution in y (in meter)

FC-ctF[0] filename, FC--ctF[0] filename

colour table in ASCII format having 5 columns: id R G B ALFA (0: transparent, 255: solid)

FC-a_srsF[0] EPSG:number, FC--a_srsF[0] EPSG:number

Override the projection for the output file

FC-dF[0] description, FC--descriptionF[0] description

Set image description

FC-vF[0], FC--verboseF[0]

verbose

EXAMPLE

We create a 5x5 raster dataset containing all 0 values except for its center (value=1)

FC
cat input.txt
0 0 0 0 0
0 0 0 0 0
0 0 1 0 0
0 0 0 0 0
0 0 0 0 0
      F[0]

We define a pixel size of 25 by 25 m and and a projected in the European projection system ETS89/LAEA. The location of the upper left pixel is also set:

FC
pkascii2img -i input.txt -o output.tif -of GTiff -a_srs epsg:3035 -dx 25 -dy 25 -ulx 1000000 -uly 1000000
      F[0]

We check the output GeoTIFF using gdalinfo:

FC
gdalinfo output.tif
Driver: GTiff/GeoTIFF
Files: output.tif
Size is 5, 5
Coordinate System is:
PROJCS["ETRS89 / LAEA Europe",
GEOGCS["ETRS89",
DATUM["European_Terrestrial_Reference_System_1989",
SPHEROID["GRS 1980",6378137,298.2572221010042,
AUTHORITY["EPSG","7019"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6258"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4258"]],
PROJECTION["Lambert_Azimuthal_Equal_Area"],
PARAMETER["latitude_of_center",52],
PARAMETER["longitude_of_center",10],
PARAMETER["false_easting",4321000],
PARAMETER["false_northing",3210000],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","3035"]]
Origin = (1000000.000000000000000,1000000.000000000000000)
Pixel Size = (25.000000000000000,-25.000000000000000)
Metadata:
AREA_OR_POINT=Area
TIFFTAG_DATETIME=2014:11:08 11:33:52
TIFFTAG_DOCUMENTNAME=output.tif
TIFFTAG_SOFTWARE=pktools 2.5.4 by Pieter Kempeneers
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 1000000.000, 1000000.000) ( 23d13' 9.44"W, 25d30'21.01"N)
Lower Left ( 1000000.000, 999875.000) ( 23d13' 8.02"W, 25d30'17.24"N)
Upper Right ( 1000125.000, 1000000.000) ( 23d13' 5.27"W, 25d30'22.72"N)
Lower Right ( 1000125.000, 999875.000) ( 23d13' 3.85"W, 25d30'18.95"N)
Center ( 1000062.500, 999937.500) ( 23d13' 6.65"W, 25d30'19.98"N)
Band 1 Block=5x5 Type=Byte, ColorInterp=Gray
gdalinfo output.tif
Driver: GTiff/GeoTIFF
Files: output.tif
Size is 5, 5
Coordinate System is:
PROJCS["ETRS89 / LAEA Europe",
GEOGCS["ETRS89",
DATUM["European_Terrestrial_Reference_System_1989",
SPHEROID["GRS 1980",6378137,298.2572221010042,
AUTHORITY["EPSG","7019"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6258"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4258"]],
PROJECTION["Lambert_Azimuthal_Equal_Area"],
PARAMETER["latitude_of_center",52],
PARAMETER["longitude_of_center",10],
PARAMETER["false_easting",4321000],
PARAMETER["false_northing",3210000],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","3035"]]
Origin = (1000000.000000000000000,1000000.000000000000000)
Pixel Size = (25.000000000000000,-25.000000000000000)
Metadata:
AREA_OR_POINT=Area
TIFFTAG_DATETIME=2014:11:08 11:33:52
TIFFTAG_DOCUMENTNAME=output.tif
TIFFTAG_SOFTWARE=pktools 2.5.4 by Pieter Kempeneers
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 1000000.000, 1000000.000) ( 23d13' 9.44"W, 25d30'21.01"N)
Lower Left ( 1000000.000, 999875.000) ( 23d13' 8.02"W, 25d30'17.24"N)
Upper Right ( 1000125.000, 1000000.000) ( 23d13' 5.27"W, 25d30'22.72"N)
Lower Right ( 1000125.000, 999875.000) ( 23d13' 3.85"W, 25d30'18.95"N)
Center ( 1000062.500, 999937.500) ( 23d13' 6.65"W, 25d30'19.98"N)
Band 1 Block=5x5 Type=Byte, ColorInterp=Gray
      F[0]

SEE ALSO

pkascii2ogr(1)

History