在MapGIS67“投影变换”子系统中的“用户文件投影转换”功能模块中, 打开坐标文本文件,设置好相关参数, 点击“投影变换”,然后点击“写到文件”, 就把转换好的坐标数据保存成文本文件了!
回答你的问题费老大劲了,opencv用的不熟
我运行输出:
row=14,col=5
14 5 248 242 234 237
gray2 crop: 220 232 219 5
中心点坐标(109,108),宽高(w=218, h=217)
[95, 134] [83, 173]
斜率k= 325
#解题思路:因为我下载的你的,黑色周围还有白色,所以先去除了周围的白色,保证只有黑色部分#对转灰度图,这样方便计算,每个像素点的值就是一个0-255的值,0为黑色,255为白色
#然后获得灰度图的row和col,与正常思维的width和height相反,row对应height,col对应width
#遍历row和col,先获得第一个白点,再获得最后一个白点,然后根据这两个坐标执行tan计算斜率k值
#读取图像
img2 = cv2imread("/img/blackWhitepng")
#转成单通道黑白图
gray = cv2cvtColor(img2, cv2COLOR_BGR2GRAY)
cv2imshow('img2', img2)
cv2imshow('gray', gray)
sp = grayshape
rows = sp[0]
cols = sp[1]
cr = 0
cl = 0
#裁剪左边和上边空白
for row in range(rows):
isBreak = False
for col in range(cols):
if gray[row, col]==0:
cr = row
cl = col
print("\nrow=%d,col=%d"%(row, col))
isBreak = True
break
#print(gray[row, col], end='')
#if(col==cols-1):
#print("\nrow=%d,col=%d-------------------------------"%(row, col))
if isBreak==True:
break
print (cr, cl, rows, cols, rows-cr, cols-cl)
gray2 = gray[cr:rows-cr, cl:cols-cl]
cv2imshow('gray2', gray2)
#裁剪右下角空白
sp = gray2shape
rows = sp[0]
cols = sp[1]
row = rows-1
print ("gray2 crop:", rows, cols, row, col)
while row!=0:
isBreak = False
col = cols - 1
while col!=0:
if gray2[row,col]==0:
isBreak = True
break
col-=1
if isBreak==True:
break
row-=1
gray3 = gray2[0:row+1, 0:col+1]
cv2imshow('gray3', gray3)
#对裁剪后的gray3求中心点坐标
sp = gray3shape
rows = sp[0]
cols = sp[1]
x = cols//2
y = rows//2
print("中心点坐标(%d,%d),宽高(w=%d, h=%d)"%(x,y,cols,rows))
#求斜率 (y2-y1)/(x2-x1)
#遍历出第一个白点和最后一个白点,做计算
x1y1=[]
x2y2=[]
#求第一个白点坐标
for row in range(rows):
isBreak = False
for col in range(cols):
if gray3[row,col]==255:
x1y1append(col)
x1y1append(row)
isBreak = True
break
if isBreak==True:
break
#求最后一个白点坐标
row = rows-1
while row!=0:
isBreak = False
col = cols - 1
while col!=0:
if gray2[row,col]==255:
x2y2append(col)
x2y2append(row)
isBreak = True
break
col-=1
if isBreak==True:
break
row-=1
print(x1y1,x2y2)
#计算斜率tan值
k = abs(x2y2[1]-x1y1[1])/abs(x2y2[0]-x1y1[0])
print("斜率k=",k)
cv2waitKey(0)
cv2destroyAllWindows()
for i in X:
for j in Y:
d=[i,j]
print(d)
运行结果:
后面没有截图,太长了,一共144行。
希望可以帮到你
以上就是关于如何用Python批量获取经纬度坐标全部的内容,包括:如何用Python批量获取经纬度坐标、python怎么在一群点集中,提取中心坐标、python从两个列表中取值,生成坐标等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)