Estou tentando traçar um mapa raster usando Pyhton. A imagem não está sendo salva, embora apareça. Além disso, os pixels não estão aparecendo no local exato. Estou obtendo a imagem conforme mostrado abaixo (o que está errado, no entanto, a rotulação lat-lon está correta.)
Abaixo está a orientação correta, que eu exportei do ArcGIS.
Há algum erro em colocar o pixel no local certo, mas não estou entendendo. Como corrigi-lo?
Aqui está o meu roteiro.
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
from osgeo import gdal
from numpy import linspace
from numpy import meshgrid
import numpy as
#################################################################################
ds = gdal.Open('Path\\To\\Raster.tif')
data = ds.ReadAsArray()
gt = ds.GetGeoTransform()
proj = ds.GetProjection()
#################################################################################
xres = gt[1]
yres = gt[5]
xmin = gt[0] + xres * 0.5
xmax = gt[0] + (xres * ds.RasterXSize) - xres * 0.5
ymin = gt[3] + (yres * ds.RasterYSize) + yres * 0.5
ymax = gt[3] - yres * 0.5
x_center=(xmin+xmax)/2
y_center=(ymin+ymax)/2
#################################################################################
fig = plt.figure(figsize=(18, 10))
ax = fig.add_subplot(111, axisbg='w', frame_on=True)
m = Basemap(llcrnrlon=xmin,llcrnrlat=ymin,urcrnrlon=xmax,urcrnrlat=ymax,
projection='tmerc', lat_0 = y_center, lon_0 = x_center)
parallels = np.arange(15.,30.,0.25)
m.drawparallels(parallels,labels=[1,0,0,0],fontsize=12, linewidth=0.4)
meridians = np.arange(70.,90.,0.25)
m.drawmeridians(meridians,labels=[0,0,0,1],fontsize=12, linewidth=0.4)
x = linspace(0, m.urcrnrx, data.shape[1])
y = linspace(0, m.urcrnry, data.shape[0])
xx, yy = meshgrid(x, y)
m.pcolormesh(xx, yy, data, cmap=plt.cm.jet)
plt.show()
plt.savefig('Path\\To\\Save_Image.png', bbox_inches='tight', pad_inches=.2, dpi=600)
imshow
vez depcolormesh
? É preciso um canto diferente como a origem da imagem por padrão.pcolormesh
porimshow
ele diz NameError: name 'imshow' não está definido Como usar o `imshow 'aqui? @countrymanRespostas:
Eu resolvi esse problema.
Eu estava lendo o arquivo de varredura em linhas de cima para baixo e plotando-o em linhas de baixo para cima.
Desde então, não posso fazer nada com plotagem, virei o array e funcionou.
Abaixo está a modificação que eu fiz,
fonte