Upgrade to Pro — share decks privately, control downloads, hide ads and more …

BETTER URBAN CHANGE DETECTION

BETTER URBAN CHANGE DETECTION

Presented by:
Frank Obusek - Hexagon Geospatial

More Decks by Texas Natural Resources Information System

Other Decks in Technology

Transcript

  1. BETTER URBAN CHANGE DETECTION Andy Zusmanis – Remote Sensing Consultant

    Frank Obusek – Applications Engineer Texas GIS Forum October 22nd, 2014
  2. 2 Why is Change Detection important? •  Change detection is

    one of the most relevant applications of geospatial technology •  Millions of dollars are invested in geospatial content (raster data and point clouds) but those data are often not utilized to their full potential •  Millions of dollars are invested in geospatial technology that can be used to perform change detection
  3. 3 Why is Change Detection important? •  Monitor urban growth

    and mitigate its environmental impacts •  Monitor tree canopy loss and focus tree planting efforts •  More efficient approach to update your GIS database •  Utilizes high resolution ortho imagery •  Utilizes elevation data •  Perfect application for point cloud data
  4. 4 Why is Change Detection important? •  The process has

    been historically complex •  Readily available true color imagery is not ideal •  Elevation data is not accurately co-registered with raster pixels •  Need a geospatial expert to perform change detection •  Because of the complexities, it is not a realistic possibility for most organizations. •  Requires outsourcing
  5. 7 El Paso Study Area – Climate, weather and tree

    mortality •  Hot desert climate (Köppen BWh) •  Hot summers with little humidity •  Mild, dry winters •  9.7 in (250 mm) rainfall per year •  Most during summer from July through September •  Winter Freeze •  February 2011 – damaging winter storm •  Feb 2nd cold episode •  High was 15 °F (−9 °C) – lowest daily maximum on record •  Low was 1 °F (−17 °C) •  Loss of desert vegetation •  Extreme Drought in 2011 •  More than 110 consecutive days without a trace of rain •  Driest on record •  Many trees died
  6. 10 Semi-Global Matching Enabled Change Detection •  4 spectral bands

    (Blue, Green, Red, NIR) •  Two study areas •  Two image dates •  August 26, 2010 •  August 30, 2012 •  Toward the end of the summer wet season A section of study area imagery displayed in ERDAS IMAGINE A section of study area DSM displayed in ERDAS IMAGINE •  SGM Digital Surface Model •  Raster product created •  First 4 bands are the orthorectified image •  SGM DSM used as the elevation source •  5th band Digital Surface Model (DSM) •  5-band Info-Stack
  7. 11 SGM Change Results – New buildings and houses 2010

    2012 Results over 2012 New commercial building in upper left (orange) and some new houses (yellow and orange) 2010 2012 Results over 2012 Tall commercial buildings have been added (orange)
  8. 12 SGM Change Results – New buildings and houses 2010

    2012 Results over 2012 New barracks completed (orange and yellow) 2010 2012 Results over 2012 Highway interchange completed
  9. 13 SGM Change Results – miscellaneous changes 2010 2012 Results

    over 2012 Quarry with significant areas of cut and fill 2010 2012 Results over 2012 Blue and magenta where buildings have been removed. Yellow and orange shows two areas of fill and a new building.
  10. 14 SGM Change Results – Removed trees 2010 2012 Results

    over 2012 Green where trees have been removed 2010 2012 Results over 2012 Many trees have been removed (in green)
  11. 15 SGM Change Results – Removed houses and trees 2010

    2012 Results over 2012 Small barracks and trees have been removed
  12. 16 Solution for Urban Change Detection: Data Prep •  Successful

    Change Detection depends on the input data •  Fly high-resolution digital cameras capable of recording 4-band CIR data (R,G,B,NIR) •  Use a photogrammetric technique called Semi-Global Matching to derive very dense point clouds from the stereo overlapping images •  Produce ortho-rectified data consisting of 5 bands of information •  Standard R, G, B, NIR •  Plus the SGM-derived elevation points as a raster •  Heights are precisely aligned with the image pixels •  Further refine the image bands •  Convert brightness values to reflectance •  Makes comparing images easier •  Calling this 5 band raster file an info-stack in this presentation
  13. 17 Point  Cloud  is  the  “Third”  Type  of  Data  (Eleva:on)

      Vector Point Measurements and Contours have been used historically to represent terrain surfaces. These are combined with break lines to create Triangulated Irregular Networks (TINs) from which surface points can be interpolated. The data representation is a sparse set of highly irregularly space {X,Y,Z} values. Raster They have been converted to gridded formats using various techniques to produce Raster datasets. Delivered as Digital Elevation Models (DEMs). The representation is a dense set of regularly spaced {Z} values. Point Cloud LiDAR data is a collection of points with attributes. The representation is a dense set of semi-regularly spaced {X,Y,Z, Attribute..} values. 17
  14. 18 Creating Point Cloud Data •  Active Laser Sensor -

    LiDAR •  Extracted from Stereo Imagery Stereo Imagery
  15. 25 Spatial Modeling •  Instead of continually reinventing the wheel,

    build your own sophisticated, reusable spatial models and save time and money every time the analyst runs the model. •  Your resident image analyst creates the model once •  Re-usable to your end-users and your customers •  Distributable and Repeatable
  16. 27 Urban Change Spatial Model •  Spatial Model Objectives • 

    Map urban change •  El Paso study area •  High growth area •  Significant tree mortality •  Focus model on changes in buildings and trees •  Re-usable •  Same model and parameters can be applied to other info-stacks •  Ease of use •  Minimize input parameters •  Provide intelligent default parameters
  17. 29 End-user Runs Urban Change Spatial Model •  A dialog

    is automatically created when the model is run •  User can change the values to fine-tune the results •  Lower elevation and area thresholds will find more areas, but will also create more false positives •  Decreasing the amount of noise filtering will speed up processing and find more areas, but will create more false positives and areas with less compact shapes •  Region of interest •  Upper right and lower left map coordinates of the region to process •  If not specified, entire overlap area is used The above dialog is automatically created when the model is run in ERDAS IMAGINE
  18. Web Processing Service (WPS) Publishing the model to a Server

    Running the model using the Geospatial Portal
  19. 32 06 Web Processing Service (WPS) •  Publish the Spatial

    Model as a Web Processing Service (OGC WPS) •  Open an ERDAS APOLLO web service in ERDAS IMAGINE •  Click the ‘Publish to ERDAS APOLLO’ icon in the Spatial Model editor •  Once published, remote managers and decision makers can execute the same process on other data and specific regions •  Don’t need domain specific expertise •  Don’t need an ERDAS IMAGINE license
  20. 33 06 Web Processing Service (WPS) •  Geospatial Portal • 

    Select the published model •  Model description and user input parameters are displayed •  Specify the before and after info- stacks as inputs •  Zoom into the viewer and specify the processing area of interest •  Keep the default or modify the values for model input parameters •  Then click execute to run the model •  Progress meter The Geospatial Portal viewer is zoomed and an area of interest box is defined (shown in red)
  21. 34 06 Web Processing Service (WPS) •  The model runs

    and creates a temporary file •  Click “Show result” and add it to the map •  Yellow and orange areas where ‘buildings’ have been added •  Blue and magenta areas where ‘buildings’ have been removed •  Green areas were trees have been removed Displaying the model results in the Geospatial Portal
  22. 35 06 Web Processing Service (WPS) •  Geospatial Portal • 

    Turn off the info-stack image •  Model results are now visible over an Open Street Maps base map •  Export as a vector dataset •  Integrate with existing GIS data Displaying the model results over Open Street Maps using the Geospatial Portal
  23. 37 Urban Change Spatial Model Description (1) •  Noise filtering

    •  Elevation differenced •  Before vegetation •  After vegetation
  24. 38 Urban Change Spatial Model Description (2) •  Removed trees

    are found •  Buildings and man-made structures are found •  Materials Removal and Displacement
  25. 39 Spatial Model Results •  Processing time •  Ran on

    a laptop system •  Area A (west El Paso) •  25 square miles (7935 rows x 8527 columns @ 1-meter) •  12 minutes to process •  Area B (east El Paso) •  23 square miles (7033 rows x 8684 columns @ 1-meter) •  10 minutes to process •  Resulting Classes •  Categories 1 and 2 (medium or large positive change) •  Areas that are higher and are not vegetation in the after image most often are new houses and buildings •  Buildings are taller and have more category 2 •  Non house or building areas include fill of the ground surface, roads and bridges •  Very little omissions in the results •  Some false positives – mainly due to bad areas in the input surface models •  Using higher resolution data most likely will eliminate most of the false positives
  26. 40 Spatial Model Results •  Result Classes continued •  Categories

    3 and 4 (medium or large negative change) •  Areas that are lower and are not vegetation in the before image most often are removed houses and buildings •  Buildings are taller and have more category 4 •  Non house or building areas are usually cuts of the ground surface •  Very little omissions in the results •  Some false positives – mainly due to bad areas in the input surface models •  Using higher resolution data most likely will eliminate most of the false positives •  Category 5 •  Negative elevation change and vegetation in the before image •  Mainly removed trees •  Some removed trees are omitted •  Reducing the elevation threshold finds more trees but increases false positives •  Using slightly higher resolution data most likely will find most removed trees
  27. 42 Conclusion (1) •  Adding Semi-Global Matching (SGM) to 4-band

    imagery greatly improves change detection capabilities. •  Many ortho providers have the SGM capability but image data distributors and their customers either do not know about it or they are not asking for it. •  SGM Requires Stereo Imagery •  SGM creates true orthos right out of the box for correlated areas.
  28. 43 Conclusions (2) •  Creating a Spatial Model enables the

    distribution of a complex solution to non-expert geospatial users. •  Publishing a Spatial Model as an OGC Web Service enables quick answers to a complex geospatial solution for a broad audience to make decisions that concern developing areas and vegetation loss.