pkascii2img(1)

NAME
pkascii2img - program to create raster image based on ascii file

SYNOPSIS
\FCpkascii2img\F[]
[\FC-i\F[] input] [\FC-o\F[] output] [\FC-ot\F[] type] [\FC-of\F[] GDALformat] [\FC-co\F[] NAME=VALUE] [\FC-dx\F[] value] [\FC-dy\F[] value] [\FC-ulx\F[] value] [\FC-uly\F[] value] [\FC-ct\F[] filename] [\FC-a_srs\F[] EPSG:number] [\FC-d\F[] 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-dx\F[] \FC-dy\F[]), upper left position (\FC-ulx\F[] \FC-uly\F[]) and the projection (\FC-a_srs\F[]). Some dataset formats can also store a description (\FC-d\F[]) and a color table (\FC-ct\F[]).

OPTIONS

\FC-i\F[] filename, \FC--input\F[] filename
    input ASCII file

\FC-o\F[] filename, \FC--output\F[] filename
    Output image file

\FC-ot\F[] type, \FC--otype\F[] 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-of\F[] GDALformat, \FC--oformat\F[] GDALformat
    Output image format (see also gdal_translate(1)). Empty string: inherit from input image

\FC-co\F[] NAME=VALUE, \FC--co\F[] NAME=VALUE
    Creation option for output file. Multiple options can be specified.

\FC-ulx\F[] value, \FC--ulx\F[] value
    Upper left x value bounding box (in geocoordinates if georef is true)

\FC-uly\F[] value, \FC--uly\F[] value
    Upper left y value bounding box (in geocoordinates if georef is true)

\FC-dx\F[] value, \FC--dx\F[] value
    Output resolution in x (in meter)

\FC-dy\F[] value, \FC--dy\F[] value
    Output resolution in y (in meter)

\FC-ct\F[] filename, \FC--ct\F[] filename
    colour table in ASCII format having 5 columns: id R G B ALFA (0: transparent, 255: solid)

\FC-a_srs\F[] EPSG:number, \FC--a_srs\F[] EPSG:number
    Override the projection for the output file

\FC-d\F[] description, \FC--description\F[] description
    Set image description

\FC-v\F[], \FC--verbose\F[]
    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[]

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[]

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[]

SEE ALSO
pkascii2ogr(1)