用arcpy做个退地类计算方法

已于 2025-01-24 11:01:31 修改
阅读量402 收藏 8
点赞数 3
分类专栏： 自然资源 文章标签： python
于 2024-11-15 12:45:47 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_15607445/article/details/143794045
版权

用arcpy做个退地类计算方法

退地类算法示意图

方法说明

本次采用arcpy方法进行计算分析，通过调用GeoProcessing达到分析的目的。

为减少GeoProcessing调用次数，大概的方式是，统一处理全覆盖的面（所有的DLTB），再分别处理非全覆盖面

为缩小分析范围，先采用SelectLayerByLocation的方式控制下分析范围，再进行分析。

具体思路如下：

  1. 先用待分析图形，选中所有的图层（DLTB,CZCDYD）
  2. 用待分析图形clip一下2023的DLTB
  3. 用2023_clip数据进行叠加所有的DLTB，生成2023_clip_union
  4. 分别用待分析图形clip一下CZCDYD
  5. 再用union的方法，将所有的分析过程数据（2023_clip_union和所有的CZCDYD_clip）进行union

剩下的就是统计了，可统计结果：

  • 2023的核定、非核定现状分析（含权属地类面积）
  • 2022的核定、非核定现状分析（含权属地类面积）
  • 2023的退地类分析

后续扩展

1、把平差的方法加进来，做好地类平差
2、可改成arctool，作为一个自定义分析方法，输出如上的结果。也可改成服务以供调用。
3、代码需要修改的地方：临时图层的命名，可以加上时间戳，防止出现重复，代码执行前做删除判断，防止重复，再代码执行后执行删除，防止数据冗余变大。

# 待分析地块，先选择下2023的现状图层，然后进行叠加分析
# 选择属于建设用地的部分，然后再与2018和2009的图层进行分析
# 选择数据其他农用地的，
import arcpy
import time

f_time=time.time()
# 地类字典
# 三调地类字典
#耕地，其他农用地，建设用地，未利用地
DL_SD_GD="('0101','0102','0103')"
DL_SD_QTNYD="('0201','0201K','0202','0202K','0203','0203K','0204','0204K','0301','0301K','0302','0302K','0303','0304','0305','0306','0307','0307K','0401','0402','0403','0403K','0404','1006','1103','1104','1104A','1104K','1107','1107A','1202','1203')"
DL_SD_JSYD="('0501','0502','0503','0504','0505','0506','0507','0508','05H1','0601','0602','0603','0604','0701','0702','0801','0802','0803','0804','0805','0806','0807','0808','0809','0810','0810A','08H1','08H2','08H2A','09','0901','0902','0903','0904','0905','0906','1001','1002','1003','1004','1005','1007','1008','1009','1109','1201')"
DL_SD_WLYD="('1101','1102','1105','1106','1108','1110','1204','1205','1206','1207')"
# 二调地类字典
DL_ED_GD="('011','012','013')"
DL_ED_QTNYD="('021','022','023','031','032','033','041','042','104','114','117','122','123')"
DL_ED_JSYD="('051','052','053','054','061','062','063','071','072','081','082','083','084','085','086','087','088','091','092','093','094','095','101','102','103','105','106','107','113','118','121','201','202','203','204','205')"
DL_ED_WLYD="('111','112','115','116','119','043','124','125','126','127')"
# 读取下待分析的地块
arcpy.env.workspace = "C:\\Users\\Administrator\\Desktop\\GHTJ.gdb"
polygon_layer = "Exp"

xz_2023_layer="TDLY_DLTB_SD_2023"
xz_2022_layer="TDLY_DLTB_SD_2022"
xz_dltb_2018="TDLY_DLTB_2018"
xz_dltb_2009="TDLY_DLTB_2009"
czc_22= "CZCDYD_SD_2022"
czc_23="CZCDYD_SD_2023"

out_xz2023="xz2023_temp"
out_result="result_temp"
out_result_23="result_23_temp"
out_result_22="result_22_temp"
out_result_X="result_x_temp"
# 清除掉临时图层
try:
    if arcpy.Exists(out_xz2023):
        arcpy.management.Delete(out_xz2023)
    if arcpy.Exists(out_result):
        arcpy.management.Delete(out_result)
    if arcpy.Exists(out_result_23):
        arcpy.management.Delete(out_result_23)
    if arcpy.Exists(out_result_22):
        arcpy.management.Delete(out_result_22)
    if arcpy.Exists(out_result_X):
        arcpy.management.Delete(out_result_X)
    print("清空临时图层")
except exec as e : 
    print(e)


def GetArea(layer_name,where_sql):
    sum_area=0
    with arcpy.da.SearchCursor(layer_name,where_clause=where_sql,field_names=['SHAPE@AREA']) as jsydCursor:
        for row in jsydCursor:
            sum_area=sum_area+row[0]
    return sum_area


