ArcGIS用Field Calculator 批量处理shp属性，添加标准图幅号（并且标记跨图幅图斑）

这个可以用arcpy做，但太麻烦了，想了个简单的法子，直接用属性计算器计算了。
我用的是10.2版本的arcgis，而且是没汉化的英文版
主要分两部分，一部分是提取图斑包络四边形顶点，另一个是计算图幅号

转地理坐标系

如果发现计算的坐标不是经纬度，那得转地理坐标系。我转的是CGCS2000，大家按照数据本来的坐标系转对应的坐标系


不想转的话也可以，直接改图层frame属性

到时候计算的时候选下面的选项：

提取图斑包络四边形顶点

方法参考于：liuxupiaoshi的博客文章
第一步、新建一个属性ID 直接从零开始赋值（等于FID就可以了）
第二步、生成包络图形

第三步、提取特征点


第四步、新建X 、 Y 属性 ，计算点坐标

注意计算的时候坐标系别搞错了

第五步、把点特征叠加到包络矩形，并作汇总
这里用的是相交工具：



对结果用提前设置好的ID作汇总


输出的是一个dbf表格

第六步、把这个表格用接回去
在原始图斑的属性表内利用ID属性进行join操作

结果如下

计算图幅号

其中，图幅号计算方法思路参考了darlun的博客文章中的代码

第一步、新建多个属性字段：
除了那个计数的属性，其他都是TEXT属性的

第二步，计算四个角点的图幅号，并用逗号隔开输出到 “总输出“ 属性
打开Field Calculator

开始写脚本

LatCharArray = ['A','B','C','D','E','F','G','H','I','J','K','L','M','N','O','P','Q','R','S','T','U','V']

#图幅比例尺  
scaleArray = [1000000,500000,250000,100000,50000,25000,10000,5000]

#图幅纬差,单位秒 
latDArray = [14400,7200,3600,1200,600,300,150,75]

#图幅经差,单位秒 
lonDArray = [21600,10800,5400,1800,900,450,225,112.5]

SCALE_1000000 = 0
SCALE_500000 = 1
SCALE_250000 = 2
SCALE_100000 = 3
SCALE_50000 = 4
SCALE_25000 = 5
SCALE_10000 = 6
SCALE_5000 = 7

#lat 为 纬度， lon为经度   1：5000的scaleID为7
def getSheetNumber(lat,lon,scaleID):
	f = lat*3600;		#以秒表示的纬度
	r = lon*3600;		#以秒表示的经度
	numString = ''
	a = int(f/(4*3600));   		#纬度序号
	b = int(r/(6*3600))+31;   	#经度序号  
	c = int(4*3600/latDArray[scaleID])-int((f%(4*3600))/latDArray[scaleID]); #在1:1000000万图幅中的序号 纬度 
	d = int((r%(6*3600))/lonDArray[scaleID])+1;   #在1:1000000万图幅中的序号  经度
			
	numString += LatCharArray[a];
	numString += ""+ str(b);


	#1:1000000万
	if(scaleID==0):
		return numString;    

	numString +=LatCharArray[scaleID];
	if c<10:         
		numString += "00"
		numString += str(c)
	elif c<100:   
		numString += "0"
		numString += str(c)
	elif c<1000:
		numString += str(c) 
	else:   
		return ""    
			      
	if d<10:        
		numString += "00"      
		numString += str(d) 
	elif d<100:      
		numString += "0" 
		numString += str(d) 
	elif d<1000:        
		numString += str(d);   
	else:     
		return ""
	return numString



def f(xmin,xmax,ymin,ymax):

	return getSheetNumber(ymax,xmax,7)+ ',' + getSheetNumber(ymax,xmin,7)+ ',' + getSheetNumber(ymin,xmax,7)+ ',' + getSheetNumber(ymin,xmin,7)

# F1 F2 是我自己写的小代码，
# 补零位
def F1(a):
	if a<10:
		strout = '00' + str(int(a))
	elif a < 100:
		strout = '0' + str(int(a))
	else:
		strout = str(int(a))
	return strout

def F2(x,y):
	#区域编号头部
	str1 = "I50H"
	#1:5000标准经差
	d1 = 1.0/60 + 52.5/3600
	#标准纬度差
	d2 = 1.0/60 + 15.0/3600
	c = 4/d2 - int(( y % 4) /d2)
	d = int((x % 6) / d1) + 1
	str2 = F1(c)
	str3 = F1(d)
	str4 = str1 + str2 + str3
	return str4

用field计算器开始运行代码，把代码粘到对应的框子里面



结果

第三步，计算重复与输出结果的中间步骤
中间输出格式为 [跨图幅数],[所有的图幅号]
代码是

#中间
def f(str0):
	str0 = str0.split(',')
	a = str0[0]
	b = str0[1]
	c = str0[2]
	d = str0[3]
	count = 1

	str1 = a
	if b!=a:
		count +=1
		str1 =str1 +  ',' + b
	if c !=a and c != b:
		count +=1
		str1 =str1 +  ',' + c
	if d != a and d != b and d != c:
		count +=1
		str1 =str1 +   ',' + d
	return str(count) +  ',' + str1

第四步，输出结果
分别提取跨图幅数和图幅号

#count
def f (a):
	return int(a[0])

#图幅号
def f (a):
	return a[2:]

结果如下

第五步，修饰
解除join

按照不同的图幅号显示颜色；欧克 完成了