我试图将坐标从经度/纬度转换为像素,但在此过程中我失去了点。将坐标转换为失去点的像素
,我正在使用的代码如下:
from matplotlib import pyplot as plt
import numpy as np
import pandas as pd
df=pd.read_csv('cords.csv')
cords=df.as_matrix(columns=['x','y'])
gt=[7.6445503225, 5.4065168747250134e-06, 0.0, 45.07436634583334, 0.0, -5.406516856707135e-06]
index=np.zeros(cords.shape)
index[:,1]=((cords[:,1] - gt[3])/gt[5]).round()
index[:,0]=((cords[:,0] - gt[0])/gt[1]).round()
index=index.astype(int)
index[:,0]=index[:,0]-min(index[:,0])+1
index[:,1]=index[:,1]-min(index[:,1])+1
row=max(index[:,1])
col=max(index[:,0])
image=np.zeros([row+1,col+1])
for i in range(0,len(index)):
image[index[i,1],index[i,0]]=255
正如你所看到的,在将纬度/经度转换为像素数字时有一些缺失点。黄色是255值,紫色是0值。 这怎么解决?
在这里,你会发现我使用cords.csv
这里你可以找到与需要被设置到每个像素的值的坐标的坐标。 cords_valus.csv
什么是'gt'在做什么? –
它是从经纬度到轴坐标的转换向量。使用GDAL获得。 –