【python】数据挖掘 实验：中国二级城市经纬度聚类分析

题目：数据仓库与挖掘实验 —— 聚类分析

1. 有如下中国城市数据，要求按照经纬度对中国二级城市聚类：
2. 首先显示城市的分布图
3. 观察设置k的合理值
4. 对数据进行聚类，并显示图形聚类结果
5. 输出每一类对应的城市
6. 完整数据在本篇博客最后

7.     北京市 北京市 北纬39.55 东经116.24

       福建省 福州 北纬26.05 东经119.18

       福建省 长乐 北纬25.58 东经119.31

       福建省 福安 北纬27.06 东经119.39

       福建省 福清 北纬25.42 东经119.23

       福建省 建瓯 北纬27.03 东经118.20

       福建省 建阳 北纬27.21 东经118.07

 

代码分析：

1.显示城市分布图

import pandas as pd
import matplotlib.pyplot as plt

inputfile = './data/secondary_province_data.xlsx' 
outputfile = './tmp2/data_type_latitude_longitude.xls'

latitude = pd.read_excel(inputfile, usecols=[3])    #北纬
longtitude = pd.read_excel(inputfile, usecols=[4])  #东经

plt.figure(figsize=(10, 10), dpi=100)
plt.scatter(longtitude, latitude, marker='*', c='red')
plt.rcParams['font.sans-serif']=['SimHei'] #用来正常显示中文标签
plt.rcParams['axes.unicode_minus']=False #用来正常显示负号
plt.xlabel("东经", fontdict={'size': 16})
plt.ylabel("北纬", fontdict={'size': 16})
plt.title("中国城市分布图", fontdict={'size': 20})
plt.show()

 

2. 聚类分析

import pandas as pd

inputfile = './data/secondary_province_data.xlsx' 
outputfile = './tmp2/data_type_latitude_longitude.xls'

k = 3 #聚类类别
iteration = 500 #聚类最大循环次数

data = pd.read_excel(inputfile, usecols=[0, 3, 4], index_col=[0])
data2 = pd.read_excel(inputfile)
# print(data2['secondary'])
# print(data2['id'])
# print(data2.columns)

data_zs = 1.0*(data - data.mean())/data.std()#数据标准化
# print(data_zs)

from sklearn.cluster import KMeans
model = KMeans(n_clusters = k, n_jobs = 1, max_iter = iteration) #分为k类
model.fit(data_zs) #开始聚类

#简单打印结果
r1 = pd.Series(model.labels_).value_counts() #统计各个类别的数目
r2 = pd.DataFrame(model.cluster_centers_) #找出聚类中心
# print(r1)
r = pd.concat([r2, r1], axis = 1) #横向连接（0是纵向），得到聚类中心对应的类别下的数目
# print(data.columns)
r.columns = list(data.columns) + [u'类别数目'] #重命名表头
print(r)

#详细输出原始数据及其类别
r = pd.concat([data, pd.Series(model.labels_, index = data.index)], axis = 1)  #详细输出每个样本对应的类别
r.columns = list(data.columns) + [u'聚类类别']#重命名表头

r.to_excel(outputfile) #保存结果
#  + [data2[u'secondary']]  
print(r)


def density_plot(data): #自定义作图函数
  import matplotlib.pyplot as plt
  plt.rcParams['font.sans-serif'] = ['SimHei'] #用来正常显示中文标签
  plt.rcParams['axes.unicode_minus'] = False #用来正常显示负号
  p = data.plot(kind='kde', linewidth = 2, subplots = True, sharex = False)
#   [p[i].set_ylabel(u'密度') for i in range(k)]
  plt.legend()
  return plt

pic_output = './tmp2/pd_' #概率密度图文件名前缀
for i in range(3): 
 density_plot(data[r[u'聚类类别']==i]).savefig(u'%s%s.png' %(pic_output, i))

 

3. 聚类结果图

inputfile = './tmp2/data_type_latitude_longitude.xls' #销量及其他属性数据
df = pd.read_excel(inputfile, usecols=[0, 1, 2, 3])
# print(df['id'])
# print(df.columns)

