如何用Python批量获取经纬度坐标

如何用Python批量获取经纬度坐标,第1张

在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从两个列表中取值,生成坐标等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!

欢迎分享,转载请注明来源:内存溢出

原文地址: http://outofmemory.cn/web/9324203.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2023-04-27
下一篇 2023-04-27

发表评论

登录后才能评论

评论列表(0条)

保存