# 先与2023年的现状数据进行叠加分析,此处选用clip方法，防止待分析地块干扰

S_2023 = arcpy.management.SelectLayerByLocation(xz_2023_layer,overlap_type="INTERSECT",select_features=polygon_layer)
S_2022 = arcpy.management.SelectLayerByLocation(xz_2022_layer,overlap_type="INTERSECT",select_features=polygon_layer)
S_2018 = arcpy.management.SelectLayerByLocation(xz_dltb_2018,overlap_type="INTERSECT",select_features=polygon_layer)
S_2009=arcpy.management.SelectLayerByLocation(xz_dltb_2009,overlap_type="INTERSECT",select_features=polygon_layer)

arcpy.analysis.Clip(S_2023,polygon_layer,out_xz2023)
# 接下来进行叠加分析 全面叠加
arcpy.analysis.Intersect([out_xz2023,S_2022,S_2018,S_2009],out_result)
# 图形先clip23年的20x，  图形再clip22年的20x，      然后进行联合分析，输入要素包括（叠加的base面数据，23的clip数据，22的clip数据）
S_CZC_2023=arcpy.management.SelectLayerByLocation(czc_23,overlap_type="INTERSECT",select_features=polygon_layer)
arcpy.analysis.Clip(S_CZC_2023,polygon_layer,out_result_23)

S_CZC_2022=arcpy.management.SelectLayerByLocation(czc_22,overlap_type="INTERSECT",select_features=polygon_layer)
arcpy.analysis.Clip(S_CZC_2022,polygon_layer,out_result_22)

# 进行union分析
arcpy.analysis.Union([out_result,out_result_23,out_result_22], out_result_X)
# 计算建设用地面积 CZCLX不为空 且 DLBM_2,DLBM_3都是建设用地的，为建设用地
jsyd_where="(CZCLX <> '' OR DLBM IN %s) AND DLBM_12 IN %s AND DLBM_12_13 IN %s" %( DL_SD_JSYD, DL_ED_JSYD,DL_ED_JSYD)
jsydzmj=GetArea(out_result_X,jsyd_where)
# print("建设用地总面积为：" + jsydzmj)

#计算09-2-1  CZCLX不为空 且 DLBM_2 是建设用地 且 DLBM_3 是其他农用地
r_09_2_1_where="(CZCLX <> '' OR DLBM IN %s) AND DLBM_12 IN %s AND DLBM_12_13 IN %s" %(DL_SD_JSYD,DL_ED_JSYD,DL_ED_QTNYD)
r_09_2_1=GetArea(out_result_X,r_09_2_1_where)
# print("09-2-1面积为：" + r_09_2_1)

#计算09-2-2  CZCLX不为空 且 DLBM_2 是建设用地 且 DLBM_3 是耕地
r_09_2_2_where="(CZCLX <> '' OR DLBM IN %s) AND DLBM_12 IN %s AND DLBM_12_13 IN %s" %(DL_SD_JSYD,DL_ED_JSYD,DL_ED_GD)
r_09_2_2=GetArea(out_result_X,r_09_2_2_where)
# print("09-2-2面积为：" + r_09_2_2)

#计算09-2-3  CZCLX不为空 且 DLBM_2 是建设用地 且 DLBM_3 是未利用地
r_09_2_3_where="(CZCLX <> '' OR DLBM IN %s) AND DLBM_12 IN %s AND DLBM_12_13 IN %s" %(DL_SD_JSYD,DL_ED_JSYD,DL_ED_WLYD)
r_09_2_3=GetArea(out_result_X,r_09_2_3_where)
# print("09-2-3面积为：" + r_09_2_3)

#计算18-1  CZCLX不为空 且 DLBM_2是其他农用地
r_18_1_where="(CZCLX <> '' OR DLBM IN %s) AND DLBM_12 IN %s " %(DL_SD_JSYD,DL_ED_QTNYD)
r_18_1=GetArea(out_result_X,r_18_1_where)
# print("18-1面积为：" + r_18_1)

#计算18-2  CZCLX不为空 且 DLBM_2是耕地
r_18_2_where="(CZCLX <> '' OR DLBM IN %s) AND DLBM_12 IN %s " %(DL_SD_JSYD,DL_ED_GD)
r_18_2=GetArea(out_result_X,r_18_2_where)
# print("18-2面积为：" + r_18_2)

#计算18-3  CZCLX不为空 且 DLBM_2是未利用地
r_18_3_where="(CZCLX <> '' OR DLBM IN %s) AND DLBM_12 IN %s " %(DL_SD_JSYD,DL_ED_WLYD)
r_18_3=GetArea(out_result_X,r_18_3_where)
# print("18-3面积为：" + r_18_3)

