I am new to GIS processing, and I am using Python 3.6 to do so.
My current task is to read in several files of geographic data (specifically tiles of elevation data from JAXA), combine them for my area of interest, and then analyze them using Python.
Unfortunately, the data for each tile from AW3D30 come as a folder of six files: three
tif files and two
txt files. I haven't found any documentation about what each of these files contains, so I don't know which one(s) I need for the elevation data (although it must be at least the DSM file). And I don't know how to properly read them in and process them together, connect the proper lon/lat data to the pixels, etc.
For example, I read the
*_AVE_DSM.tif file into Python using gdal, found out there is one band, and the band has a min and max value of
altitudeData1 = gdal.Open('N034E138/N034E138_AVE_DSM.tif') print("[ RASTER BAND COUNT ]: ", altitudeData1.RasterCount) #-> 1 print("[ MIN ] = ", altitudeData1.GetRasterBand(1).GetMinimum()) #-> None print("[ MAX ] = ", altitudeData1.GetRasterBand(1).GetMaximum()) #-> None
So I'm clearly missing something important.
My final desired output is geoPandas dataframe with a geometry column capturing the polygon of each pixel (i.e., not the center point, but the square's geometry) and a column for elevation. I am open to various workflows and intermediary steps.