x0 = df.loc[df['聚类类别']==0,['longitude']]
y0 = df.loc[df['聚类类别']==0,['latitude']]
x1 = df.loc[df['聚类类别']==1,['longitude']]
y1 = df.loc[df['聚类类别']==1,['latitude']]
x2 = df.loc[df['聚类类别']==2,['longitude']]
y2 = df.loc[df['聚类类别']==2,['latitude']]

plt.figure(figsize=(10, 10), dpi=100)
plt.scatter(x0, y0, c = "red", marker='o', label='label0')
plt.scatter(x1, y1, c = "blue", marker='*', label='label1')
plt.scatter(x2, y2, c = "green", marker='+', label='label2')
plt.xlabel("东经", fontdict={'size': 16})
plt.ylabel("北纬", fontdict={'size': 16})
plt.legend(loc=2)
plt.title("聚类结果图", fontdict={'size': 20})
plt.show()

 

4. 打印标签对应的城市

res = pd.read_excel('./tmp2/data_type_latitude_longitude.xls')
# print(res)
# print(type(res))
res.insert(1, 'secondary', data2['secondary'])
# print(res)

for i in range(k):
    print(res.loc[res['聚类类别']==i])

 

完整数据 ：

    北京市 北京市 北纬39.55 东经116.24
    福建省 福州 北纬26.05 东经119.18
    福建省 长乐 北纬25.58 东经119.31
    福建省 福安 北纬27.06 东经119.39
    福建省 福清 北纬25.42 东经119.23
    福建省 建瓯 北纬27.03 东经118.20
    福建省 建阳 北纬27.21 东经118.07
    福建省 晋江 北纬24.49 东经118.35
    福建省 龙海 北纬24.26 东经117.48
    福建省 龙岩 北纬25.06 东经117.01
    福建省 南安 北纬24.57 东经118.23
    福建省 南平 北纬26.38 东经118.10
    福建省 宁德 北纬26.39 东经119.31
    福建省 莆田 北纬24.26 东经119.01
    福建省 泉州 北纬24.56 东经118.36
    福建省 三明 北纬26.13 东经117.36
    福建省 邵武 北纬27.20 东经117.29
    福建省 石狮 北纬24.44 东经118.38
    福建省 武夷山 北纬27.46 东经118.02
    福建省 厦门 北纬24.27 东经118.06
    福建省 永安 北纬25.58 东经117.23
    福建省 漳平 北纬25.17 东经117.24
    福建省 漳州 北纬24.31 东经117.39
    甘肃省 兰州 北纬36.04 东经103.51
    甘肃省 白银 北纬36.33 东经104.12
    甘肃省 敦煌 北纬40.08 东经94.41
    甘肃省 嘉峪关 北纬39.48 东经98.14
    甘肃省 金昌 北纬38.28 东经102.10
    甘肃省 酒泉 北纬39.44 东经98.31
    甘肃省 临夏 北纬35.37 东经103.12
    甘肃省 平凉 北纬35.32 东经106.40
    甘肃省 天水 北纬34.37 东经105.42
    甘肃省 武威 北纬37.56 东经102.39
    甘肃省 西峰 北纬35.45 东经107.40
    甘肃省 玉门 北纬39.49 东经97.35
    甘肃省 张掖 北纬38.56 东经100.26
    广东省 广州 北纬23.08 东经113.14
    广东省 潮阳 北纬23.16 东经116.36
    广东省 潮州 北纬23.40 东经116.38
    广东省 澄海 北纬23.28 东经116.46
    广东省 从化 北纬23.33 东经113.33
    广东省 东莞 北纬23.02 东经113.45
    广东省 恩平 北纬22.12 东经112.19
    广东省 佛山 北纬23.02 东经113.06
    广东省 高明 北纬22.53 东经112.50
    广东省 高要 北纬23.02 东经112.26
    广东省 高州 北纬21.54 东经110.50
    广东省 鹤山 北纬22.46 东经112.57
    广东省 河源 北纬23.43 东经114.41
    广东省 花都 北纬23.23 东经113.12
    广东省 化州 北纬21.39 东经110.37
    广东省 惠阳 北纬22.48 东经114.28
    广东省 惠州 北纬23.05 东经114.22
    广东省 江门 北纬22.35 东经113.04
    广东省 揭阳 北纬22.32 东经116.21
    广东省 开平 北纬22.22 东经112.40
    广东省 乐昌 北纬25.09 东经113.21
    广东省 雷州 北纬20.54 东经110.04
    广东省 廉江 北纬21.37 东经110.17
    广东省 连州 北纬24.48 东经112.23
    广东省 罗定 北纬22.46 东经111.33
    广东省 茂名 北纬21.40 东经110.53
    广东省 梅州 北纬24.19 东经116.07
    广东省 南海 北纬23.01 东经113.09
    广东省 番禺 北纬22.57 东经113.22
    广东省 普宁 北纬23.18 东经116.10
    广东省 清远 北纬23.42 东经113.01
    广东省 三水 北纬23.10 东经112.52
    广东省 汕头 北纬23.22 东经116.41
    广东省 汕尾 北纬22.47 东经115.21
    广东省 韶关 北纬24.48 东经113.37
    广东省 深圳 北纬22.33 东经114.07
    广东省 顺德 北纬22.50 东经113.15
    广东省 四会 北纬23.21 东经112.41
    广东省 台山 北纬22.15 东经112.48
    广东省 吴川 北纬21.26 东经110.47
    广东省 新会 北纬22.32 东经113.01
    广东省 兴宁 北纬24.09 东经115.43
    广东省 阳春 北纬22.10 东经111.48
    广东省 阳江 北纬21.50 东经111.58
    广东省 英德 北纬24.10 东经113.22
    广东省 云浮 北纬22.57 东经112.02
    广东省 增城 北纬23.18 东经113.49
    广东省 湛江 北纬21.11 东经110.24
    广东省 肇庆 北纬23.03 东经112.27
    广东省 中山 北纬22.31 东经113.22
    广东省 珠海 北纬22.17 东经113.34
    海南省 海口 北纬20.02 东经110.20
    海南省 儋州 北纬19.31 东经109.34
    海南省 琼海 北纬19.14 东经110.28
    海南省 琼山 北纬19.59 东经110.21
    海南省 三亚 北纬18.14 东经109.31
    海南省 通什 北纬18.46 东经109.31
    河北省 石家庄 北纬38.02 东经114.30
    河北省 安国 北纬38.24 东经115.20
    河北省 保定 北纬38.51 东经115.30
    河北省 霸州 北纬39.06 东经116.24
    河北省 泊头 北纬38.04 东经116.34
    河北省 沧州 北纬38.18 东经116.52
    河北省 承德 北纬40.59 东经117.57
    河北省 定州 北纬38.30 东经115.00
    河北省 丰南 北纬39.34 东经118.06
    河北省 高碑店 北纬39.20 东经115.51
    河北省 蒿城 北纬38.02 东经114.50
    河北省 邯郸 北纬36.36 东经114.28
    河北省 河间 北纬38.26 东经116.05
    河北省 衡水 北纬37.44 东经115.42
    河北省 黄骅 北纬38.21 东经117.21
    河北省 晋州 北纬38.02 东经115.02
    河北省 冀州 北纬37.34 东经115.33
    河北省 廓坊 北纬39.31 东经116.42
    河北省 鹿泉 北纬38.04 东经114.19
    河北省 南宫 北纬37.22 东经115.23
    河北省 秦皇岛 北纬39.55 东经119.35
    河北省 任丘 北纬38.42 东经116.07
    河北省 三河 北纬39.58 东经117.04
    河北省 沙河 北纬36.51 东经114.30
    河北省 深州 北纬38.01 东经115.32
    河北省 唐山 北纬39.36 东经118.11
    河北省 武安 北纬36.42 东经114.11
    河北省 邢台 北纬37.04 东经114.30
    河北省 辛集 北纬37.54 东经115.12
    河北省 新乐 北纬38.20 东经114.41
    河北省 张家口 北纬40.48 东经114.53
    河北省 涿州 北纬39.29 东经115.59
    河北省 遵化 北纬40.11 东经117.58
    黑龙江省 哈尔滨 北纬45.44 东经126.36
    黑龙江省 阿城 北纬45.32 东经126.58
    黑龙江省 安达 北纬46.24 东经125.18
    黑龙江省 北安 北纬48.15 东经126.31
    黑龙江省 大庆 北纬46.36 东经125.01
    黑龙江省 富锦 北纬47.15 东经132.02
    黑龙江省 海林 北纬44.35 东经129.21
    黑龙江省 海伦 北纬47.28 东经126.57
    黑龙江省 鹤岗 北纬47.20 东经130.16
    黑龙江省 黑河 北纬50.14 东经127.29
    黑龙江省 佳木斯 北纬46.47 东经130.22
    黑龙江省 鸡西 北纬45.17 东经130.57
    黑龙江省 密山 北纬45.32 东经131.50
    黑龙江省 牡丹江 北纬44.35 东经129.36
    黑龙江省 讷河 北纬48.29 东经124.51
    黑龙江省 宁安 北纬44.21 东经129.28
    黑龙江省 齐齐哈尔 北纬47.20 东经123.57
    黑龙江省 七台河 北纬45.48 东经130.49
    黑龙江省 双城 北纬45.22 东经126.15
    黑龙江省 尚志 北纬45.14 东经127.55
    黑龙江省 双鸭山 北纬46.38 东经131.11
    黑龙江省 绥芬河 北纬44.25 东经131.11
    黑龙江省 绥化 北纬46.38 东经126.59
    黑龙江省 铁力 北纬46.59 东经128.01
    黑龙江省 同江 北纬47.39 东经132.30
    黑龙江省 五常 北纬44.55 东经127.11
    黑龙江省 五大连池 北纬48.38 东经126.07
    黑龙江省 伊春 北纬47.42 东经128.56
    黑龙江省 肇东 北纬46.04 东经125.58
    湖北省 武汉 北纬30.35 东经114.17
    湖北省 安陆 北纬31.15 东经113.41
    湖北省 当阳 北纬30.50 东经111.47
    湖北省 丹江口 北纬32.33 东经108.30
    湖北省 大冶 北纬30.06 东经114.58
    湖北省 恩施 北纬30.16 东经109.29
    湖北省 鄂州 北纬30.23 东经114.52
    湖北省 广水 北纬31.37 东经113.48
    湖北省 洪湖 北纬29.48 东经113.27
    湖北省 黄石 北纬30.12 东经115.06
    湖北省 黄州 北纬30.27 东经114.52
    湖北省 荆门 北纬31.02 东经112.12
    湖北省 荆沙 北纬30.18 东经112.16
    湖北省 老河口 北纬32.23 东经111.40
    湖北省 利川 北纬30.18 东经108.56
    湖北省 麻城 北纬31.10 东经115.01
    湖北省 浦圻 北纬29.42 东经113.51
    湖北省 潜江 北纬30.26 东经112.53
    湖北省 石首 北纬29.43 东经112.24
    湖北省 十堰 北纬32.40 东经110.47
    湖北省 随州 北纬31.42 东经113.22
    湖北省 天门 北纬60.39 东经113.10
    湖北省 武穴 北纬29.51 东经115.33
    湖北省 襄樊 北纬32.02 东经112.08
    湖北省 咸宁 北纬29.53 东经114.17
    湖北省 仙桃 北纬30.22 东经113.27
    湖北省 孝感 北纬30.56 东经113.54
    湖北省 宜昌 北纬30.42 东经111.17
    湖北省 宜城 北纬31.42 东经112.15
    湖北省 应城 北纬30.57 东经113.33
    湖北省 枣阳 北纬32.07 东经112.44
    湖北省 枝城 北纬30.23 东经111.27
    湖北省 钟祥 北纬31.10 东经112.34
    江苏省 南京 北纬32.03 东经118.46
    江苏省 常熟 北纬31.39 东经120.43
    江苏省 常州 北纬31.47 东经119.58
    江苏省 丹阳 北纬32.00 东经119.32
    江苏省 东台 北纬32.51 东经120.19
    江苏省 高邮 北纬32.47 东经119.27
    江苏省 海门 北纬31.53 东经121.09
    江苏省 淮安 北纬33.30 东经119.09
    江苏省 淮阴 北纬33.36 东经119.02
    江苏省 江都 北纬32.26 东经119.32
    江苏省 姜堰 北纬32.34 东经120.08
    江苏省 江阴 北纬31.54 东经120.17
    江苏省 靖江 北纬32.02 东经120.17
    江苏省 金坛 北纬31.46 东经119.33
    江苏省 昆山 北纬31.23 东经120.57
    江苏省 连去港 北纬34.36 东经119.10
    江苏省 溧阳 北纬31.26 东经119.29
    江苏省 南通 北纬32.01 东经120.51
    江苏省 邳州 北纬34.19 东经117.59
    江苏省 启乐 北纬31.48 东经121.39
    江苏省 如皋 北纬32.23 东经120.33
    江苏省 宿迁 北纬33.58 东经118.18
    江苏省 苏州 北纬31.19 东经120.37
    江苏省 太仓 北纬31.27 东经121.06
    江苏省 泰兴 北纬32.10 东经120.01
    江苏省 泰州 北纬32.30 东经119.54
    江苏省 通州 北纬32.05 东经121.03
    江苏省 吴江 北纬31.10 东经120.39
    江苏省 无锡 北纬31.34 东经120.18
    江苏省 兴化 北纬32.56 东经119.50
    江苏省 新沂 北纬34.22 东经118.20
    江苏省 徐州 北纬34.15 东经117.11
    江苏省 盐在 北纬33.22 东经120.08
    江苏省 扬中 北纬32.14 东经119.49
    江苏省 扬州 北纬32.23 东经119.26
    江苏省 宜兴 北纬31.21 东经119.49
    江苏省 仪征 北纬32.16 东经119.10
    江苏省 张家港 北纬31.52 东经120.32
    江苏省 镇江 北纬32.11 东经119.27
    内自治区 呼和浩特 北纬40.48 东经111.41
    内自治区 包头 北纬40.39 东经109.49
    内自治区 赤峰 北纬42.17 东经118.58
    内自治区 东胜 北纬39.48 东经109.59
    内自治区 二连浩特 北纬43.38 东经111.58
    内自治区 额尔古纳 北纬50.13 东经120.11
    内自治区 丰镇 北纬40.27 东经113.09
    内自治区 根河 北纬50.48 东经121.29
    内自治区 海拉尔 北纬49.12 东经119.39
    内自治区 霍林郭勒 北纬45.32 东经119.38
    内自治区 集宁 北纬41.02 东经113.06
    内自治区 临河 北纬40.46 东经107.22
    内自治区 满洲里 北纬49.35 东经117.23
    内自治区 通辽 北纬43.37 东经122.16
    内自治区 乌兰浩特 北纬46.03 东经122.03
    内自治区 乌海 北纬39.40 东经106.48
    内自治区 锡林浩特 北纬43.57 东经116.03
    内自治区 牙克石 北纬49.17 东经120.40
    内自治区 扎兰屯 北纬48.00 东经122.47
    宁夏自治区 银川 北纬38.27 东经106.16
    宁夏自治区 青铜峡 北纬37.56 东经105.59
    宁夏自治区 石嘴山 北纬39.02 东经106.22
    宁夏自治区 吴忠 北纬37.59 东经106.11
    青海省 西宁 北纬36.38 东经101.48
    青海省 德令哈 北纬37.22 东经97.23
    青海省 格尔木 北纬36.26 东经94.55
    山东省 济南 北纬36.40 东经117.00
    山东省 安丘 北纬36.25 东经119.12
    山东省 滨州 北纬37.22 东经118.02
    山东省 昌邑 北纬39.52 东经119.24
    山东省 德州 北纬37.26 东经116.17
    山东省 东营 北纬37.27 东经118.30
    山东省 肥城 北纬36.14 东经116.46
    山东省 高密 北纬36.22 东经119.44
    山东省 菏泽 北纬35.14 东经115.26
    山东省 胶南 北纬35.53 东经119.58
    山东省 胶州 北纬36.17 东经120.00
    山东省 即墨 北纬36.22 东经120.28
    山东省 济宁 北纬35.23 东经116.33
    山东省 莱芜 北纬36.12 东经117.40
    山东省 莱西 北纬36.52 东经120.31
    山东省 莱阳 北纬36.58 东经120.42
    山东省 莱州 北纬37.10 东经119.57
    山东省 乐陵 北纬37.44 东经117.12
    山东省 聊城 北纬36.26 东经115.57
    山东省 临清 北纬36.51 东经115.42
    山东省 临沂 北纬35.03 东经118.20
    山东省 龙口 北纬37.39 东经120.21
    山东省 蓬莱 北纬37.48 东经120.45
    山东省 平度 北纬36.47 东经119.58
    山东省 青岛 北纬36.03 东经120.18
    山东省 青州 北纬36.42 东经118.28
    山东省 曲阜 北纬35.36 东经116.58
    山东省 日照 北纬35.23 东经119.32
    山东省 荣成 北纬37.10 东经122.25
    山东省 乳山 北纬36.54 东经121.31
    山东省 寿光 北纬36.53 东经118.44
    山东省 泰安 北纬36.11 东经117.08
    山东省 滕州 北纬35.06 东经117.09
    山东省 潍坊 北纬36.43 东经119.06
    山东省 威海 北纬37.31 东经122.07
    山东省 文登 北纬37.12 东经122.03
    山东省 新泰 北纬35.54 东经117.45
    山东省 烟台 北纬37.32 东经121.24
    山东省 兖州 北纬35.32 东经116.49
    山东省 禹城 北纬36.56 东经116.39
    山东省 枣庄 北纬34.52 东经117.33
    山东省 章丘 北纬36.43 东经117.32
    山东省 招远 北纬37.21 东经120.23
    山东省 诸城 北纬35.59 东经119.24
    山东省 淄博 北纬36.48 东经118.03
    山东省 邹城 北纬35.24 东经116.58
    陕西省 西安 北纬34.17 东经108.57
    陕西省 安康 北纬32.41 东经109.01
    陕西省 宝鸡 北纬34.22 东经107.09
    陕西省 韩城 北纬35.28 东经110.27
    陕西省 汉中 北纬33.04 东经107.01
    陕西省 华阴 北纬34.34 东经110.05
    陕西省 商州 北纬33.52 东经109.57
    陕西省 铜川 北纬35.06 东经109.07
    陕西省 渭南 北纬34.30 东经109.30
    陕西省 咸阳 北纬34.20 东经108.43
    陕西省 兴平 北纬34.18 东经108.29
    陕西省 延安 北纬36.35 东经109.28
    陕西省 榆林 北纬38.18 东经109.47
    上海市 上海市 北纬31.14 东经121.29
    四川省 成都 北纬30.40 东经104.04
    四川省 巴中 北纬31.51 东经106.43
    四川省 崇州 北纬30.39 东经103.40
    四川省 达川 北纬31.14 东经107.29
    四川省 德阳 北纬31.09 东经104.22
    四川省 都江堰 北纬31.01 东经103.37
    四川省 峨眉山 北纬29.36 东经103.29
    四川省 涪陵 北纬29.42 东经107.22
    四川省 广汉 北纬30.58 东经104.15
    四川省 广元 北纬32.28 东经105.51
    四川省 华蓥 北纬30.26 东经106.44
    四川省 简阳 北纬30.24 东经104.32
    四川省 江油 北纬31.48 东经104.42
    四川省 阆中 北纬31.36 东经105.58
    四川省 乐山 北纬29.36 东经103.44
    四川省 泸州 北纬28.54 东经105.24
    四川省 绵阳 北纬31.30 东经104.42
    四川省 南充 北纬30.49 东经106.04
    四川省 内江 北纬29.36 东经105.02
    四川省 攀枝花 北纬26.34 东经101.43
    四川省 彭州 北纬30.59 东经103.57
    四川省 邛崃 北纬30.26 东经103.28
    四川省 遂宁 北纬30.31 东经105.33
    四川省 万县 北纬30.50 东经108.21
    四川省 万源 北纬32.03 东经108.03
    四川省 西昌 北纬27.54 东经102.16
    四川省 雅安 北纬29.59 东经102.59
    四川省 宜宾 北纬28.47 东经104.34
    四川省 自贡 北纬29.23 东经104.46
    四川省 资阳 北纬30.09 东经104.38
    台湾省 台北市 北纬25.03 东经121.30
    天津市 天津市 北纬39.02 东经117.12
    西藏自治区 拉萨 北纬29.39 东经91.08
    西藏自治区 日喀则 北纬29.16 东经88.51
    香港 香港市 北纬21.23 东经115.12
    新疆自治区 乌鲁木齐 北纬43.45 东经87.36
    新疆自治区 阿克苏 北纬41.09 东经80.19
    新疆自治区 阿勒泰 北纬47.50 东经88.12
    新疆自治区 阿图什 北纬39.42 东经76.08
    新疆自治区 博乐 北纬44.57 东经82.08
    新疆自治区 昌吉 北纬44.02 东经87.18
    新疆自治区 阜康 北纬44.09 东经87.58
    新疆自治区 哈密 北纬42.50 东经93.28
    新疆自治区 和田 北纬37.09 东经79.55
    新疆自治区 克拉玛依 北纬45.36 东经84.51
    新疆自治区 喀什 北纬39.30 东经75.59
    新疆自治区 库尔勒 北纬41.46 东经86.07
    新疆自治区 奎屯 北纬44.27 东经84.56
    新疆自治区 石河子 北纬44.18 东经86.00
    新疆自治区 塔城 北纬46.46 东经82.59
    新疆自治区 吐鲁番 北纬42.54 东经89.11
    新疆自治区 伊宁 北纬43.55 东经81.20
    云南省 昆明 北纬25.04 东经102.42
    云南省 保山 北纬25.08 东经99.10
    云南省 楚雄 北纬25.01 东经101.32
    云南省 大理 北纬25.34 东经100.13
    云南省 东川 北纬26.06 东经103.12
    云南省 个旧 北纬23.21 东经103.09
    云南省 景洪 北纬22.01 东经100.48
    云南省 开远 北纬23.43 东经103.13
    云南省 曲靖 北纬25.30 东经103.48
    云南省 瑞丽 北纬24.00 东经97.50
    云南省 思茅 北纬22.48 东经100.58
    云南省 畹町 北纬24.06 东经98.04
    云南省 宣威 北纬26.13 东经104.06
    云南省 玉溪 北纬24.22 东经102.32
    云南省 昭通 北纬27.20 东经103.42
    浙江省 杭州 北纬30.16 东经120.10
    浙江省 慈溪 北纬30.11 东经121.15
    浙江省 东阳 北纬29.16 东经120.14
    浙江省 奉化 北纬29.39 东经121.24
    浙江省 富阳 北纬30.03 东经119.57
    浙江省 海宁 北纬30.32 东经120.42
    浙江省 湖州 北纬30.52 东经120.06
    浙江省 建德 北纬29.29 东经119.16
    浙江省 江山 北纬28.45 东经118.37
    浙江省 嘉兴 北纬30.46 东经120.45
    浙江省 金华 北纬29.07 东经119.39
    浙江省 兰溪 北纬29.12 东经119.28
    浙江省 临海 北纬28.51 东经121.08
    浙江省 丽水 北纬28.27 东经119.54
    浙江省 龙泉 北纬28.04 东经119.08
    浙江省 宁波 北纬29.52 东经121.33
    浙江省 平湖 北纬30.42 东经121.01
    浙江省 衢州 北纬28.58 东经118.52
    浙江省 瑞安 北纬27.48 东经120.38
    浙江省 上虞 北纬30.01 东经120.52
    浙江省 绍兴 北纬30.00 东经120.34
    浙江省 台州 北纬28.41 东经121.27
    浙江省 桐乡 北纬30.38 东经120.32
    浙江省 温岭 北纬28.22 东经121.21
    浙江省 温州 北纬28.01 东经120.39
    浙江省 萧山 北纬30.09 东经120.16
    浙江省 义乌 北纬29.18 东经120.04
    浙江省 乐清 北纬28.08 东经120.58
    浙江省 余杭 北纬30.26 东经120.18
    浙江省 余姚 北纬30.02 东经121.10
    浙江省 永康 北纬29.54 东经120.01
    浙江省 舟山 北纬30.01 东经122.06
    浙江省 诸暨 北纬29.43 东经120.14
    重庆市 重庆市 北纬29.35 东经106.33
    重庆市 合川市 北纬30.02 东经106.15
    重庆市 江津市 北纬29.18 东经106.16
    重庆市 南川市 北纬29.10 东经107.05
    重庆市 永川市 北纬29.23 东经105.53