#计算1  CZCLX为空，DLBM是其他农用地
r_1_where="CZCLX = '' AND DLBM IN %s " %(DL_SD_QTNYD)
r_1=GetArea(out_result_X,r_1_where)
# print("1面积为：" + r_1)

#计算2  CZCLX为空，DLBM是耕地
r_2_where="CZCLX = '' AND DLBM IN %s " %(DL_SD_GD)
r_2=GetArea(out_result_X,r_2_where)
# print("2面积为：" + r_2)

#计算3  CZCLX为空，DLBM是未利用地
r_3_where="CZCLX = '' AND DLBM IN %s " %(DL_SD_WLYD)
r_3=GetArea(out_result_X,r_3_where)
# print("3面积为：" + r_3)

# 计算6 CZCLX为空，CZCLX_1为空，DLBM_1为耕地
r_6_where="CZCLX = '' AND CZCLX_1 = '' AND DLBM IN %s AND DLBM_1 IN %s " %(DL_SD_QTNYD,DL_SD_GD)
r_6=GetArea(out_result_X,r_6_where)
# print("6面积为：" + r_6)

#计算9-1-2 CZCLX为空，CZCLX_1不为空，DLBM_3为耕地
r_9_1_2_where="CZCLX = '' AND (CZCLX_1 = '' OR DLBM_1 IN %s ) AND DLBM IN %s AND DLBM_1 IN %s AND DLBM_12_13 IN %s " %(DL_SD_JSYD,DL_SD_QTNYD,DL_SD_JSYD,DL_ED_GD)
r_9_1_2=GetArea(out_result_X,r_9_1_2_where)
# print("9-1-2面积为：" + r_9_1_2)
print("退地类算法")
print("合计其他农用地面积：%s" %(r_1+r_18_1+r_09_2_1-r_6-r_9_1_2) ) 
print("耕地面积为：%s" %(r_2+r_6+r_9_1_2+r_18_2+r_09_2_2) ) 
print("未利用地面积：%s" %(r_3+r_18_3+r_09_2_3) ) 
print("建设用地总面积为：" + str(jsydzmj))

print ("总用地面积"+str(r_1+r_18_1+r_09_2_1-r_6-r_9_1_2+   r_2+r_6+r_9_1_2+r_18_2+r_09_2_2+    r_3+r_18_3+r_09_2_3+   jsydzmj))

a_time=time.time()
print ("一共消耗%s秒" % str( a_time-f_time))

print("三调分析")

ed_qtnyd_sql="CZCLX = '' AND DLBM IN %s" %(DL_SD_QTNYD)
ed_qtnyd=GetArea(out_result_X,ed_qtnyd_sql)
print ("其他农用地总面积为：%s" %str(ed_qtnyd))

ed_gd_sql="CZCLX = '' AND DLBM IN %s" %(DL_SD_GD)
ed_gd=GetArea(out_result_X,ed_gd_sql)
print ("耕地总面积为：%s" %str(ed_gd))     

ed_jsyd_sql="(CZCLX <> '' OR DLBM IN %s)" %(DL_SD_JSYD)
ed_jsyd=GetArea(out_result_X,ed_jsyd_sql)
print ("建设用地总面积为：%s" %str(ed_jsyd))

ed_wlyd_sql="CZCLX = '' AND DLBM IN %s" %(DL_SD_WLYD)
ed_wlyd=GetArea(out_result_X,ed_wlyd_sql)
print ("未利用地总面积为：%s" % str(ed_wlyd))   

a_time=time.time()
print ("一共消耗%s秒" % str( a_time-f_time))

以上代码存在一些问题，如DLBM字段在多个图层中都有出现，在进行相交等运算的时候，字段就变得不可控了，对以上代码进行了优化，将一些重复性字段进行了处理，在统计的时候能过够更加明确字段信息；同时增加了叠加报批成果分析的内容。当然这也牺牲了一些运算时间，相较以上代码，对大宗地块运算来讲，大概牺牲了10秒。以下代码的运行环境是在arcgis pro 的python环境下运行的，arcgis10版本的需要将arcpy.management.CalculateField方法的参数进行修改才可，代码如下：

# 待分析地块，先选择下2023的现状图层，然后进行叠加分析
# 选择属于建设用地的部分，然后再与2018和2009的图层进行分析
# 选择数据其他农用地的，
import arcpy
import time

