notebooks/Add Offset Metadata Fields.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Add Offset Metadata fields\n",
    "\n",
    "This notebook is designed for use in ArcGIS Pro 2.9+.  See [this topic](https://pro.arcgis.com/en/pro-app/2.8/arcpy/get-started/pro-notebooks.htm) For information on how to use Jupyter notebooks in ArcGIS Pro.\n",
    "\n",
    "- Prior to running this script you should ensure that the [Add GPS Metadata Fields](https://pro.arcgis.com/en/pro-app/latest/tool-reference/data-management/add-gps-metadata-fields.htm) Geoprocessing tool has been run on the `input_feature_class`. \n",
    "- Use this script to update an existing point feature class to support offset workflows\n",
    "- This adds a series of fields and necessary coded value domains\n",
    "\n",
    "### Steps\n",
    "1. Update the `input_feaature_class` variable\n",
    "2. Run all cells"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 70,
   "metadata": {},
   "outputs": [],
   "source": [
    "# update to location of line or polygon feature class\n",
    "input_feature_class = r\"C:\\Users\\doug-m\\Documents\\ArcGIS\\Projects\\MyProject\\MyGDB.gdb\\water_valves\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 71,
   "metadata": {},
   "outputs": [],
   "source": [
    "import arcpy\n",
    "import os\n",
    "\n",
    "def get_geodatabase_path(input_feature_class):\n",
    "    \"\"\"\n",
    "    Gets the parent geodatabase of the layer\n",
    "    :param input_feature_class: (string) The feature layer to get the parent database of\n",
    "    :return: (string) The path to the geodatabase\n",
    "    \"\"\"\n",
    "    workspace = os.path.dirname(arcpy.Describe(input_feature_class).catalogPath)\n",
    "    if [any(ext) for ext in ('.gdb', '.mdb', '.sde') if ext in os.path.splitext(workspace)]:\n",
    "        return workspace\n",
    "    else:\n",
    "        return os.path.dirname(workspace)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 72,
   "metadata": {},
   "outputs": [],
   "source": [
    "def check_and_create_domains(geodatabase):\n",
    "    \"\"\"\n",
    "    Checks if the Fix Type domain already exists.\n",
    "    If the domains do not exist, they are created\n",
    "    :param geodatabase: (string) the path to the geodatabase to check\n",
    "    :return:\n",
    "    \"\"\"\n",
    "    domains = arcpy.da.ListDomains(geodatabase)\n",
    "    domain_names = [domain.name.lower() for domain in domains]\n",
    "    if 'esri_fix_type_domain' in domain_names:\n",
    "        for domain in domains:\n",
    "            if domain.name.lower() == 'esri_fix_type_domain':\n",
    "                # check if cvs 0,1,2,4,5 are in the codedValues\n",
    "                values = [cv for cv in domain.codedValues]\n",
    "                if not {0, 1, 2, 4, 5}.issubset(values):\n",
    "                    return \"esri_fix_type_domain is missing a coded value pair.\"\n",
    "    else:\n",
    "        # Add the domain and values\n",
    "      \n",
    "        arcpy.CreateDomain_management(in_workspace=geodatabase,\n",
    "                                      domain_name=\"esri_fix_type_domain\",\n",
    "                                      domain_description=\"Fix Type\",\n",
    "                                      field_type=\"SHORT\",\n",
    "                                      domain_type=\"CODED\",\n",
    "                                      split_policy=\"DEFAULT\",\n",
    "                                      merge_policy=\"DEFAULT\")\n",
    "\n",
    "        arcpy.AddCodedValueToDomain_management(in_workspace=geodatabase,\n",
    "                                               domain_name=\"esri_fix_type_domain\",\n",
    "                                               code=\"0\",\n",
    "                                               code_description=\"Fix not valid\")\n",
    "        arcpy.AddCodedValueToDomain_management(in_workspace=geodatabase,\n",
    "                                               domain_name=\"esri_fix_type_domain\",\n",
    "                                               code=\"1\",\n",
    "                                               code_description=\"GPS\")\n",
    "        arcpy.AddCodedValueToDomain_management(in_workspace=geodatabase,\n",
    "                                               domain_name=\"esri_fix_type_domain\",\n",
    "                                               code=\"2\",\n",
    "                                               code_description=\"Differential GPS\")\n",
    "        arcpy.AddCodedValueToDomain_management(in_workspace=geodatabase,\n",
    "                                               domain_name=\"esri_fix_type_domain\",\n",
    "                                               code=\"4\",\n",
    "                                               code_description=\"RTK Fixed\")\n",
    "        arcpy.AddCodedValueToDomain_management(in_workspace=geodatabase,\n",
    "                                               domain_name=\"esri_fix_type_domain\",\n",
    "                                               code=\"5\",\n",
    "                                               code_description=\"RTK Float\")\n",
    "        \n",
    "    if 'esri_offset_method_domain' in domain_names:\n",
    "        for domain in domains:\n",
    "            if domain.name.lower() == 'esri_offset_method_domain':\n",
    "                # check if offset methods are in the codedValues\n",
    "                values = [cv for cv in domain.codedValues]\n",
    "                if not {'azimuth-azimuth', \n",
    "                        'distance-azimuth', \n",
    "                        'distance-distance', \n",
    "                        'distance-angle', \n",
    "                        'azimuth-azimuth-backsight', \n",
    "                        'distance-azimuth-backsight',\n",
    "                        'distance-distance-backsight',\n",
    "                        'user-defined'}.issubset(values):\n",
    "                    return \"esri_offset_method_domain is missing a coded value pair.\"\n",
    "    else:\n",
    "        # Add the domain and values\n",
    "      \n",
    "        arcpy.CreateDomain_management(in_workspace=geodatabase,\n",
    "                                      domain_name=\"esri_offset_method_domain\",\n",
    "                                      domain_description=\"Offset Method Domain\",\n",
    "                                      field_type=\"TEXT\",\n",
    "                                      domain_type=\"CODED\",\n",
    "                                      split_policy=\"DEFAULT\",\n",
    "                                      merge_policy=\"DEFAULT\")\n",
    "\n",
    "        arcpy.AddCodedValueToDomain_management(in_workspace=geodatabase,\n",
    "                                               domain_name=\"esri_offset_method_domain\",\n",
    "                                               code=\"azimuth-azimuth\",\n",
    "                                               code_description=\"Azimuth-Azimuth\")\n",
    "        arcpy.AddCodedValueToDomain_management(in_workspace=geodatabase,\n",
    "                                               domain_name=\"esri_offset_method_domain\",\n",
    "                                               code=\"distance-angle\",\n",
    "                                               code_description=\"Distance-Angle\")\n",
    "        arcpy.AddCodedValueToDomain_management(in_workspace=geodatabase,\n",
    "                                               domain_name=\"esri_offset_method_domain\",\n",
    "                                               code=\"distance-azimuth\",\n",
    "                                               code_description=\"Distance-Azimuth\")\n",
    "        arcpy.AddCodedValueToDomain_management(in_workspace=geodatabase,\n",
    "                                               domain_name=\"esri_offset_method_domain\",\n",
    "                                               code=\"distance-distance\",\n",
    "                                               code_description=\"Distance-Distance\")\n",
    "        arcpy.AddCodedValueToDomain_management(in_workspace=geodatabase,\n",
    "                                               domain_name=\"esri_offset_method_domain\",\n",
    "                                               code=\"azimuth-azimuth-backsight\",\n",
    "                                               code_description=\"Azimuth-Azimuth-Backsight\")\n",
    "        arcpy.AddCodedValueToDomain_management(in_workspace=geodatabase,\n",
    "                                               domain_name=\"esri_offset_method_domain\",\n",
    "                                               code=\"distance-azimuth-backsight\",\n",
    "                                               code_description=\"Distance-Azimuth-Backsight\")\n",
    "        arcpy.AddCodedValueToDomain_management(in_workspace=geodatabase,\n",
    "                                               domain_name=\"esri_offset_method_domain\",\n",
    "                                               code=\"distance-distance-backsight\",\n",
    "                                               code_description=\"Distance-Distance-Backsight\")\n",
    "        arcpy.AddCodedValueToDomain_management(in_workspace=geodatabase,\n",
    "                                               domain_name=\"esri_offset_method_domain\",\n",
    "                                               code=\"user-defined\",\n",
    "                                               code_description=\"User-defined\")        "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "def add_offset_fields(feature_class):\n",
    "    \"\"\"\n",
    "    This adds specific fields required for laser offset workflows\n",
    "        This will report errors if:\n",
    "            1) Any of the fields already exist\n",
    "            2) The input layer is not a point layer\n",
    "            3) The layer is not found\n",
    "            4) The layer is a shapefile\n",
    "    \"\"\"\n",
    "\n",
    "    if not arcpy.Exists(feature_class):\n",
    "        return \"Feature layer: {} not found!\".format(feature_class)\n",
    "    if arcpy.Describe(feature_class).shapeType != \"Point\":\n",
    "        return \"Feature layer: {} is not a point layer.\".format(feature_class)\n",
    "    if arcpy.Describe(feature_class).dataType == \"ShapeFile\":\n",
    "        return \"Shapefiles are not supported.\"\n",
    "\n",
    "    # Check the domains to see if they exist and are valid\n",
    "    # will update if necessary\n",
    "    geodatabase = get_geodatabase_path(feature_class)\n",
    "    check_and_create_domains(geodatabase)    \n",
    "\n",
    "    # Add GNSS metadata fields\n",
    "    existing_fields = [field.name.lower() for field in arcpy.ListFields(feature_class)]\n",
    "\n",
    "    if 'esrisnsr_offset_method' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrisnsr_offset_method',\n",
    "                                  field_type=\"TEXT\",\n",
    "                                  field_alias='Offset method',\n",
    "                                  field_is_nullable=\"NULLABLE\",\n",
    "                                  field_domain=\"esri_offset_method_domain\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrisnsr_offset_device' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrisnsr_offset_device',\n",
    "                                  field_type=\"TEXT\",\n",
    "                                  field_alias='Rangefinder name',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrignss_antenna_height' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_antenna_height',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='GNSS Antenna height (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrisnsr_offset_laser_height' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrisnsr_offset_laser_height',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='Rangefinder height (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrisnsr_offset_laser_target_ht' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrisnsr_offset_laser_target_ht',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='Laser target height (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrisnsr_offset_mag_var' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrisnsr_offset_mag_var',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='Rangefinder Magnetic variation (°)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "        \n",
    "    if 'esrignss_cp1_latitude' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_latitude',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP1 Latitude',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrignss_cp1_longitude' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_longitude',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP1 Longitude',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )       \n",
    "\n",
    "    if 'esrignss_cp1_altitude' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_altitude',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP1 Ellipsoidal height (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrignss_cp1_h_rms' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_h_rms',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP1 Horizontal accuracy (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )     \n",
    "\n",
    "    if 'esrignss_cp1_v_rms' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_v_rms',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP1 Vertical accuracy (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrignss_cp1_fixdatetime' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_fixdatetime',\n",
    "                                  field_type=\"DATE\",\n",
    "                                  field_alias='CP1 Fix time',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )     \n",
    "\n",
    "    if 'esrignss_cp1_fixtype' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_fixtype',\n",
    "                                  field_type=\"SHORT\",\n",
    "                                  field_alias='CP1 Fix type',\n",
    "                                  field_is_nullable=\"NULLABLE\",\n",
    "                                  field_domain=\"esri_fix_type_domain\"\n",
    "                                  )     \n",
    "\n",
    "    if 'esrignss_cp1_corr_age' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_corr_age',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP1 Correction age',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrignss_cp1_stationid' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_stationid',\n",
    "                                  field_type=\"SHORT\",\n",
    "                                  field_alias='CP1 Station ID',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )    \n",
    "\n",
    "    if 'esrignss_cp1_numsats' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_numsats',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP1 Number of satellites',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrignss_cp1_pdop' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_pdop',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP1 PDOP',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )    \n",
    "\n",
    "    if 'esrignss_cp1_hdop' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_hdop',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP1 HDOP',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )    \n",
    "        \n",
    "    if 'esrignss_cp1_vdop' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_vdop',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP1 VDOP',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )        \n",
    "\n",
    "    if 'esrignss_cp1_avg_h_rms' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_avg_h_rms',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP1 Average horizontal accuracy (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )    \n",
    "\n",
    "    if 'esrignss_cp1_avg_v_rms' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_avg_v_rms',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP1 Average vertical accuracy (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )    \n",
    "\n",
    "    if 'esrignss_cp1_avg_positions' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_avg_positions',\n",
    "                                  field_type=\"SHORT\",\n",
    "                                  field_alias='CP1 Averaged positions',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrignss_cp1_h_stddev' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp1_h_stddev',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP1 Standard deviation (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )    \n",
    "\n",
    "    if 'esrisnsr_cp1_slope_dist' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrisnsr_cp1_slope_dist',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP1 Laser slope distance (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )    \n",
    "        \n",
    "    if 'esrisnsr_cp1_inclination' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrisnsr_cp1_inclination',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP1 Laser inclination (°)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )    \n",
    "        \n",
    "    if 'esrisnsr_cp1_azimuth' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrisnsr_cp1_azimuth',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP1 Laser azimuth (°)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )        \n",
    "\n",
    "    if 'esrignss_cp2_latitude' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_latitude',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP2 Latitude',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrignss_cp2_longitude' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_longitude',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP2 Longitude',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )       \n",
    "\n",
    "    if 'esrignss_cp2_altitude' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_altitude',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP2 Ellipsoidal height (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrignss_cp2_h_rms' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_h_rms',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP2 Horizontal accuracy (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )     \n",
    "\n",
    "    if 'esrignss_cp2_v_rms' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_v_rms',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP2 Vertical accuracy (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrignss_cp2_fixdatetime' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_fixdatetime',\n",
    "                                  field_type=\"DATE\",\n",
    "                                  field_alias='CP2 Fix time',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )     \n",
    "\n",
    "    if 'esrignss_cp2_fixtype' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_fixtype',\n",
    "                                  field_type=\"SHORT\",\n",
    "                                  field_alias='CP2 Fix type',\n",
    "                                  field_is_nullable=\"NULLABLE\",\n",
    "                                  field_domain=\"esri_fix_type_domain\"\n",
    "                                  )     \n",
    "\n",
    "    if 'esrignss_cp2_corr_age' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_corr_age',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP2 Correction age',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrignss_cp2_stationid' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_stationid',\n",
    "                                  field_type=\"SHORT\",\n",
    "                                  field_alias='CP2 Station ID',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )    \n",
    "\n",
    "    if 'esrignss_cp2_numsats' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_numsats',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP2 Number of satellites',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrignss_cp2_pdop' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_pdop',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP2 PDOP',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )    \n",
    "\n",
    "    if 'esrignss_cp2_hdop' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_hdop',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP2 HDOP',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )    \n",
    "        \n",
    "    if 'esrignss_cp2_vdop' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_vdop',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP2 VDOP',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )        \n",
    "\n",
    "    if 'esrignss_cp2_avg_h_rms' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_avg_h_rms',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP2 Average horizontal accuracy (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )    \n",
    "\n",
    "    if 'esrignss_cp2_avg_v_rms' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_avg_v_rms',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP2 Average vertical accuracy (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )    \n",
    "\n",
    "    if 'esrignss_cp2_avg_positions' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_avg_positions',\n",
    "                                  field_type=\"SHORT\",\n",
    "                                  field_alias='CP2 Averaged positions',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    if 'esrignss_cp2_h_stddev' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrignss_cp2_h_stddev',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP2 Standard deviation (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )    \n",
    "\n",
    "    if 'esrisnsr_cp2_slope_dist' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrisnsr_cp2_slope_dist',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP2 Laser slope distance (m)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )    \n",
    "        \n",
    "    if 'esrisnsr_cp2_inclination' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrisnsr_cp2_inclination',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP2 Laser inclination (°)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )    \n",
    "        \n",
    "    if 'esrisnsr_cp2_azimuth' not in existing_fields:\n",
    "        arcpy.AddField_management(feature_class,\n",
    "                                  'esrisnsr_cp2_azimuth',\n",
    "                                  field_type=\"DOUBLE\",\n",
    "                                  field_alias='CP2 Laser azimuth (°)',\n",
    "                                  field_is_nullable=\"NULLABLE\"\n",
    "                                  )\n",
    "\n",
    "    # Update GNSS metadata fields with Domains\n",
    "\n",
    "    domain_fields = [field for field in arcpy.ListFields(feature_class) if field.name == 'esrisnsr_offset_method']\n",
    "\n",
    "    for field in domain_fields:\n",
    "        if field.name == 'esrisnsr_offset_method' and not field.domain:\n",
    "            arcpy.AssignDomainToField_management(feature_class, field, 'esri_offset_method_domain')\n",
    "\n",
    "    domain_fields = [field for field in arcpy.ListFields(feature_class) if field.name == 'esrignss_cp1_fixtype']\n",
    "\n",
    "    for field in domain_fields:\n",
    "        if field.name == 'esrignss_cp1_fixtype' and not field.domain:\n",
    "            arcpy.AssignDomainToField_management(feature_class, field, 'esri_fix_type_domain')\n",
    "\n",
    "    domain_fields = [field for field in arcpy.ListFields(feature_class) if field.name == 'esrignss_cp2_fixtype']\n",
    "\n",
    "    for field in domain_fields:\n",
    "        if field.name == 'esrignss_cp2_fixtype' and not field.domain:\n",
    "            arcpy.AssignDomainToField_management(feature_class, field, 'esri_fix_type_domain')\n",
    "\n",
    "    return \"Successfully added offset metadata fields and domains.\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "add_offset_fields(input_feature_class)"
   ]
  }
 ],
 "metadata": {
  "interpreter": {
   "hash": "2f7d5e7f44707441fdf59f9ea9e91257aa419b0f47527e9654c7357783801722"
  },
  "kernelspec": {
   "display_name": "ArcGISPro",
   "language": "Python",
   "name": "python3"
  },
  "language_info": {
   "file_extension": ".py",
   "name": "python",
   "version": "3.7.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}