percelenErfpacht.py

import arcpy
import sys
sys.path.append('.')

import pandas as pd
from xlsxwriter.workbook import Workbook
import matplotlib.colors as mcolors



arcpy.env.overwriteOutput = True
arcpy.env.workspace = r'C:\Users\vince\Documents\ArcGIS\Projects\percelen erfpacht januari 2023\percelen erfpacht januari 2023.gdb'

excel = r'C:\Users\vince\Documents\ArcGIS\Projects\percelen erfpacht januari 2023\output\uitvoer_percelen_januari2023.xlsx'

meetpunten = "meetpunten_update_24_jan"
meetpunt_veld = "Werknummer"
percelen = "meetpunten_percelen_allen_update_24_jan"
velden = ["KADTOTAAL","ZR_OMSCHRIJVING","NAAM","VOORN","VOORL","VOORV","TITEL","STRAAT","HUISNR","POSTCODE","WOONPLAATS","aantal_go","go_nummers"]
velden_eigenaren = ["ZR_OMSCHRIJVING","NAAM","VOORN","VOORL","VOORV","TITEL","STRAAT","HUISNR","POSTCODE","WOONPLAATS"]
array = arcpy.da.FeatureClassToNumPyArray(percelen, velden)
df = pd.DataFrame(array)


def add_measurements():
    # meetpunten toevoegen aan dataset
    join_cursor = arcpy.da.UpdateCursor(percelen,velden)
    for jRow in join_cursor:
        meetpunt_nummers = []
        perceel_nummer = jRow[0]
        where = '"' + "KADTOTAAL" + '" = ' + "'"+perceel_nummer+"'"
        temp_perceel= "temp_perceel"
        arcpy.Select_analysis(percelen, temp_perceel, where)
        arcpy.MakeFeatureLayer_management(meetpunten, "temp_meetpunten") 
        arcpy.SelectLayerByLocation_management("temp_meetpunten", "INTERSECT", "temp_perceel", "", "NEW_SELECTION", "NOT_INVERT")
        arcpy.CopyFeatures_management("temp_meetpunten", "temp_meetpunten_output")
        tempCursor = arcpy.da.SearchCursor("temp_meetpunten_output",meetpunt_veld)
        for tRow in tempCursor:
            meetpunt_nummers.append("{}, ".format(tRow[0]))
        
        del tempCursor
        alle_meetpunten = ' '.join([str(elem) for elem in meetpunt_nummers])
        alle_meetpunten = alle_meetpunten[:len(alle_meetpunten)-2]
        jRow[12] = alle_meetpunten
        join_cursor.updateRow(jRow)
        print (perceel_nummer)


def writer():
    array = arcpy.da.FeatureClassToNumPyArray(percelen, velden)
    df = pd.DataFrame(array)
    # sorteren
    df["volgorde"] = ""
    df["Soorten pachters"] = ""

    for index, row in df.iterrows():
        if row['ZR_OMSCHRIJVING'] == "Eigendom (recht van)":
            # row["volgorde"] = 0
            df.at[index, 'volgorde'] = 0

        elif "Erfpacht" in row['ZR_OMSCHRIJVING']:
            #   row["volgorde"] = 1
            df.at[index, 'volgorde'] = 1


        else:
            # row["volgorde"] = 2
            df.at[index, 'volgorde'] = 2

        
    df = df.sort_values(by=['volgorde'])



    # df_selected = df[(df["ZR_OMSCHRIJVING"] == "Eigendom (recht van)") | (df["ZR_OMSCHRIJVING"] == "Erfpacht (recht van)") |(df["ZR_OMSCHRIJVING"] == "Erfpacht en Opstal (recht van)")]
    df_selected = df




    # opbouw xlsx
    workbook = Workbook(excel)
    worksheet = workbook.add_worksheet()

    # stijl toevoegen voor headers
    bold = workbook.add_format({'bold': True})
    bold.set_bg_color("#25A4C6")
    colors = list(mcolors.cnames.values())



    # schrijf default kolomnamen
    worksheet.write(0, 0, "Perceelnummer", bold)
    worksheet.write(0, 1, "Aantal eigenaren/pachters", bold)
    worksheet.write(0, 2, "Soorten pachters", bold)
    worksheet.write(0, 3, "Aantal GO-locaties", bold)
    worksheet.write(0, 4, "GO-nummers", bold)




    # sorteer op perceelnummer (KADTOTAAL)
    groep_percelen = df_selected.groupby('KADTOTAAL')

    rijnummer = 1

    for name, groep in groep_percelen:


        # erfpachters = False

        

        # check of erpachters aanwezig, anders skippen
        # for index, row in groep.iterrows():
        #     if "Erfpacht" in row["ZR_OMSCHRIJVING"]:
        #         erfpachters = True

        # check of pachters aanwezig, anders skippen

        pachters = True
        pachter_values = ["Erfpacht","Gebruik","Opstal","Zakelijk"]
        pachter_list = []


        for index, row in groep.iterrows():
            for item in pachter_values:
                if item in row["ZR_OMSCHRIJVING"]:
                    if item in pachter_list:
                        pass
                    else:
                        pachter_list.append(item)
                    
                    pachters = True

        
        pachter_list_string = [str(i) for i in pachter_list]
        pachter_list_string = ','.join(pachter_list_string)

        
        # if erfpachters == True:
        if pachters == True:
            
            kolomnummer = 5

            # per perceel een nieuwe rij aanmaken
            worksheet.write(rijnummer, 0, name)
            

            # per rij aantal erfpachters weergeven in kolom
            groupmembers = 0
            for index, row in groep.iterrows():

            
                groupmembers += 1
                
                # write extra header
                for veld in velden_eigenaren:
                    

                    if row[veld] == "None":
                        fill_value = ''
                    else:
                        fill_value = row[veld]
                    

                
                    worksheet.write(0, kolomnummer, "{} {}".format(veld, groupmembers), bold)
                    worksheet.write(rijnummer, kolomnummer, fill_value)

                    kolomnummer += 1

            worksheet.write(rijnummer, 1, groupmembers)
            worksheet.write(rijnummer, 2, pachter_list_string)
            worksheet.write(rijnummer, 3, row['aantal_go'])
            worksheet.write(rijnummer, 4, row['go_nummers'])
        


            rijnummer += 1
        else:
            pass

            

    workbook.close()







# add_measurements()
writer()