f_time=time.time()
# 地类字典
# 三调地类字典
#耕地，其他农用地，建设用地，未利用地
DL_SD_GD="('0101','0102','0103')"
DL_SD_QTNYD="('0201','0201K','0202','0202K','0203','0203K','0204','0204K','0301','0301K','0302','0302K','0303','0304','0305','0306','0307','0307K','0401','0402','0403','0403K','0404','1006','1103','1104','1104A','1104K','1107','1107A','1202','1203')"
DL_SD_JSYD="('0501','0502','0503','0504','0505','0506','0507','0508','05H1','0601','0602','0603','0604','0701','0702','0801','0802','0803','0804','0805','0806','0807','0808','0809','0810','0810A','08H1','08H2','08H2A','09','0901','0902','0903','0904','0905','0906','1001','1002','1003','1004','1005','1007','1008','1009','1109','1201')"
DL_SD_WLYD="('1101','1102','1105','1106','1108','1110','1204','1205','1206','1207')"
# 二调地类字典
DL_ED_GD="('011','012','013')"
DL_ED_QTNYD="('021','022','023','031','032','033','041','042','104','114','117','122','123')"
DL_ED_JSYD="('051','052','053','054','061','062','063','071','072','081','082','083','084','085','086','087','088','091','092','093','094','095','101','102','103','105','106','107','113','118','121','201','202','203','204','205')"
DL_ED_WLYD="('111','112','115','116','119','043','124','125','126','127')"

# 读取下待分析的地块
arcpy.env.workspace = "C:\\Users\\Administrator\\Desktop\\GHTJ.gdb"
polygon_layer = "Exp"

# 建设用地报批成果层
# BP_CG_layer="CG_JSYDBP"

BP_CG_layer="CG_JSYDBP_TEMP"

xz_2023_layer="TDLY_DLTB_SD_2023"
xz_2022_layer="TDLY_DLTB_SD_2022"
# xz_2023_layer="TDLY_DLTB_SD_2022"
# xz_2022_layer="TDLY_DLTB_SD_2021"
xz_dltb_2018="TDLY_DLTB_2018"
xz_dltb_2009="TDLY_DLTB_2009"

czc_22= "CZCDYD_SD_2022"
czc_23="CZCDYD_SD_2023"

# czc_23="CZCDYD_SD_2022"
# czc_22= "CZCDYD_SD_2021"


out_xz2023="xz2023_temp"
out_result="result_temp"
out_result_23="result_23_temp"
out_result_22="result_22_temp"
out_result_X="result_x_temp"
# 清除掉临时图层
try:
    if arcpy.Exists(out_xz2023):
        arcpy.management.Delete(out_xz2023)
    if arcpy.Exists(out_result):
        arcpy.management.Delete(out_result)
    if arcpy.Exists(out_result_23):
        arcpy.management.Delete(out_result_23)
    if arcpy.Exists(out_result_22):
        arcpy.management.Delete(out_result_22)
    if arcpy.Exists(out_result_X):
        arcpy.management.Delete(out_result_X)
        
    if arcpy.Exists("TmpAnaLayer"):
        arcpy.management.Delete("TmpAnaLayer")
    if arcpy.Exists("YBP_TEMP"):
        arcpy.management.Delete("YBP_TEMP")
    if arcpy.Exists("Intersect_2022"):
        arcpy.management.Delete("Intersect_2022")
    if arcpy.Exists("Intersect_2018"):
        arcpy.management.Delete("Intersect_2018")
    if arcpy.Exists("Intersect_2009"):
        arcpy.management.Delete("Intersect_2009")
        
        
    print("清空临时图层")
except exec as e : 
    print(e)

# 通用方法--------根据筛选条件计算面积
def GetArea(layer_name,where_sql=""):
    sum_area=0
    with arcpy.da.SearchCursor(layer_name,where_clause=where_sql,field_names=['SHAPE@AREA']) as jsydCursor:
        for row in jsydCursor:
            sum_area=sum_area+row[0]
    return sum_area

# 通用方法--------对DLTB图层进行处理
def HandleDLTB(dltb_layer,suffix):
    arcpy.management.AddField(dltb_layer,"DLBM"+suffix,field_type="TEXT",field_length=5)
    arcpy.management.AddField(dltb_layer,"QSXZ"+suffix,field_type="TEXT",field_length=3)
    arcpy.management.CalculateField(dltb_layer,"DLBM"+suffix,"!DLBM!",'PYTHON3')
    arcpy.management.CalculateField(dltb_layer,"QSXZ"+suffix,"!QSXZ!",'PYTHON3')
    # 删除其他字段
    arcpy.management.DeleteField(dltb_layer,["DLBM"+suffix,"QSXZ"+suffix],method="KEEP_FIELDS")
    
# 通用方法--------对CZCDYD图层进行处理
def HandleCZCDYD(dltb_layer,suffix):
    arcpy.management.AddField(dltb_layer,"CZCLX"+suffix,field_type="TEXT",field_length=4)
    arcpy.management.CalculateField(dltb_layer,"CZCLX"+suffix,"!CZCLX!",'PYTHON3')
    # 删除其他字段
    arcpy.management.DeleteField(dltb_layer,["CZCLX"+suffix],method="KEEP_FIELDS")
    

