利用Arcpy从shape files采点,批量

问题描述

M个study area, 每一个study area有多种类别,每一种类别有N个polygons,要求从每个study area的每一类中采样得到N个点。

解决办法

  • 采样可以用arcMap
  • 批量处理可以用python遍历文件夹

相关

  • kml to shapefile: ogrmerge.py
  • setup python and ArcGIS
  • 可以直接在arcmap的python窗口里面load代码然后运行
  • 如果用pycharm运行python,需要设置project interpreter

Sample code

1 convert kml files to shapefiles

#!/bin/bash
set -e

#the dir of the kml files
Dir='/data/qiu/LCZ_wudpat_train/raw/*/'
files=$(find $Dir -name '*.km*')

for f in $files; do
    echo $f

    city0=${f%/*}
    #cut f: delete the first (starting from right) / and everything rightwards.
    # echo $city0
    city=${city0##*/}
  #cut city0: delete the first (starting from left) / and everything leftwards.
    echo $city

    rm -r -f  $city0/$city'_shapeFile'
    rm -r -f  /data/qiu/LCZ_wudpat_train/shapeFile/$city

    ogrmerge.py -o /data/qiu/LCZ_wudpat_train/shapeFile/$city $f -nln $city'_'{LAYER_NAME}
    #the output can be one file or multiple files for multiple subdir in the kml file

done;

# Ref:
# https://gdal.org/programs/ogrmerge.html

2 sampling with arcpy

import os
import arcpy
import numpy as np
import shutil
import glob

#top_folder = r"F:\0Data4Paper\LCZ_wudpat_train\test"
#'folder that contains multiple folders for study areas'
top_folder = r"F:\0Data4Paper\LCZ_wudpat_train\shapeFile"

#'tmp dir for intermediate results'
intermFolder = r"F:\0Data4Paper\LCZ_wudpat_train\tmp"

for path, dirs, files in os.walk(top_folder):
    for d in dirs:

        #print(path, dirs, files)
        sh_path = os.path.join(top_folder, d)
        #print(sh_path)

        #'''clear existing files'''
        saveDisFolder=os.path.join(r'F:\0Data4Paper\LCZ_wudpat_train\ds', d)
        print(saveDisFolder)
        if os.path.exists(saveDisFolder):
            shutil.rmtree(saveDisFolder)
        if not os.path.exists(saveDisFolder):
            os.mkdir(saveDisFolder)

        pointsFolder = os.path.join(r"F:\0Data4Paper\LCZ_wudpat_train\points", d)
        print(pointsFolder)
        if os.path.exists(pointsFolder):
            shutil.rmtree(pointsFolder)
        if not os.path.exists(pointsFolder):
            os.mkdir(pointsFolder)

        if os.path.exists(intermFolder):
            shutil.rmtree(intermFolder)
        if not os.path.exists(intermFolder):
            os.mkdir(intermFolder)

        #for each class:
        for i in np.arange(1,18):

            c=i
            if i>10:
                i=i+90

            file=os.path.join(sh_path, d+'_'+str(i) + "." + 'shp')
            print(file)
            if os.path.isfile(file):

                #dissolve polygons into one and save temporally
                saveDis=os.path.join(saveDisFolder, d+'_'+str(i))
                arcpy.Dissolve_management(file, saveDis, "", "", "", "")

                #sample points from the dissolved polygons
                arcpy.CreateRandomPoints_management(intermFolder, d+'_'+str(i), saveDis+ "." + 'shp', "", 30, "100 Meter", "MULTIPOINT")

                savePoints=os.path.join(intermFolder, d+'_'+str(i))

                #add a field to assign classes
                arcpy.AddField_management(savePoints+ "." + 'shp', "class", "SHORT")
                arcpy.CalculateField_management(savePoints+ "." + 'shp', "class", int(c))
            else:
                print('maybe this class not exist!')


        fileList = glob.glob(intermFolder+ "\*.shp")
        # #print(fileList)
        #'merge points into one file for each study area'
        arcpy.Merge_management(fileList, os.path.join(pointsFolder, d), "")
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 221,331评论 6 515
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 94,372评论 3 398
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 167,755评论 0 360
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 59,528评论 1 296
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 68,526评论 6 397
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 52,166评论 1 308
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,768评论 3 421
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,664评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 46,205评论 1 319
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 38,290评论 3 340
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 40,435评论 1 352
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 36,126评论 5 349
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,804评论 3 333
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 32,276评论 0 23
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,393评论 1 272
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,818评论 3 376
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 45,442评论 2 359