Visualize text file with location info and intensity

Question

I have an ascii text file containing location data (column 9-lat and 10-long) and intensity(column 20)

  200010  207  020311    40658.5  406593 52 344927.31 7100203.50  -26.2078720  127.4491855 345060.64 7100369.14   26.4  650.3  628.0 55471.293    20.168 55648.817 55637.523  -146.062

the text file has many lines 10k+

I am trying to export this ascii text into a GeoTIFF using GDAL, but not sure how to proceed.

Ideas?

Answer 1

With GDAL, you would need to create a VRT file. This needs to be done by hand for what you want, but it's not difficult and is explained here. But for completeness, this is what you need to do:

Using a text editor, copy this snippet into it and save it as whatever.vrt The fields should be self explanitory, but note that without column headings in the CSV file (that is the source ASCII file you have) the fields you use must be "field_n" where n=1 is the first column.

<OGRVRTDataSource>
    <OGRVRTLayer name="test">
        <SrcDataSource>test.csv</SrcDataSource>
        <GeometryType>wkbPoint</GeometryType>
        <LayerSRS>WGS84</LayerSRS>
        <GeometryField encoding="PointFromColumns" x="field_10" y="field_9" z="field_20"/>
    </OGRVRTLayer>
</OGRVRTDataSource>

The next step is to run gdal_grid which will interpolate between points if you ask it to. The minimum you need on the command line is:

gdal_grid -ot Float32 -l test test.vrt test.tif

Where -l is the name of the layer as specified in the OGRVRTLayer tag in the VRT. There are many options you can add to gdal_grid to get the results you want, so look here for them all, as well as some extra stuff about VRT files.

If you're not confident with the command line, the latest QGIS provides the gdal_grid function as an option in the Raster menu, but the building of the VRT part has to be done by hand because the QGIS one only deals with rasters.

Answer 2

GDAL doesn't have it's own viewer as far as I know. But there are many GIS using GDAL, for example QGIS.

How about using QGIS Delimited Text Layer plugin? It will load every line as a point with all columns as attributes. You will want to add a header to the text file before you import it.