# 第一步，先与报批进行相交，计算出相交的面积
# 第二步，地块扣除相交部分，得到要进行分析的图形
s_bp = arcpy.management.SelectLayerByLocation(BP_CG_layer,overlap_type="INTERSECT",select_features=polygon_layer)
arcpy.analysis.Intersect([polygon_layer,s_bp],out_feature_class="YBP_TEMP",join_attributes="ONLY_FID")
ybp_area=GetArea("YBP_TEMP")
print("与报批相交的面积为："+str(ybp_area))
if (ybp_area==0):
    analysis_layer=polygon_layer
else:
    arcpy.analysis.Erase(polygon_layer,"YBP_TEMP","TmpAnaLayer")
    analysis_layer="TmpAnaLayer"




# 先与2023年的现状数据进行叠加分析,此处选用clip方法，防止待分析地块干扰
S_2023 = arcpy.management.SelectLayerByLocation(xz_2023_layer,overlap_type="INTERSECT",select_features=analysis_layer)
S_2022 = arcpy.management.SelectLayerByLocation(xz_2022_layer,overlap_type="INTERSECT",select_features=analysis_layer)
arcpy.CopyFeatures_management(S_2022,"Intersect_2022")
# 对图层字段进行处理
HandleDLTB("Intersect_2022","2022")

S_2018 = arcpy.management.SelectLayerByLocation(xz_dltb_2018,overlap_type="INTERSECT",select_features=analysis_layer)
arcpy.CopyFeatures_management(S_2018,"Intersect_2018")
HandleDLTB("Intersect_2018","2018")

S_2009=arcpy.management.SelectLayerByLocation(xz_dltb_2009,overlap_type="INTERSECT",select_features=analysis_layer)
arcpy.CopyFeatures_management(S_2009,"Intersect_2009")
HandleDLTB("Intersect_2009","2009")

arcpy.analysis.Clip(S_2023,analysis_layer,out_xz2023)



# 接下来进行叠加分析 全面叠加
arcpy.analysis.Intersect([out_xz2023,"Intersect_2022","Intersect_2018","Intersect_2009"],out_result)

# 图形先clip23年的20x，  图形再clip22年的20x，      然后进行联合分析，输入要素包括（叠加的base面数据，23的clip数据，22的clip数据）
S_CZC_2023=arcpy.management.SelectLayerByLocation(czc_23,overlap_type="INTERSECT",select_features=analysis_layer)
arcpy.analysis.Clip(S_CZC_2023,analysis_layer,out_result_23)

S_CZC_2022=arcpy.management.SelectLayerByLocation(czc_22,overlap_type="INTERSECT",select_features=analysis_layer)
arcpy.analysis.Clip(S_CZC_2022,analysis_layer,out_result_22)

# 字段要重命名一下，以免出现无法确认字段来源的情况，删除无效的字段
HandleCZCDYD(out_result_22,'2022')

# 进行union分析
arcpy.analysis.Union([out_result,out_result_23,out_result_22], out_result_X)



# 计算建设用地面积 CZCLX不为空 且 DLBM2018,DLBM2009都是建设用地的，为建设用地
jsyd_where="(CZCLX <> '' OR DLBM IN %s) AND DLBM2018 IN %s AND DLBM2009 IN %s" %( DL_SD_JSYD, DL_ED_JSYD,DL_ED_JSYD)
jsydzmj=GetArea(out_result_X,jsyd_where) + ybp_area
# print("建设用地总面积为：" + jsydzmj)

#计算09-2-1  CZCLX不为空 且 DLBM2018 是建设用地 且 DLBM2009 是其他农用地
r_09_2_1_where="(CZCLX <> '' OR DLBM IN %s) AND DLBM2018 IN %s AND DLBM2009 IN %s" %(DL_SD_JSYD,DL_ED_JSYD,DL_ED_QTNYD)
r_09_2_1=GetArea(out_result_X,r_09_2_1_where)
# print("09-2-1面积为：" + r_09_2_1)

#计算09-2-2  CZCLX不为空 且 DLBM2018 是建设用地 且 DLBM2009 是耕地
r_09_2_2_where="(CZCLX <> '' OR DLBM IN %s) AND DLBM2018 IN %s AND DLBM2009 IN %s" %(DL_SD_JSYD,DL_ED_JSYD,DL_ED_GD)
r_09_2_2=GetArea(out_result_X,r_09_2_2_where)
# print("09-2-2面积为：" + r_09_2_2)

