# -*- coding: utf-8 -*- #by gisoracle 2021.01.23import arcpy import math #保存数据def Save(polygon,rows):row = rows.newRow()#row.setValue(shapefieldname,featureList)row.shape=polygonrows.insertRow(row)del rowx=500000 y=0 r=200.0 pi=3.1415926 outFeature=arcpy.GetParameterAsText(0) array = arcpy.Array() #point = arcpy.Point()for i in range(0,180):x1=x+r*math.cos(i*pi/180)Z=r*math.sin(i*pi/180)point...
年前写了一个用渔网工具制作图幅接图表的文章，链接在这里：使用ArcMap做一个1:5000标准分幅图并编号本文提供一个使用ArcMap利用接图表图斑裁切一幅影像为多幅的方法。第一步，将接图表拆分，使之每一个图斑为一个要素类（Shapefile文件）拆分方法很多，提供两个参考：1、可以使用Python脚本2、可以使用ArcToolbox——分析工具——split拆分的结果是每个影像框图斑为一个Shapefile文件，所有Shapefile文件生成在一个文件夹下。第二...
本文翻译自：Qingkai‘s Blog当使用python的Basemap库绘制地图时，选择一个漂亮的底图会为图片增色不少，但是使用map.bluemarble()、map.etopo()或者map.shadedrelief()等函数时，由于分辨率的原因，将图片缩小到较小范围会使得底图非常模糊。其实，创建一个高分辨率的底图的方式时使用arcgisimage方法。您可以在这里看到Basemap的详细说明。我将在下边的例子当中展示几个非常漂亮的例子。为了使用arcgisimage方法，你需要为不同的...
最近要处理个简单数据，一个字段中为文本类型，包含各种描述。要求是包含平方米的数值提取出来，变成数值，如果包含多个，则把各个值累加起来。比如 字段值为 “非法占用100平方米” 处理后结果为 100字段值为 “工棚3间” 处理结果为0字段值为 “工棚3间，瓦屋2间，非法占用100平方米” 处理结果100字段值为 “工棚100平方米，瓦屋100平方米” 处理结果为200 由于一直写.net，上述功能在.net中做起来并不麻烦。但是觉得应该是字段...
arcgis python 刷新arcpy.RefreshActiveView() 刷新地图和布局窗口arcpy.RefreshTOC() 刷新内容列表arcpy.RefreshCatalog(r"F:\tknew106\tool") 刷新目录树??原文：https://blog.51cto.com/u_12139363/3023493
# -*- coding: cp936 -*-import arcpy import os import ylpy import random def main():num=ylpy.getCount(inFeature)n=int(num*mscale/100)ylpy.clearselect(inFeature)A=random.sample(range(num),n)FID=ylpy.getOIDField(inFeature)sql=str(FID)+" in ("for i in range(0,n):sql=sql+str(A[i])if i<n-1:sql=sql+","sql=sql+")"arcpy.Select_analysis(inFeature,outFeature,sql)inFeature = arcpy.GetParameterAsText(0) # msca...
属性说明area 面要素的面积值；对于所有其他要素类型为空centroid 如果质心位于要素之内或要素之上则为真；否则返回标注点（返回一个点对象）范围 返回一个范围对象firstPoint 要素的第一个坐标；返回一个点对象hullRectangle 凸包矩形的坐标对isMultipart 如果此几何的部分数大于一，则为真labelPoint 标注所在的点；返回一个点对象（标注点始终位于要素之内或要素之上）。 lastPoint 要素的最后一个坐标；返回一个点对象长度 线状...
一、程序代码：import arcpy arcpy.CheckOutExtension("spatial") arcpy.gp.overwriteOutput=1 arcpy.env.workspace = "F:\\Modis_16\\1Moasic" rasters = arcpy.ListRasters("*", "tif") mask= "F:\\Vegetation Change\\Data\\Bound\\bound_Buffer_Polygon.shp" for raster in rasters:print(raster)out= "F:\\Vegetation Change\\Data\\GIMMS Data\\new\\"+"ma_"+raster[0:8]arcpy.gp.ExtractByMask_sa(raster, mask, out)print(...
一、引言在arcgis打开一个图层的属性表，可以对属性表的某个字段进行计算，但是在平常一般都是使用arcgis提供的字段计算器的界面进行傻瓜式的简答的赋值操作，并没有使用到脚本对字段值进行逻辑的操作。由于最近一直在学python脚本，刚好又碰上一好基友需要我的助攻（使用arcgis制图），这就用上了。本以为能够轻松搞定的，没想到搬石头砸脚了，下面就来说我是如何被砸脚的吧。 二、问题描述：将test字段中值为“湖南”的变为“湖南...
import arcpydesc = arcpy.Describe(r‘C:\Users\dell\Documents\ArcGIS\default.gdb‘) print desc.release地理数据库版本值ArcGIS 版本2.2.09.22.3.09.3、9.3.13.0.010.0、10.1、10.2、10.2.1、10.2.2、10.3、10.3.1、10.4、10.4.1、10.5、10.5.1 ArcGIS Pro 1.0、1.1、1.2、1.3、1.4ArcGIS Python查看geodatabase（.gdb）版本标签：gis odata head 不能 10.2.1 keyword des dbf geo 本文系统来源：https://ww...
Python和ArcGIS自动化制图完全指南（三）：分配模板 前言：在《指南》第二章中，我们制作了3种不同尺寸的 mxd 模板，就是为了适配不同的制图单位，比如适配自贡市的 mxd 模板就不适用于广元市（广元市面积大）。而如何让计算机给不同制图单位分配大小合适的模板就是本章需要解决的问题 。 文章目录 Python和ArcGIS自动化制图完全指南（三）：分配模板1.模板信息标准化2.最小边界几何3.更新**MBG_Width** 和 **MBG_Length** 字段4.更...
1 模型的创建与保存 可以在ARCGIS当中通过右键单击现有工具箱并选择，新建 > 模型，来创建一个新模型。这样便可在该工具箱中创建一个具有默认名称的模型，并打开该模型以执行编辑。 2 模型基本元素 2.1 工具 ARCGIS当中的数据处理工具可以导入到模型当中，可以在模型当中设置各个工具所需的输入输出变量。 note： 1，右键->参数，可将值设置为参数，当模型保存后打开时，会在对话框当中显示该参数 2，输出变量不会自动添加到地图文...
# -*- coding: utf-8 -*- import arcpy import sys def updatebsm(inFeature,inField):i=1cursor=arcpy.UpdateCursor(inFeature)for row in cursor:row.setValue(inField,i)i=i+1cursor.updateRow(row)del cursordef updatebsm2(inFeature,inField):i=1cursor=arcpy.da.UpdateCursor(inFeature,inField)for row in cursor:row[0]=ii=i+1cursor.updateRow(row)del cursor inFeature = arcpy.GetParameterAsText(0) # inField = arcp...
1.模型生成的二值预测图如图所示（已将原始数据的坐标系通过gdal转到预测图中）2.采用如下代码将其转化为面矢量 1 def polygonize(inImgPath, outImgPath):2 inImg = CXImage()3 inImg.Open(inImgPath)4 srcband = inImg.proDataset.GetRasterBand(1)5 prj = osr.SpatialReference()6 prj.ImportFromWkt(inImg.proDataset.GetProjection())7 maskband = srcband.GetMaskBand()8 9 driver = ogr.GetDri...
Python和ArcGIS自动化制图完全指南（四）：自动制图 前言：在完成了《指南》第二章和第三章后，终于来到了制图。运用 python 和 arcpy 自动化完成 mxd 模板图层定义查询语句的自动更新、图框的自动居中、比例尺的校正、另存mxd等操作。 关键字：python、arcpy、mxd PAGESIZE 字段是 定义查询 语句中使用的字段。作用是给每个 mxd 模板分配任务。 文章目录 Python和ArcGIS自动化制图完全指南（四）：自动制图1.限制 MXD 模板2.遍历...