--- name: "pkascii2img" section: "1" platform: "linux" package: "pktools" version: "2.6.7.6+ds-6+b3" release: "Debian 13 (Trixie)" platforms: ["linux"] url: "https://new.carta.tech/pkascii2img" --- # pkascii2img(1) program to create raster image based on ascii file ## NAME pkascii2img - program to create raster image based on ascii file ## SYNOPSIS FC**pkascii2img**F\[0\] \[FC**-i**F\[0\] *input*\] \[FC**-o**F\[0\] *output*\] \[FC**-ot**F\[0\] *type*\] \[FC**-of**F\[0\] *GDALformat*\] \[FC**-co**F\[0\] *NAME=VALUE*\] \[FC**-dx**F\[0\] *value*\] \[FC**-dy**F\[0\] *value*\] \[FC**-ulx**F\[0\] *value*\] \[FC**-uly**F\[0\] *value*\] \[FC**-ct**F\[0\] *filename*\] \[FC**-a\_srs**F\[0\] *EPSG:number*\] \[FC**-d**F\[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**-dx**F\[0\] FC**-dy**F\[0\]), upper left position (FC**-ulx**F\[0\] FC**-uly**F\[0\]) and the projection (FC**-a\_srs**F\[0\]). Some dataset formats can also store a description (FC**-d**F\[0\]) and a color table (FC**-ct**F\[0\]). ## OPTIONS - FC**-i**F\[0\] *filename*, FC**--input**F\[0\] *filename* input ASCII file - FC**-o**F\[0\] *filename*, FC**--output**F\[0\] *filename* Output image file - FC**-ot**F\[0\] *type*, FC**--otype**F\[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**-of**F\[0\] *GDALformat*, FC**--oformat**F\[0\] *GDALformat* Output image format (see also [**gdal\_translate**(1)](/gdal_translate/1.md)). Empty string: inherit from input image - FC**-co**F\[0\] *NAME=VALUE*, FC**--co**F\[0\] *NAME=VALUE* Creation option for output file. Multiple options can be specified. - FC**-ulx**F\[0\] *value*, FC**--ulx**F\[0\] *value* Upper left x value bounding box (in geocoordinates if georef is true) - FC**-uly**F\[0\] *value*, FC**--uly**F\[0\] *value* Upper left y value bounding box (in geocoordinates if georef is true) - FC**-dx**F\[0\] *value*, FC**--dx**F\[0\] *value* Output resolution in x (in meter) - FC**-dy**F\[0\] *value*, FC**--dy**F\[0\] *value* Output resolution in y (in meter) - FC**-ct**F\[0\] *filename*, FC**--ct**F\[0\] *filename* colour table in ASCII format having 5 columns: id R G B ALFA (0: transparent, 255: solid) - FC**-a\_srs**F\[0\] *EPSG:number*, FC**--a\_srs**F\[0\] *EPSG:number* Override the projection for the output file - FC**-d**F\[0\] *description*, FC**--description**F\[0\] *description* Set image description - FC**-v**F\[0\], FC**--verbose**F\[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)](/pkascii2ogr/1.md)