Построение файла csv в матрице высот (geotiff) - PullRequest
Новая
       54

Построение файла csv в матрице высот (geotiff)

0 голосов
/ 28 мая 2020

У меня есть файл csv с именем "Coordinates.csv", который содержит широту и долготу точки. Я хочу построить его на DEM (geotiff), а также у меня есть еще один файл с несколькими CSV, который содержит данные о величине и векторе. Я хочу нанести его на точную точку в DEM (geotiff) после нанесения широты и долготы в Coordinates.csv. Вы можете помочь мне? Я застрял в построении Coordinates.csv. Спасибо!

Вот мой пример кода: 

from osgeo import gdal, gdalconst
import numpy as np
import matplotlib 
import matplotlib.pyplot as plt
import elevation
import affine
import cartopy.crs as ccrs

filename = "Sample-30m-DEM.tif"
gdal_data = gdal.Open(filename)
gdal_band = gdal_data.GetRasterBand(1)
nodatava1 = gdal_band.GetNoDataValue()

# convert to a numpy array
data_array = gdal_data.ReadAsArray().astype(np.float)
data_array


# replace missing values if necessary
if np.any(data_array == nodatava1):
    data_array[data_array == nodatava1] = np.nan


# Plot out data with Matplotlob's 'contour'
fig = plt.figure(figsize = (12, 8))
ax = fig.add_subplot(111)
plt.contour(data_array, cmap=plt.cm.terrain,
            levels = list(range(0, 2500, 100)))


ax.set(xlabel='Longitude (Degrees)', ylabel='Latitude
(Degrees)', title='Sample Plot')
ax.xaxis.grid(color='gray', linestyle='dashed')
ax.yaxis.grid(color='gray', linestyle='dashed')

cbar = plt.colorbar()

plt.gca().set_aspect('equal', adjustable='box')

plt.show()
...