Model and Methodology for Creating PAG Annual Existing Land Use Layer - 2008 Developed by: Pima Association of Governments Data Services Section (Kristen Zimmerman, Mike Gutierrez, Brooke Pope) Planning Technologies, LLC (Mike Corlett) Date: 07-29-2008 Summary The PAG 2008 existing land use shapefile represents July 1 land uses and residential densities in Pima County. There are approximately 50,000 polygons in this dataset. It builds on the USE code collected by the Pima County Assessor’s office, aggregates geographically continuous parcels with the same land use code (e.g., individual residential lots are aggregated into one ‘neighborhood’), removes gaps in the parcel base by ‘growing’ surrounding land uses, and estimates mid-year residential density using address level building permits issued by Pima County, the City of Tucson, the City of South Tucson, the Town of Oro Valley, the Town of Marana, and the Town of Sahuarita and Census 2000 data. Methodology Geocode Residential Building Permits 1. Create Address Locator using adcounty downloaded at same date of paregion (e.g., ~July 1, 2008). 2. Create Address Locator using adparcel downloaded at same date of paregion. 3. Add unique ID field to original permit table and make a copy of it. 4. Delete all fields except unique ID and address if necessary. 5. Geocode to adcounty 6. Select and export all features with Score >74. Name file with adcounty. 7. Switch Selection. 8. Export selection to dbf. 9. Geocode dbf to adparcel. 10. Export all features with Score > 74. Name file with adparcel. 11. Switch selection to see if record count of un-geocoded addresses is acceptable. If not geocode to paregion either with an address locator or manually. If it’s a limited number sometime it’s easier by hand. 12. Merge adcounty and adparcel exports. Add addresses that matched paregion if applicable. 13. Join back to original table. 14. Check for duplicate IDs and compare totals to original table. Note: Do not geocode to stnetall or any other feature that would create points in the right of way. This is too inaccurate and does not allow location selection of paregion for the purposes of recalculating PAGLUCODE to work properly. All residential building permits should fall within a parcel. Prep Input Shapefiles - All input shapefiles are current as of July 1, 2008. paregion Part I 1. Use the same exact file for land use inventories, update in model. 2. Remove unnecessary fields from paregion. This is most easily done by turning the unwanted fields off and exporting to a new file. Designate and keep this file unchanged and in a safe place because it is the same file to be used for the land use inventories, update in model. 3. Join the PAGLUCODE table to the exported paregion based on the USE field. The PAGLUCODE field must be an integer type. Note: this table can change over time. No USE value is equal to PAGLUCODE 999 Part II * Some building permits fall into parcels that do not have a residential PAGLUCODE. This occurs because of paregion update lag (see important note). This code must be changed to account for these new building activities. 1. Add EXLU with the final Residential building permits to ArcMap. 2. Select and create 2 layers from permit points by PAGTYPE = MF and PAGTYPE = (SFR, MH, TWN). Name layers accordingly. 3. Select by location EXLU that contain MF. Create layer from EXLU selection. 4. Select all non-residential PAGLUCODE values and recalculate PAGLUCODE to 165 5. Repeat process for SFR,MH, TWN and convert PAGLUCODE to 101 Important Note: The selection tool is finicky with large datasets. Make sure areas of 10-15 acres which are obviously non-residential features are not selected when recalculating PAGLUCODEs. Bad selections are also due to paregion update lag on new subdivisions. Do not over-write code on National Parks Forest, Monuments, rights of way, etc. This is why it is important to use the same exact paregion file. pima_all 1. Add field PAGLUCODE and calculate field = 617 PAG major roads inventory 1. Add field GRIDCODE that is equivalent to PAGLUCODE PAG land use inventories 1. Put in a file geodatabase in a feature dataset called Input. 2. Create two other feature datasets called Output and Intermediate. Build the Current Year EXLU Features (using ArcGIS Model Builder) This first step creates a single feature class of all the land use inventories. They must be added in a particular order so more detailed coverages do not over-write less detailed coverages. For example casinos should be updated after Indian Reservations because they are on reservations. The Update tool does not create overlapping topology. It performs erases and appends in one step. 1. Use the Update tool (required ArcInfo license) to “merge or append” the below listed inventories in order from A to AA. Name the output with its associated letter for convenience and save under feature dataset Output. A) ASLD Grazing B) Indian Reservations C) Military D) Open Space E) Post High School schools F) Airports G) MTLUS industrial H) MTLUS office I) MTLUS streetside commercial J) MTLUS shopping centers K) Casinos L) Hospitals M) Lodging facilities N) K-12 schools O) Active landfills P) Eldercare facilities Q) Prisons R) Dedicated government S) Special events T) Transit centers U) Mines V) Utilities W) RV/Mobile Home parks X) Multifamily residential Y) Washes Z) CAP AA) Railroads BB) Paregion (same as used for inventories) CC) pima_all (Note: See Attachment A for inventory description and detailed methodology on how each inventory dataset was constructed.) 6. Make sure there are no features without a PAGLUCODE. If there is, populate it with correct code or delete it if appropriate. 7. Clip CC to pima_all and name EXLU 8. Input EXLU into ET Geowizards tool: Explode Multi Part Features, name as EXLU_2008 Create new model of raster tools. Spatial Analyst is needed for this step. 9. Convert EXLU to grid with the Feature to Raster tool. Field is PAGLUCODE. Name output GRID and put in working geodatabase. Output cell size is 10. 10. Use the Shrink tool on GRID. Name output SHRINK and put in working geodatabase. Number of cells is 75. Zone values are 617 and 999. 11. Weird but necessary step. Use the Resample tool to create SHRINK2 that has a cell size of 5. 12. Convert SHRINK2 to polygon using Raster to Polygon tool. Check the simplify polygons box. Name output EXLU_YEAR (ex: EXLU_2008) in Output feature dataset. 13. Run Update tool on EXLU_YEAR with major roads inventory. 14. Run Check Geometry tool. 15. Run Repair Geometry tool if necessary. 16. Add field LUCode that is equivalent to the GRIDCODE field. Sweet, the Existing Land Use layer is ready!! Attribute the Current Year EXLU Features with Current Year Residential Density (using two Python Scripts) Part I This tool is used to assign dwelling units (DUs), households (HHs), and household population to a land use theme for the BASE YEAR, typically the year 2000 for which a census block group databases is available. The tool does NOT affect the input land use theme -- all results will appear in the COPY, which is written to the same directory as the source. The tool REQUIRES there to be an input existing land use theme containing polygons with appropriate land use codes AND a census block group theme source. The tool will overlay the two creating an identity theme so that we know for each census block group which residential polygons are resident in it. The tool then assigns census data to those polygons, apportioning on the basis of area. However, the method also attempts to assign DUs and HHs to each residential polygon in a way so as to prevent very large residential polygons from absorbing everything in the block group at the expense of other small residential polys. Part II This tool will UPDATE residential data in a land use shape file that you have already created with census 2000 base data, using a point file of lagged residential building permits as a source. The update can be applied periodically over time and as needed to keep the land use database current. DUs appearing in the permits file are added to those already in the land use theme. HHs and population are then estimated for them, using occupancy rates and household size averages that appear in the census block shape file. Attachment A – Land Use Inventory Construction Methodology Land Use Inventories Used in 2008 Existing Land Use Dataset (Note: Every attempt was made to update and accurately portray the actual physical extent of each feature and the PAGLUCODE attribute. Other attribute data are less reliable.) Layer Description Source Feature creation method notes Airports Publicly accessible airports (TIA, Ryan, Marana) previous year's inventory plus any new select paregion that "have their centroid in" old inventory, create centroid from old inventory; join selection to centroid to populate name & address; merge individual shapes then explode ASLD Grazing AZ State trust owned lands that are leased for grazing Most recently available AGIC datasets: "alott" and "own" own is CATEGORY = State Trust; project datasets to match paregion; select parcels and copy, paste and trim to appropriate boundaries to create paregion rectified boundaries, PAGLUCODE = 750 for Agriculture and Ranching; dissolve by PAGLUCODE; explode multi part features Dedicated Gov't Lands owned by local, state, and federal governments PCLIS bldg_gov, firestat, library, police, postoffi, cotparc, parcpima, related websites combine bldg_gov, firestat, library, police, postoffi tables to create one large table to geocode to adcounty; make sure records have correct PAGLUCODE, name address, check cotparc, parcpima and websites to find new records to add, or update cot with parc then update with new combined final shape; attribute with new PAGLUCODE; dissolve by PAGLUCODE then explode OR do select by location paregion "that have their centroid in" old inventory and find new records to add from datasets and websites Indian Reservations Lands held in trust by the Pascua Yaqui and TON Tribes PCLIS limjuris select paregion that "have their centroid in" exported Indian boundaries from limjuris, copy paste and trim remaining parcels for complete coverage; dissolve, explode; PAGLUCODE = 821 Landfills Active landfills previous year's inventory, PCLIS lfil_ex select paregion that "have their centroid in" old inventory, create centroid from old inventory; join selection with centroid to populate name & address; merge individual shapes then explode, check jurisdiction websites, and lfil_ex "status = open" for additional records Mines Active mineral and sand/gravel mines previous year's inventory, PCLIS mines and mines and mines_act select paregion that "have their centroid in" old inventory, create centroid from old inventory; join selection with centroid to populate name & address; copy, paste and trim for complete coverage, check for new records in updated county datasets merge individual shapes then explode, Open Space Lands dedicated to active and passive open space PCLIS preserve, park_all, golf select paregion that "have their centroid in" old inventory, create centroid from old inventory; join selection with centroid to populate name & address; copy, paste and trim for complete coverage, check for new records in updated county datasets, merge individual shapes then explode, attribute PAGLUCODE by dataset Layer Description Source Feature creation method notes Prisons Detention facilities (all ages) previous year's inventory select paregion that "have their centroid in" old inventory, create centroid from old inventory; join selection with centroid to populate name & address; copy, paste and trim for complete coverage, check for possible new records, merge individual shapes then explode, attribute PAGLUCODE Special Events Special events and tourist destinations (e.g., TCC, Rodeo Grounds) previous year's inventory, county and tucson websites, metro atlas Downtown (Map C) select paregion that "have their centroid in" old inventory, create centroid from old inventory; join selection with centroid to populate name & address; copy, paste and trim for complete coverage, check for possible new records, merge individual shapes then explode, attribute PAGLUCODE Transit Centers SunTran transit centers (Rondstat, Tohono T'dai, Laos) previous year's inventory select paregion that "have their centroid in" old inventory, create centroid from old inventory; join selection with centroid to populate name & address; copy, paste and trim for complete coverage, check for possible new records, merge individual shapes then explode, attribute PAGLUCODE Multifamily Residential attached dwelling units (dpx, tpx, qpx, twn, condo, apt, assisted living, boarding house) previous year's inventory, plus any new (sources for new include building permits, MTLUS, appropriate pargion use code) select paregion that "have their centroid in" old inventory, create centroid from old inventory; join selection with centroid to populate table; copy, paste and trim for complete coverage, check for possible new records, merge individual shapes then explode, attribute PAGLUCODE = 165 Lodging hotels, motels, bed/breakfast, resorts previous year's inventory, plus any new (sources for new include building permits, MTLUS, appropriate pargion use code) select paregion that "have their centroid in" old inventory, create centroid from old inventory; join selection with centroid to populate table; copy, paste and trim for complete coverage, check for possible new records, merge individual shapes then explode, attribute PAGLUCODE Major Roads estimated width of major roads PCLIS stnetgdb Query out major roads (MAJOR ROAD = Y), calculate total number of lanes per record (if Direction_of_travel = "Both Directions", total number of lanes = 2*LANE_CAT, if DIRECTION_OF_TRAVEL = 'One Way' or 'One-way reverse', total number of lanes = LANE_CAT. Buffer each feature accordin to number of lanes. (2 lanes= 26 ft buffer, 3 lanes = 32 ft buffer, 4 lanes = 38 ft buffer, 6 lanes = 68 ft buffer). Dissolve on STREET_NAME. PAGLCODE = 616, except I-10, I- 19, State Roads =615. Clean topology (overlaps), dissolve again, clean gaps including highway interchanges, dissolve again. RV/MH Parks Recreational vehicle/mobile home parks previous year's inventory, plus any new (sources fo new include appropriate paregion use code) select paregion that "have their centroid in" old inventory, create centroid from old inventory; join selection with centroid to populate table; copy, paste and trim for complete coverage, check for possible new records, merge individual shapes then explode, attribute PAGLUCODE = 161 Major Washes major water courses PCLIS wash_maj delete all records where WASH_NAME = null. PAGLUCODE = 730. Layer Description Source Feature creation method notes K-12 Schools private, charter, and public k-12 schools previous year's inventory, plus any new (sources for new include Tucson Book of Lists, AZ Daily Star school survey, AZ Dept of Education enrollment report, PCLIS schools) select paregion that "have their centroid in" old inventory, create centroid from old inventory; join selection with centroid to populate table; copy, paste and trim for complete coverage, check for possible new records, merge individual shapes then explode MTLUS_IND MTLUS Industrial facilities, >5000 sq ft MTLUS 3Q report, Table 6.1A select paregion that "have their centroid in" old inventory, create centroid from old inventory; join selection with centroid to populate table; copy, paste and trim for complete coverage, check for possible new records, update vacancy/occupancy attributes, merge individual shapes then explode MTLUS_STSIDECOM MTLUS Strip commercial located along selected corridors MTLUS, original source date was 2006. select paregion that "have their centroid in" old inventory, create centroid from old inventory; join selection with centroid to populate table; copy, paste and trim for complete coverage, check for possible new records (if updated list is obtained), update vacancy/occupancy attributes, merge individual shapes then explode MTLUS_OFFICE MTLUS Office facilities, >5000 sq ft MTLUS 3Q report, Table 5.1a select paregion that "have their centroid in" old inventory, create centroid from old inventory; join selection with centroid to populate table; copy, paste and trim for complete coverage, check for possible new records, update vacancy/occupancy attributes, merge individual shapes then explode MTLUS_SHOPCTR MTLUS Shopping centers, >25,000 sqft MTLUS 2Q report, Table 8.1 select paregion that "have their centroid in" old inventory, create centroid from old inventory; join selection with centroid to populate table; copy, paste and trim for complete coverage, check for possible new records, update vacancy/occupancy attributes, merge individual shapes then explode Cap Canal Above ground Central AZ Project Canal PCLIS cap Digitized physical extent of canal based on PCLIS cap feature overlaid on orthophoto. Geometry does not change from year to year. Casinos Actively operated casinos Tribal/casino websites Digitized physical extent of casino boundary based on orthophoto. Add/edit as changes occur. Railroad Lines Major railroads and spurs PCLIS railroad Digitized physical extent of lines based on PCLIS railroad feature overlaid on orthophoto. Geometry does not change from year to year. 2