python經緯度

『壹』 怎麼用python進行經緯度間距離的計算

以前有過類似的問題 可以參考下：

http://blog.csdn.net/avgworkld/article/details/77989578

『貳』 python如何根據經緯度輸出最近距離點

但是這樣是一連串的數字不是我們想要的結果，我們可以利用time模塊的格式化時間的方法來處理：
time.localtime(time.time())
用time.localtime()方法，作用是格式化時間戳為本地的時間。
輸出的結果是：
time.struct_time(tm_year=2010, tm_mon=7, tm_mday=19, tm_hour=22, tm_min=33, tm_sec=39, tm_wday=0, tm_yday=200, tm_isdst=0)

『叄』 python將高斯坐標轉換經緯度 經緯度坐標與高斯坐標的轉換代碼

#網上搜來的

# 高斯坐標轉經緯度演算法 # B=大地坐標X # C=大地坐標Y # IsSix=6度帶或3度帶

import math
def GetLatLon2(B, C,IsSix):
#帶號
D = math.trunc( C/ 1000000)

#中央經線(單位：弧度)
K = 0
if IsSix:
K = D * 6 - 3 #6度帶計算
else:
K = D * 3 #3度帶計算
L = B/(6378245*(1-0.006693421623)*1.0050517739)
M = L +(0.00506237764 * math.sin(2*L)/2-0.00001062451*math.sin(4*L)/4+0.0000002081*math.sin(6*L)/6)/1.0050517739
N = L +(0.00506237764 * math.sin(2*M)/2-0.00001062451*math.sin(4*M)/4+0.0000002081*math.sin(6*M)/6)/1.0050517739
O = L +(0.00506237764 * math.sin(2*N)/2-0.00001062451*math.sin(4*N)/4+0.0000002081*math.sin(6*N)/6)/1.0050517739
P = L +(0.00506237764 * math.sin(2*O)/2-0.00001062451*math.sin(4*O)/4+0.0000002081*math.sin(6*O)/6)/1.0050517739
Q = L +(0.00506237764 * math.sin(2*P)/2-0.00001062451*math.sin(4*P)/4+0.0000002081*math.sin(6*P)/6)/1.0050517739
R = L +(0.00506237764 * math.sin(2*Q)/2-0.00001062451*math.sin(4*Q)/4+0.0000002081*math.sin(6*Q)/6)/1.0050517739
S = math.tan(R)
T = 0.006738525415*(math.cos(R))**2
U = 6378245/math.sqrt(1-0.006693421623*(math.sin(R))**2)
V = 6378245*(1-0.006693421623)/(math.sqrt((1-0.006693421623*(math.sin(R))**2)))**3
W = 5+3*S**2+T-9*T*S**2
X = 61+90*S**2+45*S**4
Y = 1+2*S**2+T**2
Z = 5+28*S**2+24*S**4+6*T+8*T*S**2
Lat= (180/math.pi)*(R-(C-D*1000000-500000)**2*S/(2*V*U)+(C-D*1000000-500000)**4*W/(24*U**3*V)-(C-D*1000000-500000)**6*X/(7200*U**5*V))
Lon= (180/math.pi)*(C-D*1000000-500000)*(1-(C-D*1000000-500000)**2*Y/(6*U**2)+(C-D*1000000-500000)**4*Z/(120*U**4))/(U*math.cos(P))
Lat = Lat
Lon = K + Lon
return (Lon, Lat)

『肆』 利用python編程，經緯度的十進制與度分秒的轉換

