日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 人文社科 > 生活经验 >内容正文

生活经验

python爬取晋江_[Arcpy] 爬取高德地图面状地物(半自动化)—获取任意地区的行政边界...

發(fā)布時間:2023/11/27 生活经验 44 豆豆
生活随笔 收集整理的這篇文章主要介紹了 python爬取晋江_[Arcpy] 爬取高德地图面状地物(半自动化)—获取任意地区的行政边界... 小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.

# -*- coding: utf-8 -*-

# @Time : 2019/07/30

# @Author : PasserQi

# @Email : passerqi@gmail.com

# @File : param_poiinfo

# @Software: PyCharm

# @Version :

# @Desc : 解析高德地圖poiInfo

# 例子:https://www.amap.com/service/poiInfo?query_type=IDQ&pagesize=20&pagenum=1&qii=true&cluster_state=5&need_utd=true&utd_sceneid=1000&div=PC1000&addr_poi_merge=true&is_classify=true&zoom=12&id=B0FFFRKVWU&city=110105&geoobj=117.983877%7C24.408902%7C118.346426%7C24.578237&keywords=%E4%B8%89%E7%8E%AF

import os

import sys

import arcpy, arcgisscripting

from Tools import coordinate_conversion

# coordinate_conversion為坐標轉換庫

# 在工程目錄下創(chuàng)建Tools文件夾,再創(chuàng)建一個coordinate_conversion.py文件,將下列鏈接的代碼放入

# 鏈接:https://blog.csdn.net/summer_dew/article/details/80723434

# python2編碼問題

defaultencoding = 'utf-8'

if sys.getdefaultencoding() != defaultencoding:

reload(sys)

sys.setdefaultencoding(defaultencoding)

# 參數(shù)修改

out_dir = r"C:\Users\PasserQi\Desktop" #保存路徑

out_name = u"晉江市.shp" #保存文件名

# 從F12控制臺中復制的經緯度坐標串

shp_str = "118.586651;24.887213;118.589856;24.879356;118.593491;24.875708;118.599304;24.872889;118.605415;24.871446;118.624269;24.869941;118.636227;24.866015;118.643555;24.860675;118.645981;24.853015;118.651331;24.845679;118.672849;24.849566;118.683505;24.853791;118.687718;24.856569;118.691014;24.860458;118.694682;24.866351;118.695950;24.867216;118.698993;24.867843;118.702429;24.865521;118.701906;24.860221;118.698566;24.852632;118.699479;24.847657;118.705192;24.844329;118.717886;24.841981;118.729753;24.841322;118.732901;24.838998;118.734929;24.826089;118.730478;24.821229;118.721094;24.815886;118.662294;24.802783;118.654385;24.786973;118.634813;24.777266;118.633370;24.780991;118.631945;24.782315;118.626722;24.780853;118.626513;24.777269;118.630588;24.766748;118.630501;24.760373;118.629336;24.757815;118.624140;24.756222;118.620786;24.753992;118.611424;24.753472;118.601873;24.745861;118.590506;24.745574;118.589022;24.744357;118.590419;24.738893;118.587148;24.734097;118.584997;24.732579;118.584364;24.727241;118.585118;24.724764;118.594070;24.720963;118.597384;24.715788;118.598217;24.710946;118.607938;24.705602;118.614243;24.705638;118.618699;24.710424;118.627161;24.711369;118.629420;24.712478;118.636956;24.719601;118.646632;24.718483;118.649363;24.720204;118.653519;24.719240;118.654841;24.716864;118.654058;24.713976;118.654771;24.711012;118.656701;24.707392;118.662572;24.703136;118.661858;24.699239;118.659753;24.695504;118.659457;24.692357;118.661979;24.686969;118.662231;24.684593;118.670172;24.679759;118.671894;24.675294;118.677487;24.667206;118.687014;24.660570;118.694317;24.654833;118.699377;24.648414;118.696846;24.628668;118.689942;24.584811;118.670581;24.547085;118.606412;24.507367;118.565712;24.500086;118.527168;24.513751;118.486898;24.542900;118.467133;24.559559;118.463662;24.564981;118.444639;24.614950;118.439675;24.627817;118.438801;24.643452;118.448609;24.659979;118.448791;24.664517;118.446418;24.674074;118.446705;24.691943;118.449859;24.705397;118.448084;24.708736;118.446169;24.708721;118.444366;24.712447;118.445809;24.715681;118.443717;24.721596;118.440145;24.725632;118.438648;24.729647;118.435093;24.730647;118.437170;24.736516;118.437231;24.742043;118.434615;24.753379;118.432986;24.755409;118.426753;24.756416;118.424586;24.759202;118.427950;24.764948;118.433767;24.768907;118.433802;24.771921;118.431609;24.775570;118.427769;24.777823;118.420381;24.776585;118.414336;24.776898;118.412149;24.778426;118.411074;24.780699;118.414839;24.784620;118.425490;24.790356;118.430258;24.799103;118.436013;24.795953;118.437443;24.796188;118.441377;24.802283;118.442132;24.813739;118.441283;24.818974;118.441708;24.825343;118.437048;24.835288;118.442298;24.839163;118.440141;24.839870;118.436788;24.843957;118.437802;24.845896;118.437629;24.849291;118.439076;24.850142;118.441346;24.849188;118.443772;24.849774;118.444967;24.849479;118.447947;24.845168;118.452840;24.847323;118.458156;24.844471;118.460537;24.845787;118.461628;24.849809;118.467324;24.849775;118.471089;24.852175;118.477666;24.852824;118.479613;24.854082;118.481949;24.858032;118.488213;24.858919;118.488723;24.860474;118.488109;24.862494;118.480886;24.870891;118.478991;24.875961;118.478948;24.879302;118.482427;24.883961;118.482609;24.890224;118.485014;24.894003;118.483821;24.900109;118.484652;24.902581;118.494151;24.903732;118.498484;24.898901;118.498735;24.896760;118.502143;24.894579;118.503934;24.895435;118.504790;24.898209;118.505058;24.896756;118.509504;24.893287;118.515083;24.890659;118.518396;24.887982;118.522523;24.879981;118.524487;24.879354;118.530449;24.884382;118.533408;24.884942;118.540254;24.881551;118.541543;24.879434;118.541499;24.876131;118.537327;24.870327;118.537214;24.867560;118.549264;24.867025;118.551264;24.869404;118.552001;24.873515;118.561944;24.883271;118.569222;24.888781;118.578600;24.888129;118.585676;24.889032;118.586651;24.887213"

out_path = os.path.join(out_dir, out_name)

gp = arcgisscripting.create() #GP

spat_ref = "4326" #坐標系

gp.CreateFeatureClass_management(out_dir, out_name, "POLYGON", "", "", "", spat_ref)

cur = gp.InsertCursor(out_path)

newRow = cur.newRow()

XYsStr = shp_str

XYarray = gp.CreateObject("array")

XYList = XYsStr.split(';')

for i in range(0, len(XYList), 2):

XY = [0.0 , 0.0]

XY[0], XY[1] = float(XYList[i]), float(XYList[i+1])

point = gp.CreateObject("point")

point.X, point.Y = coordinate_conversion.gcj02towgs84(XY[0], XY[1])

XYarray.add(point)

newRow.setValue("Shape", XYarray)

cur.InsertRow(newRow)

del cur,newRow

總結

以上是生活随笔為你收集整理的python爬取晋江_[Arcpy] 爬取高德地图面状地物(半自动化)—获取任意地区的行政边界...的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內容還不錯,歡迎將生活随笔推薦給好友。