Python melaksanakan membaca maklumat ketinggian dalam DEM menggunakan koordinat titik latitud dan longitud (contoh terperinci)

Artikel ini membawakan anda isu berkaitan tentang cara menggunakan 10 baris kod untuk membaca maklumat ketinggian dalam DEM menggunakan koordinat titik latitud dan longitud Idea utama ialah menggunakan GDAL untuk membaca data DEM, mendapatkan jalur pertama dan membaca Ambilnya sebagai tatasusunan, kemudian hitung nombor baris dan lajur menggunakan longitud dan latitud, dan gunakan nombor baris dan lajur untuk membaca maklumat ketinggian saya harap ia akan membantu semua orang.

Pembelajaran yang disyorkan: tutorial video mysql

Syarat pelaksanaan: 1. Perpustakaan yang diperlukan ialah GDAL; 2.DEMdata;

1 Idea

(1) Gunakan GDAL untuk membaca data DEM, mula-mula dapatkan jalur pertama dan baca sebagai tatasusunan, Fungsi-fungsi yang digunakan ialah:

ds.GetRasterBand(1)和band.ReadAsArray()；

(2) Kira nombor baris dan lajur mengikut longitud dan latitud

(3) Maklumat ketinggian boleh dibaca terus mengikut baris dan lajur; nombor.

2. Kod

Kod dan ulasan adalah seperti berikut:

from osgeo import gdal

gdal.UseExceptions()
#以湖北DEM数据为例
ds = gdal.Open('./data/hubei_wgs84.tif')

band = ds.GetRasterBand(1)
elevation = band.ReadAsArray()
nrows, ncols = elevation.shape

x0, dx, dxdy, y0, dydx, dy = ds.GetGeoTransform()
#指定经纬度点坐标
latitude,longitude=31.15,111.24
#根据经纬度计算行列号，dx=dy为分辨率，不相等的时候(y0-latitude)/dx改为(y0-latitude)/-dy
new_ncols,new_nrows=int((y0-latitude)/dx),int((longitude-x0)/dx)

#根据行列号读取并打印输出指定坐标点高程
print(elevation[new_ncols][new_nrows])

3 Keputusan

932, iaitu ketinggian 932m.

4. Pengesahan Arcgis

Gunakan alat Pergi Ke XY, alat paling kanan bar alat

Masukkan koordinat latitud dan longitud untuk mencari titik:

Tambah titik pada peta

Gunakan Identiti untuk melihat atribut titik:

Hasilnya adalah konsisten.

Pembelajaran yang disyorkan: tutorial video mysql

Atas ialah kandungan terperinci Python melaksanakan membaca maklumat ketinggian dalam DEM menggunakan koordinat titik latitud dan longitud (contoh terperinci). Untuk maklumat lanjut, sila ikut artikel berkaitan lain di laman web China PHP!