file=input("請輸入文件位置：")
f=open(file)
k=open("轉換後.txt",'a')
foriinf:
line=f.readline()
_1=line.find("°")
=line[:_1]
fen_1=line.find("′")
fen=line[_1+1:fen_1]
miao_1=line.find('''″''')
miao=line[fen_1+1:miao_1]
k.write(line.strip()+"----------"+str(int()+int(fen)/60+int(miao)/3600)[:12])
k.write('
')
f.close()
k.close()

『伍』 使用Python實時將gps返回的經緯度轉化為圖片

主要就是做了兩件事情：
1.生成一張有文本信息的JPG圖片
2.寫入EXIF信息

生成照片需要PIL和libjpeg

import Imageimport ImageDrawimport ImageFontdef create_pic(path, text=[], type='jpeg'):
img = Image.new("RGB", (2448, 3264), '#37b6ce')#顏色和大小
draw = ImageDraw.Draw(img)
font = ImageFont.truetype('ziti.ttf', 120)#ttf是字體，120是字型大小
for h in range(0, len(text)):#多行文本
draw.text((256, 256 + 120 * h), text[h], font=font)
img.save(path, type)#保存
# img.show()

讀寫EXIF信息需要pyexiv2，獲取google的經緯度需要geopy
順便說下經緯度的表示：
一般exif里看到的都是這樣的57°55'56.6"，是度，分，秒這么展示的，google上獲取來的是十進制的57.9323888888888
所以需要轉換一下再寫進去
公式：57°55'56.6" =57+55/60+56.6/3600=57.9323888888888

etemplate=pyexiv2.ImageMetadata('IMG_4408.JPG')template.read()#exif信息很多，所以找個真正手機拍攝的照片當模版googlev3=geopy.GoogleV3()place,gps=googlev3.geocode(location)#獲取gps信息，location寫地名，比如『北京王府井』，偶爾會被牆，最好掛個代理defset_exif(path,date_time=None,gps=()):
"""
datetime=2014:10:0412:41:38
geo=(lat=39.12315,lng=115.12231)
"""
metadata=pyexiv2.ImageMetadata(path)
metadata.read()
forkintemplate.exif_keys:
metadata[k]=pyexiv2.ExifTag(k,template[k].value)
ifnotdate_time:
date_str=pyexiv2.utils.exif(date_time)
metadata['Exif.Photo.DateTimeOriginal']=date_str
metadata['Exif.Photo.DateTimeDigitized']=date_str
metadata['Exif.Image.DateTime']=date_str
iflen(geo)>0:
c_lat=decimal2coordinate(geo[0],['S','N'])
c_lng=decimal2coordinate(geo[1],['W','E'])
metadata["Exif.GPSInfo.GPSLatitude"]=coordinate2rational(c_lat[0],c_lat[1],c_lat[2])
metadata["Exif.GPSInfo.GPSLatitudeRef"]=c_lat[3]
metadata["Exif.GPSInfo.GPSLongitude"]=coordinate2rational(c_lng[0],c_lng[1],c_lng[2])
metadata["Exif.GPSInfo.GPSLongitudeRef"]=c_lng[3]
else:
metadata._delete_exif_tag("Exif.GPSInfo.GPSLatitude")
metadata._delete_exif_tag("Exif.GPSInfo.GPSLatitudeRef")
metadata._delete_exif_tag("Exif.GPSInfo.GPSLongitude")
metadata._delete_exif_tag("Exif.GPSInfo.GPSLongitudeRef")
metadata.write()defdecimal2coordinate(value,loc):
"""
loc=lat=>["S","N"],lng=>["W","E"]
retrunD,M,S,locate
"""
ifvalue<0:
loc_value=loc[0]
elifvalue>0:
loc_value=loc[1]
else:
loc_value=""
abs_value=abs(value)
deg=int(abs_value)
t1=(abs_value-deg)*60
min=int(t1)
sec=round((t1-min)*60,5)
return(deg,min,sec,loc_value)defcoordinate2rational(D,M,S):
return(fractions.Fraction(D,1),fractions.Fraction(int((M+S/60)*100),100),fractions.Fraction(0,1))

『陸』 如何用Python批量獲取經緯度坐標

使用Python程序，直接嵌入即可。
需要注意的是網路API獲取的是墨卡托坐標，而實際使用的是WGS84坐標。代碼已添加轉換，只要設定語言一致即可實現。

『柒』 python 十進制經緯度計算10米為多少度

text = '這是一個經度：120°30′30″，這是一個緯度：
30°12′34″'
text.replace(/(\-?\d{1,3})°(\d{1,3})′(\d{1,3})″/g, function (match, de, mi, se) {
return Number(de) + Number(mi)/60 + Number(se)/3600
})

『捌』 Python 如何判斷 經緯度信息

你要說明白點。比方一個地方。和一個模塊。和自己的一個自定

『玖』 python如何講兩列經緯度變成一列的坐標

為了迎合大量數據處理的需求，以java為基礎的大數據開始了。Hadoop為一批數據處理，發展以java為基礎的架構關鍵，相對於其他處理工具，Hadoop慢許多，但是無比的准確可被後端資料庫分析廣泛使用，和Hive搭配的很好。