#计算09-2-3  CZCLX不为空 且 DLBM2018 是建设用地 且 DLBM2009 是未利用地
r_09_2_3_where="(CZCLX <> '' OR DLBM IN %s) AND DLBM2018 IN %s AND DLBM2009 IN %s" %(DL_SD_JSYD,DL_ED_JSYD,DL_ED_WLYD)
r_09_2_3=GetArea(out_result_X,r_09_2_3_where)
# print("09-2-3面积为：" + r_09_2_3)

#计算18-1  CZCLX不为空 且 DLBM2018是其他农用地
r_18_1_where="(CZCLX <> '' OR DLBM IN %s) AND DLBM2018 IN %s " %(DL_SD_JSYD,DL_ED_QTNYD)
r_18_1=GetArea(out_result_X,r_18_1_where)
# print("18-1面积为：" + r_18_1)

#计算18-2  CZCLX不为空 且 DLBM2018是耕地
r_18_2_where="(CZCLX <> '' OR DLBM IN %s) AND DLBM2018 IN %s " %(DL_SD_JSYD,DL_ED_GD)
r_18_2=GetArea(out_result_X,r_18_2_where)
# print("18-2面积为：" + r_18_2)

#计算18-3  CZCLX不为空 且 DLBM2018是未利用地
r_18_3_where="(CZCLX <> '' OR DLBM IN %s) AND DLBM2018 IN %s " %(DL_SD_JSYD,DL_ED_WLYD)
r_18_3=GetArea(out_result_X,r_18_3_where)
# print("18-3面积为：" + r_18_3)

#计算1  CZCLX为空，DLBM是其他农用地
r_1_where="CZCLX = '' AND DLBM IN %s " %(DL_SD_QTNYD)
r_1=GetArea(out_result_X,r_1_where)
# print("1面积为：" + r_1)

#计算2  CZCLX为空，DLBM是耕地
r_2_where="CZCLX = '' AND DLBM IN %s " %(DL_SD_GD)
r_2=GetArea(out_result_X,r_2_where)
# print("2面积为：" + r_2)

#计算3  CZCLX为空，DLBM是未利用地
r_3_where="CZCLX = '' AND DLBM IN %s " %(DL_SD_WLYD)
r_3=GetArea(out_result_X,r_3_where)
# print("3面积为：" + r_3)

# 计算6 CZCLX为空，CZCLX_1为空，DLBM_1为耕地
r_6_where="CZCLX = '' AND CZCLX2022 = '' AND DLBM IN %s AND DLBM2022 IN %s " %(DL_SD_QTNYD,DL_SD_GD)
r_6=GetArea(out_result_X,r_6_where)
# print("6面积为：" + r_6)

#计算9-1-2 CZCLX为空，CZCLX_1不为空，DLBM_3为耕地
r_9_1_2_where="CZCLX = '' AND (CZCLX2022 = '' OR DLBM2022 IN %s ) AND DLBM IN %s AND DLBM2022 IN %s AND DLBM2009 IN %s " %(DL_SD_JSYD,DL_SD_QTNYD,DL_SD_JSYD,DL_ED_GD)
r_9_1_2=GetArea(out_result_X,r_9_1_2_where)
# print("9-1-2面积为：" + r_9_1_2)
print("退地类算法")
print("耕地面积为：%s" %(r_2+r_6+r_9_1_2+r_18_2+r_09_2_2) ) 
print("合计其他农用地面积：%s" %(r_1+r_18_1+r_09_2_1-r_6-r_9_1_2) ) 

print("建设用地总面积为：" + str(jsydzmj))
print("未利用地面积：%s" %(r_3+r_18_3+r_09_2_3) ) 
print ("总用地面积"+str(r_1+r_18_1+r_09_2_1-r_6-r_9_1_2+   r_2+r_6+r_9_1_2+r_18_2+r_09_2_2+    r_3+r_18_3+r_09_2_3+   jsydzmj ))
a_time=time.time()
print ("一共消耗%s秒" % str( a_time-f_time))


print("三调分析")

ed_gd_sql="CZCLX = '' AND DLBM IN %s" %(DL_SD_GD)
ed_gd=GetArea(out_result_X,ed_gd_sql)
print ("耕地总面积为：%s" %str(ed_gd))    

ed_qtnyd_sql="CZCLX = '' AND DLBM IN %s" %(DL_SD_QTNYD)
ed_qtnyd=GetArea(out_result_X,ed_qtnyd_sql)
print ("其他农用地总面积为：%s" %str(ed_qtnyd))

ed_jsyd_sql="(CZCLX <> '' OR DLBM IN %s)" %(DL_SD_JSYD)
ed_jsyd=GetArea(out_result_X,ed_jsyd_sql)
print ("建设用地总面积为：%s" %str(ed_jsyd+ybp_area))

ed_wlyd_sql="CZCLX = '' AND DLBM IN %s" %(DL_SD_WLYD)
ed_wlyd=GetArea(out_result_X,ed_wlyd_sql)
print ("未利用地总面积为：%s" % str(ed_wlyd))   

a_time=time.time()
print ("一共消耗%s秒" % str( a_time-f_time))
学习后的第一个arcpy实用工具出炉--三调统计口径下的地类图斑面积统计
nxy_918的专栏
04-19 1236
真正实现三调统计口径下不同行政区（县市区、乡镇）的地类图斑面积统计。
有源淹没分析arcgis_基于DEM的有源洪水淹没分析算法
weixin_42184548的博客
12-23 3771
成静摘 要：近年以来，基于数字高程模型(Digital Elevation Model，DEM)，将水动力模型与地理信息系统技术相结合，来计算、显示洪水淹没范围是水利应用领域的研究热点之一。洪水淹没模拟分为有源淹没和无源淹没两种情形，针对有源淹没递归算法计算效率低和稳定性不高的缺点，本文基于“栈”数据结构的特点，提出了一种基于淹没区缓冲堆栈的有源洪水淹没分析种子蔓延算法，实现了更为高效的有源洪水淹...
ARCPY使用手册ARCPY使用手册
08-11
ARCPY使用手册
建设用地审查报批地类认定工具-地类追溯（退地类工具）
xyt1172228177的博客
03-27 4561
建设用地审查报批地类追溯
【ArcGIS】确定生态退耕的地块
dession_Wu的博客
12-27 464
分析饮用水源保护区范围内耕地分布情况（总面积、斑块数、平均斑块大小、水田旱地比例、坡度情况）。确定需要进行生态退耕的地块：（1）一级、二级保护区范围内，某一个完整的耕地地块如有一半以上面积坡度大于25度，则认为该地块需要进行生态退耕；（2）一级、二级保护区范围内，距离村庄超过300米的耕地地块，需进行生态退耕。
第三次国土调查数据库字母缩写代表含义
地信网（地理信息网|测绘小站）
04-25 4717
CCWJQ拆除未尽区，CJDCQ村界调查区，CJDCQJX村界调查区界线，CLKZD测量控制点，CSKFBJ城市开发边界，CZCDYD城镇村等用地，DGX等高线，DLTB地类图斑，DZGY地质公园，FJMSQ风景名胜区，GCZJD高程注记点，GDDB耕地等别，GFBQ光伏板区，GJGY国家公园，JZKZD数字正射影像图纠正控制点，KFYQ开发园区，LMFW路面范围，LSYD临时用地，PZWJSTD批准为建设用地，QTJZKFQ其他禁止开发区，SCZZBHQ水产种植资源保护区，SDGY湿地公园
国土调查数据库字母缩写各自代表的意思？
weixin_48478010的博客
07-03 1万+
CCWJQ拆除未尽区，CJDCQ村界调查区，CJDCQJX村界调查区界线......
国土调查数据库标准 缩写含义
xzy123138的博客
07-22 2413
JZKZD数字正射影像图纠正控制点。SCZZBHQ水产种植资源保护区。CJDCQJX村界调查区界线。PZWJSTD批准为建设用地。QTJZKFQ其他禁止开发区。YJJBNTTB永久基本农田。ZRYCBHQ自然遗产保护区。ZYXMYD重要项目用地。#国土调查数据库字母缩写。CSKFBJ城市开发边界。CZCDYD城镇村等用地。STBHHX生态保护红线。YYSSYD饮用水水源地。CCWJQ拆除未尽区。CJDCQ村界调查区。CLKZD测量控制点。FJMSQ风景名胜区。GCZJD高程注记点。WJMHD无居民海岛。
MySQL存储引擎总结
An342647823的专栏
04-09 1053
转载：http://www.jb51.net/article/55849.htm 前言 在数据库中存的就是一张张有着千丝万缕关系的表，所以表设计的好坏，将直接影响着整个数据库。而在设计表的时候，我们都会关注一个问题，使用什么存储引擎。等一下，存储引擎？什么是存储引擎？ 什么是存储引擎？ 关系数据库表是用于存储和组织信息的数据结构，可以将表理解为由行和列组成的表格，类似于
ArcGIS学习：用Python实现地图四色填充
05-13
##### 3.3 使用Python编写工具计算每个省份的颜色 在ArcGIS中使用Python编写脚本实现四色填充的过程如下： 1. **导入必要的模块**：首先需要导入ArcPy模块，以便能够调用ArcGIS的功能。 2. **定义参数**：定义需要...
VIC 模型_参数准备-气候驱动
UCAS_Leon的博客
03-16 423
基流也叫底水,是河道中常年存在的那部分径流,基流系数是基流占河流径流的比重。 流量历时曲线(flow-duration curve)，是反映流量在某一时段内(或年内某一季节、一年)超过某一数值持续天数的一种统计特性曲线。 参数获取 植被参数： 中科院数据分类不能满足要求，因此需要重分类： 根据冬天的 NDVI 数据，以 1：10 万林地做掩膜，首先区分常绿和落叶；根据 NDVI 和其它气 候指标...
ARCGIS自定义工具箱：打造个性化土地利用分析功能
首先，概述了土地利用分析的理论基础，包括分类标准、常用分析方法，以及ARCGIS在此领域的平台优势和核心功能。接着，详细阐述了构建ARCGIS自定义工具箱的步骤，包括创建、管理、模型构建、参数设置和性能优化等方面...
三调图斑、地类图斑字段
weixin_45864817的博客
03-22 2万+
标题三调包含哪些图斑？每个图斑代表什么含义？ 图斑 含义 CCWJQ 拆除未尽区 CJDCQ 村界调查区 CJDCQJX 村界调查区界线 CLKZD 测量控制点 CSKFBJ 城市开发边界 CZCDYD 城镇村等用地 DGX 等高线 DLTB 地类图斑 DZGY 地质公园 FJMSQ 风景名胜区 GCZJD 高程注记点 GDDB 耕地等别 GFBQ 光伏板区 GJGY 国家公园 JZKZD 数字正射影像图纠正控制点 KFYQ 开发园区 LMFW 路面范围 LSYD 临时用地 PZWJSTD 批准为建设用地
三调数据库标准和输出表格
流浪的地信人
05-10 2万+
TD 中华人民共和国国土资源行业标准 TD/T 1016—2017 土地利用数据库标准 Standard For Land Use Database 20XX-XX-XX 发布 20XX-XX-XX 实施 中华人民共和国国土资源部发布 f 前 言 II 引 言 III 1范围 1 2规范性引用文件 1 3术语和定义 1 4数据库内容和要素分类编码 ,2 5数据库结构定义 5 6数据交换文件命名规则 25 7数据交换内容与格式 ，.霆* . … . 28 8 元数据 28 -* t —1— 刖 s 本标准按照
【arcmap】SelectLayerByLocation_management和SelectLayerByAttribute_management
weixin_43958438的博客
09-07 1067
arcpy.SelectLayerByLocation_management 和 arcpy.SelectLayerByAttribute_management 是 ArcGIS 中两个常用的选择要素的方法。
三调数据库及DLTB各个字段含义
热门推荐
qq_45697428的博客
03-31 8万+
三调包含哪些图斑？每个图斑代表什么含义？ 图斑 含义 CCWJQ 拆除未尽区 CJDCQ 村界调查区 CJDCQJX 村界调查区界线 CLKZD 测量控制点 CSKFBJ 城市开发边界 CZCDYD 城镇村等用地 DGX 等高线 DLTB 地类图斑 DZGY 地质公园 FJMSQ 风景名胜区 GCZJD 高程注记点 GDDB 耕地等别 GFBQ 光伏板区 GJGY 国家公园 JZKZD 数字正射影像图纠正控制点 KFYQ 开发园区
使用ArcPy进行图斑内洞检查
a_dev的专栏
05-19 1190
文章目录需求思路实现 需求 有一面状地类图斑要素类，需要检测出所有位于要素面范围内的孔洞，根据其面积大小进行不同的处理（比如，大于50亩需要人工复核，小于50亩报错）。 思路 第一种思路：逐个遍历，检测面中是否包含内洞。具体方式是将其转为GeometryCollection对象，判断是否具有面积为负值的多边形，再判断面积。此种方式执行效率较低，速度慢，尤其是图斑数量多而孔洞少的时候。 第二种思路...
Python学习笔记（三）（程序流程控制）
最新发布
wsys_yysn的博客
04-24 422
条件语句：通过实现分支逻辑。循环语句for遍历序列，while满足条件时循环。流程控制break终止循环，continue跳过当前迭代，pass占位。异常处理保证程序健壮性。掌握这些流程控制工具，可以编写更灵活、高效的 Python 程序！
2025-04-22 李沐深度学习5 —— 线性回归
zheliku的博客
04-22 1229
​ 在美国买房时，买家需根据房屋信息（如卧室数量、卫生间数量、面积等）预测合理的成交价。输入数据：房屋特征（如X1=卧室数X2=卫生间数X3=面积输出目标：预测成交价Y。实际挑战：卖家的标价和网站估价（如Redfin）仅为参考，最终需通过竞价决定成交价，因此准确预测至关重要。​ 计算输入特征X\mathbf{X}X和模型权重w\mathbf{w}w的矩阵-向量乘法后加上偏置bbb。​ 注意，上面的XwXw是一个向量，而bbb是一个标量。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值